Search

The Dark Energy Spectroscopic Instrument (DESI) survey is mapping the large-scale distribution of millions of Emission Line Galaxies (ELGs) over vast cosmic volumes to measure the growth history of the Universe. However, compared to Luminous Red Galaxies (LRGs), very little is known about the connection of ELGs with the underlying matter field. In this paper, we employ a novel theoretical model, SHAMe-SF, to infer the connection between ELGs and their host dark matter subhaloes. SHAMe-SF is a version of subhalo abundance matching that incorporates prescriptions for multiple processes, including star formation, tidal stripping, environmental correlations, and quenching. We analyse the public measurements of the projected and redshift-space ELGs correlation functions at $z=1.0$ and $z=1.3$ from DESI One Percent data release, which we fit over a broad range of scales $r \in [0.1, 30]/h^{-1}$Mpc to within the statistical uncertainties of the data. We also validate the inference pipeline using two mock DESI ELG catalogues built from hydrodynamical (TNG300) and semi-analytical galaxy formation models (\texttt{L-Galaxies}). SHAMe-SF is able to reproduce the clustering of DESI-ELGs and the mock DESI samples within statistical uncertainties. We infer that DESI ELGs typically reside in haloes of $\sim 10^{11.8}h^{-1}$M$_{\odot}$ when they are central, and $\sim 10^{12.5}h^{-1}$M$_{\odot}$ when they are a satellite, which occurs in $\sim$30 \% of the cases. In addition, compared to the distribution of dark matter within halos, satellite ELGs preferentially reside both in the outskirts and inside haloes, and have a net infall velocity towards the centre. Finally, our results show evidence of assembly bias and conformity., Comment: 23 pages, 17 figures. To be submitted to A&A

Context. Mock galaxy catalogues are essential for correctly interpreting current and future generations of galaxy surveys. Despite their significance in galaxy formation and cosmology, little to no work has been done to validate the predictions of these mocks for high-order clustering statistics. Aims. We compare the predicting power of the latest generation of empirical models used in the creation of mock galaxy catalogues: a 13-parameter Halo Occupation Distribution (HOD) and an extension of the SubHalo Abundance Matching technique (SHAMe). Methods. We build GalaxyEmu-Planck, an emulator that makes precise predictions for the two-point correlation function, galaxy-galaxy lensing (restricted to distances greater than 1 $h^{-1} {\rm Mpc}$ to avoid baryonic effects), and other high-order statistics resulting from the evaluation of SHAMe and HOD models. Results. We evaluate the precision of GalaxyEmu-Planck using two galaxy samples extracted from the FLAMINGO hydrodynamical simulation that mimic the properties of DESI-BGS and BOSS galaxies, finding that the emulator reproduces all the predicted statistics precisely. The HOD showed comparable performance when fitting galaxy clustering and galaxy-galaxy lensing. In contrast, the SHAMe model showed better predictions for higher-order statistics, especially regarding the galaxy assembly bias. We also tested the performance of the models after removing some of their extensions, finding that we can withdraw two of the HOD parameters without a loss of performance. Conclusions. The results of this paper validate the current generation of empirical models as a way to reproduce galaxy clustering, galaxy-galaxy lensing and other high-order statistics. The excellent performance of the SHAMe model with a small number of free parameters suggests that it is a valid method to extract cosmological constraints from galaxy clustering., Comment: 22 pages, 12 figures, submitted to A&A. Comments welcomed!

We study the impact of black hole nuclear activity on both the global and radial star formation rate (SFR) profiles in X-ray-selected active galactic nuclei (AGN) in the field of miniJPAS, the precursor of the much wider J-PAS project. Our sample includes 32 AGN with z < 0.3 detected via the XMM-Newton and Chandra surveys. For comparison, we assembled a control sample of 71 star-forming (SF) galaxies with similar magnitudes, sizes, and redshifts. To derive the global properties of both the AGN and the control SF sample, we used CIGALE to fit the spectral energy distributions derived from the 56 narrowband and 4 broadband filters from miniJPAS. We find that AGN tend to reside in more massive galaxies than their SF counterparts. After matching samples based on stellar mass and comparing their SFRs and specific SFRs (sSFRs), no significant differences appear. This suggests that the presence of AGN does not strongly influence overall star formation. However, when we used miniJPAS as an integral field unit (IFU) to dissect galaxies along their position angle, a different picture emerges. We find that AGN tend to be more centrally concentrated in mass with respect to SF galaxies. Moreover, we find a suppression of the sSFR up to 1Re and then an enhancement beyond 1Re , strongly contrasting with the decreasing radial profile of sSFRs in SF galaxies. This could point to an inside-out quenching of AGN host galaxies. These findings suggest that the reason we do not see differences on a global scale is because star formation is suppressed in the central regions and enhanced in the outer regions of AGN host galaxies. While limited in terms of sample size, this work highlights the potential of the upcoming J-PAS as a wide-field low-resolution IFU for thousands of nearby galaxies and AGN., Comment: 21 pages, 20 figures, Accepted for Publication in A&A

We present a list of quasar candidates including photometric redshift estimates from the miniJPAS Data Release constructed using SQUEzE. This work is based on machine-learning classification of photometric data of quasar candidates using SQUEzE. It has the advantage that its classification procedure can be explained to some extent, making it less of a `black box' when compared with other classifiers. Another key advantage is that using user-defined metrics means the user has more control over the classification. While SQUEzE was designed for spectroscopic data, here we adapt it for multi-band photometric data, i.e. we treat multiple narrow-band filters as very low-resolution spectra. We train our models using specialized mocks from Queiroz et al. (2022). We estimate our redshift precision using the normalized median absolute deviation, $\sigma_{\rm NMAD}$ applied to our test sample. Our test sample returns an $f_1$ score (effectively the purity and completeness) of 0.49 for quasars down to magnitude $r=24.3$ with $z\geq2.1$ and 0.24 for quasars with $z<2.1$. For high-z quasars, this goes up to 0.9 for $r<21.0$. We present two catalogues of quasar candidates including redshift estimates: 301 from point-like sources and 1049 when also including extended sources. We discuss the impact of including extended sources in our predictions (they are not included in the mocks), as well as the impact of changing the noise model of the mocks. We also give an explanation of SQUEzE reasoning. Our estimates for the redshift precision using the test sample indicate a $\sigma_{NMAD}=0.92\%$ for the entire sample, reduced to 0.81\% for $r<22.5$ and 0.74\% for $r<21.3$. Spectroscopic follow-up of the candidates is required in order to confirm the validity of our findings., Comment: Accepted in A&A 24 pages, 24 figures, 7 tables

We present the first measurements of Lyman-$\alpha$ (Ly$\alpha$) forest correlations using early data from the Dark Energy Spectroscopic Instrument (DESI). We measure the auto-correlation of Ly$\alpha$ absorption using 88,509 quasars at $z>2$, and its cross-correlation with quasars using a further 147,899 tracer quasars at $z\gtrsim1.77$. Then, we fit these correlations using a 13-parameter model based on linear perturbation theory and find that it provides a good description of the data across a broad range of scales. We detect the BAO peak with a signal-to-noise ratio of $3.8\sigma$, and show that our measurements of the auto- and cross-correlations are fully-consistent with previous measurements by the Extended Baryon Oscillation Spectroscopic Survey (eBOSS). Even though we only use here a small fraction of the final DESI dataset, our uncertainties are only a factor of 1.7 larger than those from the final eBOSS measurement. We validate the existing analysis methods of Ly$\alpha$ correlations in preparation for making a robust measurement of the BAO scale with the first year of DESI data.

We present the one-dimensional Lyman-$\alpha$ forest power spectrum measurement using the first data provided by the Dark Energy Spectroscopic Instrument (DESI). The data sample comprises $26,330$ quasar spectra, at redshift $z > 2.1$, contained in the DESI Early Data Release and the first two months of the main survey. We employ a Fast Fourier Transform (FFT) estimator and compare the resulting power spectrum to an alternative likelihood-based method in a companion paper. We investigate methodological and instrumental contaminants associated to the new DESI instrument, applying techniques similar to previous Sloan Digital Sky Survey (SDSS) measurements. We use synthetic data based on log-normal approximation to validate and correct our measurement. We compare our resulting power spectrum with previous SDSS and high-resolution measurements. With relatively small number statistics, we successfully perform the FFT measurement, which is already competitive in terms of the scale range. At the end of the DESI survey, we expect a five times larger Lyman-$\alpha$ forest sample than SDSS, providing an unprecedented precise one-dimensional power spectrum measurement., Comment: 23 pages, 23 figures, Journal version

The Lyman-$\alpha$ forest offers a unique avenue for studying the distribution of matter in the high redshift universe and extracting precise constraints on the nature of dark matter, neutrino masses, and other $\Lambda$CDM extensions. However, interpreting this observable requires accurate modelling of the thermal and ionisation state of the intergalactic medium, and therefore resorting to computationally expensive hydrodynamical simulations. In this work, we build a neural network that serves as a surrogate model for rapid predictions of the one-dimensional \lya flux power spectrum ($P_{\rm 1D}$), thereby making Bayesian inference feasible for this observable. Our emulation technique is based on modelling $P_{\rm 1D}$ as a function of the slope and amplitude of the linear matter power spectrum rather than as a function of cosmological parameters. We show that our emulator achieves sub-percent precision across the full range of scales ($k_{\parallel }=0.1$ to 4Mpc$^{-1}$) and redshifts ($z=2$ to 4.5) considered, and also for three $\Lambda$CDM extensions not included in the training set: massive neutrinos, running of the spectral index, and curvature. Furthermore, we show that it performs at the 1% level for ionisation and thermal histories not present in the training set and performs at the percent level when emulating down to $k_{\parallel}$=8Mpc$^{-1}$. These results affirm the efficacy of our emulation strategy in providing accurate predictions even for cosmologies and reionisation histories that were not explicitly incorporated during the training phase, and we expect it to play a critical role in the cosmological analysis of the DESI survey., Comment: 15 pages, 17 figures

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

This paper is part of large effort within the J-PAS collaboration that aims to classify point-like sources in miniJPAS, which were observed in 60 optical bands over $\sim$ 1 deg$^2$ in the AEGIS field. We developed two algorithms based on artificial neural networks (ANN) to classify objects into four categories: stars, galaxies, quasars at low redshift ($z < 2.1)$, and quasars at high redshift ($z \geq 2.1$). As inputs, we used miniJPAS fluxes for one of the classifiers (ANN$_1$) and colours for the other (ANN$_2$). The ANNs were trained and tested using mock data in the first place. We studied the effect of augmenting the training set by creating hybrid objects, which combines fluxes from stars, galaxies, and quasars. Nevertheless, the augmentation processing did not improve the score of the ANN. We also evaluated the performance of the classifiers in a small subset of the SDSS DR12Q superset observed by miniJPAS. In the mock test set, the f1-score for quasars at high redshift with the ANN$_1$ (ANN$_2$) are $0.99$ ($0.99$), $0.93$ ($0.92$), and $0.63$ ($0.57$) for $17 < r \leq 20$, $20 < r \leq 22.5$, and $22.5 < r \leq 23.6$, respectively, where $r$ is the J-PAS rSDSS band. In the case of low-redshift quasars, galaxies, and stars, we reached $0.97$ ($0.97$), $0.82$ ($0.79$), and $0.61$ ($0.58$); $0.94$ ($0.94$), $0.90$ ($0.89$), and $0.81$ ($0.80$); and $1.0$ ($1.0$), $0.96$ ($0.94$), and $0.70$ ($0.52$) in the same r bins. In the SDSS DR12Q superset miniJPAS sample, the weighted f1-score reaches 0.87 (0.88) for objects that are mostly within $20 < r \leq 22.5$. Finally, we estimate the number of point-like sources that are quasars, galaxies, and stars in miniJPAS.

Astrophysical surveys rely heavily on the classification of sources as stars, galaxies or quasars from multi-band photometry. Surveys in narrow-band filters allow for greater discriminatory power, but the variety of different types and redshifts of the objects present a challenge to standard template-based methods. In this work, which is part of larger effort that aims at building a catalogue of quasars from the miniJPAS survey, we present a Machine Learning-based method that employs Convolutional Neural Networks (CNNs) to classify point-like sources including the information in the measurement errors. We validate our methods using data from the miniJPAS survey, a proof-of-concept project of the J-PAS collaboration covering $\sim$ 1 deg$^2$ of the northern sky using the 56 narrow-band filters of the J-PAS survey. Due to the scarcity of real data, we trained our algorithms using mocks that were purpose-built to reproduce the distributions of different types of objects that we expect to find in the miniJPAS survey, as well as the properties of the real observations in terms of signal and noise. We compare the performance of the CNNs with other well-established Machine Learning classification methods based on decision trees, finding that the CNNs improve the classification when the measurement errors are provided as inputs. The predicted distribution of objects in miniJPAS is consistent with the putative luminosity functions of stars, quasars and unresolved galaxies. Our results are a proof-of-concept for the idea that the J-PAS survey will be able to detect unprecedented numbers of quasars with high confidence., Comment: 16 pages, 15 figures, published by MNRAS

Satellite galaxies undergo a variety of physical processes when they are accreted by groups and clusters, often resulting in the loss of baryonic and dark matter (DM) mass. In this work, we evaluate the predictions from the IllustrisTNG hydrodynamical simulation regarding the evolution of the matter content of satellites, focusing on a population that are accreted at $z>1$ and retain their identity as satellites down to $z=0$. At fixed host halo mass, the amount of DM and stellar mass stripped depends mostly on the pericentric distance, $d_{\rm peri}$, here normalised by host halo virial radius. The closest encounters result in significant loss of DM, with subhaloes retaining between 20 and a few per cent of their $z=1$ mass. At fixed $d_{\rm peri}$, DM mass stripping seems more severe in lower mass haloes. Conversely, the average satellite in higher mass haloes has its stellar mass growth halted earlier, having lost a higher fraction of stellar mass by $z=0$. We also show that mass stripping has a strong impact on the quenched fractions. The IllustrisTNG boxes are qualitatively consistent in these predictions, with quantitative differences mostly originating from the distinct subhalo mass ranges covered by the boxes. Finally, we have identified DM-deficient systems in all TNG boxes. These objects are preferentially found in massive clusters ($M_{\rm host } \gtrsim 10^{13}$ M$_\odot$), had very close encounters with their central galaxies ($d_{\rm peri}\simeq0.05\, R_{\rm vir}$), and were accreted at high redshift ($z_{\rm infall} \gtrsim 1.4$), reinforcing the notion that tidal stripping is responsible for their remarkable lack of DM., Comment: 16 pages, 12 figures. Submitted to MNRAS, comments welcome

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

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

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

We present Diffstar, a smooth parametric model for the in-situ star formation history (SFH) of galaxies. Diffstar is distinct from conventional SFH models that are used to interpret the spectral energy distribution (SED) of an observed galaxy, because our model is parametrized directly in terms of basic features of galaxy formation physics. The Diffstar model assumes that star formation is fueled by the accretion of gas into the dark matter halo of the galaxy, and at the foundation of Diffstar is a parametric model for halo mass assembly, Diffmah. We include parametrized ingredients for the fraction of accreted gas that is eventually transformed into stars, $\epsilon_{\rm ms},$ and for the timescale over which this transformation occurs, $\tau_{\rm cons};$ some galaxies in Diffstar experience a quenching event at time $t_{\rm q},$ and may subsequently experience rejuvenated star formation. We fit the SFHs of galaxies predicted by the IllustrisTNG (TNG) and UniverseMachine (UM) simulations with the Diffstar parameterization, and show that our model is sufficiently flexible to describe the average stellar mass histories of galaxies in both simulations with an accuracy of $\sim0.1$ dex across most of cosmic time. We use Diffstar to compare TNG to UM in common physical terms, finding that: (i) star formation in UM is less efficient and burstier relative to TNG; (ii) galaxies in UM have longer gas consumption timescales, $\tau_{\rm cons}$, relative to TNG; (iii) rejuvenated star formation is ubiquitous in UM, whereas quenched TNG galaxies rarely experience sustained rejuvenation; and (iv) in both simulations, the distributions of $\epsilon_{\rm ms}$, $\tau_{\rm cons}$, and $t_{\rm q}$ share a common characteristic dependence upon halo mass, and present significant correlations with halo assembly history. [Abridged], Comment: 26 pages, 21 figures

The Javalambre-Physics of the Accelerating Universe Astrophysical Survey (J-PAS) is expected to map thousands of square degrees of the northern sky with 56 narrowband filters in the upcoming years. This will make J-PAS a very competitive and unbiased emission line survey compared to spectroscopic or narrowband surveys with fewer filters. The miniJPAS survey covered 1 deg$^2$, and it used the same photometric system as J-PAS, but the observations were carried out with the pathfinder J-PAS camera. In this work, we identify and characterize the sample of emission line galaxies (ELGs) from miniJPAS with a redshift lower than $0.35$. Using a method based on artificial neural networks, we detect the ELG population and measure the equivalent width and flux of the $H\alpha$, $H\beta$, [OIII], and [NII] emission lines. We explore the ionization mechanism using the diagrams [OIII]/H$\beta$ versus [NII]/H$\alpha$ (BPT) and EW(H$\alpha$) versus [NII]/H$\alpha$ (WHAN). We identify 1787 ELGs ($83$%) from the parent sample (2154 galaxies) in the AEGIS field. For the galaxies with reliable EW values that can be placed in the WHAN diagram (2000 galaxies in total), we obtained that $72.8 \pm 0.4$%, $17.7 \pm 0.4$% , and $9.4 \pm 0.2$% are star-forming (SF), active galactic nucleus (Seyfert), and quiescent galaxies, respectively. Based on the flux of $H\alpha$ we find that the star formation main sequence is described as $\log$ SFR $[M_\mathrm{\odot} \mathrm{yr}^{-1}] = 0.90^{+ 0.02}_{-0.02} \log M_{\star} [M_\mathrm{\odot}] -8.85^{+ 0.19}_{-0.20}$ and has an intrinsic scatter of $0.20^{+ 0.01}_{-0.01}$. The cosmic evolution of the SFR density ($\rho_{\text{SFR}}$) is derived at three redshift bins: $0 < z \leq 0.15$, $0.15 < z \leq 0.25$, and $0.25 < z \leq 0.35$, which agrees with previous results that were based on measurements of the $H\alpha$ emission line., Comment: 22 pages, 19 figures

We compute the linear order, general relativistic corrections to angular redshift fluctuations (ARF), a new cosmological observable built upon density-weighted two-dimensional (2D) maps of galaxy redshifts. We start with an existing approach for galaxy/source counts developed in the Newtonian gauge, and generalize it to ARF, modifying for this purpose a standard Boltzmann code. Our calculations allow us identifying the velocity terms as the leading corrections on large scales, emphasizing the sensitivity of ARF to peculiar, cosmological velocity fields. Just like for standard 2D clustering, the impact of gravitational lensing on ARF is dominant on small angular scales and for wide redshift shells, while the signatures associated to gravitational potentials are extremely small and hardly detectable. The ARF also present interesting correlation properties to anisotropies of the Cosmic Microwave Background (CMB): they are highly correlated to CMB lensing potential fluctuations, while also exhibiting a significant (S/N$\sim 4$-$5$) {\em anti-}correlation with the Integrated Sachs-Wolfe effect (ISW). This negative ARF$\times$ISW signal is quite complementary to the standard 2D clustering$\times$ISW correlation, since the former appears mostly at higher redshift ($z\sim 2$) than the latter ($z\lesssim 1)$, and the combination of the two observables significantly increases the $\chi^2$ statistics testing the null (no ISW) hypothesis. We conclude that ARF constitute a novel, alternative, and potentially powerful tool to constrain the nature of Dark Energy component that gives rise to the ISW., Comment: 25 pages, 7 figures, comments welcome!

In this series of papers, we employ several machine learning (ML) methods to classify the point-like sources from the miniJPAS catalogue, and identify quasar candidates. Since no representative sample of spectroscopically confirmed sources exists at present to train these ML algorithms, we rely on mock catalogues. In this first paper we develop a pipeline to compute synthetic photometry of quasars, galaxies and stars using spectra of objects targeted as quasars in the Sloan Digital Sky Survey. To match the same depths and signal-to-noise ratio distributions in all bands expected for miniJPAS point sources in the range $17.5\leq r<24$, we augment our sample of available spectra by shifting the original $r$-band magnitude distributions towards the faint end, ensure that the relative incidence rates of the different objects are distributed according to their respective luminosity functions, and perform a thorough modeling of the noise distribution in each filter, by sampling the flux variance either from Gaussian realizations with given widths, or from combinations of Gaussian functions. Finally, we also add in the mocks the patterns of non-detections which are present in all real observations. Although the mock catalogues presented in this work are a first step towards simulated data sets that match the properties of the miniJPAS observations, these mocks can be adapted to serve the purposes of other photometric surveys., Comment: 20 pages, 18 figures, submitted to MNRAS

Models of stellar population synthesis (SPS) are the fundamental tool that relates the physical properties of a galaxy to its spectral energy distribution (SED). In this paper, we present DSPS: a python package for stellar population synthesis. All of the functionality in DSPS is implemented natively in the JAX library for automatic differentiation, and so our predictions for galaxy photometry are fully differentiable, and directly inherit the performance benefits of JAX, including portability onto GPUs. DSPS also implements several novel features, such as i) a flexible empirical model for stellar metallicity that incorporates correlations with stellar age, and ii) support for the diffstar model that provides a physically-motivated connection between the star formation history of a galaxy (SFH) and the mass assembly of its underlying dark matter halo. We detail a set of theoretical techniques for using autodiff to calculate gradients of predictions for galaxy SEDs with respect to SPS parameters that control a range of physical effects, including SFH, stellar metallicity, nebular emission, and dust attenuation. When forward modeling the colors of a synthetic galaxy population, we find that DSPS can provide a factor of 5 speedup over standard SPS codes on a CPU, and a factor of 300-400 on a modern GPU. When coupled with gradient-based techniques for optimization and inference, DSPS makes it practical to conduct expansive likelihood analyses of simulation-based models of the galaxy--halo connection that fully forward model galaxy spectra and photometry., Comment: 13 pages, 2 appendices. Version accepted for publication in MNRAS

Photometric redshift estimation algorithms are often based on representative data from observational campaigns. Data-driven methods of this type are subject to a number of potential deficiencies, such as sample bias and incompleteness. Motivated by these considerations, we propose using physically motivated synthetic spectral energy distributions in redshift estimation. In addition, the synthetic data would have to span a domain in colour-redshift space concordant with that of the targeted observational surveys. With a matched distribution and realistically modelled synthetic data in hand, a suitable regression algorithm can be appropriately trained; we use a mixture density network for this purpose. We also perform a zero-point re-calibration to reduce the systematic differences between noise-free synthetic data and the (unavoidably) noisy observational data sets. This new redshift estimation framework, SYTH-Z, demonstrates superior accuracy over a wide range of redshifts compared to baseline models trained on observational data alone. Approaches using realistic synthetic data sets can therefore greatly mitigate the reliance on expensive spectroscopic follow-up for the next generation of photometric surveys., Comment: 14 pages, 8 figures

Precise measurements of black hole masses are essential to understanding the coevolution of these sources and their host galaxies. We develop a novel approach for computing black hole virial masses using measurements of continuum luminosities and emission line widths from partially overlapping, narrow-band observations of quasars; we refer to this technique as single-epoch photometry. This novel method relies on forward-modelling quasar observations for estimating emission line widths, which enables unbiased measurements even for lines coarsely resolved by narrow-band data. We assess the performance of this technique using quasars from the Sloan Digital Sky Survey (SDSS) observed by the miniJPAS survey, a proof-of-concept project of the Javalambre Physics of the Accelerating Universe Astrophysical Survey (J-PAS) collaboration covering $\simeq1\,\mathrm{deg}^2$ of the northern sky using the 56 J-PAS narrow-band filters. We find remarkable agreement between black hole masses from single-epoch SDSS spectra and single-epoch miniJPAS photometry, with no systematic difference between these and a scatter ranging from 0.4 to 0.07 dex for masses from $\log(M_\mathrm{BH})\simeq8$ to 9.75, respectively. Reverberation mapping studies show that single-epoch masses present approximately 0.4 dex precision, letting us conclude that our novel technique delivers black hole masses with only mildly lower precision than single-epoch spectroscopy. The J-PAS survey will soon start observing thousands of square degrees without any source preselection other than the photometric depth in the detection band, and thus single-epoch photometry has the potential to provide details on the physical properties of quasar populations that do not satisfy the preselection criteria of previous spectroscopic surveys., Comment: 16 pages, 10 figures, accepted by A&A

Among the properties shaping the light of a galaxy, the star formation history (SFH) is one of the most challenging to model due to the variety of correlated physical processes regulating star formation. In this work, we leverage the stellar population synthesis model FSPS, together with SFHs predicted by the hydrodynamical simulation IllustrisTNG and the empirical model UNIVERSEMACHINE, to study the impact of star formation variability on galaxy colours. We start by introducing a model-independent metric to quantify the burstiness of a galaxy formation model, and we use this metric to demonstrate that UNIVERSEMACHINE predicts SFHs with more burstiness relative to IllustrisTNG. Using this metric and principal component analysis, we construct families of SFH models with adjustable variability, and we show that the precision of broad-band optical and near-infrared colours degrades as the level of unresolved short-term variability increases. We use the same technique to demonstrate that variability in metallicity and dust attenuation presents a practically negligible impact on colours relative to star formation variability. We additionally provide a model-independent fitting function capturing how the level of unresolved star formation variability translates into imprecision in predictions for galaxy colours; our fitting function can be used to determine the minimal SFH model that reproduces colours with some target precision. Finally, we show that modelling the colours of individual galaxies with percent-level precision demands resorting to complex SFH models, while producing precise colours for galaxy populations can be achieved using models with just a few degrees of freedom., Comment: 18 pages, 15 figures, accepted by MNRAS

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

We present a new empirical model for the mass assembly of dark matter halos. We approximate the growth of individual halos as a simple power-law function of time, where the power-law index smoothly decreases as the halo transitions from the fast-accretion regime at early times, to the slow-accretion regime at late times. Using large samples of halo merger trees taken from high-resolution cosmological simulations, we demonstrate that our 3-parameter model, Diffmah, can approximate halo growth with a typical accuracy of 0.1 dex for t > 1 Gyr for all halos of present-day mass greater than 10^11Msun, including subhalos and host halos in gravity-only simulations, as well as in the TNG hydrodynamical simulation. We additionally present a new model for the assembly of halo populations, DiffmahPop, which not only reproduces average mass growth across time, but also faithfully captures the diversity with which halos assemble their mass. Our python implementation is based on the autodiff library JAX, and so our model self-consistently captures the mean and variance of halo mass accretion rate across cosmic time. We show that the connection between halo assembly and the large-scale density field, known as halo assembly bias, is accurately captured by Diffmah, and that residual errors in our approximations to halo assembly history exhibit a negligible residual correlation with the density field. Our publicly available source code can be used to generate Monte Carlo realizations of cosmologically representative halo histories; our differentiable implementation facilitates the incorporation of our model into existing analytical halo model frameworks., Comment: Main body: 11 pages, 8 figures. Appendices: 11 pages, 11 figures. Update to the version originally accepted by the Open Journal of Astrophysics; revisions limited to textual changes in the nomenclature used to refer to the model

Halo assembly bias is the secondary dependence of the clustering of dark-matter haloes on their assembly histories at fixed halo mass. This established dependence is expected to manifest itself on the clustering of the galaxy population, a potential effect commonly known as galaxy assembly bias. Using the IllustrisTNG300 magnetohydrodynamical simulation, we analyse the dependence of the properties and clustering of galaxies on the shape of the specific mass accretion history of their hosting haloes (sMAH). We first show that several halo and galaxy properties strongly correlate with the slope of the sMAH ($\beta$) at fixed halo mass. Namely, haloes with increasingly steeper $\beta$ increment their halo masses faster at early times, and their hosted galaxies present larger stellar-to-halo mass ratios, lose their gas faster, reach the peak of their star formation histories at higher redshift, and become quenched earlier. We also demonstrate that $\beta$ is more directly connected to these key galaxy formation properties than other broadly employed halo proxies, such as formation time. Finally, we measure the secondary dependence of galaxy clustering on $\beta$ at fixed halo mass as a function of redshift. By tracing back the evolution of individual haloes, we show that the amplitude of the galaxy assembly bias signal for the progenitors of $z=0$ galaxies increases with redshift, reaching a factor of 2 at $z = 1$ for haloes of $M_\mathrm{halo}=10^{11.5}-10^{12}$ $h^{-1}\mathrm{M}_\odot$. The measurement of the evolution of assembly bias along the merger tree provides a new theoretical perspective to the study of secondary bias. Our findings, which show a tight relationship between halo accretion and both the clustering and the observational properties of the galaxy population, have also important implications for the generation of mock catalogues for upcoming cosmological surveys., Comment: 10 pages, 6 figures. Submitted to MNRAS

Handling big data has largely been a major bottleneck in traditional statistical models. Consequently, when accurate point prediction is the primary target, machine learning models are often preferred over their statistical counterparts for bigger problems. But full probabilistic statistical models often outperform other models in quantifying uncertainties associated with model predictions. We develop a data-driven statistical modeling framework that combines the uncertainties from an ensemble of statistical models learned on smaller subsets of data carefully chosen to account for imbalances in the input space. We demonstrate this method on a photometric redshift estimation problem in cosmology, which seeks to infer a distribution of the redshift -- the stretching effect in observing the light of far-away galaxies -- given multivariate color information observed for an object in the sky. Our proposed method performs balanced partitioning, graph-based data subsampling across the partitions, and training of an ensemble of Gaussian process models.

Fluctuations in sky maps of the galaxy redshifts, dubbed as angular redshift fluctuations (ARF), contain precise information about the growth rate of structures and the nature of gravity in the Universe. Specifically, ARF constrain the combination of cosmological parameters $H/H_0\,f\sigma_8(z)$, while being an intrinsically tomographic probe and largely insensitive to many observational systematic errors, all this without requiring the assumption of any redshift-to-distance relation under a given fiducial cosmology. We present the first cosmological constraints derived from ARF by using BOSS LOWZ+CMASS DR12 galaxy samples, obtaining 7\%-accurate constraints on $H/H_0 f\sigma_8(z)$ at more than 20 redshifts over the range $z \in [0.26,0.72]$. Our best-fitting value is $10\%$ larger, but compatible at the $1.4\sigma$ level, than the $\Lambda$CDM expectation set by {\it Planck} observations of the Cosmic Microwave Background (CMB) radiation. Our tomographic measurements, combined with these CMB data, provides one of the strongest constraints on the gravity index $\gamma$, $\gamma=0.44^{+0.09}_{-0.07}$, which lies within $2\sigma$ from the prediction of General Relativity ($\gamma_{\rm GR}\simeq 0.55$)., Comment: 5 pages, 3 figures, matches accepted version in MNRAS letters

We propose the use of angular fluctuations in the galaxy redshift field as a new way to extract cosmological information in the Universe. This new probe $\delta z (\hat{n})$ consists on the statistics of sky maps built by projecting redshifts under a Gaussian window of width $\sigma_z$ centred upon a redshift $z_{\rm obs}$, and weighted by the galaxy density field. We compute the angular power spectrum of the $\delta z (\hat{n})$ field in both numerical simulations and in linear perturbation theory. From these we find that the $\delta z (\hat{n})$ field: {\it (i)} is sensitive to the underlying density and peculiar velocity fields; {\it (ii)} is highly correlated, at the $\gtrsim 60\,\%$ level, to the line-of-sight projected peculiar velocity field; {\it (iii)} for narrow windows $(\sigma_z < 0.03$), it is almost completely uncorrelated to the projected galaxy angular density field under the same redshift window; and {\it (iv)} it is largely unaffected by multiplicative and additive systematic errors on the observed number of galaxies that are redshift-independent over $\sim\sigma_z$. We conclude that $\delta z (\hat{n})$ is a simple and robust tomographic measure of the cosmic density and velocity fields, complementary to angular clustering, that will contribute to more complete exploitations of current and upcoming galaxy redshift surveys., Comment: 6 pages, 3 figures, matches accepted version in MNRAS letters, companion paper of Chaves-Montero et al., arXiv:1911.10690

A complete census of baryons in the late universe is a long-standing challenge due to the intermediate temperate and rarefied character of the majority of cosmic gas. To gain insight into this problem, we extract measurements of the kinematic Sunyaev-Zel'dovich (kSZ) effect from the cross-correlation of angular redshift fluctuations maps, which contain precise information about the cosmic density and velocity fields, and CMB maps high-pass filtered using aperture photometry; we refer to this technique as ARF-kSZ tomography. Remarkably, we detect significant cross-correlation for a wide range of redshifts and filter apertures using 6dF galaxies, BOSS galaxies, and SDSS quasars as tracers, yielding a 11 sigma detection of the kSZ effect. We then leverage these measurements to set constraints on the location, density, and abundance of gas inducing the kSZ effect, finding that this gas resides outside dark matter haloes, presents densities ranging from 10 to 250 times the cosmic average, and comprises half of cosmic baryons. Taken together, these findings indicate that ARF-kSZ tomography provides a nearly complete census of intergalactic gas from z=0 to 5., Comment: 17 pages, 14 figures, accepted for publication in MNRAS

The spectral energy distribution of a galaxy emerges from the complex interplay of many physical ingredients, including its star formation history (SFH), metallicity evolution, and dust properties. Using GALAXPY, a new galaxy spectral prediction tool, and SFHs predicted by the empirical model UNIVERSEMACHINE and the cosmological hydrodynamical simulation IllustrisTNG, we isolate the influence of SFH on optical and near-infrared colours from 3200 to 10800 \r{A} at z=0. By carrying out a principal component analysis, we show that physically-motivated SFH variations modify galaxy colours along a single direction in colour space: the SFH-direction. We find that the projection of a galaxy's present-day colours onto the SFH-direction is almost completely regulated by the fraction of stellar mass that the galaxy formed over the last billion years. Together with cosmic downsizing, this results in galaxies becoming redder as their host halo mass increases. We additionally study the change in galaxy colours due to variations in metallicity, dust attenuation, and nebular emission lines, finding that these properties vary broad-band colours along distinct directions in colour space relative to the SFH-direction. Finally, we show that the colours of low-redshift SDSS galaxies span an ellipsoid with significant extent along two independent dimensions, and that the SFH-direction is well-aligned with the major axis of this ellipsoid. Our analysis supports the conclusion that variations in star formation history are the dominant influence on present-day galaxy colours, and that the nature of this influence is strikingly simple., Comment: 17 pages, 12 figures, accepted by MNRAS

The expansion of the Universe is understood to have accelerated during two epochs: in its very first moments during a period of Inflation and much more recently, at $z < 1$, when Dark Energy is hypothesized to drive cosmic acceleration. The undiscovered mechanisms behind these two epochs represent some of the most important open problems in fundamental physics. The large cosmological volume at $2 < z < 5$, together with the ability to efficiently target high-$z$ galaxies with known techniques, enables large gains in the study of Inflation and Dark Energy. A future spectroscopic survey can test the Gaussianity of the initial conditions up to a factor of ~50 better than our current bounds, crossing the crucial theoretical threshold of $\sigma(f_{NL}^{\rm local})$ of order unity that separates single field and multi-field models. Simultaneously, it can measure the fraction of Dark Energy at the percent level up to $z = 5$, thus serving as an unprecedented test of the standard model and opening up a tremendous discovery space., Comment: Science white paper submitted to the Astro2020 Decadal Survey

Our current understanding of the Universe is established through the pristine measurements of structure in the cosmic microwave background (CMB) and the distribution and shapes of galaxies tracing the large scale structure (LSS) of the Universe. One key ingredient that underlies cosmological observables is that the field that sources the observed structure is assumed to be initially Gaussian with high precision. Nevertheless, a minimal deviation from Gaussianityis perhaps the most robust theoretical prediction of models that explain the observed Universe; itis necessarily present even in the simplest scenarios. In addition, most inflationary models produce far higher levels of non-Gaussianity. Since non-Gaussianity directly probes the dynamics in the early Universe, a detection would present a monumental discovery in cosmology, providing clues about physics at energy scales as high as the GUT scale., Comment: 5 pages + references; Submitted to the Astro2020 call for science white papers. This version: fixed author list

SPHEREx is a proposed NASA MIDEX mission selected for Phase A study. SPHEREx would carry out the first all-sky spectral survey in the near infrared. At the end of its two-year mission, SPHEREx would obtain 0.75-to-5$\mu$m spectra of every 6.2 arcsec pixel on the sky, with spectral resolution R>35 and a 5-$\sigma$ sensitivity AB$>$19 per spectral/spatial resolution element. More details concerning SPHEREx are available at http://spherex.caltech.edu. The SPHEREx team has proposed three specific science investigations to be carried out with this unique data set: cosmic inflation, interstellar and circumstellar ices, and the extra-galactic background light. Though these three themes are undoubtedly compelling, they are far from exhausting the scientific output of SPHEREx. Indeed, SPHEREx would create a unique all-sky spectral database including spectra of very large numbers of astronomical and solar system targets, including both extended and diffuse sources. These spectra would enable a wide variety of investigations, and the SPHEREx team is dedicated to making the data available to the community to enable these investigations, which we refer to as Legacy Science. To that end, we have sponsored two workshops for the general scientific community to identify the most interesting Legacy Science themes and to ensure that the SPHEREx data products are responsive to their needs. In February of 2016, some 50 scientists from all fields met in Pasadena to develop these themes and to understand their implications for the SPHEREx mission. The 2016 workshop highlighted many synergies between SPHEREx and other contemporaneous astronomical missions, facilities, and databases. Consequently, in January 2018 we convened a second workshop at the Center for Astrophysics in Cambridge to focus specifically on these synergies. This white paper reports on the results of the 2018 SPHEREx workshop., Comment: 50 pages, 24 figures, more details at http://spherex.caltech.edu

We present ELDAR, a new method that exploits the potential of medium- and narrow-band filter surveys to securely identify active galactic nuclei (AGN) and determine their redshifts. Our methodology improves on traditional approaches by looking for AGN emission lines expected to be identified against the continuum, thanks to the width of the filters. To assess its performance, we apply ELDAR to the data of the ALHAMBRA survey, which covered an effective area of $2.38\,{\rm deg}^2$ with 20 contiguous medium-band optical filters down to F814W$\simeq 24.5$. Using two different configurations of ELDAR in which we require the detection of at least 2 and 3 emission lines, respectively, we extract two catalogues of type-I AGN. The first is composed of 585 sources ($79\,\%$ of them spectroscopically-unknown) down to F814W$=22.5$ at $z_{\rm phot}>1$, which corresponds to a surface density of $209\,{\rm deg}^{-2}$. In the second, the 494 selected sources ($83\,\%$ of them spectroscopically-unknown) reach F814W$=23$ at $z_{\rm phot}>1.5$, for a corresponding number density of $176\,{\rm deg}^{-2}$. Then, using samples of spectroscopically-known AGN in the ALHAMBRA fields, for the two catalogues we estimate a completeness of $73\,\%$ and $67\,\%$, and a redshift precision of $1.01\,\%$ and $0.86\,\%$ (with outliers fractions of $8.1\,\%$ and $5.8\,\%$). At $z>2$, where our selection performs best, we reach $85\,\%$ and $77\,\%$ completeness and we find no contamination from galaxies., Comment: 23 pages, 18 figures, submitted to MNRAS

In the upcoming era of high-precision galaxy surveys, it becomes necessary to understand the impact of redshift uncertainties on cosmological observables. In this paper we explore the effect of sub-percent photometric redshift errors (photo-$z$ errors) on galaxy clustering and baryonic acoustic oscillations (BAO). Using analytic expressions and results from $1\,000$ $N$-body simulations, we show how photo-$z$ errors modify the amplitude of moments of the 2D power spectrum, their variances, the amplitude of BAO, and the cosmological information in them. We find that: a) photo-$z$ errors suppress the clustering on small scales, increasing the relative importance of shot noise, and thus reducing the interval of scales available for BAO analyses; b) photo-$z$ errors decrease the smearing of BAO due to non-linear redshift-space distortions (RSD) by giving less weight to line-of-sight modes; and c) photo-$z$ errors (and small-scale RSD) induce a scale dependence on the information encoded in the BAO scale, and that reduces the constraining power on the Hubble parameter. Using these findings, we propose a template that extracts unbiased cosmological information from samples with photo-$z$ errors with respect to cases without them. Finally, we provide analytic expressions to forecast the precision in measuring the BAO scale, showing that spectro-photometric surveys will measure the expansion history of the Universe with a precision competitive to that of spectroscopic surveys., Comment: 19 pages, 15 figures, resubmission. Included full analysis of power spectrum moments l=0, 2, and 4

Subhalo abundance matching (SHAM) is a widely-used method to connect galaxies with dark matter structures in numerical simulations. SHAM predictions agree remarkably well with observations, yet they still lack strong theoretical support. We examine the performance, implementation, and assumptions of SHAM using the EAGLE project simulations. We find that $V_{\rm relax}$, the highest value of the circular velocity attained by a subhalo while it satisfies a relaxation criterion, is the subhalo property that correlates most strongly with galaxy stellar mass ($M_{\rm star}$). Using this parameter in SHAM, we retrieve the real-space clustering of EAGLE to within our statistical uncertainties on scales greater than $2$ Mpc for galaxies with $8.77<\log_{10}(M_{\rm star}[M_\odot])<10.77$. Conversely, clustering is overestimated by $30\%$ on scales below $2$ Mpc for galaxies with $8.77<\log_{10}(M_{\rm star}[M_\odot])<9.77$ because SHAM slightly overpredicts the fraction of satellites in massive haloes compared to EAGLE. The agreement is even better in redshift-space, where the clustering is recovered to within our statistical uncertainties for all masses and separations. Additionally, we analyse the dependence of galaxy clustering on properties other than halo mass, i.e. the assembly bias. We demonstrate assembly bias alters the clustering in EAGLE by $20\%$ and that $V_{\rm relax}$ captures its effect to within $15\%$. We trace small differences in the clustering to the failure of SHAM as typically implemented, i.e. the $M_{\rm star}$ assigned to a subhalo does not depend on i) its host halo mass, ii) whether it is a central or a satellite. In EAGLE we find that these assumptions are not completely satisfied., Comment: 20 pages, 14 figures