We present constraints on the timescale of star formation variability and the correlation between star formation and host halo accretion histories of star-forming (SF) central galaxies from the scatter of the stellar-to-halo mass relation (SHMR). SF galaxies are found to have a tight relationship between their star formation rates and stellar masses on the so-called "star-forming sequence" (SFS), which characterizes their SF histories and $M_*$ growths. Meanwhile, observed SHMR constraints connect $M_*$ growth to halo accretion history. Combining these observed trends with a cosmological $N$-body simulation, we present flexible models that track the SF, $M_*$, and halo accretion histories of SF central galaxies at $z<1$ while reproducing the observed stellar mass function and SFS in SDSS DR7. We find that the scatter in SHMR at $M_h=10^{12}M_\odot$, $\sigma_{M_*|M_h=10^{12}M_\odot}$, is sensitive to the timescale of star formation variability, $t_{\rm duty}$, and the correlation coefficient, $r$, between SF and halo accretion histories: shorter $t_{\rm duty}$ and higher $r$ both result in tighter $\sigma_{M_*|M_h=10^{12}M_\odot}$. To reproduce a constant 0.2 dex scatter over $z=1-0$, our models require $t_{\rm duty}\leq1.5$ Gyr for $r=0.99$ or $r>0.6$ for $t_{\rm duty}=0.1$ Gyr. For $r\sim0.6$, as found in the literature, $t_{\rm duty}<0.2$ Gyr is necessary. Meanwhile, to reproduce the tightening of $\sigma_{M_*|M_h=10^{12}M_\odot}=0.35$ to 0.2 dex from $z=1-0$ in hydrodynamical simulations, our models require $t_{\rm duty}=0.1$ Gyr for $r>0.5$. Although, the lack of consensus on $\sigma_{M_*|M_h=10^{12}M_\odot}$ at $M_h=10^{12}M_\odot$ and at $z=1$ from observations and galaxy formation models remains the main bottleneck in precisely constraining $r$ and $t_{\rm duty}$, we demonstrate that SHMR can be used to constrain the SF and host halo accretion histories of SF central galaxies.

Central galaxies make up the majority of the galaxy population, including the majority of the quiescent population at M∗ > 1010M⊙. Thus, the mechanism(s) responsible for quenching central galaxies play a crucial role in galaxy evolution as whole. We combine a high-resolution cosmological N-body simulation with observed evolutionary trends of the "star formation main sequence," quiescent fraction, and stellar mass function at z < 1 to construct a model that statistically tracks the star formation histories and quenching of central galaxies. Comparing this model to the distribution of central galaxy star formation rates in a group catalog of the SDSS Data Release 7, we constrain the timescales over which physical processes cease star formation in central galaxies. Over the stellar mass range 109.5-1011M⊙ we infer quenching e-folding times that span 1.5-0.5 Gyr with more massive central galaxies quenching faster. For M∗ = 1010.5M⊙, this implies a total migration time of ∼4 Gyr from the star formation main sequence to quiescence. Compared to satellites, central galaxies take ∼2 Gyr longer to quench their star formation, suggesting that different mechanisms are responsible for quenching centrals versus satellites. Finally, the central galaxy quenching timescale we infer provides key constraints for proposed star formation quenching mechanisms. Our timescale is generally consistent with gas depletion timescales predicted by quenching through strangulation. However, the exact physical mechanism(s) responsible for this remain unclear.

We measure how the properties of star-forming central galaxies correlate with large-scale environment, d, measured on 10 h-1 Mpc scales. We use galaxy group catalogues to isolate a robust sample of central galaxies with high purity and completeness. The galaxy properties we investigate are star formation rate (SFR), exponential disc scale length Rexp, and Sersic index of the galaxy light profile, nS. We find that, at all stellar masses, there is an inverse correlation between SFR and δ, meaning that above-average star-forming centrals live in underdense regions. For nS and Rexp, there is no correlation with δ at M* ≲ 1010.5 M⊙, but at higher masses there are positive correlations; a weak correlation with Rexp and a strong correlation with nS. These data are evidence of assembly bias within the star-forming population. The results for SFR are consistent with a model in which SFR correlates with present-day halo accretion rate, Mh. In thismodel, galaxies are assigned to haloes using the abundance-matching ansatz, which maps galaxy stellar mass onto halo mass. At fixed halo mass, SFR is then assigned to galaxies using the same approach, but Mh is used to map onto SFR. The best-fitting model requires some scatter in the Mh-SFR relation. The Rexp and nS measurements are consistent with a model in which both of these quantities are correlated with the spin parameter of the halo, λ. Halo spin does not correlate with δ at low halo masses, but for higher mass haloes, high-spin haloes live in higher density environments at fixed Mh. Put together with the earlier instalments of this series, these data demonstrate that quenching processes have limited correlation with halo formation history, but the growth of active galaxies, as well as other detailed galaxies properties, are influenced by the details of halo assembly.

Using group catalogues from the Sloan Digital Sky Survey (SDSS) Data Release 7, we measure galactic conformity in the local universe. We measure the quenched fraction of neighbour galaxies around isolated primary galaxies, dividing the isolated sample into star-forming and quiescent objects. We restrict our measurements to scales > 1 Mpc to probe the correlations between halo formation histories. Over the stellar mass range 109.7 ≥ M*/M⊙ = 1010.9, we find minimal evidence for conformity. We further compare these data to predictions of the halo age-matching model, in which the oldest galaxies are associated with the oldest haloes. For models with strong correlations between halo and stellar age, the conformity is too large to be consistent with the data. Weaker implementations of the age-matching model would not produce a detectable signal in SDSS data. We reproduce the results of Kauffmann et al., in which the star formation rates of neighbour galaxies are reduced around primary galaxies when the primaries are low star formers. However, we find this result is mainly driven by contamination in the isolation criterion; when removing the small fraction of satellite galaxies in the sample, the conformity signal largely goes away. Lastly, we show that small conformity signals, i.e. 2-5 per cent differences in the quenched fractions of neighbour galaxies, can be produced by mechanisms other than halo assembly bias. For example, if passive galaxies occupy more massive haloes than star-forming galaxies of the same stellar mass, a conformity signal that is consistent with recent measurements from PRIMUS (Berti et al.) can be produced.

We investigate the effects of galaxy environment on the evolution of the quiescent fraction (fQ) from z = 0.8 to 0.0 using spectroscopic redshifts and multi-wavelength imaging data from the PRIsm MUlti-object Survey (PRIMUS) and the Sloan Digital Sky Survey (SDSS). Our stellar mass limited galaxy sample consists of ∼14,000 PRIMUS galaxies within z = 0.2-0.8 and ∼64,000 SDSS galaxies within z = 0.05-0.12. We classify the galaxies as quiescent or star-forming (SF) based on an evolving specific star formation cut, and as low or high density environments based on fixed cylindrical aperture environment measurements on a volume-limited environment defining population. For quiescent and SF galaxies in low or high density environments, we examine the evolution of their stellar mass function (SMF). Then using the SMFs we compute fQ (M∗) and quantify its evolution within our redshift range. We find that the quiescent fraction is higher at higher masses and in denser environments. The quiescent fraction rises with cosmic time for all masses and environments. At a fiducial mass of 1010.5 M⊙, from z ∼ 0.7 to 0.1, the quiescent fraction rises by 15% at the lowest environments and by 25% at the highest environments we measure. These results suggest that for a minority of galaxies their cessation of star formation is due to external influences on them. In other words, in the recent universe a substantial fraction of the galaxies that cease forming stars do so due to internal processes.

The PRObabilistic Value-added Bright Galaxy Survey (PROVABGS) catalog will provide measurements of galaxy properties, such as stellar mass (M *), star formation rate (SFR), stellar metallicity (Z), and stellar age (t age), for >10 million galaxies of the Dark Energy Spectroscopic Instrument (DESI) Bright Galaxy Survey. Full posterior distributions of the galaxy properties will be inferred using state-of-the-art Bayesian spectral energy distribution (SED) modeling of DESI spectroscopy and Legacy Surveys photometry. In this work, we present the SED model, the neural emulator for the model, and the Bayesian inference framework of PROVABGS. Furthermore, we apply the PROVABGS SED modeling on realistic synthetic DESI spectra and photometry, constructed using the L-Galaxies semi-analytic model. We compare the inferred galaxy properties to the true values of the simulation using a hierarchical Bayesian framework to quantify accuracy and precision. Overall, we accurately infer the true M *, SFR, Z, and t age of the simulated galaxies. However, the priors on galaxy properties induced by the SED model have a significant impact on the posteriors, which we characterize in detail. This work also demonstrates that a joint analysis of spectra and photometry significantly improves the constraints on galaxy properties over photometry alone and is necessary to mitigate the impact of the priors. With the methodology presented and validated in this work, PROVABGS will maximize information extracted from DESI observations and extend current galaxy studies to new regimes and unlock cutting-edge probabilistic analyses. https://github.com/changhoonhahn/provabgs/

We present the steps taken to produce a reliable and complete input galaxy catalogue for the Dark Energy Spectroscopic Instrument (DESI) Bright Galaxy Survey (BGS) using the photometric Legacy Survey DR8 DECam. We analyse some of the main issues faced in the selection of targets for the DESI BGS, such as star-galaxy separation, contamination by fragmented stars and bright galaxies. Our pipeline utilizes a new way to select BGS galaxies using Gaia photometry and we implement geometrical and photometric masks that reduce the number of spurious objects. The resulting catalogue is cross-matched with the Galaxy And Mass Assembly (GAMA) survey to assess the completeness of the galaxy catalogue and the performance of the target selection. We also validate the clustering of the sources in our BGS catalogue by comparing with mock catalogues and the Sloan Digital Sky Survey (SDSS) data. Finally, the robustness of the BGS selection criteria is assessed by quantifying the dependence of the target galaxy density on imaging and other properties. The largest systematic correlation we find is a 7 per cent suppression of the target density in regions of high stellar density.

Abstract: The Dark Energy Spectroscopic Instrument (DESI) will execute a nearly magnitude-limited survey of low redshift galaxies (0.05 ≤ z ≤ 0.4, median z ≈ 0.2). Clustering analyses of this Bright Galaxy Survey (BGS) will yield the most precise measurements to date of baryon acoustic oscillations and redshift-space distortions at low redshift. DESI BGS will comprise two target classes: (i) BRIGHT (r < 19.5 mag), and (ii) FAINT (19.5 < r < 20 mag). Here we present a summary of the star-galaxy separation, and different photometric and geometrical masks, used in BGS to reduce the number of spurious targets. The selection results in a total density of ∼800 objects deg−2 for the BRIGHT and ∼600 objects deg−2 for the FAINT selections. A full characterization of the BGS selection can be found in Ruiz-Macias et al.

We present an extended validation of semi-analytical, semi-empirical covariance matrices for the two-point correlation function (2PCF) on simulated catalogs representative of luminous red galaxies (LRGs) data collected during the initial 2 months of operations of the Stage-IV ground-based Dark Energy Spectroscopic Instrument (DESI). We run the pipeline on multiple effective Zel'dovich (EZ) mock galaxy catalogs with the corresponding cuts applied and compare the results with the mock sample covariance to assess the accuracy and its fluctuations. We propose an extension of the previously developed formalism for catalogs processed with standard reconstruction algorithms. We consider methods for comparing covariance matrices in detail, highlighting their interpretation and statistical properties caused by sample variance, in particular, non-trivial expectation values of certain metrics even when the external covariance estimate is perfect. With improved mocks and validation techniques, we confirm a good agreement between our predictions and sample covariance. This allows one to generate covariance matrices for comparable data sets without the need to create numerous mock galaxy catalogs with matching clustering, only requiring 2PCF measurements from the data itself. The code used in this paper is publicly available at https://github.com/oliverphilcox/RascalC.

In this investigation, we leverage the combination of the Dark Energy Spectroscopic Instrument (DESI) Legacy Imaging Surveys Data Release 9, Survey Validation 3, and Year 1 data sets to estimate the conditional luminosity functions and conditional stellar mass functions (CLFs and CSMFs) of galaxies across various halo mass bins and redshift ranges. To support our analysis, we utilize a realistic DESI mock galaxy redshift survey (MGRS) generated from a high-resolution Jiutian simulation. An extended halo-based group finder is applied to both MGRS catalogs and DESI observation. By comparing the r- and z-band luminosity functions (LFs) and stellar mass functions (SMFs) derived using both photometric and spectroscopic data, we quantified the impact of photometric redshift (photo-z) errors on the galaxy LFs and SMFs, especially in the low-redshift bin at the low-luminosity/mass end. By conducting prior evaluations of the group finder using MGRS, we successfully obtain a set of CLF and CSMF measurements from observational data. We find that at low redshift, the faint-end slopes of CLFs and CSMFs below ∼109 h −2 L ⊙ (or h −2 M ⊙) evince a compelling concordance with the subhalo mass functions. After correcting the cosmic variance effect of our local Universe following Chen et al., the faint-end slopes of the LFs/SMFs turn out to also be in good agreement with the slope of the halo mass function.