Baryon Oscillation Spectroscopic Survey CMASS Research Articles

Galaxy clustering analysis with SimBIG and the wavelet scattering transform

The non-Gaussian spatial distribution of galaxies traces the large-scale structure of the Universe and therefore constitutes a prime observable to constrain cosmological parameters. We conduct Bayesian inference of the ΛCDM parameters Ωm, Ωb, h, ns, and σ8 from the Baryon Oscillation Spectroscopic Survey CMASS galaxy sample by combining the wavelet scattering transform (WST) with a simulation-based inference approach enabled by the SimBIG forward model. We design a set of reduced WST statistics that leverage symmetries of redshift-space data. Posterior distributions are estimated with a conditional normalizing flow trained on 20,000 simulated SimBIG galaxy catalogs with survey realism. We assess the accuracy of the posterior estimates using simulation-based calibration and quantify generalization and robustness to the change of forward model using a suite of 2000 test simulations. When probing scales down to kmax=0.5 h/Mpc, we are able to derive accurate posterior estimates that are robust to the change of forward model for all parameters, except σ8. We mitigate the robustness issues with σ8 by removing the WST coefficients that probe scales smaller than k∼0.3 h/Mpc. Applied to the Baryon Oscillation Spectroscopic Survey CMASS sample, our WST analysis yields seemingly improved constraints obtained from a standard perturbation-theory-based power spectrum analysis with kmax=0.25 h/Mpc for all parameters except h. However, we still raise concerns on these results. The observational predictions significantly vary across different normalizing flow architectures, which we interpret as a form of model misspecification. This highlights a key challenge for forward modeling approaches when using summary statistics that are sensitive to detailed model-specific or observational imprints on galaxy clustering. Published by the American Physical Society 2024

Evidence for baryon acoustic oscillations from galaxy–ellipticity correlations

The baryon acoustic oscillation feature in the clustering of galaxies or quasars provides a ‘standard ruler’ for distance measurements in cosmology. In this work, we report a 2–3σ signal of the baryon acoustic oscillation dip feature in the galaxy density–ellipticity cross-correlation functions using the spectroscopic sample of the Baryon Oscillation Spectroscopic Survey CMASS, combined with the deep Dark Energy Spectroscopic Instrument Legacy Imaging Surveys for precise galaxy shape measurements. We measure the galaxy–ellipticity correlation functions and model them using the linear alignment model. We constrain the distance DV/rd to redshift 0.57 to a precision of 3–5%, depending on the details of modelling. The galaxy–ellipticity measurement reduces the uncertainty of distance measurement by ~10% on top of that derived from the galaxy–galaxy correlation. More importantly, for future large and deep galaxy surveys, the independent galaxy–ellipticity measurements can help sort out the systematics in the baryon acoustic oscillation studies.

Atacama Cosmology Telescope: Modeling the gas thermodynamics in BOSS CMASS galaxies from kinematic and thermal Sunyaev-Zel’dovich measurements

The thermal and kinematic Sunyaev-Zel'dovich effects (tSZ, kSZ) probe the thermodynamic properties of the circumgalactic and intracluster medium (CGM and ICM) of galaxies, groups, and clusters, since they are proportional, respectively, to the integrated electron pressure and momentum along the line-of-sight. We present constraints on the gas thermodynamics of CMASS galaxies in the Baryon Oscillation Spectroscopic Survey (BOSS) using new measurements of the kSZ and tSZ signals obtained in a companion paper. Combining kSZ and tSZ measurements, we measure within our model the amplitude of energy injection $\epsilon M_\star c^2$, where $M_\star$ is the stellar mass, to be $\epsilon=(40\pm9)\times10^{-6}$, and the amplitude of the non-thermal pressure profile to be $\alpha_{\rm Nth}<0.2$ (2$\sigma$), indicating that less than 20% of the total pressure within the virial radius is due to a non-thermal component. We estimate the effects of including baryons in the modeling of weak-lensing galaxy cross-correlation measurements using the best-fit density profile from the kSZ measurement. Our estimate reduces the difference between the original theoretical model and the weak-lensing galaxy cross-correlation measurements in arXiv:1611.08606 by half but does not fully reconcile it. Comparing the tSZ measurements to cosmological simulations, we find that simulations underestimate the CGM pressure at large radii while they fare better in comparison with the kSZ measurements. This suggests that the energy injected via feedback models in the simulations that we compared against does not sufficiently heat the gas at these radii. We do not find significant disagreement at smaller radii. These measurements provide novel tests of current and future simulations. This work demonstrates the power of joint, high signal-to-noise kSZ and tSZ observations, upon which future cross-correlation studies will improve.

Can assembly bias explain the lensing amplitude of the BOSS CMASS sample in a Planck cosmology?

ABSTRACT In this paper, we investigate whether galaxy assembly bias can reconcile the 20–40 ${{\ \rm per\ cent}}$ disagreement between the observed galaxy projected clustering signal and the galaxy–galaxy lensing signal in the Baryon Oscillation Spectroscopic Survey CMASS galaxy sample. We use the suite of abacuscosmos lambda cold dark matter simulations at Planck best-fitting cosmology and two flexible implementations of extended halo occupation distribution (HOD) models that incorporate galaxy assembly bias to build forward models and produce joint fits of the observed galaxy clustering signal and the galaxy–galaxy lensing signal. We find that our models using the standard HODs without any assembly bias generalizations continue to show a 20–40 ${{\ \rm per\ cent}}$ overprediction of the observed galaxy–galaxy lensing signal. We find that our implementations of galaxy assembly bias do not reconcile the two measurements at Planck best-fitting cosmology. In fact, despite incorporating galaxy assembly bias, the satellite distribution parameter, and the satellite velocity bias parameter into our extended HOD model, our fits still strongly suggest a $\sim \! 34{{\ \rm per\ cent}}$ discrepancy between the observed projected clustering and galaxy–galaxy lensing measurements. It remains to be seen whether a combination of other galaxy assembly bias models, alternative cosmological parameters, or baryonic effects can explain the amplitude difference between the two signals.

The environmental dependence of the baryon acoustic peak in the Baryon Oscillation Spectroscopic Survey CMASS sample

The environmental dependence of galaxy clustering encodes information about the physical processes governing the growth of cosmic structure. We analyze the baryon acoustic peak as a function of environment in the galaxy correlation function of the Baryon Oscillation Spectroscopic Survey CMASS sample. Dividing the sample into three subsets by smoothed local overdensity, we detect acoustic peaks in the six separate auto-correlation and cross-correlation functions of the sub-samples. Fitting models to these correlation functions, calibrated by mock galaxy and dark matter catalogues, we find that the inferred distance scale is independent of environment, and consistent with the result of analyzing the combined sample. The shape of the baryon acoustic feature, and the accuracy of density-field reconstruction in the Zeldovich approximation, varies with environment. By up-weighting underdense regions and down-weighting overdense regions in their contribution to the full-sample correlation function, by up to 50%, we achieve a fractional improvement of a few per cent in the precision of baryon acoustic oscillation fits to the CMASS data and mock catalogues: the scatter in the preferred-scale fits to the ensemble of mocks improves from 1.45% to 1.34% (pre-reconstruction) and 1.03% to 1.00% (post-reconstruction). These results are consistent with the notion that the acoustic peak is sharper in underdense environments.

THE EVOLUTION OF POST-STARBURST GALAXIES FROM z ∼ 1 TO THE PRESENT

ABSTRACT Post-starburst galaxies are in the transitional stage between blue, star-forming galaxies and red, quiescent galaxies and therefore hold important clues for our understanding of galaxy evolution. In this paper, we systematically searched for and identified a large sample of post-starburst galaxies from the spectroscopic data set of the Sloan Digital Sky Survey (SDSS) Data Release 9. In total, we found more than 6000 objects with redshifts between z ∼ 0.05 and z ∼ 1.3, making this the largest sample of post-starburst galaxies in the literature. We calculated the luminosity function of the post-starburst galaxies using two uniformly selected subsamples: the SDSS main galaxy sample and the Baryon Oscillation Spectroscopic Survey CMASS sample. The luminosity functions are reasonably fit by half-Gaussian functions. The peak magnitudes shift as a function of redshift from M ∼ −23.5 at z ∼ 0.8 to M ∼ −20.3 at z ∼ 0.1. This is consistent with the downsizing trend, whereby more massive galaxies form earlier than low-mass galaxies. We compared the mass of the post-starburst stellar population found in our sample to the decline of the global star formation rate and found that only a small amount (∼1%) of all star formation quenching in the redshift range z = 0.2–0.7 results in post-starburst galaxies in the luminosity range our sample is sensitive to. Therefore, luminous post-starburst galaxies are not the place where most of the decline in the star formation rate of the universe is happening.

TheXMMCluster Survey: the halo occupation number of BOSS galaxies in X-ray clusters

We present a direct measurement of the mean halo occupation distribution (HOD) of galaxies taken from the eleventh data release (DR11) of the Sloan Digital Sky Survey-III Baryon Oscillation Spectroscopic Survey (BOSS). The HOD of BOSS low-redshift (LOWZ: $0.2 < z < 0.4$) and Constant-Mass (CMASS: $0.43 <z <0.7$) galaxies is inferred via their association with the dark-matter halos of 174 X-ray-selected galaxy clusters drawn from the XMM Cluster Survey (XCS). Halo masses are determined for each galaxy cluster based on X-ray temperature measurements, and range between ${\rm log_{10}} (M_{180}/M_{\odot}) = 13-15$. Our directly measured HODs are consistent with the HOD-model fits inferred via the galaxy-clustering analyses of Parejko et al. for the BOSS LOWZ sample and White et al. for the BOSS CMASS sample. Under the simplifying assumption that the other parameters that describe the HOD hold the values measured by these authors, we have determined a best-fit alpha-index of 0.91$\pm$0.08 and $1.27^{+0.03}_{-0.04}$ for the CMASS and LOWZ HOD, respectively. These alpha-index values are consistent with those measured by White et al. and Parejko et al. In summary, our study provides independent support for the HOD models assumed during the development of the BOSS mock-galaxy catalogues that have subsequently been used to derive BOSS cosmological constraints.

The clustering of galaxies in the SDSS-III Baryon Oscillation Spectroscopic Survey: single-probe measurements from CMASS anisotropic galaxy clustering

With the largest spectroscopic galaxy survey volume drawn from the SDSS-III Baryon Oscillation Spectroscopic Survey (BOSS), we can extract cosmological constraints from the measurements of redshift and geometric distortions at quasi-linear scales (e.g. above 50 $h^{-1}$Mpc). We analyze the broad-range shape of the monopole and quadrupole correlation functions of the BOSS Data Release 12 (DR12) CMASS galaxy sample, at the effective redshift $z=0.59$, to obtain constraints on the Hubble expansion rate $H(z)$, the angular-diameter distance $D_A(z)$, the normalized growth rate $f(z)\sigma_8(z)$, and the physical matter density $\Omega_mh^2$. We obtain robust measurements by including a polynomial as the model for the systematic errors, and find it works very well against the systematic effects, e.g., ones induced by stars and seeing. We provide accurate measurements $\{D_A(0.59)r_{s,fid}/r_s$ $\rm Mpc$, $H(0.59)r_s/r_{s,fid}$ $km s^{-1} Mpc^{-1}$, $f(0.59)\sigma_8(0.59)$, $\Omega_m h^2\}$ = $\{1427\pm26$, $97.3\pm3.3$, $0.488 \pm 0.060$, $0.135\pm0.016\}$, where $r_s$ is the comoving sound horizon at the drag epoch and $r_{s,fid}=147.66$ Mpc is the sound scale of the fiducial cosmology used in this study. The parameters which are not well constrained by our galaxy clustering analysis are marginalized over with wide flat priors. Since no priors from other data sets, e.g., cosmic microwave background (CMB), are adopted and no dark energy models are assumed, our results from BOSS CMASS galaxy clustering alone may be combined with other data sets, i.e., CMB, SNe, lensing or other galaxy clustering data to constrain the parameters of a given cosmological model. The uncertainty on the dark energy equation of state parameter, $w$, from CMB+CMASS is about 8 per cent. The uncertainty on the curvature fraction, $\Omega_k$, is 0.3 per cent. We do not find deviation from flat $\Lambda$CDM.

Modeling the reconstructed BAO in Fourier space

The density field reconstruction technique, which was developed to partially reverse the nonlinear degradation of the Baryon Acoustic Oscillation (BAO) feature in the galaxy redshift surveys, has been successful in substantially improving the cosmology constraints from recent galaxy surveys such as Baryon Oscillation Spectroscopic Survey (BOSS). We estimate the efficiency of the reconstruction method as a function of various reconstruction details. To directly quantify the BAO information in nonlinear density fields before and after reconstruction, we calculate the cross-correlations (i.e., propagators) of the pre(post)-reconstructed density field with the initial linear field using a mock galaxy sample that is designed to mimic the clustering of the BOSS CMASS galaxies. The results directly provide the BAO damping as a function of wavenumber that can be implemented into the Fisher matrix analysis. We focus on investigating the dependence of the propagator on a choice of smoothing filters and on two major different conventions of the redshift-space density field reconstruction that have been used in literature. By estimating the BAO signal-to-noise for each case, we predict constraints on the angular diameter distance and Hubble parameter using the Fisher matrix analysis. We thus determine an optimal Gaussian smoothing filter scale for the signal-to-noise level of the BOSS CMASS. We also present appropriate BAO fitting models for different reconstruction methods based on the first and second order Lagrangian perturbation theory in Fourier space. Using the mock data, we show that the modified BAO fitting model can substantially improve the accuracy of the BAO position in the best fits as well as the goodness of the fits.

Dive in the cosmic web: voids with Delaunay triangulation from discrete matter tracer distributions

We present a novel parameter-free cosmological void finder (dive, Delaunay TrIangulation Void findEr) based on Delaunay Triangulation (DT), which efficiently computes the empty spheres constrained by a discrete set of tracers. We define the spheres as DT voids, and describe their properties, including a universal density profile together with an intrinsic scatter. We apply this technique on 100 halo catalogues with volumes of 2.5 h−1Gpc side each, with a bias and number density similar to the Baryon Oscillation Spectroscopic Survey CMASS luminous red galaxies, performed with the patchy code. Our results show that there are two main species of DT voids, which can be characterized by the radius: they have different responses to halo redshift space distortions, to number density of tracers, and reside in different dark matter environments. Based on dynamical arguments using the tidal field tensor, we demonstrate that large DT voids are hosted in expanding regions, whereas the haloes used to construct them reside in collapsing ones. Our approach is therefore able to efficiently determine the troughs of the density field from galaxy surveys, and can be used to study their clustering. We further study the power spectra of DT voids, and find that the bias of the two populations are different, demonstrating that the small DT voids are essentially tracers of groups of haloes.

Photometric redshifts and model spectral energy distributions of galaxies from the SDSS-III BOSS DR10 data

We construct a set of model spectra specifically designed to match the colours of the Baryon Oscillation Spectroscopic Survey CMASS galaxies and to be used with photometric redshift template fitting techniques. As a basis, we use a set of spectral energy distributions(SEDs) of single and composite stellar population models. These models cannot describe well the whole colour range populated by the CMASS galaxies at all redshifts, wherefore we modify them by multiplying the SEDs with λ −β for λ > λ i for different values of λ i and β. When fitting these SEDs to the colours of the CMASS sample, with a burst and dust components in superposition, we can recreate the location in colour spaces inhabited by the CMASS galaxies. From the best-fitting models, we select a small subset in a two-dimensional plane, whereto the galaxies were mapped by a self-organizing map. These models are used for the estimation of photometric redshifts with a Bayesian template fitting code. The photometric redshifts with the novel templates have a very small outlier rate of 0.22 per cent, a low bias 〈Δz/(1 + z)〉 = 2.0 × 10 −3 , and scatter of σ 68 = 0.026 in the rest frame. Using our models, the galaxy colours are reproduced to a better extent with the photometric redshifts of this work than with photometric redshifts of Sloan Digital Sky Survey.

The information content of anisotropic Baryon Acoustic Oscillation scale measurements

Anisotropic measurements of the Baryon Acoustic Oscillation (BAO) feature within a galaxy survey enable joint inference about the Hubble parameter H(z) and angular diameter distance DA(z). These measurements are typically obtained from moments of the measured 2-point clustering statistics, with respect to the cosine of the angle to the line of sight mu. The position of the BAO features in each moment depends on a combination of DA(z) and H(z), and measuring the positions in two or more moments breaks this parameter degeneracy. We derive analytic formulae for the parameter combinations measured from moments given by Legendre polynomials, power laws and top-hat Wedges in μ, showing explicitly what is being measured by each in real-space for both the correlation function and power spectrum, and in redshift-space for the power spectrum. The large volume covered by modern galaxy samples means that the correlation function can be well approximated as having no correlations at different mu on the BAO scale, and that the errors on this scale are approximately independent of mu. Using these approximations, we derive the information content of various moments. We show that measurements made using either the monopole and quadrupole, or the monopole and μ2 power-law moment, are optimal for anisotropic BAO measurements, in that they contain all of the available information using two moments, the minimal number required to measure both H(z) and DA(z). We test our predictions using 600 mock galaxy samples, matched to the SDSS-III Baryon Oscillation Spectroscopic Survey CMASS sample, finding a good match to our analytic predictions. Our results should enable the optimal extraction of information from future galaxy surveys such as eBOSS, DESI and Euclid.

The clustering of galaxies in the SDSS-III Baryon Oscillation Spectroscopic Survey: measuring H(z) and DA(z) at z = 0.57 with clustering wedges

We analyze the 2D correlation function of the SDSS-III Baryon Oscillation Spectroscopic Survey (BOSS) CMASS sample of massive galaxies of the ninth data release to measure cosmic expansion H and the angular diameter distance D_A at a mean redshift of <z> = 0.57. We apply, for the first time, a new correlation function technique called clustering wedges. Using a physically motivated model, the anisotropic baryonic acoustic feature in the galaxy sample is detected at a significance level of 4.7 sigma compared to a featureless model. The baryonic acoustic feature is used to obtain model independent constraints cz/H/r_s = 12.28 +- 0.82 (6.7 per-cent accuracy) and D_A/r_s = 9.05 +- 0.27 (3.0 per-cent) with a correlation coefficient of -0.5, where r_s is the sound horizon scale at the end of the baryonic drag era. We conduct thorough tests on the data and 600 simulated realizations, finding robustness of the results regardless of the details of the analysis method. Combining with r_s constraints from the Cosmic Microwave Background we obtain H(0.57) = 90.8 +- 6.2 kms-1Mpc-1 and D_A(0.57) = 1386 +- 45 Mpc. We use simulations to forecast results of the final BOSS CMASS data set. We apply the reconstruction technique on the simulations demonstrating that the sharpening of the anisotropic baryonic acoustic feature should improve the detection as well as tighten constraints of H and D_A by 30 per-cent on average.

The clustering of galaxies in the SDSS-III Baryon Oscillation Spectroscopic Survey: single-probe measurements and the strong power of f(z) 8(z) on constraining dark energy

We present measurements of the anisotropic galaxy clustering from the Data Release 9 (DR9) CMASS sample of the SDSS-III Baryon Oscillation Spectroscopic Survey (BOSS). We analyze the broad-range shape of the monopole and quadrupole correlation functions to obtain constraints, at the effective redshift $z=0.57$ of the sample, on the Hubble expansion rate $H(z)$, the angular-diameter distance $D_A(z)$, the normalized growth rate $f(z)\sigma_8(z)$, the physical matter density $\Omega_m h^2$, and the biased amplitude of matter fluctuation b\sigma_8(z). We obtain {$H(0.57)$, $D_A(0.57)$, $f(0.57)\sigma_8(0.57)$, $\Omega_m h^2$, $b\sigma_8(0.57)$} = {$87.6_{-6.8}^{+6.7}$, $1396\pm73$, $0.126_{-0.010}^{+0.008}$, $1.19\pm0.14$, $0.428\pm0.066$} and their covariance matrix as well. The parameters which are not well constrained by our of galaxy clustering analysis are marginalized over with wide flat priors. Since no priors from other data sets (i.e., CMB) are adopted and no dark energy models are assumed, our results from BOSS CMASS galaxy clustering alone may be combined with other data sets, i.e. CMB, SNe, lensing or other galaxy clustering data to constrain the parameters of a given cosmological model. We show that the major power on constraining dark energy from the anisotropic galaxy clustering signal, as compared to the angular-averaged one (monopole), arises from including the normalized growth rate $f(z)\sigma_8(z)$. In the case of the wCDM cosmological model our single-probe CMASS constraints, combined with CMB (WMAP9+SPT), yield a value for the dark energy equation of state parameter of $w=-0.90\pm0.11$. Therefore, it is important to include $f(z)\sigma_8(z)$ while investigating the nature of dark energy with current and upcoming large-scale galaxy surveys.

The clustering of galaxies in the SDSS-III Baryon Oscillation Spectroscopic Survey: the low-redshift sample

We report on the small scale (0.5<r<40h^-1 Mpc) clustering of 78895 massive (M*~10^11.3M_sun) galaxies at 0.2<z<0.4 from the first two years of data from the Baryon Oscillation Spectroscopic Survey (BOSS), to be released as part of SDSS Data Release 9 (DR9). We describe the sample selection, basic properties of the galaxies, and caveats for working with the data. We calculate the real- and redshift-space two-point correlation functions of these galaxies, fit these measurements using Halo Occupation Distribution (HOD) modeling within dark matter cosmological simulations, and estimate the errors using mock catalogs. These galaxies lie in massive halos, with a mean halo mass of 5.2x10^13 h^-1 M_sun, a large scale bias of ~2.0, and a satellite fraction of 12+/-2%. Thus, these galaxies occupy halos with average masses in between those of the higher redshift BOSS CMASS sample and the original SDSS I/II LRG sample.