Showing 1–50 of 196 results for author: Padilla, N

ODIN: Identifying Protoclusters and Cosmic Filaments Traced by Ly$α$-emitting Galaxies

Authors: Vandana Ramakrishnan, Kyoung-Soo Lee, Maria Celeste Artale, Eric Gawiser. Yujin Yang, Changbom Park, Robin Ciardullo, Lucia Guaita, Sang Hyeok Im, Seongjae Kim, Ankit Kumar, Jaehyun Lee, Seong-Kook Lee, Byeongha Moon, Nelson Padilla, Alexandra Pope, Roxana Popescu, Hyunmi Song, Paulina Troncoso, Francisco Valdes, Ann Zabludoff

Abstract: To understand the formation and evolution of massive cosmic structures, studying them at high redshift, in the epoch when they formed the majority of their mass is essential. The One-hundred-deg$^2$ DECam Imaging in Narrowbands (ODIN) survey is undertaking the widest-area narrowband program to date, to use Ly$α$-emitting galaxies (LAEs) to trace the large-scale structure (LSS) of the Universe at t… ▽ More To understand the formation and evolution of massive cosmic structures, studying them at high redshift, in the epoch when they formed the majority of their mass is essential. The One-hundred-deg$^2$ DECam Imaging in Narrowbands (ODIN) survey is undertaking the widest-area narrowband program to date, to use Ly$α$-emitting galaxies (LAEs) to trace the large-scale structure (LSS) of the Universe at three cosmic epochs. In this work, we present results at $z$ = 3.1 based on early ODIN data in the COSMOS field. We identify and characterize protoclusters and cosmic filaments using multiple methods and discuss their strengths and weaknesses. We then compare our observations against the IllustrisTNG suite of cosmological hydrodynamical simulations. The two are in excellent agreement, with a similar number and angular size of structures identified above a specified density threshold. We are able to recover the simulated protoclusters with $\log$(M$_{z=0}$/$M_\odot$) $\gtrsim$ 14.4 in $\sim$ 60\% of the cases. With these objects we show that the descendant masses of the protoclusters in our sample can be estimated purely based on our 2D measurements, finding a median $z$ = 0 mass of $\sim10^{14.5}$M$_\odot$. The lack of information on the radial extent of each protocluster introduces a $\sim$0.4~dex uncertainty in its descendant mass. Finally, we show that the recovery of the cosmic web in the vicinity of protoclusters is both efficient and accurate. The similarity of our observations and the simulations imply that our structure selection is likewise robust and efficient, demonstrating that LAEs are reliable tracers of the LSS. △ Less

Submitted 12 June, 2024; originally announced June 2024.

Comments: 26 pages, 18 figures; submitted to ApJ

Impact of Primordial Black Hole Dark Matter on Gas Properties at Very High Redshift: A Semi-Analytical Model

Authors: C. Casanueva-Villarreal, P. B. Tissera, N. Padilla, B. Liu, V. Bromm, S. Pedrosa, L. Bignone, R. Dominguez-Tenreiro

Abstract: Context. Primordial black holes (PBHs) have been proposed as potential candidates for dark matter (DM) and have garnered significant attention in recent years. Aims. Our objective is to delve into the distinct impact of PBHs on gas properties and their potential role in shaping the cosmic structure. Specifically, we aim to analyze the evolving gas properties while considering the presence of accre… ▽ More Context. Primordial black holes (PBHs) have been proposed as potential candidates for dark matter (DM) and have garnered significant attention in recent years. Aims. Our objective is to delve into the distinct impact of PBHs on gas properties and their potential role in shaping the cosmic structure. Specifically, we aim to analyze the evolving gas properties while considering the presence of accreting PBHs with varying monochromatic masses and in different quantities. By studying the feedback effects produced by this accretion, our final goal is to assess the plausibility of PBHs as candidates for DM. Methods. We develop a semi-analytical model which works on top of the CIELO hydrodynamical simulation around $z\sim23$. This model enables a comprehensive analysis of the evolution of gas properties influenced by PBHs. Our focus lies on the temperature and hydrogen abundances, placing specific emphasis on the region closest to the halo center. We explore PBH masses of $1$, $33$, and $100~\Msun$, located within mass windows where a substantial fraction of DM could exist in the form of PBHs. We investigate various DM fractions composed of these PBHs ($f_{\rm{PBH}}>10^{-4}$). Results. Our findings suggest that the existence of PBHs with masses of $1~\Msun$ and fractions greater than or equal to approximately $10^{-2}$ would be ruled out due to the significant changes induced in gas properties. The same applies to PBHs with a mass of $33~\Msun$ and $100~\Msun$ and fractions greater than approximately $10^{-3}$. These effects are particularly pronounced in the region nearest to the halo center, potentially leading to delayed galaxy formation within haloes. △ Less

Submitted 3 May, 2024; originally announced May 2024.

Structure formation with primordial black holes to alleviate early star formation tension revealed by JWST

Authors: P. E. Colazo, F. Stasyszyn, N. Padilla

Abstract: This Letter explores the potential role of primordial black holes (PBHs) to address cosmological tensions as the presence of more massive than expected galaxies at high redshifts, as indicated by recent James Webb Space Telescope observations. Motivated by inflation models that enhance the power at scales beyond the observable range that produce PBHs with Schechter-like mass functions, we aim to e… ▽ More This Letter explores the potential role of primordial black holes (PBHs) to address cosmological tensions as the presence of more massive than expected galaxies at high redshifts, as indicated by recent James Webb Space Telescope observations. Motivated by inflation models that enhance the power at scales beyond the observable range that produce PBHs with Schechter-like mass functions, we aim to explain the excess of high redshift galaxies via a modification of the $Λ$ cold dark matter power spectrum that consists in adding (i) a blue spectral index $n_b$ at $k_{\text{piv}}=10/$Mpc and (ii) Poisson and isocurvature contributions from massive PBHs that only make up $0.5\%$ of the dark matter.We simulated these models using the SWIFT code and find an increased abundance of high redshift galaxies in simulations that include PBHs. We compared these models to estimates from James Webb Space Telescope observations.Unlike the $Λ$ cold dark matter model, the inclusion of PBHs allowed us to reproduce the the observations with reasonable values for the star formation efficiency. Furthermore, the power spectra we adopted potentially produce PBHs that can serve as seeds for supermassive black holes with masses $7.57 \times 10^4 M_{\odot}$. △ Less

Submitted 19 April, 2024; originally announced April 2024.

Comments: 5 pages, 2 figures, Accepted for publication in A&A

ODIN: Improved Narrowband Ly$α$ Emitter Selection Techniques for $z$ = 2.4, 3.1, and 4.5

Authors: Nicole M. Firestone, Eric Gawiser, Vandana Ramakrishnan, Kyoung-Soo Lee, Francisco Valdes, Changbom Park, Yujin Yang, Robin Ciardullo, María Celeste Artale, Barbara Benda, Adam Broussard, Lana Eid, Rameen Farooq, Caryl Gronwall, Lucia Guaita, Stephen Gwyn, Ho Seong Hwang, Sang Hyeok Im, Woong-Seob Jeong, Shreya Karthikeyan, Dustin Lang, Byeongha Moon, Nelson Padilla, Marcin Sawicki, Eunsuk Seo , et al. (3 additional authors not shown)

Abstract: Lyman-Alpha Emitting galaxies (LAEs) are typically young, low-mass, star-forming galaxies with little extinction from interstellar dust. Their low dust attenuation allows their Ly$α$ emission to shine brightly in spectroscopic and photometric observations, providing an observational window into the high-redshift universe. Narrowband surveys reveal large, uniform samples of LAEs at specific redshif… ▽ More Lyman-Alpha Emitting galaxies (LAEs) are typically young, low-mass, star-forming galaxies with little extinction from interstellar dust. Their low dust attenuation allows their Ly$α$ emission to shine brightly in spectroscopic and photometric observations, providing an observational window into the high-redshift universe. Narrowband surveys reveal large, uniform samples of LAEs at specific redshifts that probe large scale structure and the temporal evolution of galaxy properties. The One-hundred-deg$^2$ DECam Imaging in Narrowbands (ODIN) utilizes three custom-made narrowband filters on the Dark Energy Camera (DECam) to discover LAEs at three equally spaced periods in cosmological history. In this paper, we introduce the hybrid-weighted double-broadband continuum estimation technique, which yields improved estimation of Ly$α$ equivalent widths. Using this method, we discover 6339, 6056, and 4225 LAE candidates at $z =$ 2.4, 3.1, and 4.5 in the extended COSMOS field ($\sim$9 deg$^2$). We find that [O II] emitters are a minimal contaminant in our LAE samples, but that interloping Green Pea-like [O III] emitters are important for our redshift 4.5 sample. We introduce an innovative method for identifying [O II] and [O III] emitters via a combination of narrowband excess and galaxy colors, enabling their study as separate classes of objects. We present scaled median stacked SEDs for each galaxy sample, revealing the overall success of our selection methods. We also calculate rest-frame Ly$α$ equivalent widths for our LAE samples and find that the EW distributions are best fit by exponential functions with scale lengths of $w_0$ = 55 $\pm$ 1, 65 $\pm$ 1, and 62 $\pm$ 1 Angstroms, respectively. △ Less

Submitted 26 December, 2023; originally announced December 2023.

Comments: Submitted to ApJ, 23 pages, 14 figures, comments welcome

Galaxy evolution in modified gravity simulations: using passive galaxies to constrain gravity with upcoming surveys

Authors: Diego Pallero, Facundo A. Gómez, Nelson D. Padilla, Yara L. Jaffé, Carlton M. Baugh, Baojiu Li, César Hernández-Aguayo, Christian Arnold

Abstract: We present a quantitative analysis of the properties of galaxies and structures evolving in universes dominated by different modified gravitational models, including two variants of the f(R)-gravity (F) and two of the Dvali-Gabdadze-Poratti (N) braneworld model, which respectively feature the chameleon and Vainshtein screening mechanisms. Using the Simulation HYdrodynamics BeyONd Einstein (SHYBONE… ▽ More We present a quantitative analysis of the properties of galaxies and structures evolving in universes dominated by different modified gravitational models, including two variants of the f(R)-gravity (F) and two of the Dvali-Gabdadze-Poratti (N) braneworld model, which respectively feature the chameleon and Vainshtein screening mechanisms. Using the Simulation HYdrodynamics BeyONd Einstein (SHYBONE) cosmological hydrodynamical full-physics simulations suite, we study the departures in the properties of galaxies residing in different environments with respect to the standard model (GR). Using two different criteria to compare, we find that structures formed within modified gravity tend to show a denser gas density profile than their GR counterparts. Within the different modified gravity models, N1 and F5 gravity models show greater departures from the standard model, with gas density profiles $ρ_{\rm IGM} \geq 30\%$ denser in the outskirts for the N1 model, and in the inner parts for the F5 model. Additionally, we find that haloes evolving in MG universes show, in general, larger quenched fractions than GR, reaching up to $20\%$ larger quenching fractions in F5 regardless of the stellar mass of the galaxy. With respect to the other models, F6, N1 and N5 show slightly larger quenched fractions, but no strong differences can be found. These results directly impact the colour distribution of galaxies, making them in MG models redder and older than their GR counterparts. Like GR, once the environment starts to play a role, galaxies rapidly get quenched and the differences between models vanish. △ Less

Submitted 3 October, 2023; originally announced October 2023.

Comments: 17 pages, 12 figures, submitted to MNRAS

Cosmological constraints from density-split clustering in the BOSS CMASS galaxy sample

Authors: Enrique Paillas, Carolina Cuesta-Lazaro, Will J. Percival, Seshadri Nadathur, Yan-Chuan Cai, Sihan Yuan, Florian Beutler, Arnaud de Mattia, Daniel Eisenstein, Daniel Forero-Sanchez, Nelson Padilla, Mathilde Pinon, Vanina Ruhlmann-Kleider, Ariel G. Sánchez, Georgios Valogiannis, Pauline Zarrouk

Abstract: We present a clustering analysis of the BOSS DR12 CMASS galaxy sample, combining measurements of the galaxy two-point correlation function and density-split clustering down to a scale of $1\,h^{-1}{\rm Mpc}$. Our theoretical framework is based on emulators trained on high-fidelity mock galaxy catalogues that forward model the cosmological dependence of the clustering statistics within an extended-… ▽ More We present a clustering analysis of the BOSS DR12 CMASS galaxy sample, combining measurements of the galaxy two-point correlation function and density-split clustering down to a scale of $1\,h^{-1}{\rm Mpc}$. Our theoretical framework is based on emulators trained on high-fidelity mock galaxy catalogues that forward model the cosmological dependence of the clustering statistics within an extended-$Λ$CDM framework, including redshift-space and Alcock-Paczynski distortions. Our base-$Λ$CDM analysis finds $ω_{\rm cdm} = 0.1201\pm 0.0022$, $σ_8 = 0.792\pm 0.034$, and $n_s = 0.970\pm 0.018$, corresponding to $fσ_8 = 0.462\pm 0.020$ at $z \approx 0.525$, which is in agreement with Planck 2018 predictions and various clustering studies in the literature. We test single-parameter extensions to base-$Λ$CDM, varying the running of the spectral index, the dark energy equation of state, and the density of massless relic neutrinos, finding no compelling evidence for deviations from the base model. We model the galaxy-halo connection using a halo occupation distribution framework, finding signatures of environment-based assembly bias in the data. We validate our pipeline against mock catalogues that match the clustering and selection properties of CMASS, showing that we can recover unbiased cosmological constraints even with a volume 84 times larger than the one used in this study. △ Less

Submitted 2 October, 2023; v1 submitted 28 September, 2023; originally announced September 2023.

Comments: Submitted to MNRAS. Source code to generate the figures available in the captions. Updated to add missing references and fix legend of Fig. 6

SUNBIRD: A simulation-based model for full-shape density-split clustering

Authors: Carolina Cuesta-Lazaro, Enrique Paillas, Sihan Yuan, Yan-Chuan Cai, Seshadri Nadathur, Will J. Percival, Florian Beutler, Arnaud de Mattia, Daniel Eisenstein, Daniel Forero-Sanchez, Nelson Padilla, Mathilde Pinon, Vanina Ruhlmann-Kleider, Ariel G. Sánchez, Georgios Valogiannis, Pauline Zarrouk

Abstract: Combining galaxy clustering information from regions of different environmental densities can help break cosmological parameter degeneracies and access non-Gaussian information from the density field that is not readily captured by the standard two-point correlation function (2PCF) analyses. However, modelling these density-dependent statistics down to the non-linear regime has so far remained cha… ▽ More Combining galaxy clustering information from regions of different environmental densities can help break cosmological parameter degeneracies and access non-Gaussian information from the density field that is not readily captured by the standard two-point correlation function (2PCF) analyses. However, modelling these density-dependent statistics down to the non-linear regime has so far remained challenging. We present a simulation-based model that is able to capture the cosmological dependence of the full shape of the density-split clustering (DSC) statistics down to intra-halo scales. Our models are based on neural-network emulators that are trained on high-fidelity mock galaxy catalogues within an extended-$Λ$CDM framework, incorporating the effects of redshift-space, Alcock-Paczynski distortions and models of the halo-galaxy connection. Our models reach sub-percent level accuracy down to $1\,h^{-1}{\rm Mpc}$ and are robust against different choices of galaxy-halo connection modelling. When combined with the galaxy 2PCF, DSC can tighten the constraints on $ω_{\rm cdm}$, $σ_8$, and $n_s$ by factors of 2.9, 1.9, and 2.1, respectively, compared to a 2PCF-only analysis. DSC additionally puts strong constraints on environment-based assembly bias parameters. Our code is made publicly available on Github. △ Less

Submitted 2 October, 2023; v1 submitted 28 September, 2023; originally announced September 2023.

Comments: Submitted to MNRAS. Source code to generate the figures available in the captions. Updated to add missing references

Gravitational to Coulomb force ratio and the origin of the Cosmic magnetic field

Authors: Nelson D. Padilla, Juan Racker, Ignacio J. Araya, Federico Stasyszyn

Abstract: Aims: The origin of the seeds of galactic magnetic fields is a subject that remains under debate; here we will explore a new possible mechanism that could provide such seed from charge asymmetries in slowly rotating protogalaxies. Methods: We use current knowledge of galaxy formation and evolution to estimate that a charge imbalance of $1$ every $\sim10^{38\pm 5}$ charge carriers in slowly rotatin… ▽ More Aims: The origin of the seeds of galactic magnetic fields is a subject that remains under debate; here we will explore a new possible mechanism that could provide such seed from charge asymmetries in slowly rotating protogalaxies. Methods: We use current knowledge of galaxy formation and evolution to estimate that a charge imbalance of $1$ every $\sim10^{38\pm 5}$ charge carriers in slowly rotating protogalaxies can provide an alternative mechanism to produce such seeds. Results: The gravitationally induced charge imbalance within plasmas in hydrostatic equilibrium lies in the upper range of this estimation, and is higher than possible surviving primordial charge asymmetries. We conclude that this known, but unappreciated thermodynamic property of plasmas in gravitational fields, may have important cosmological consequences. △ Less

Submitted 22 April, 2024; v1 submitted 26 September, 2023; originally announced September 2023.

Comments: letter, comments welcome

The electric charge of black holes within galaxies

Authors: Nelson D. Padilla, Ignacio J. Araya, Federico Stasyszyn

Abstract: We present improved estimates of the electric charge that black holes could hold when these are embedded in the ionised plasma within galaxies. We have implemented the spontaneous emission of charges of opposite sign to that of the black hole via athermal Hawking evaporation, including its dependence on black hole spin, and we have estimated the equilibrium charge that arises as this charge loss i… ▽ More We present improved estimates of the electric charge that black holes could hold when these are embedded in the ionised plasma within galaxies. We have implemented the spontaneous emission of charges of opposite sign to that of the black hole via athermal Hawking evaporation, including its dependence on black hole spin, and we have estimated the equilibrium charge that arises as this charge loss is balanced by the continuous accretion of charges from the surrounding plasma. The resulting charge can be several orders of magnitude lower than previously estimated upper limits, but it can surpass the pair production limit noted by Gibbons (1974) by a margin that increases with the amplitude of the black hole spin and the density of the plasma. We also implement a calculation for the net charge of galaxies due to their stellar black holes and to primordial black holes in the case these make up a fraction of the dark matter. We find that these charges lie within the range that would produce adequate magnetic field seeds at the onset of galaxy formation. △ Less

Submitted 16 April, 2024; v1 submitted 26 September, 2023; originally announced September 2023.

Comments: Accepted for publication in JCAP

The One-hundred-deg^2 DECam Imaging in Narrowbands (ODIN): Survey Design and Science Goals

Authors: Kyoung-Soo Lee, Eric Gawiser, Changbom Park, Yujin Yang, Francisco Valdes, Dustin Lang, Vandana Ramakrishnan, Byeongha Moon, Nicole Firestone, Stephen Appleby, Maria Celeste Artale, Moira Andrews, Franz E. Bauer, Barbara Benda, Adam Broussard, Yi-Kuan Chiang, Robin Ciardullo, Arjun Dey, Rameen Farooq, Caryl Gronwall, Lucia Guaita, Yun Huang, Ho Seong Hwang, Sanghyeok Im, Woong-Seob Jeong , et al. (17 additional authors not shown)

Abstract: We describe the survey design and science goals for ODIN (One-hundred-deg^2 DECam Imaging in Narrowbands), a NOIRLab survey using the Dark Energy Camera (DECam) to obtain deep (AB~25.7) narrow-band images over an unprecedented area of sky. The three custom-built narrow-band filters, N419, N501, and N673, have central wavelengths of 419, 501, and 673 nm and respective full-widthat-half-maxima of 7.… ▽ More We describe the survey design and science goals for ODIN (One-hundred-deg^2 DECam Imaging in Narrowbands), a NOIRLab survey using the Dark Energy Camera (DECam) to obtain deep (AB~25.7) narrow-band images over an unprecedented area of sky. The three custom-built narrow-band filters, N419, N501, and N673, have central wavelengths of 419, 501, and 673 nm and respective full-widthat-half-maxima of 7.2, 7.4, and 9.8 nm, corresponding to Lya at z=2.4, 3.1, and 4.5 and cosmic times of 2.8, 2.1, and 1.4 Gyr, respectively. When combined with even deeper, public broad-band data from Hyper Suprime-Cam, DECam, and in the future, LSST, the ODIN narrow-band images will enable the selection of over 100,000 Lya-emitting (LAE) galaxies at these epochs. ODIN-selected LAEs will identify protoclusters as galaxy overdensities, and the deep narrow-band images enable detection of highly extended Lya blobs (LABs). Primary science goals include measuring the clustering strength and dark matter halo connection of LAEs, LABs, and protoclusters, and their respective relationship to filaments in the cosmic web. The three epochs allow the redshift evolution of these properties to be determined during the period known as Cosmic Noon, where star formation was at its peak. The two narrow-band filter wavelengths are designed to enable interloper rejection and further scientific studies by revealing [O II] and [O III] at z=0.34, Lya and He II 1640 at z=3.1, and Lyman continuum plus Lya at z=4.5. Ancillary science includes similar studies of the lower-redshift emission-line galaxy samples and investigations of nearby star-forming galaxies resolved into numerous [O III] and [S II] emitting regions. △ Less

Submitted 18 September, 2023; originally announced September 2023.

Comments: 21 pages, 8 figures, submitted to Astrophysical Journal

A new test of gravity -- II: Application of marked correlation functions to luminous red galaxy samples

Authors: Joaquin Armijo, Carlton M. Baugh, Peder Norberg, Nelson D. Padilla

Abstract: We apply the marked correlation function test proposed by Armijo et al. (Paper I) to samples of luminous red galaxies (LRGs) from the final data release of the Sloan Digital Sky Survey (SDSS) III. The test assigns a density-dependent mark to galaxies in the estimation of the projected marked correlation function. Two gravity models are compared: general relativity (GR) and $f(R)$ gravity. We build… ▽ More We apply the marked correlation function test proposed by Armijo et al. (Paper I) to samples of luminous red galaxies (LRGs) from the final data release of the Sloan Digital Sky Survey (SDSS) III. The test assigns a density-dependent mark to galaxies in the estimation of the projected marked correlation function. Two gravity models are compared: general relativity (GR) and $f(R)$ gravity. We build mock catalogues which, by construction, reproduce the measured galaxy number density and two-point correlation function of the LRG samples, using the halo occupation distribution model (HOD). A range of HOD models give acceptable fits to the observational constraints, and this uncertainty is fed through to the error in the predicted marked correlation functions. The uncertainty from the HOD modelling is comparable to the sample variance for the SDSS-III LRG samples. Our analysis shows that current galaxy catalogues are too small for the test to distinguish a popular $f(R)$ model from GR. However, upcoming surveys with a better measured galaxy number density and smaller errors on the two-point correlation function, or a better understanding of galaxy formation, may allow our method to distinguish between viable gravity models. △ Less

Submitted 9 March, 2024; v1 submitted 18 September, 2023; originally announced September 2023.

Comments: 7 pages, 4 Figures. Published by the Monthly Notices of the Royal Astronomical Society journal

Journal ref: Monthly Notices of the Royal Astronomical Society, Volume 528, Issue 4, March 2024, Pages 6631

Not Hydro: Using Neural Networks to estimate galaxy properties on a Dark-Matter-Only simulation

Authors: Cristian Hernández Cuevas, Roberto E. González, Nelson D. Padilla

Abstract: Using data from TNG300-2, we train a neural network (NN) to recreate the stellar mass ($M^*$) and star formation rate (SFR) of central galaxies in a dark-matter-only simulation. We consider 12 input properties from the halo and sub-halo hosting the galaxy and the near environment. $M^*$ predictions are robust, but the machine does not fully reproduce its scatter. The same happens for SFR, but the… ▽ More Using data from TNG300-2, we train a neural network (NN) to recreate the stellar mass ($M^*$) and star formation rate (SFR) of central galaxies in a dark-matter-only simulation. We consider 12 input properties from the halo and sub-halo hosting the galaxy and the near environment. $M^*$ predictions are robust, but the machine does not fully reproduce its scatter. The same happens for SFR, but the predictions are not as good as for $M^*$. We chained neural networks, improving the predictions on SFR to some extent. For SFR, we time-averaged this value between $z=0$ and $z=0.1$, which improved results for $z=0$. Predictions of both variables have trouble reproducing values at lower and higher ends. We also study the impact of each input variable in the performance of the predictions using a leave-one-covariate-out approach, which led to insights about the physical and statistical relation between input variables. In terms of metrics, our machine outperforms similar studies, but the main discoveries in this work are not linked with the quality of the predictions themselves, but to how the predictions relate to the input variables. We find that previously studied relations between physical variables are meaningful to the machine. We also find that some merger tree properties strongly impact the performance of the machine. %We highlight the value of machine learning (ML) methods in helping understand the information contained in different variables, since with its help we were able to obtain useful insights resulting from studying the impact of input variables on the resulting behaviour of galaxy properties. We conclude that ML models are useful tools to understand the significance of physical different properties and their impact on target characteristics, as well as strong candidates for potential simulation methods. △ Less

Submitted 24 July, 2023; originally announced July 2023.

Comments: 17 pages, 16 figures, to be published in MNRAS

Active galactic nuclei and gravitational redshifts

Authors: N. D. Padilla, S. Carneiro, J. Chaves-Montero, C. J. Donzelli, C. Pigozzo, P. Colazo, J. S. Alcaniz

Abstract: Context: Gravitational redshift is a classical effect of Einstein's General Relativity, already measured in stars, quasars and clusters of galaxies. Aims: We here aim to identify the signature of gravitational redshift in the emission lines of active galaxies due to supermassive black holes, and compare to what is found for inactive galaxies. Methods: Using the virial theorem, we estimate gravitat… ▽ More Context: Gravitational redshift is a classical effect of Einstein's General Relativity, already measured in stars, quasars and clusters of galaxies. Aims: We here aim to identify the signature of gravitational redshift in the emission lines of active galaxies due to supermassive black holes, and compare to what is found for inactive galaxies. Methods: Using the virial theorem, we estimate gravitational redshifts for quasars from the 14th data release of the Sloan Digital Sky Survey, and compare these with measured ones from the difference between the redshifts of emission lines of Sydney Australian Astronomical Observatory Multi-object Integral Field (SAMI) galaxies in central and outer annuli of their integral field spectra. Results: Firstly, from the full width at half maximum of $H_β$ lines of 57 Seyfert type I galaxies of the AGN Black Hole Mass Database, we derive a median gravitational redshift $z_g = 1.18 \times 10^{-4}$. Expanding this analysis to 86755 quasars from DR14 of SDSS we have a median value $z_g = 1.52 \times 10^{-4}$. Then, by comparing the redshifts of $34$ lines measured at central and outer regions of LINER galaxies in the SAMI survey we obtain $z_g = (0.68 \pm 0.09) \times 10^{-4}$, which increases to $z_g = (1.0 \pm 0.1) \times 10^{-4}$ when using $H_α$ and $H_β$ lines. These numbers are compatible with central black holes of $\approx 10^9$ solar masses and broad line regions of $\approx 1$pc. For non-AGN galaxies the gravitational redshift is compatible with zero. △ Less

Submitted 20 December, 2023; v1 submitted 24 April, 2023; originally announced April 2023.

Comments: Accepted for publication in A&A, 7 pages, 6 figures

Journal ref: Astronomy & Astrophysics 683, A120 (2024)

A new test of gravity -- I: Introduction to the method

Authors: Joaquin Armijo, Carlton M. Baugh, Peder Norberg, Nelson D. Padilla

Abstract: We introduce a new scheme based on the marked correlation function to probe gravity using the large-scale structure of the Universe. We illustrate our approach by applying it to simulations of the metric-variation $f(R)$ modified gravity theory and general relativity (GR). The modifications to the equations in $f(R)$ gravity lead to changes in the environment of large-scale structures that could,… ▽ More We introduce a new scheme based on the marked correlation function to probe gravity using the large-scale structure of the Universe. We illustrate our approach by applying it to simulations of the metric-variation $f(R)$ modified gravity theory and general relativity (GR). The modifications to the equations in $f(R)$ gravity lead to changes in the environment of large-scale structures that could, in principle, be used to distinguish this model from GR. Applying the Monte Carlo Markov Chain algorithm, we use the observed number density and two-point clustering to fix the halo occupation distribution (HOD) model parameters and build mock galaxy catalogues from both simulations. To generate a mark for galaxies when computing the marked correlation function we estimate the local density using a Voronoi tessellation. Our approach allows us to isolate the contribution to the uncertainty in the predicted marked correlation function that arises from the range of viable HOD model parameters, in addition to the sample variance error for a single set of HOD parameters. This is critical for assessing the discriminatory power of the method. In a companion paper we apply our new scheme to a current large-scale structure survey. △ Less

Submitted 9 March, 2024; v1 submitted 12 April, 2023; originally announced April 2023.

Comments: 11 pages, 7 figures. Accepted for publication in the Monthly Notcies of the Royal Astronomical Society

Redshift evolution of the dark matter haloes shapes

Authors: Pedro Cataldi, Susana Pedrosa, Patricia Tissera, Maria Celeste Artale, Nelson Padilla, Rosa Dominguez-Tenreiro, Lucas Bignone, Roberto Gonzalez, Leonardo Pellizza

Abstract: In this work, we aim at investigating the morphology evolution of Milky Way mass-like dark matter haloes selected from the CIELO and IllustrisTNG Projects. The connection between halo shapes and their environment has been studied in previous works at z=0 but their connection remains yet to be fully understood. We focus on the evolution across cosmic time of the halo shapes and the relation with th… ▽ More In this work, we aim at investigating the morphology evolution of Milky Way mass-like dark matter haloes selected from the CIELO and IllustrisTNG Projects. The connection between halo shapes and their environment has been studied in previous works at z=0 but their connection remains yet to be fully understood. We focus on the evolution across cosmic time of the halo shapes and the relation with the infalling material, using hydrodynamical simulations. Our findings show that haloes tend to be more triaxial at earlier times as a consequence of stronger accretion in the direction of the filaments. As the haloes evolve towards a dominant isotropic accretion mode and relaxation, their shape at 20 percent of the virial mass becomes more spherical. In agreement with previous results, baryons have an important effect within the inner regions of the haloes, driving them from triaxial to rounder shapes. We also find a correlation between the strength of the quadrupole infalling mode and the degree of ellipticity of the haloes: as the filament strength decreases steadily with redshift, the haloes became more spherical and less elliptical. △ Less

Submitted 17 February, 2023; originally announced February 2023.

ODIN: Where Do Lyman-alpha Blobs Live? Contextualizing Blob Environments within the Large-Scale Structure

Authors: Vandana Ramakrishnan, Byeongha Moon, Sang Hyeok Im, Rameen Farooq, Kyoung-Soo Lee, Eric Gawiser, Yujin Yang, Changbom Park, Ho Seong Hwang, Francisco Valdes, Maria Celeste Artale, Robin Ciardullo, Arjun Dey, Caryl Gronwall, Lucia Guaita, Woong-Seob Jeong, Nelson Padilla, Akriti Singh, Ann Zabludoff

Abstract: While many Lyman-alpha Blobs (LABs) are found in and around several well-known protoclusters at high redshift, how they trace the underlying large-scale structure is still poorly understood. In this work, we utilize 5,352 Lyman-alpha emitters (LAEs) and 129 LABs at z=3.1 identified over a $\sim$ 9.5 sq. degree area in early data from the ongoing One-hundred-deg$^2$ DECam Imaging in Narrowbands (OD… ▽ More While many Lyman-alpha Blobs (LABs) are found in and around several well-known protoclusters at high redshift, how they trace the underlying large-scale structure is still poorly understood. In this work, we utilize 5,352 Lyman-alpha emitters (LAEs) and 129 LABs at z=3.1 identified over a $\sim$ 9.5 sq. degree area in early data from the ongoing One-hundred-deg$^2$ DECam Imaging in Narrowbands (ODIN) survey to investigate this question. Using LAEs as tracers of the underlying matter distribution, we identify overdense structures as galaxy groups, protoclusters, and filaments of the cosmic web. We find that LABs preferentially reside in regions of higher-than-average density and are located in closer proximity to overdense structures, which represent the sites of protoclusters and their substructures. Moreover, protoclusters hosting one or more LABs tend to have a higher descendant mass than those which do not. Blobs are also strongly associated with filaments of the cosmic web, with $\sim$ 70% of the population being within a projected distance of 2.4 pMpc from a filament. We show that the proximity of LABs to protoclusters is naturally explained by their association with filaments as large cosmic structures are where many filaments converge. The contiguous wide-field coverage of the ODIN survey allows us for the first time to firmly establish a connection between LABs as a population and their environment. △ Less

Submitted 7 July, 2023; v1 submitted 15 February, 2023; originally announced February 2023.

Comments: 25 pages, 17 figures; accepted for publication in ApJ

Guess the cheese flavour by the size of its holes: A cosmological test using the abundance of Popcorn voids

Authors: Dante J. Paz, Carlos M. Correa, Sebastián R. Gualpa, Andres N. Ruiz, Carlos S. Bederián, R. Dario Graña, Nelson D. Padilla

Abstract: We present a new definition of cosmic void and a publicly available code with the algorithm that implements it. Underdense regions are defined as free-form objects, called popcorn voids, made from the union of spheres of maximum volume with a given joint integrated underdensity contrast.The method is inspired by the excursion-set theory and consequently no rescaling processing is needed, the remov… ▽ More We present a new definition of cosmic void and a publicly available code with the algorithm that implements it. Underdense regions are defined as free-form objects, called popcorn voids, made from the union of spheres of maximum volume with a given joint integrated underdensity contrast.The method is inspired by the excursion-set theory and consequently no rescaling processing is needed, the removal of overlapping voids and objects with sizes below the shot noise threshold is inherent in the algorithm. The abundance of popcorn voids in the matter field can be fitted using the excursion-set theory provided the relationship between the linear density contrast of the barrier and the threshold used in void identification is modified relative to the spherical evolution model. We also analysed the abundance of voids in biased tracer samples in redshift space. We show how the void abundance can be used to measure the geometric distortions due to the assumed fiducial cosmology, in a test similar to an Alcock-Paczyński test. Using the formalism derived from previous works, we show how to correct the abundance of popcorn voids for redshift-space distortion effects. Using this treatment, in combination with the excursion-set theory, we demonstrate the feasibility of void abundance measurements as cosmological probes. We obtain unbiased estimates of the target parameters, albeit with large degeneracies in the parameter space. Therefore, we conclude that the proposed test in combination with other cosmological probes has potential to improve current cosmological parameter constraints. △ Less

Submitted 3 April, 2023; v1 submitted 13 December, 2022; originally announced December 2022.

Comments: Updated manuscript sent to the MNRAS after referee report: 16 pages, 8 figures. Corrections were made to Fig. 4, some related conclusions were modified. The main conclusions remain unchanged

Extended primordial black hole mass functions with a spike

Authors: Juan Magaña, Marco San Martín, Joaquín Sureda, Marcelo E. Rubio, Ignacio J. Araya, Nelson D. Padilla

Abstract: We introduce a modification of the Press-Schechter formalism aimed to derive general mass functions for primordial black holes (PBHs). In this case, we start from primordial power spectra (PPS) which include a monochromatic spike, typical of ultra slow-roll inflation models. We consider the PBH formation as being associated to the amplitude of the spike on top of the linear energy density fluctuat… ▽ More We introduce a modification of the Press-Schechter formalism aimed to derive general mass functions for primordial black holes (PBHs). In this case, we start from primordial power spectra (PPS) which include a monochromatic spike, typical of ultra slow-roll inflation models. We consider the PBH formation as being associated to the amplitude of the spike on top of the linear energy density fluctuations, coming from a PPS with a blue index. By modelling the spike with a log-normal function, we study the properties of the resulting mass function spikes, and compare these to the underlying extended mass distributions. When the spike is at PBH masses which are much lower than the exponential cutoff of the extended distribution, very little mass density is held by the PBHs within the spike, and it is not ideal to apply the Press-Schechter formalism in this case as the resulting characteristic overdensity is too different from the threshold for collapse. It is more appropriate to do so when the spike mass is similar to, or larger than the cutoff mass. Additionally, it can hold a similar mass density as the extended part. Such particular mass functions also contain large numbers of small PBHs, especially if stable PBH relics are considered, and they can provide $\sim 1000M_\odot$ seeds for the supermassive black holes at the centres of present-day galaxies. The constraints on the fraction of dark matter in PBHs for monochromatic mass functions are somewhat relaxed when there is an additional underlying extended distribution of masses. △ Less

Submitted 8 January, 2023; v1 submitted 27 July, 2022; originally announced July 2022.

Comments: 14 pages, 4 figures, 1 table, refereces added, discussion extended and conclusion unchanged. Accepted for publication in MNRAS

Fingerprints of modified gravity on galaxies in voids

Authors: Pedro Cataldi, Susana Pedrosa, Nelson Padilla, Susana Landau, Christian Arnold, Baojiu Li

Abstract: We search for detectable signatures of f(R) gravity and its chameleon screening mechanism in the baryonic and dark matter (DM) properties of simulated void galaxies. The enhancement of the gravitational acceleration can have a meaningful impact on the scaling relations as well as on the halo morphology. The galaxy rotational velocity field (calculated with the velocity of the gas disc and the acce… ▽ More We search for detectable signatures of f(R) gravity and its chameleon screening mechanism in the baryonic and dark matter (DM) properties of simulated void galaxies. The enhancement of the gravitational acceleration can have a meaningful impact on the scaling relations as well as on the halo morphology. The galaxy rotational velocity field (calculated with the velocity of the gas disc and the acceleration fields) deviates from the typical values of the Tully-Fisher Relation (TFR) in GR. For a given stellar mass, f(R) gravity tends to produce greater maximum velocities. On the other hand, the mass in haloes in f(R) gravity is more concentrated than their counterparts in GR. This trend changes when the concentration is calculated with the dynamical density profile, which takes into account the unscreened outer regions of the halo. Stellar discs interact with the overall potential well in the central regions, modifying the morphology of the screening regions and reshaping them. We find a trend for galaxies with a more dominant stellar disc to deviate further from round screening regions. We find that small haloes are less triaxial and more round in f(R) than their GR counterparts. The difference between halo morphology becomes smaller in f(R) haloes whose inner regions are screened. These results suggest possible observables that could unveil modified gravity effects on galaxies in voids in future cosmological tests of gravity. △ Less

Submitted 26 July, 2022; originally announced July 2022.

Comments: 19 pages, 10 figures

Dark matter from primordial black holes would hold charge

Authors: Ignacio J. Araya, Nelson D. Padilla, Marcelo E. Rubio, Joaquín Sureda, Juan Magaña, Loreto Osorio

Abstract: We explore the possibility that primordial black holes (PBHs) contain electric charge down to the present day. We find that PBHs should hold a non-zero net charge at their formation, due to either Poisson fluctuations at horizon crossing or high-energy particle collisions. Although initial charge configurations are subject to fast discharge processes through particle accretion or quantum particle… ▽ More We explore the possibility that primordial black holes (PBHs) contain electric charge down to the present day. We find that PBHs should hold a non-zero net charge at their formation, due to either Poisson fluctuations at horizon crossing or high-energy particle collisions. Although initial charge configurations are subject to fast discharge processes through particle accretion or quantum particle emission, we show that maximally rotating PBHs could produce magnetic fields able to shield them from discharge. Moreover, given that electrons are the lightest and fastest charge carriers, we show that the plasma within virialised dark matter haloes can endow PBHs with net negative charge. We report charge-to-mass ratios between $10^{-31}\,C/\mbox{kg}$ and $10^{-15}\,C/\mbox{kg}$ for PBHs within the mass windows that allow them to constitute the entirety of the dark matter in the Universe. △ Less

Submitted 7 February, 2023; v1 submitted 12 July, 2022; originally announced July 2022.

Comments: 25 pages, 6 figures, accepted in JCAP

Journal ref: JCAP02(2023)030

A deep learning approach to halo merger tree construction

Authors: Sandra Robles, Jonathan S. Gómez, Adín Ramírez Rivera, Nelson D. Padilla, Diego Dujovne

Abstract: A key ingredient for semi-analytic models (SAMs) of galaxy formation is the mass assembly history of haloes, encoded in a tree structure. The most commonly used method to construct halo merger histories is based on the outcomes of high-resolution, computationally intensive N-body simulations. We show that machine learning (ML) techniques, in particular Generative Adversarial Networks (GANs), are a… ▽ More A key ingredient for semi-analytic models (SAMs) of galaxy formation is the mass assembly history of haloes, encoded in a tree structure. The most commonly used method to construct halo merger histories is based on the outcomes of high-resolution, computationally intensive N-body simulations. We show that machine learning (ML) techniques, in particular Generative Adversarial Networks (GANs), are a promising new tool to tackle this problem with a modest computational cost and retaining the best features of merger trees from simulations. We train our GAN model with a limited sample of merger trees from the Evolution and Assembly of GaLaxies and their Environments (EAGLE) simulation suite, constructed using two halo finders-tree builder algorithms: SUBFIND-D-TREES and ROCKSTAR-ConsistentTrees. Our GAN model successfully learns to generate well-constructed merger tree structures with high temporal resolution, and to reproduce the statistical features of the sample of merger trees used for training, when considering up to three variables in the training process. These inputs, whose representations are also learned by our GAN model, are mass of the halo progenitors and the final descendant, progenitor type (main halo or satellite) and distance of a progenitor to that in the main branch. The inclusion of the latter two inputs greatly improves the final learned representation of the halo mass growth history, especially for SUBFIND-like ML trees. When comparing equally sized samples of ML merger trees with those of the EAGLE simulation, we find better agreement for SUBFIND-like ML trees. Finally, our GAN-based framework can be utilised to construct merger histories of low- and intermediate-mass haloes, the most abundant in cosmological simulations. △ Less

Submitted 27 June, 2022; v1 submitted 31 May, 2022; originally announced May 2022.

Comments: 17 pages, 12 figures, 3 tables, 2 appendices. Minor editorial improvements, matches published version

Report number: KCL-PH-TH/2022-33

Journal ref: MNRAS 514, 3692-3708 (2022)

Satellite galaxies in groups in the CIELO Project I. Gas removal from galaxies and its re-distribution in the intragroup medium

Authors: S. Rodríguez, D. Garcia Lambas, N. D. Padilla, P. Tissera, L. Bignone, R. Dominguez-Tenreiro, R. Gonzalez, S. Pedrosa

Abstract: We study the impact of the environment on galaxies as they fall in and orbit in the potential well of a Local Group (LG) analogue, following them with high cadence. The analysis is performed on eight disc satellite galaxies from the CIELO suite of hydrodynamical simulations. All galaxies have stellar masses within the range $[10^{8.1} - 10^{9.56}] M_{\sun} $h$^{-1}$. We measure tidal torques, ram… ▽ More We study the impact of the environment on galaxies as they fall in and orbit in the potential well of a Local Group (LG) analogue, following them with high cadence. The analysis is performed on eight disc satellite galaxies from the CIELO suite of hydrodynamical simulations. All galaxies have stellar masses within the range $[10^{8.1} - 10^{9.56}] M_{\sun} $h$^{-1}$. We measure tidal torques, ram pressure and specific star formation rates (sSFR) as a function of time, and correlate them with the amount of gas lost by satellites along their orbits. Stronger removal episodes occur when the disc plane is oriented perpendicular to the direction of motion. More than one peripassage is required to significantly modify the orientations of the discs with respect to the orbital plane. The gas removed during the interaction with the central galaxies may be also found opposite to the direction of motion, depending on the orbital configuration. Satellites are not totally quenched when the galaxies reach their first peripassage, and continue forming about $10\%$ of the final stellar mass after this event. The fraction of removed gas is found to be the product of the joint action of tidal torque and ram pressure, which can also trigger new star formation activity and subsequent supernova feedback. △ Less

Submitted 13 May, 2022; originally announced May 2022.

Comments: 17 pages, 13 figures, accepted for publication in MNRAS

Revealing the properties of void galaxies and their assembly using the EAGLE simulation

Authors: Yetli Rosas-Guevara, Patricia Tissera, Claudia del P. Lagos, Enrique Paillas, Nelson Padilla

Abstract: We explore the properties of central galaxies living in voids using the EAGLE cosmological hydrodynamic simulations. Based on the minimum void-centric distance, we define four galaxy samples: inner void, outer void, wall, and skeleton. We find that inner void galaxies with host halo masses $<10^{12}M_\odot$ have lower stellar mass and stellar mass fractions than those in denser environments, and t… ▽ More We explore the properties of central galaxies living in voids using the EAGLE cosmological hydrodynamic simulations. Based on the minimum void-centric distance, we define four galaxy samples: inner void, outer void, wall, and skeleton. We find that inner void galaxies with host halo masses $<10^{12}M_\odot$ have lower stellar mass and stellar mass fractions than those in denser environments, and the fraction of galaxies with star formation (SF) activity and atomic hydrogen (HI) gas decreases with increasing void-centric distance, in agreement with observations. To mitigate the influence of stellar (halo) mass, we compare inner void galaxies to subsamples of fixed stellar (halo) mass. Compared to denser environments, inner void galaxies with $M_{*}= 10^{[9.0-9.5]}M_\odot$ have comparable SF activity and HI gas fractions, but the lowest quenched galaxy fraction. Inner void galaxies with $M_{*}= 10^{[9.5-10.5]}M_\odot$ have the lowest HI gas fraction, the highest quenched fraction and the lowest gas metallicities. On the other hand, inner void galaxies with $M_{*}>10^{10.5}M_\odot$ have comparable SF activity and HI gas fractions to their analogues in denser environments. They retain the highest metallicity gas that might be linked to physical processes that act with lower efficiency in underdense regions, such as AGN feedback. Furthermore, inner void galaxies have the lowest fraction of positive gas-phase metallicity gradients, which are typically associated with external processes or feedback events, suggesting they have more quiet merger histories than galaxies in denser environments. Our findings shed light on how galaxies are influenced by their large-scale environment. △ Less

Submitted 8 September, 2022; v1 submitted 9 April, 2022; originally announced April 2022.

Comments: 20 pages,16 figures, revised version with a discussion section and edition in the text. Accepted to MNRAS

The Seventeenth Data Release of the Sloan Digital Sky Surveys: Complete Release of MaNGA, MaStar and APOGEE-2 Data

Authors: Abdurro'uf, Katherine Accetta, Conny Aerts, Victor Silva Aguirre, Romina Ahumada, Nikhil Ajgaonkar, N. Filiz Ak, Shadab Alam, Carlos Allende Prieto, Andres Almeida, Friedrich Anders, Scott F. Anderson, Brett H. Andrews, Borja Anguiano, Erik Aquino-Ortiz, Alfonso Aragon-Salamanca, Maria Argudo-Fernandez, Metin Ata, Marie Aubert, Vladimir Avila-Reese, Carles Badenes, Rodolfo H. Barba, Kat Barger, Jorge K. Barrera-Ballesteros, Rachael L. Beaton , et al. (316 additional authors not shown)

Abstract: This paper documents the seventeenth data release (DR17) from the Sloan Digital Sky Surveys; the fifth and final release from the fourth phase (SDSS-IV). DR17 contains the complete release of the Mapping Nearby Galaxies at Apache Point Observatory (MaNGA) survey, which reached its goal of surveying over 10,000 nearby galaxies. The complete release of the MaNGA Stellar Library (MaStar) accompanies… ▽ More This paper documents the seventeenth data release (DR17) from the Sloan Digital Sky Surveys; the fifth and final release from the fourth phase (SDSS-IV). DR17 contains the complete release of the Mapping Nearby Galaxies at Apache Point Observatory (MaNGA) survey, which reached its goal of surveying over 10,000 nearby galaxies. The complete release of the MaNGA Stellar Library (MaStar) accompanies this data, providing observations of almost 30,000 stars through the MaNGA instrument during bright time. DR17 also contains the complete release of the Apache Point Observatory Galactic Evolution Experiment 2 (APOGEE-2) survey which publicly releases infra-red spectra of over 650,000 stars. The main sample from the Extended Baryon Oscillation Spectroscopic Survey (eBOSS), as well as the sub-survey Time Domain Spectroscopic Survey (TDSS) data were fully released in DR16. New single-fiber optical spectroscopy released in DR17 is from the SPectroscipic IDentification of ERosita Survey (SPIDERS) sub-survey and the eBOSS-RM program. Along with the primary data sets, DR17 includes 25 new or updated Value Added Catalogs (VACs). This paper concludes the release of SDSS-IV survey data. SDSS continues into its fifth phase with observations already underway for the Milky Way Mapper (MWM), Local Volume Mapper (LVM) and Black Hole Mapper (BHM) surveys. △ Less

Submitted 13 January, 2022; v1 submitted 3 December, 2021; originally announced December 2021.

Comments: 40 pages, 8 figures, 6 tables. In press at ApJSS (arxiv v2 corrects some minor typos and updates references)

Making use of sub-resolution halos in N-body simulations

Authors: Joaquin Armijo, Carlton M. Baugh, Nelson D. Padilla, Peder Norberg, Christian Arnold

Abstract: Conservative mass limits are often imposed on the dark matter halo catalogues extracted from N-body simulations. By comparing simulations with different mass resolutions, at $z=0$ we find that even for halos resolved by 100 particles, the lower resolution simulation predicts a cumulative halo abundance that is 5 per cent lower than in the higher resolution simulation. We propose a simple weighting… ▽ More Conservative mass limits are often imposed on the dark matter halo catalogues extracted from N-body simulations. By comparing simulations with different mass resolutions, at $z=0$ we find that even for halos resolved by 100 particles, the lower resolution simulation predicts a cumulative halo abundance that is 5 per cent lower than in the higher resolution simulation. We propose a simple weighting scheme to utilise the halos that are usually regarded as being `sub-resolution'. With the scheme, we are able to use halos which contain only 11 particles to reproduce the clustering measured in the higher resolution simulation to within 5 per cent on scales down to $2 h^{-1}$ Mpc, thereby extending the useful halo resolution by a factor of ten below the mass at which the mass functions in the two simulations first start to deviate. The performance of the method is slightly worse at higher redshift. Our method allows a simulation to be used to probe a wider parameter space in clustering studies, for example, in a halo occupation distribution analysis. This reduces the cost of generating many simulations to estimate the covariance matrix on measurements or using a larger volume simulation to make large-scale clustering predictions. △ Less

Submitted 23 November, 2021; v1 submitted 22 November, 2021; originally announced November 2021.

Comments: 5 pages, 3 figures. Accepted for publication in MNRAS

The origin of star-gas misalignments in simulated galaxies

Authors: Catalina I. Casanueva, Claudia del P. Lagos, Nelson D. Padilla, Thomas A. Davison

Abstract: We study the origin of misalignments between the stellar and star-forming gas components of simulated galaxies in the EAGLE simulations. We focus on galaxies with stellar masses $\geq 10^9$ M$_\odot$ at 0$\leq$z$\leq$1. We compare the frequency of misalignments with observational results from the SAMI survey and find that overall, EAGLE can reproduce the incidence of misalignments in the field and… ▽ More We study the origin of misalignments between the stellar and star-forming gas components of simulated galaxies in the EAGLE simulations. We focus on galaxies with stellar masses $\geq 10^9$ M$_\odot$ at 0$\leq$z$\leq$1. We compare the frequency of misalignments with observational results from the SAMI survey and find that overall, EAGLE can reproduce the incidence of misalignments in the field and clusters, as well as the dependence on stellar mass and optical colour within the uncertainties. We study the dependence on kinematic misalignments with internal galaxy properties and different processes related to galaxy mergers and sudden changes in stellar and star-forming gas mass. We found that despite the environment being relevant in setting the conditions to misalign the star-forming gas, the internal galaxy properties play a crucial role in determining whether the gas quickly aligns with the stellar component or not. Hence, galaxies that are more triaxial and more dispersion dominated display more misalignments because they are inefficient at realigning the star-forming gas towards the stellar angular momentum vector. △ Less

Submitted 2 March, 2022; v1 submitted 1 October, 2021; originally announced October 2021.

Evolution mapping: a new approach to describe matter clustering in the non-linear regime

Authors: Ariel G. Sanchez, Andrés N. Ruiz, Jenny Gonzalez Jara, Nelson D. Padilla

Abstract: We present a new approach to describe statistics of the non-linear matter density field that exploits a degeneracy in the impact of different cosmological parameters on the linear dimensionless matter power spectrum, $Δ^2_{\rm L}(k)$. We classify all cosmological parameters into two groups, shape parameters, which determine the shape of $Δ^2_{\rm L}(k)$, and evolution parameters, which only affect… ▽ More We present a new approach to describe statistics of the non-linear matter density field that exploits a degeneracy in the impact of different cosmological parameters on the linear dimensionless matter power spectrum, $Δ^2_{\rm L}(k)$. We classify all cosmological parameters into two groups, shape parameters, which determine the shape of $Δ^2_{\rm L}(k)$, and evolution parameters, which only affect its amplitude at any given redshift. With this definition, the time evolution of $Δ^2_{\rm L}(k)$ in models with identical shape parameters but different evolution parameters can be mapped from one to the other by relabelling the redshifts that correspond to the same clustering amplitude, which we characterize by the linear mass fluctuation in spheres of radius $12\,{\rm Mpc}$, $σ_{12}(z)$. We use N-body simulations to show that the same evolution mapping relation gives a good description of the non-linear power spectrum, the halo mass function, or the full density field. The deviations from the exact degeneracy are the result of the different structure formation histories experienced by each model to reach the same clustering amplitude and can be accurately described in terms of differences in the suppression factor $g(a) = D(a)/a$. These relations can be used to drastically reduce the number of parameters required to describe the cosmology dependence of the power spectrum. We show how this can help to speed up the inference of parameter constraints from cosmological observations. We also present a new design of an emulator of the non-linear power spectrum whose predictions can be adapted to an arbitrary choice of evolution parameters and redshift. △ Less

Submitted 15 June, 2022; v1 submitted 28 August, 2021; originally announced August 2021.

Comments: 13 pages, 11 figures, replaced to match version accepted by MNRAS

Redshift-space effects in voids and their impact on cosmological tests. Part II: the void-galaxy cross-correlation function

Authors: Carlos M. Correa, Dante J. Paz, Nelson D. Padilla, Ariel G. Sánchez, Andrés N. Ruiz, Raúl E. Angulo

Abstract: This is the second part of a thorough investigation of the redshift-space effects that affect void properties and the impact they have on cosmological tests. Here, we focus on the void-galaxy cross-correlation function, specifically, on the projected versions that we developed in a previous work. The pillar of the analysis is the one-to-one relationship between real and redshift-space voids above… ▽ More This is the second part of a thorough investigation of the redshift-space effects that affect void properties and the impact they have on cosmological tests. Here, we focus on the void-galaxy cross-correlation function, specifically, on the projected versions that we developed in a previous work. The pillar of the analysis is the one-to-one relationship between real and redshift-space voids above the shot-noise level identified with a spherical void finder. Under this mapping, void properties are affected by three effects: (i) a systematic expansion as a consequence of the distortions induced by galaxy dynamics, (ii) the Alcock-Paczynski volume effect, which manifests as an overall expansion or contraction depending on the fiducial cosmology, and (iii) a systematic off-centring along the line of sight as a consequence of the distortions induced by void dynamics. We found that correlations are also affected by an additional source of distortions: the ellipticity of voids. This is the first time that distortions due to the off-centring and ellipticity effects are detected and quantified. With a simplified test, we verified that the Gaussian streaming model is still robust provided all these effects are taken into account, laying the foundations for improvements in current models in order to obtain unbiased cosmological constraints from spectroscopic surveys. Besides this practical importance, this analysis also encodes key information about the structure and dynamics of the Universe at the largest scales. Furthermore, some of the effects constitute cosmological probes by themselves, as is the case of the void ellipticity. △ Less

Submitted 19 January, 2022; v1 submitted 2 July, 2021; originally announced July 2021.

Comments: 14 pages, 9 figures, published in MNRAS, accepted after minor comments

Report number: MNRAS, Volume 509, Issue 2, pp.1871-1884, January 2022

Halo Merger Tree Comparison: Impact on Galaxy Formation Models

Authors: Jonathan S. Gómez, Nelson D. Padilla, John C. Helly, Cedric G. Lacey, Carlton M. Baugh, Claudia del P. Lagos

Abstract: We examine the effect of using different halo finders and merger tree building algorithms on galaxy properties predicted using the GALFORM semi-analytical model run on a high resolution, large volume dark matter simulation. The halo finders/tree builders HBT, ROCKSTAR, SUBFIND and VELOCIRAPTOR differ in their definitions of halo mass, on whether only spatial or phase-space information is used, and… ▽ More We examine the effect of using different halo finders and merger tree building algorithms on galaxy properties predicted using the GALFORM semi-analytical model run on a high resolution, large volume dark matter simulation. The halo finders/tree builders HBT, ROCKSTAR, SUBFIND and VELOCIRAPTOR differ in their definitions of halo mass, on whether only spatial or phase-space information is used, and in how they distinguish satellite and main haloes; all of these features have some impact on the model galaxies, even after the trees are post-processed and homogenised by GALFORM. The stellar mass function is insensitive to the halo and merger tree finder adopted. However, we find that the number of central and satellite galaxies in GALFORM does depend slightly on the halo finder/tree builder. The number of galaxies without resolved subhaloes depends strongly on the tree builder, with VELOCIRAPTOR, a phase-space finder, showing the largest population of such galaxies. The distributions of stellar masses, cold and hot gas masses, and star formation rates agree well between different halo finders/tree builders. However, because VELOCIRAPTOR has more early progenitor haloes, with these trees GALFORM produces slightly higher star formation rate densities at high redshift, smaller galaxy sizes, and larger stellar masses for the spheroid component. Since in all cases these differences are small we conclude that, when all of the trees are processed so that the main progenitor mass increases monotonically, the predicted GALFORM galaxy populations are stable and consistent for these four halo finders/tree builders. △ Less

Submitted 23 June, 2021; originally announced June 2021.

Comments: 20 pages, 16 figures

Simultaneous Estimation of Large-Scale Structure and Milky Way Dust Extinction from Galaxy Surveys

Authors: Matías Bravo, Eric Gawiser, Nelson D. Padilla, Joseph DeRose, Risa H. Wechsler, The LSST Dark Energy Science Collaboration

Abstract: The high cosmological precision offered by the next generation of galaxy surveys hinges on improved corrections for Galactic dust extinction. We explore the possibility of estimating both the dust extinction and large-scale structure from a single photometric galaxy survey, making use of the predictable manner in which Milky Way dust affects the measured brightness and colors of galaxies in a give… ▽ More The high cosmological precision offered by the next generation of galaxy surveys hinges on improved corrections for Galactic dust extinction. We explore the possibility of estimating both the dust extinction and large-scale structure from a single photometric galaxy survey, making use of the predictable manner in which Milky Way dust affects the measured brightness and colors of galaxies in a given sky location in several redshift bins. To test our method, we use a synthetic catalog from a cosmological simulation designed to model the Vera C. Rubin Observatory Legacy Survey of Space and Time. At high Galactic latitude ($|b|\gtrsim20^\circ$) and a resolution of $1^\circ$ ($7'$), we predict the uncertainty in the measurement of dust extinction, $E(B-V)$, to be $0.005\ \mathrm{mag}$ ($0.015\ \mathrm{mag}$). This is similar to the uncertainty of existing dust maps, illustrating the feasibility of our method. Simultaneous estimation of large-scale structure is predicted to recover the galaxy overdensity $δ$ with a precision of $\sim0.01$ ($\sim0.05$) at $1^\circ$ ($7'$) resolution. We also introduce a Bayesian formalism that combines prior information from existing dust maps with the likelihood of Galactic dust extinction determined from the excursion of observed galaxy properties. △ Less

Submitted 16 June, 2021; originally announced June 2021.

Comments: 22 pages, 11 figures. Submitted to ApJ

Gravitational redshifting of galaxies in the SPIDERS cluster catalogue

Authors: C. T. Mpetha, C. A. Collins, N. Clerc, A. Finoguenov, J. A. Peacock, J. Comparat, D. Schneider, R. Capasso, S. Damsted, K. Furnell, A. Merloni, N. D. Padilla, A. Saro

Abstract: Data from the SPectroscopic IDentification of ERosita Sources (SPIDERS) are searched for a detection of the gravitational redshifting of light from $\sim\!20\,000$ galaxies in $\sim\!2500$ galaxy clusters using three definitions of the cluster centre: its Brightest Cluster Galaxy (BCG), the redMaPPer identified Central Galaxy (CG), or the peak of X-ray emission. Distributions of velocity offsets b… ▽ More Data from the SPectroscopic IDentification of ERosita Sources (SPIDERS) are searched for a detection of the gravitational redshifting of light from $\sim\!20\,000$ galaxies in $\sim\!2500$ galaxy clusters using three definitions of the cluster centre: its Brightest Cluster Galaxy (BCG), the redMaPPer identified Central Galaxy (CG), or the peak of X-ray emission. Distributions of velocity offsets between galaxies and their host cluster's centre, found using observed redshifts, are created. The quantity $\hatΔ$, the average of the radial velocity difference between the cluster members and the cluster systemic velocity, reveals information on the size of a combination of effects on the observed redshift, dominated by gravitational redshifting. The change of $\hatΔ$ with radial distance is predicted for SPIDERS galaxies in General Relativity (GR), and $f(R)$ gravity, and compared to the observations. The values of $\hatΔ=-13.5\pm4.7$ km s$^{-1}$, $\hatΔ=-12.5\pm5.1$ km s$^{-1}$, and $\hatΔ=-18.6\pm4.8$ km s$^{-1}$ for the BCG, X-ray and CG cases respectively broadly agree with the literature. There is no significant preference of one gravity theory over another, but all cases give a clear detection ($>2.5σ$) of $\hatΔ$. The BCG centroid is deemed to be the most robust method in this analysis, due to no well defined central redshift when using an X-ray centroid, and CGs identified by redMaPPer with no associated spectroscopic redshift. For future gravitational redshift studies, an order of magnitude more galaxies, $\sim\!500\,000$, will be required-a possible feat with the forthcoming Vera C. Rubin Observatory, Euclid and eROSITA. △ Less

Submitted 23 February, 2021; v1 submitted 22 February, 2021; originally announced February 2021.

Comments: Accepted for publication in MNRAS. 11 pages, 10 figures

Redshift-space distortions with split densities

Authors: Enrique Paillas, Yan-Chuan Cai, Nelson Padilla, Ariel Sánchez

Abstract: Accurate modelling of redshift-space distortions (RSD) is challenging in the non-linear regime for two-point statistics e.g. the two-point correlation function (2PCF). We take a different perspective to split the galaxy density field according to the local density, and cross-correlate those densities with the entire galaxy field. Using mock galaxies, we demonstrate that combining a series of cross… ▽ More Accurate modelling of redshift-space distortions (RSD) is challenging in the non-linear regime for two-point statistics e.g. the two-point correlation function (2PCF). We take a different perspective to split the galaxy density field according to the local density, and cross-correlate those densities with the entire galaxy field. Using mock galaxies, we demonstrate that combining a series of cross-correlation functions (CCFs) offers improvements over the 2PCF as follows: 1. The distribution of peculiar velocities in each split density is nearly Gaussian. This allows the Gaussian streaming model for RSD to perform accurately within the statistical errors of a ($1.5\,h^{-1}$Gpc)$^3$ volume for almost all scales and all split densities. 2. The PDF of the density field at small scales is non-Gaussian, but the CCFs of split densities capture the non-Gaussianity, leading to improved cosmological constraints over the 2PCF. We can obtain unbiased constraints on the growth parameter $fσ_{12}$ at the per-cent level, and Alcock-Paczynski (AP) parameters at the sub-per-cent level with the minimal scale of $15\,h^{-1}{\rm Mpc}$. This is a $\sim$30 per cent and $\sim$6 times improvement over the 2PCF, respectively. The diverse and steep slopes of the CCFs at small scales are likely to be responsible for the improved constraints of AP parameters. 3. Baryon acoustic oscillations (BAO) are contained in all CCFs of split densities. Including BAO scales helps to break the degeneracy between the line-of-sight and transverse AP parameters, allowing independent constraints on them. We discuss and compare models for RSD around spherical densities. △ Less

Submitted 2 June, 2021; v1 submitted 24 January, 2021; originally announced January 2021.

Comments: 23 pages, 16 figures, MNRAS accepted version after peer review, minor comments

SPIDERS: An Overview of The Largest Catalogue of Spectroscopically Confirmed X-ray Galaxy Clusters

Authors: C. C. Kirkpatrick, N. Clerc, A. Finoguenov, S. Damsted, J. Ider Chitham, A. E. Kukkola, A. Gueguen, K. Furnell, E. Rykoff, J. Comparat, A. Saro, R. Capasso, N. Padilla, G. Erfanianfar, G. A. Mamon, C. Collins, A. Merloni, J. R. Brownstein, D. P. Schneider

Abstract: SPIDERS is the spectroscopic follow-up effort of the Sloan Digital Sky Survey IV (SDSS-IV) project for the identification of X-ray selected galaxy clusters. We present our catalogue of 2740 visually inspected galaxy clusters as a part of the SDSS Data Release 16 (DR16). Here we detail the target selection, our methods for validation of the candidate clusters, performance of the survey, the constru… ▽ More SPIDERS is the spectroscopic follow-up effort of the Sloan Digital Sky Survey IV (SDSS-IV) project for the identification of X-ray selected galaxy clusters. We present our catalogue of 2740 visually inspected galaxy clusters as a part of the SDSS Data Release 16 (DR16). Here we detail the target selection, our methods for validation of the candidate clusters, performance of the survey, the construction of the final sample, and a full description of what is found in the catalogue. Of the sample, the median number of members per cluster is approximately 10, with 818 having 15 or greater. We find that we are capable of validating over 99% of clusters when 5 redshifts are obtained below $z<0.3$ and when 9 redshifts are obtained above $z>0.3$. We discuss the improvements of this catalogue's identification of cluster using 33,340 redshifts, with $Δz_{\rm{phot}} / Δz_{\rm{spec}} \sim 100$, over other photometric and spectroscopic surveys, as well as present an update to previous ($σ- L_{X}$) and ($σ- λ$) relations. Finally, we present our cosmological constraints derived using the velocity dispersion function. △ Less

Submitted 12 January, 2021; originally announced January 2021.

Comments: 16 pages, 14 figures, Accepted for publication in MNRAS

Magnetic field generation from primordial black hole distributions

Authors: Ignacio J. Araya, Marcelo E. Rubio, Marco San Martin, Federico A. Stasyszyn, Nelson D. Padilla, Juan Magana, Joaquin Sureda

Abstract: We introduce a statistical method for estimating magnetic field fluctuations generated from primordial black hole (PBH) populations. To that end, we consider monochromatic and extended Press-Schechter PBH mass functions, such that each constituent is capable of producing its own magnetic field due to some given physical mechanism. Assuming linear correlation between magnetic field fluctuations and… ▽ More We introduce a statistical method for estimating magnetic field fluctuations generated from primordial black hole (PBH) populations. To that end, we consider monochromatic and extended Press-Schechter PBH mass functions, such that each constituent is capable of producing its own magnetic field due to some given physical mechanism. Assuming linear correlation between magnetic field fluctuations and matter over-densities, our estimates depend on the mass function, the physical field generation mechanism by each PBH constituent, and the characteristic PBH separation. After computing the power spectrum of magnetic field fluctuations, we apply our formalism to study the plausibility that two particular field generation mechanisms could have given rise to the expected seed fields according to current observational constraints. The first mechanism is the Biermann battery and the second one is due to the accretion of magnetic monopoles at PBH formation, constituting magnetic PBHs. Our results show that, for monochromatic distributions, it does not seem to be possible to generate sufficiently intense seed fields in any of the two field generation mechanisms. For extended distributions, it is also not possible to generate the required seed field by only assuming a Biermann battery mechanism. In fact, we report an average seed field by this mechanism of about 10^{-47} G, at z = 20. For the case of magnetic monopoles we instead assume that the seed values from the literature are achieved and calculate the necessary number density of monopoles. In this case we obtain values that are below the upper limits from current constraints. △ Less

Submitted 8 March, 2021; v1 submitted 17 December, 2020; originally announced December 2020.

Comments: 15 pages, 5 figures, accepted for publication in MNRAS

Too dense to go through: The importance of low-mass clusters for satellite quenching

Authors: Diego Pallero, Facundo A. Gómez, Nelson D. Padilla, Yannick M. Bahé, Cristian A. Vega-Martínez, S. Torres-Flores

Abstract: We study the evolution of satellite galaxies in clusters of the \textsc{c-eagle} simulations, a suite of 30 high-resolution cosmological hydrodynamical zoom-in simulations based on the \textsc{eagle} code. We find that the majority of galaxies that are quenched at $z=0$ ($\gtrsim$ 80$\%$) reached this state in a dense environment (log$_{10}$M$_{200}$[M$_{\odot}$]$\geq$13.5). At low redshift, regar… ▽ More We study the evolution of satellite galaxies in clusters of the \textsc{c-eagle} simulations, a suite of 30 high-resolution cosmological hydrodynamical zoom-in simulations based on the \textsc{eagle} code. We find that the majority of galaxies that are quenched at $z=0$ ($\gtrsim$ 80$\%$) reached this state in a dense environment (log$_{10}$M$_{200}$[M$_{\odot}$]$\geq$13.5). At low redshift, regardless of the final cluster mass, galaxies appear to reach their quenching state in low-mass clusters. Moreover, galaxies quenched inside the cluster that they reside in at $z=0$ are the dominant population in low-mass clusters, while galaxies quenched in a different halo dominate in the most massive clusters. When looking at clusters at $z>0.5$, their in situ quenched population dominates at all cluster masses. This suggests that galaxies are quenched inside the first cluster they fall into. After galaxies cross the cluster's $r_{200}$ they rapidly become quenched ($\lesssim$ 1Gyr). Just a small fraction of galaxies ($\lesssim 15\%$) is capable of retaining their gas for a longer period of time, but after 4Gyr, almost all galaxies are quenched. This phenomenon is related to ram pressure stripping and is produced when the density of the intracluster medium reaches a threshold of $ρ_{\rm ICM}$ $\sim 3 \times 10 ^{-5}$ n$_{\rm H}$ (cm$^{-3}$). These results suggest that galaxies start a rapid-quenching phase shortly after their first infall inside $r_{200}$ and that, by the time they reach $r_{500}$, most of them are already quenched. △ Less

Submitted 3 October, 2023; v1 submitted 15 December, 2020; originally announced December 2020.

Comments: 14 pages, 8 figures, 2 appendices

Journal ref: Monthly Notices of the Royal Astronomical Society, Volume 511, Issue 3, pp.3210-3227, April 2022

Clustering of CODEX clusters

Authors: Valtteri Lindholm, Alexis Finoguenov, Johan Comparat, Charles C. Kirkpatrick, Eli Rykoff, Nicolas Clerc, Chris Collins, Sanna Damsted, Jacob Ider Chitham, Nelson Padilla

Abstract: Aims. We analyze the autocorrelation function of a large contiguous sample of galaxy clusters, the Constrain Dark Energy with X-ray (CODEX) sample, in which we take particular care of cluster definition. These clusters were X-ray selected using the RASS survey and then identified as galaxy clusters using the code redMaPPer run on the photometry of the SDSS. We develop methods for precisely account… ▽ More Aims. We analyze the autocorrelation function of a large contiguous sample of galaxy clusters, the Constrain Dark Energy with X-ray (CODEX) sample, in which we take particular care of cluster definition. These clusters were X-ray selected using the RASS survey and then identified as galaxy clusters using the code redMaPPer run on the photometry of the SDSS. We develop methods for precisely accounting for the sample selection effects on the clustering and demonstrate their robustness using numerical simulations. Methods. Using the clean CODEX sample, which was obtained by applying a redshift-dependent richness selection, we computed the two-point autocorrelation function of galaxy clusters in the $0.1<z<0.3$ and $0.3<z<0.5$ redshift bins. We compared the bias in the measured correlation function with values obtained in numerical simulations using a similar cluster mass range. Results. By fitting a power law, we measured a correlation length $r_0=18.7 \pm 1.1$ and slope $γ=1.98 \pm 0.14$ for the correlation function in the full redshift range. By fixing the other cosmological parameters to their WMAP9 values, we reproduced the observed shape of the correlation function under the following cosmological conditions: $Ω_{m_0}=0.22^{+0.04}_{-0.03}$ and $S_8=σ_8 (Ω_{m_0} /0.3)^{0.5}=0.85^{+0.10}_{-0.08}$ with estimated additional systematic errors of $σ_{Ω_{m_0}} = 0.02$ and $σ_{S_8} = 0.20$. We illustrate the complementarity of clustering constraints by combining them with CODEX cosmological constraints based on the X-ray luminosity function, deriving $Ω_{m_0} = 0.25 \pm 0.01$ and $σ_8 = 0.81^{+0.01}_{-0.02}$ with an estimated additional systematic error of $σ_{Ω_{m_0}} = 0.07$ and $σ_{σ_8} = 0.04$. The mass calibration and statistical quality of the mass tracers are the dominant source of uncertainty. △ Less

Submitted 30 November, 2020; originally announced December 2020.

Comments: 13 pages, 15 figures

Journal ref: A&A 646, A8 (2021)

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

Power spectrum of density fluctuations, halo abundances and clustering with primordial black holes

Authors: Nelson Padilla, Juan Magana, Joaquin Sureda, Ignacio Araya

Abstract: We study the effect of dark matter (DM) being encapsulated in primordial black holes (PBHs) on the power spectrum of density fluctuations $P(k)$; we also look at its effect on the abundance of haloes and their clustering. We allow the growth of Poisson fluctuations since matter and radiation equality and study both monochromatic and extended PBH mass distributions. We present updated monochromatic… ▽ More We study the effect of dark matter (DM) being encapsulated in primordial black holes (PBHs) on the power spectrum of density fluctuations $P(k)$; we also look at its effect on the abundance of haloes and their clustering. We allow the growth of Poisson fluctuations since matter and radiation equality and study both monochromatic and extended PBH mass distributions. We present updated monochromatic black hole mass constraints by demanding $<10\%$ deviations from the $Λ$ cold dark matter power spectrum at a scale of $k=1$hMpc$^{-1}$. Our results show that PBHs with masses $>10^4$h$^{-1}M_\odot$ are excluded from conforming all of the DM in the Universe. We also apply this condition to our extended Press-Schechter (PS) mass functions, and find that the Poisson power is scale dependent even before applying evolution. We find that characteristic masses $M^*\leq10^2 $h$^{-1}M_\odot$ are allowed, {leaving only two characteristic PBH mass windows of PS mass functions when combining with previous constraints, at $M^*\sim10^2$h$^{-1}M_\odot$ and $\sim10^{-8}$h$^{-1}M_\odot$ where all of the DM can be in PBHs. The resulting DM halo mass functions within these windows are similar} to those resulting from cold dark matter made of fundamental particles. However, as soon as the parameters produce unrealistic $P(k)$, the resulting halo mass functions and their bias as a function of halo mass deviate strongly from the behaviour measured in the real Universe. △ Less

Submitted 19 April, 2021; v1 submitted 13 October, 2020; originally announced October 2020.

Comments: 19 pages, 9 figures, version accepted for publication in MNRAS

Press-Schechter primordial black hole mass functions and their observational constraints

Authors: Joaquin Sureda, Juan Magana, Ignacio J. Araya, Nelson D. Padilla

Abstract: We present a modification of the Press-Schechter (PS) formalism to derive general mass functions for primordial black holes (PBHs), considering their formation as being associated to the amplitude of linear energy density fluctuations. To accommodate a wide range of physical relations between the linear and non-linear conditions for collapse, we introduce an additional parameter to the PS mechanis… ▽ More We present a modification of the Press-Schechter (PS) formalism to derive general mass functions for primordial black holes (PBHs), considering their formation as being associated to the amplitude of linear energy density fluctuations. To accommodate a wide range of physical relations between the linear and non-linear conditions for collapse, we introduce an additional parameter to the PS mechanism, and that the collapse occurs at either a given cosmic time, or as fluctuations enter the horizon. We study the case where fluctuations obey Gaussian statistics and follow a primordial power spectrum of broken power-law form with a blue spectral index for small scales. We use the observed abundance of super-massive black holes (SMBH) to constrain the extended mass functions taking into account dynamical friction. We further constrain the modified PS by developing a method for converting existing constraints on the PBH mass fraction, derived assuming monochromatic mass distributions for PBHs, into constraints applicable for extended PBH mass functions. We find that when considering well established monochromatic constraints there are regions in parameter space where all the dark matter can be made of PBHs. Of special interest is the region for the characteristic mass of the distribution ~10^2 M_Sun, for a wide range of blue spectral indices in the scenario where PBHs form as they enter the horizon, where the linear threshold for collapse is of the order of the typical overdensities, as this is close to the black hole masses detected by LIGO which are difficult to explain by stellar collapse. △ Less

Submitted 24 August, 2021; v1 submitted 21 August, 2020; originally announced August 2020.

Comments: 25 pages, 9 figures, version accepted for publication in MNRAS

Redshift-space effects in voids and their impact on cosmological tests. Part I: the void size function

Authors: Carlos M. Correa, Dante J. Paz, Ariel G. Sánchez, Andrés N. Ruiz, Nelson D. Padilla, Raúl E. Angulo

Abstract: Voids are promising cosmological probes. Nevertheless, every cosmological test based on voids must necessarily employ methods to identify them in redshift space. Therefore, redshift-space distortions (RSD) and the Alcock-Paczynski effect (AP) have an impact on the void identification process itself generating distortion patterns in observations. Using a spherical void finder, we developed a statis… ▽ More Voids are promising cosmological probes. Nevertheless, every cosmological test based on voids must necessarily employ methods to identify them in redshift space. Therefore, redshift-space distortions (RSD) and the Alcock-Paczynski effect (AP) have an impact on the void identification process itself generating distortion patterns in observations. Using a spherical void finder, we developed a statistical and theoretical framework to describe physically the connection between the identification in real and redshift space. We found that redshift-space voids above the shot noise level have a unique real-space counterpart spanning the same region of space, they are systematically bigger and their centres are preferentially shifted along the line of sight. The expansion effect is a by-product of RSD induced by tracer dynamics at scales around the void radius, whereas the off-centring effect constitutes a different class of RSD induced at larger scales by the global dynamics of the whole region containing the void. The volume of voids is also altered by the fiducial cosmology assumed to measure distances, this is the AP change of volume. These three systematics have an impact on cosmological statistics. In this work, we focus on the void size function. We developed a theoretical framework to model these effects and tested it with a numerical simulation, recovering the statistical properties of the abundance of voids in real space. This description depends strongly on cosmology. Hence, we lay the foundations for improvements in current models of the abundance of voids in order to obtain unbiased cosmological constraints from redshift surveys. △ Less

Submitted 15 October, 2020; v1 submitted 23 July, 2020; originally announced July 2020.

Comments: 16 pages, 12 figures, accepted in MNRAS

The Completed SDSS-IV extended Baryon Oscillation Spectroscopic Survey: exploring the Halo Occupation Distribution model for Emission Line Galaxies

Authors: Santiago Avila, Violeta Gonzalez-Perez, Faizan G. Mohammad, Arnaud de Mattia, Cheng Zhao, Anand Raichoor, Amelie Tamone, Shadab Alam, Julian Bautista, Davide Bianchi, Etienne Burtin, Michael J. Chapman, Chia-Hsun Chuang, Johan Comparat, Kyle Dawson, Thomas Divers, Helion du Mas des Bourboux, Hector Gil-Marin, Eva-Maria Mueller, Salman Habib, Katrin Heitmann, Vanina Ruhlmann-Kleider, Nelson Padilla, Will J. Percival, Ashley J. Ross , et al. (3 additional authors not shown)

Abstract: We study the modelling of the Halo Occupation Distribution (HOD) for the eBOSS DR16 Emission Line Galaxies (ELGs). Motivated by previous theoretical and observational studies, we consider different physical effects that can change how ELGs populate haloes. We explore the shape of the average HOD, the fraction of satellite galaxies, their probability distribution function (PDF), and their density a… ▽ More We study the modelling of the Halo Occupation Distribution (HOD) for the eBOSS DR16 Emission Line Galaxies (ELGs). Motivated by previous theoretical and observational studies, we consider different physical effects that can change how ELGs populate haloes. We explore the shape of the average HOD, the fraction of satellite galaxies, their probability distribution function (PDF), and their density and velocity profiles. Our baseline HOD shape was fitted to a semi-analytical model of galaxy formation and evolution, with a decaying occupation of central ELGs at high halo masses. We consider Poisson and sub/super-Poissonian PDFs for satellite assignment. We model both NFW and particle profiles for satellite positions, also allowing for decreased concentrations. We model velocities with the virial theorem and particle velocity distributions. Additionally, we introduce a velocity bias and a net infall velocity. We study how these choices impact the clustering statistics while keeping the number density and bias fixed to that from eBOSS ELGs. The projected correlation function, $w_p$, captures most of the effects from the PDF and satellites profile. The quadrupole, $ξ_2$, captures most of the effects coming from the velocity profile. We find that the impact of the mean HOD shape is subdominant relative to the rest of choices. We fit the clustering of the eBOSS DR16 ELG data under different combinations of the above assumptions. The catalogues presented here have been analysed in companion papers, showing that eBOSS RSD+BAO measurements are insensitive to the details of galaxy physics considered here. These catalogues are made publicly available. △ Less

Submitted 27 November, 2020; v1 submitted 17 July, 2020; originally announced July 2020.

Comments: Data available here: http://popia.ft.uam.es/eBOSS_ELG_OR_mocks. A description of eBOSS and links to all associated publications can be found here: https://www.sdss.org/surveys/eboss/ ; 24 pages, 17 Figures; Published in MNRAS 25 Sep 2020

Journal ref: Monthly Notices of the Royal Astronomical Society, Volume 499, Issue 4, pp.5486-5507 (2020)

SPIDERS: overview of the X-ray galaxy cluster follow-up and the final spectroscopic data release

Authors: N. Clerc, C. C. Kirkpatrick, A. Finoguenov, R. Capasso, J. Comparat, S. Damsted, K. Furnell, A. E. Kukkola, J. Ider Chitham, A. Merloni, M. Salvato, A. Gueguen, T. Dwelly, C. Collins, A. Saro, G. Erfanianfar, D. P. Schneider, J. Brownstein, G. A. Mamon, N. Padilla, E. Jullo, D. Bizyaev

Abstract: SPIDERS (The SPectroscopic IDentification of eROSITA Sources) is a large spectroscopic programme for X-ray selected galaxy clusters as part of the Sloan Digital Sky Survey-IV (SDSS-IV). We describe the final dataset in the context of SDSS Data Release 16 (DR16): the survey overall characteristics, final targeting strategies, achieved completeness and spectral quality, with special emphasis on its… ▽ More SPIDERS (The SPectroscopic IDentification of eROSITA Sources) is a large spectroscopic programme for X-ray selected galaxy clusters as part of the Sloan Digital Sky Survey-IV (SDSS-IV). We describe the final dataset in the context of SDSS Data Release 16 (DR16): the survey overall characteristics, final targeting strategies, achieved completeness and spectral quality, with special emphasis on its use as a galaxy cluster sample for cosmology applications. SPIDERS now consists of about 27,000 new optical spectra of galaxies selected within 4,000 photometric red sequences, each associated with an X-ray source. The excellent spectrograph efficiency and a robust analysis pipeline yield a spectroscopic redshift measurement success rate exceeding 98%, with a median velocity accuracy of 20 km s$^{-1}$ (at $z=0.2$). Using the catalogue of 2,740 X-ray galaxy clusters confirmed with DR16 spectroscopy, we reveal the three-dimensional map of the galaxy cluster distribution in the observable Universe up to $z\sim0.6$. We highlight the homogeneity of the member galaxy spectra among distinct regions of the galaxy cluster phase space. Aided by accurate spectroscopic redshifts and by a model of the sample selection effects, we compute the galaxy cluster X-ray luminosity function and we present its lack of evolution up to $z=0.6$. Finally we discuss the prospects of forthcoming large multiplexed spectroscopic programmes dedicated to follow up the next generation of all-sky X-ray source catalogues. △ Less

Submitted 10 July, 2020; originally announced July 2020.

Comments: 19 pages, 21 figures. Accepted for publication in MNRAS

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.

The manifestation of secondary halo bias on the galaxy population from IllustrisTNG300

Authors: Antonio D. Montero-Dorta, M. Celeste Artale, L. Raul Abramo, Beatriz Tucci, Nelson Padilla, Gabriela Sato-Polito, Iván Lacerna, Facundo Rodriguez, Raul E. Angulo

Abstract: We use the improved IllustrisTNG300 magneto-hydrodynamical cosmological simulation to revisit the effect that secondary halo bias has on the clustering of the central galaxy population. With a side length of 205 $h^{-1}$Mpc and significant improvements on the sub-grid model with respect to the previous Illustris boxes, IllustrisTNG300 allows us to explore the dependencies of galaxy clustering over… ▽ More We use the improved IllustrisTNG300 magneto-hydrodynamical cosmological simulation to revisit the effect that secondary halo bias has on the clustering of the central galaxy population. With a side length of 205 $h^{-1}$Mpc and significant improvements on the sub-grid model with respect to the previous Illustris boxes, IllustrisTNG300 allows us to explore the dependencies of galaxy clustering over a large cosmological volume and wide halo-mass range. We show, at high statistical significance, that the halo assembly bias signal (i.e., the secondary dependence of halo bias on halo formation redshift) manifests itself on the clustering of the central galaxy population when this is split by stellar mass, colour, specific star formation rate, and surface density. A significant detection is also obtained for galaxy size: at fixed halo mass, larger central galaxies are more tightly clustered than smaller central galaxies in haloes of mass M$_{\rm vir} \lesssim 10^{12.5}$ $h^{-1}$M$_{\odot}$. This effect, however, seems to be uncorrelated with halo formation time, unlike the rest of the secondary dependencies analysed. We also explore the transmission of the halo spin bias signal, i.e., the secondary dependence of halo bias on halo spin. Although galaxy spin retains little information about the total spin of the halo, the correlation is enough to produce a significant galaxy spin bias signal. We discuss possible ways to probe the spin bias effects with observations. △ Less

Submitted 6 January, 2020; originally announced January 2020.

Comments: 16 pages, 10 figures, submitted to MNRAS

The Sixteenth Data Release of the Sloan Digital Sky Surveys: First Release from the APOGEE-2 Southern Survey and Full Release of eBOSS Spectra

Authors: Romina Ahumada, Carlos Allende Prieto, Andres Almeida, Friedrich Anders, Scott F. Anderson, Brett H. Andrews, Borja Anguiano, Riccardo Arcodia, Eric Armengaud, Marie Aubert, Santiago Avila, Vladimir Avila-Reese, Carles Badenes, Christophe Balland, Kat Barger, Jorge K. Barrera-Ballesteros, Sarbani Basu, Julian Bautista, Rachael L. Beaton, Timothy C. Beers, B. Izamar T. Benavides, Chad F. Bender, Mariangela Bernardi, Matthew Bershady, Florian Beutler , et al. (289 additional authors not shown)

Abstract: This paper documents the sixteenth data release (DR16) from the Sloan Digital Sky Surveys; the fourth and penultimate from the fourth phase (SDSS-IV). This is the first release of data from the southern hemisphere survey of the Apache Point Observatory Galactic Evolution Experiment 2 (APOGEE-2); new data from APOGEE-2 North are also included. DR16 is also notable as the final data release for the… ▽ More This paper documents the sixteenth data release (DR16) from the Sloan Digital Sky Surveys; the fourth and penultimate from the fourth phase (SDSS-IV). This is the first release of data from the southern hemisphere survey of the Apache Point Observatory Galactic Evolution Experiment 2 (APOGEE-2); new data from APOGEE-2 North are also included. DR16 is also notable as the final data release for the main cosmological program of the Extended Baryon Oscillation Spectroscopic Survey (eBOSS), and all raw and reduced spectra from that project are released here. DR16 also includes all the data from the Time Domain Spectroscopic Survey (TDSS) and new data from the SPectroscopic IDentification of ERosita Survey (SPIDERS) programs, both of which were co-observed on eBOSS plates. DR16 has no new data from the Mapping Nearby Galaxies at Apache Point Observatory (MaNGA) survey (or the MaNGA Stellar Library "MaStar"). We also preview future SDSS-V operations (due to start in 2020), and summarize plans for the final SDSS-IV data release (DR17). △ Less

Submitted 11 May, 2020; v1 submitted 5 December, 2019; originally announced December 2019.

Comments: DR16 release: Monday Dec 9th 2019. This is the alphabetical order SDSS-IV collaboration data release paper. 25 pages, 6 figures, accepted by ApJS on 11th May 2020. Minor changes clarify or improve text and figures relative to v1

Following the crumbs: Statistical effects of Ram Pressure in Galaxies

Authors: S. Rodríguez, D. Garcia Lambas, N. D. Padilla, P. Troncoso-Iribarren

Abstract: We analyse the presence of dust around galaxy group members through the reddening of background quasars. By taking into account quasar colour and their dependence on redshift and angular position, we derive mean quasar colours excess in projected regions around member galaxies and infer the associated dust mass. For disc-like galaxies perpendicular to the plane of the sky, and at group-centric dis… ▽ More We analyse the presence of dust around galaxy group members through the reddening of background quasars. By taking into account quasar colour and their dependence on redshift and angular position, we derive mean quasar colours excess in projected regions around member galaxies and infer the associated dust mass. For disc-like galaxies perpendicular to the plane of the sky, and at group-centric distances of the order of the virial radius, thus likely to reside in the infall regions of groups, we find systematic colour excess values $e \sim 0.009 \pm 0.004$ for $g-r$ colour. Under the hypothesis of Milky Way dust properties we derive dust masses of $5.8 \pm 2.5 \cdot 10^8 M_{\odot}/h$, implying that a large fraction of dust is being stripped from galaxies in their path to groups. We also studied the photometry of member galaxies to derive a colour asymmetry relative to the group centre direction from a given galaxy. We conclude that the regions of galaxies facing the centre are bluer, consistent with the effects of gas compression and star-formation. We also combine these two procedures finding that galaxies with a small colour asymmetry show the largest amounts of dust towards the external regions compared to a control sample. We conclude that dust removal is very efficient in galaxies on infall. The fact that galaxies redder towards groups centres are associated to the strongest reddening of background quasars suggest that gas removal induced by ram pressure stripping plays a key role in galaxy evolution and dust content. △ Less

Submitted 6 November, 2019; originally announced November 2019.

Comments: 7 pages, 7 figures, submitted to MNRAS

Massive Low Surface Brightness Galaxies in the EAGLE Simulation

Authors: Andrea Kulier, Gaspar Galaz, Nelson D. Padilla, James W. Trayford

Abstract: We investigate the formation and properties of low surface brightness galaxies (LSBGs) with $M_{*} > 10^{9.5} \mathrm{M_{\odot}}$ in the EAGLE hydrodynamical cosmological simulation. Galaxy surface brightness depends on a combination of stellar mass surface density and mass-to-light ratio ($M/L$), such that low surface brightness is strongly correlated with both galaxy angular momentum (low surfac… ▽ More We investigate the formation and properties of low surface brightness galaxies (LSBGs) with $M_{*} > 10^{9.5} \mathrm{M_{\odot}}$ in the EAGLE hydrodynamical cosmological simulation. Galaxy surface brightness depends on a combination of stellar mass surface density and mass-to-light ratio ($M/L$), such that low surface brightness is strongly correlated with both galaxy angular momentum (low surface density) and low specific star formation rate (high $M/L$). This drives most of the other observed correlations between surface brightness and galaxy properties, such as the fact that most LSBGs have low metallicity. We find that LSBGs are more isolated than high surface brightness galaxies (HSBGs), in agreement with observations, but that this trend is driven entirely by the fact that LSBGs are unlikely to be close-in satellites. The majority of LSBGs are consistent with a formation scenario in which the galaxies with the highest angular momentum are those that formed most of their stars recently from a gas reservoir co-rotating with a high-spin dark matter halo. However, the most extended LSBG disks in EAGLE, which are comparable in size to observed giant LSBGs, are built up via mergers. These galaxies are found to inhabit dark matter halos with a higher spin in their inner regions ($<0.1r_{200c}$), even when excluding the effects of baryonic physics by considering matching halos from a dark matter only simulation with identical initial conditions. △ Less

Submitted 2 September, 2020; v1 submitted 11 October, 2019; originally announced October 2019.

Comments: Version accepted to MNRAS. 20 pages, 13 figures, plus appendix with 2 figures

Journal ref: Monthly Notices of the Royal Astronomical Society, Volume 496, Issue 3, August 2020, Pages 3996-4016

[OII] emitters in MultiDark-Galaxies and DEEP2

Authors: G. Favole, V. Gonzalez-Perez, A. Orsi, D. Stoppacher, J. Comparat, S. A. Cora, C. A. Vega-Martinez, A. R. H. Stevens, C. Maraston, D. Croton, A. Knebe, A. J. Benson, A. D. Montero-Dorta, N. Padilla, F. Prada, D. Thomas

Abstract: We use three semi-analytic models (SAMs) of galaxy formation and evolution, run on the same 1$h^{-1}$Gpc MultiDark Planck2 cosmological simulation, to investigate the properties of [OII] emission line galaxies in the redshift range $0.6<z<1.2$. We compare model predictions with different observational data sets, including DEEP2--Firefly galaxies with absolute magnitudes. We estimate the [OII] lumi… ▽ More We use three semi-analytic models (SAMs) of galaxy formation and evolution, run on the same 1$h^{-1}$Gpc MultiDark Planck2 cosmological simulation, to investigate the properties of [OII] emission line galaxies in the redshift range $0.6<z<1.2$. We compare model predictions with different observational data sets, including DEEP2--Firefly galaxies with absolute magnitudes. We estimate the [OII] luminosity, L[OII], using simple relations derived both from the models and observations and also using a public code. This code ideally uses as input instantaneous star formation rates (SFRs), which are only provided by one of the SAMs under consideration. We use this SAM to study the feasibility of inferring galaxies' L[OII] for models that only provide average SFRs. We find that the post-processing computation of L[OII] from average SFRs is accurate for model galaxies with dust attenuated L[OII]$\lesssim10^{42.2}$erg s$^{-1}$ ($<5\%$ discrepancy). We also explore how to derive the [OII] luminosity from simple relations using global properties usually output by SAMs. Besides the SFR, the model L[OII] is best correlated with the observed-frame $u$ and $g$ broad-band magnitudes. These correlations have coefficients (r-values) above 0.64 and a dispersion that varies with L[OII]. We use these correlations and an observational one based on SFR and metallicity to derive L[OII]. These relations result in [OII] luminosity functions and halo occupation distributions with shapes that vary depending on both the model and the method used. Nevertheless, for all the considered models, the amplitude of the clustering at scales above 1$h^{-1}$Mpc remains unchanged independently of the method used to derive L[OII]. △ Less

Submitted 12 February, 2021; v1 submitted 15 August, 2019; originally announced August 2019.

Comments: 23 pages, 19 figures, published in MNRAS. Data available at: http://popia.ft.uam.es/MultiDarkEmissionLines/. Dust attenuation code: https://github.com/gfavole/dust

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

A Halo Merger Tree Generation and Evaluation Framework

Authors: Sandra Robles, Jonathan S. Gómez, Adín Ramírez Rivera, Jenny A. González, Nelson D. Padilla, Diego Dujovne

Abstract: Semi-analytic models are best suited to compare galaxy formation and evolution theories with observations. These models rely heavily on halo merger trees, and their realistic features (i.e., no drastic changes on halo mass or jumps on physical locations). Our aim is to provide a new framework for halo merger tree generation that takes advantage of the results of large volume simulations, with a mo… ▽ More Semi-analytic models are best suited to compare galaxy formation and evolution theories with observations. These models rely heavily on halo merger trees, and their realistic features (i.e., no drastic changes on halo mass or jumps on physical locations). Our aim is to provide a new framework for halo merger tree generation that takes advantage of the results of large volume simulations, with a modest computational cost. We treat halo merger tree construction as a matrix generation problem, and propose a Generative Adversarial Network that learns to generate realistic halo merger trees. We evaluate our proposal on merger trees from the EAGLE simulation suite, and show the quality of the generated trees. △ Less

Submitted 21 June, 2019; originally announced June 2019.

Comments: 11 pages, 7 figures, 2 tables, 3 appendices. Presented at the ICML 2019 Workshop on Theoretical Physics for Deep Learning