For the next generation of spectroscopic galaxy surveys, it is important to forecast their performances and to accurately interpret their large data sets. For this purpose, it is necessary to consistently simulate different populations of galaxies, in particular Emission Line Galaxies (ELGs), less used in the past for cosmological purposes. In this work, we further the forward modeling approach presented in Fagioli et al. 2018, by extending the spectra simulator Uspec to model galaxies of different kinds with improved parameters from Tortorelli et al. 2020. Furthermore, we improve the modeling of the selection function by using the image simulator Ufig. We apply this to the Sloan Digital Sky Survey (SDSS), and simulate $\sim157,000$ multi-band images. We pre-process and analyse them to apply cuts for target selection, and finally simulate SDSS/BOSS DR14 galaxy spectra. We compute photometric, astrometric and spectroscopic properties for red and blue, real and simulated galaxies, finding very good agreement. We compare the statistical properties of the samples by decomposing them with Principal Component Analysis (PCA). We find very good agreement for red galaxies and a good, but less pronounced one, for blue galaxies, as expected given the known difficulty of simulating those. Finally, we derive stellar population properties, mass-to-light ratios, ages and metallicities, for all samples, finding again very good agreement. This shows how this method can be used not only to forecast cosmology surveys, but it is also able to provide insights into studies of galaxy formation and evolution.

In any imaging survey, measuring accurately the astronomical background light is crucial to obtain good photometry. This paper introduces BKGnet, a deep neural network to predict the background and its associated error. BKGnet has been developed for data from the Physics of the Accelerating Universe Survey (PAUS), an imaging survey using a 40 narrow-band filter camera (PAUCam). Images obtained with PAUCam are affected by scattered light: an optical effect consisting of light multiply that deposits energy in specific detector regions contaminating the science measurements. Fortunately, scattered light is not a random effect, but it can be predicted and corrected for. We have found that BKGnet background predictions are very robust to distorting effects, while still being statistically accurate. On average, the use of BKGnet improves the photometric flux measurements by 7% and up to 20% at the bright end. BKGnet also removes a systematic trend in the background error estimation with magnitude in the i-band that is present with the current PAU data management method. With BKGnet, we reduce the photometric redshift outlier rate

The galaxy Luminosity Function (LF) is a key observable for galaxy formation, evolution studies and for cosmology. In this work, we propose a novel technique to forward model wide-field galaxy surveys using the fast image simulator UFig and measure the LF of galaxies in the B-band. We use Approximate Bayesian Computation (ABC) to constrain the galaxy population model parameters of the simulations and match data from the Canada-France-Hawaii Telescope Legacy Survey (CFHTLS). We define a number of distance metrics between the simulated and the survey data. By exploring the parameter space of the galaxy population model through ABC to find the set of parameters that minimize these distance metrics, we obtain constraints on the LFs of blue and red galaxies as a function of redshift. We find that $\mathrm{M^*}$ fades by $\Delta \mathrm{M}^*_{\mathrm{0.1-1.0,b}} = 0.37 \pm 0.34$ and $\Delta \mathrm{M}^*_{\mathrm{0.1-1.0,r}} = 0.58 \pm 0.42$ magnitudes between redshift $\mathrm{z = 1}$ and $\mathrm{z = 0.1}$ for blue and red galaxies, respectively. We also find that $\phi^*$ for blue galaxies stays roughly constant between redshift $\mathrm{z = 0.1}$ and $\mathrm{z=1}$, while for red galaxies it decreases of $\sim 90\%$. We compare our results to other measurements, finding good agreement at all redshifts, for both blue and red galaxies. To further test our results, we compare the redshift distributions for survey and simulated data. We use the spectroscopic redshift distribution from the VIMOS Public Extragalactic Redshift Survey (VIPERS) and we apply the same selection in colours and magnitudes on our simulations. We find a good agreement between the survey and the simulated redshift distributions. We provide best-fit values and uncertainties for the parameters of the LF. This work offers excellent prospects for measuring other galaxy population properties as a function of redshift using ABC.

The PAU Survey (PAUS) is an innovative photometric survey with 40 narrow bands at the William Herschel Telescope (WHT). The narrow bands are spaced at 100\AA\ intervals covering the range 4500\AA\ to 8500\AA\ and, in combination with standard broad bands, enable excellent redshift precision. This paper describes the technique, galaxy templates and additional photometric calibration used to determine early photometric redshifts from PAUS. Using BCNz2, a new photometric redshift code developed for this purpose, we characterise the photometric redshift performance using PAUS data on the COSMOS field. Comparison to secure spectra from zCOSMOS DR3 shows that PAUS achieves $\sigma_{68} /(1+z) = 0.0037$ to $i_{\mathrm{AB}} < 22.5$ when selecting the best 50% of the sources based on a photometric redshift quality cut. Furthermore, a higher photo-z precision ($\sigma_{68}/(1+z) \sim 0.001$) is obtained for a bright and high quality selection, which is driven by the identification of emission lines. We conclude that PAUS meets its design goals, opening up a hitherto uncharted regime of deep, wide, and dense galaxy survey with precise redshifts that will provide unique insights into the formation, evolution and clustering of galaxies, as well as their intrinsic alignments.

Classification of stars and galaxies is a well-known astronomical problem that has been treated using different approaches, most of them relying on morphological information. In this paper, we tackle this issue using the low-resolution spectra from narrow band photometry, provided by the PAUS (Physics of the Accelerating Universe) survey. We find that, with the photometric fluxes from the 40 narrow band filters and without including morphological information, it is possible to separate stars and galaxies to very high precision, 98.4% purity with a completeness of 98.8% for objects brighter than I = 22.5. This precision is obtained with a Convolutional Neural Network as a classification algorithm, applied to the objects' spectra. We have also applied the method to the ALHAMBRA photometric survey and we provide an updated classification for its Gold sample.

The Physics of the Accelerating Universe (PAU) Survey is an international project for the study of cosmological parameters associated with Dark Energy. PAU's 18-CCD camera (PAUCam), installed at the prime focus of the William Herschel Telescope at the Roque de los Muchachos Observatory (La Palma, Canary Islands), scans part of the northern sky, to collect low resolution spectral information of millions of galaxies with its unique set of 40 narrow-band filters in the optical range from 450 nm to 850 nm, and a set of 6 standard broad band filters. The PAU data management (PAUdm) team is in charge of treating the data, including data transfer from the observatory to the PAU Survey data center, hosted at Port d'Informaci\'o Cient\'ifica (PIC). PAUdm is also in charge of the storage, data reduction and, finally, of making the results available to the scientific community. We describe the technical solutions adopted to cover different aspects of the PAU Survey data management, from the computing infrastructure to support the operations, to the software tools and web services for the data process orchestration and exploration. In particular we will focus on the PAU database, developed for the coordination of the different PAUdm tasks, and to preserve and guarantee the consistency of data and metadata.

Galaxy spectra are essential to probe the spatial distribution of galaxies in our Universe. To better interpret current and future spectroscopic galaxy redshift surveys, it is important to be able to simulate these data sets. We describe Uspec, a forward modeling tool to generate galaxy spectra taking into account some intrinsic galaxy properties as well as instrumental responses of a given telescope. The model for the intrinsic properties of the galaxy population, i.e., the luminosity functions, and size and spectral coefficients distribu- tions, was developed in an earlier work for broad-band imaging surveys [1], and we now aim to test the model further using spectroscopic data. We apply Uspec to the SDSS/CMASS sample of Luminous Red Galaxies (LRGs). We construct selection cuts that match those used to build this LRG sample, which we then apply to data and simulations in the same way. The resulting real and simulated average spectra show a good statistical agreement overall, with residual differences likely coming from a bluer galaxy population of the simulated sam- ple. We also do not explore the impact of non-solar element ratios in our simulations. For a quantitative comparison, we perform Principal Component Analysis (PCA) of the sets of spectra. By comparing the PCs constructed from simulations and data, we find good agree- ment for all components. The distributions of the eigencoefficients also show an appreciable overlap. We are therefore able to properly simulate the LRG sample taking into account the SDSS/BOSS instrumental responses. The differences between the two samples can be ascribed to the intrinsic properties of the simulated galaxy population, which can be reduced by further improvements of our modelling method in the future. We discuss how these results can be useful for the forward modeling of upcoming large spectroscopic surveys.

We present a mock catalogue for the Physics of the Accelerating Universe Survey (PAUS) and use it to quantify the competitiveness of the narrow band imaging for measuring spectral features and galaxy clustering. The mock agrees with observed number count and redshift distribution data. We demonstrate the importance of including emission lines in the narrow band fluxes. We show that PAUCam has sufficient resolution to measure the strength of the 4000\AA{} break to the nominal PAUS depth. We predict the evolution of a narrow band luminosity function and show how this can be affected by the OII emission line. We introduce new rest frame broad bands (UV and blue) that can be derived directly from the narrow band fluxes. We use these bands along with D4000 and redshift to define galaxy samples and provide predictions for galaxy clustering measurements. We show that systematic errors in the recovery of the projected clustering due to photometric redshift errors in PAUS are significantly smaller than the expected statistical errors. The galaxy clustering on two halo scales can be recovered quantatively without correction, and all qualitative trends seen in the one halo term are recovered. In this analysis mixing between samples reduces the expected contrast between the one halo clustering of red and blue galaxies and demonstrates the importance of a mock catalogue for interpreting galaxy clustering results. The mock catalogue is available on request at https://cosmohub.pic.es/home.

Weak gravitational lensing is a powerful probe of the dark sector, once measurement systematic errors can be controlled. In Refregier & Amara (2014), a calibration method based on forward modeling, called MCCL, was proposed. This relies on fast image simulations (e.g., UFig; Berge et al. 2013) that capture the key features of galaxy populations and measurement effects. The MCCL approach has been used in Herbel et al. (2017) to determine the redshift distribution of cosmological galaxy samples and, in the process, the authors derived a model for the galaxy population mainly based on broad-band photometry. Here, we test this model by forward modeling the 40 narrow-band photometry given by the novel PAU Survey (PAUS). For this purpose, we apply the same forced photometric pipeline on data and simulations using Source Extractor (Bertin & Arnouts 1996). The image simulation scheme performance is assessed at the image and at the catalogues level. We find good agreement for the distribution of pixel values, the magnitudes, in the magnitude-size relation and the interband correlations. A principal component analysis is then performed, in order to derive a global comparison of the narrow-band photometry between the data and the simulations. We use a `mixing' matrix to quantify the agreement between the observed and simulated sets of Principal Components (PCs). We find good agreement, especially for the first three most significant PCs. We also compare the coefficients of the PCs decomposition. While there are slight differences for some coefficients, we find that the distributions are in good agreement. Together, our results show that the galaxy population model derived from broad-band photometry is in good overall agreement with the PAUS data. This offers good prospect for incorporating spectral information to the galaxy model by adjusting it to the PAUS narrow-band data using forward modeling.

We analyse the Kormendy relations (KRs) of the two Frontier Fields clusters, Abell S1063, at z = 0.348, and MACS J1149.5+2223, at z = 0.542, exploiting very deep Hubble Space Telescope photometry and VLT/MUSE integral field spectroscopy. With this novel dataset, we are able to investigate how the KR parameters depend on the cluster galaxy sample selection and how this affects studies of galaxy evolution based on the KR. We define and compare four different galaxy samples according to: (a) S\'ersic indices: early-type ('ETG'), (b) visual inspection: 'ellipticals', (c) colours: 'red', (d) spectral properties: 'passive'. The classification is performed for a complete sample of galaxies with m$_{\textit{F814W}} \le$ 22.5 ABmag (M$_{*}$ $\gtrsim 10^{10.0}$ M$_{\odot}$). To derive robust galaxy structural parameters, we use two methods: (1) an iterative estimate of structural parameters using images of increasing size, in order to deal with closely separated galaxies and (2) different background estimations, to deal with the Intracluster light contamination. The comparison between the KRs obtained from the different samples suggests that the sample selection could affect the estimate of the best-fitting KR parameters. The KR built with ETGs is fully consistent with the one obtained for ellipticals and passive. On the other hand, the KR slope built on the red sample is only marginally consistent with those obtained with the other samples. We also release the photometric catalogue with structural parameters for the galaxies included in the present analysis.