Search

In this paper, we set out to construct a set of reference mock galaxy redshift surveys (MGRSs) for the future Chinese Space-station Survey Telescope (CSST) observation, where subsequent survey selection effects can be added and evaluated. This set of MGRSs is generated using the dark matter subhalos extracted from a high-resolution Jiutian $N$-body simulation of the standard $\Lambda$CDM cosmogony with $\Omega_m=0.3111$, $\Omega_{\Lambda}=0.6889$, and $\sigma_8=0.8102$. The simulation has a boxsize of $1~h^{-1} {\rm Gpc}$, and consists of $6144^3$ particles with mass resolution $3.723 \times 10^{8} h^{-1} M_\odot $. In order to take into account the effect of redshift evolution, we first use all 128 snapshots in the Jiutian simulation to generate a light-cone halo/subhalo catalog. Next, galaxy luminosities are assigned to the main and subhalo populations using the subhalo abundance matching (SHAM) method with the DESI $z$-band luminosity functions at different redshifts. Multi-band photometries, as well as images, are then assigned to each mock galaxy using a 3-dimensional parameter space nearest neighbor sampling of the DESI LS observational galaxies and groups. Finally, the CSST and DESI LS survey geometry and magnitude limit cuts are applied to generate the required MGRSs. As we have checked, this set of MGRSs can generally reproduce the observed galaxy luminosity/mass functions within 0.1 dex for galaxies with $L > 10^8 L_\odot$ (or $M_* > 10^{8.5} M_\odot$) and within 1-$\sigma$ level for galaxies with $L < 10^8L_\odot$ (or $M_* < 10^{8.5} M_\odot$). Together with the CSST slitless spectra and redshifts for our DESI LS seed galaxies that are under construction, we will set out to test various slitless observational selection effects in subsequent probes., Comment: 13 pages, 9 figures, accepted for publication in MNRAS

The small-scale linear information in galaxy samples typically lost during non-linear growth can be restored to a certain level by the density field reconstruction, which has been demonstrated for improving the precision of the baryon acoustic oscillations (BAO) measurements. As proposed in the literature, a joint analysis of the power spectrum before and after the reconstruction enables an efficient extraction of information carried by high-order statistics. However, the statistics of the post-reconstruction density field are difficult to model. In this work, we circumvent this issue by developing an accurate emulator for the pre-reconstructed, post-reconstructed, and cross power spectra ($P_{\rm pre}$, $P_{\rm post}$, $P_{\rm cross}$) up to $k=0.5~h~{\rm Mpc^{-1}}$ based on the \textsc{Dark Quest} N-body simulations. The accuracy of the emulator is at percent level, namely, the error of the emulated monopole and quadrupole of the power spectra is less than $1\%$ and $10\%$ of the ground truth, respectively. A fit to an example power spectra using the emulator shows that the constraints on cosmological parameters get largely improved using $P_{\rm pre}$+$P_{\rm post}$+$P_{\rm cross}$ with $k_{\rm max}=0.25~h~{\rm Mpc^{-1}}$, compared to that derived from $P_{\rm pre}$ alone, namely, the constraints on ($\Omega_m$, $H_0$, $\sigma_8$) are tightened by $\sim41 \%-55\%$, and the uncertainties of the derived BAO and RSD parameters ($\alpha_{\perp}$, $\alpha_{||}$, $f\sigma_8$) shrink by $\sim 28\%-54\%$, respectively. This highlights the complementarity among $P_{\rm pre}$, $P_{\rm post}$ and $P_{\rm cross}$, which demonstrates the efficiency and practicability of a joint $P_{\rm pre}$, $P_{\rm post}$ and $P_{\rm cross}$ analysis for cosmological implications., Comment: 19 pages, 11 figures, 2 tables; accepted for publication in ApJ

We measure and analyze galaxy clustering and the dependence on luminosity, color, age, stellar mass and specific star formation rate using Baryon Oscillation Spectroscopic Survey (BOSS) galaxies at $0.48

We present the Aemulus $\nu$ simulations: a suite of 150 $(1.05 h^{-1}\rm Gpc)^3$ $N$-body simulations with a mass resolution of $3.51\times 10^{10} \frac{\Omega_{cb}}{0.3} ~ h^{-1} M_{\odot}$ in a $w\nu$CDM cosmological parameter space. The simulations have been explicitly designed to span a broad range in $\sigma_8$ to facilitate investigations of tension between large scale structure and cosmic microwave background cosmological probes. Neutrinos are treated as a second particle species to ensure accuracy to $0.5\, \rm eV$, the maximum neutrino mass that we have simulated. By employing Zel'dovich control variates, we increase the effective volume of our simulations by factors of $10-10^5$ depending on the statistic in question. As a first application of these simulations, we build new hybrid effective field theory and matter power spectrum surrogate models, demonstrating that they achieve $\le 1\%$ accuracy for $k\le 1\, h\,\rm Mpc^{-1}$ and $0\le z \le 3$, and $\le 2\%$ accuracy for $k\le 4\, h\,\rm Mpc^{-1}$ for the matter power spectrum. We publicly release the trained surrogate models, and estimates of the surrogate model errors in the hope that they will be broadly applicable to a range of cosmological analyses for many years to come., Comment: 37 pages, 15 figures, matching version accepted by JCAP

The source of the tension between local SN Ia based Hubble constant measurements and those from the CMB or BAO+BBN measurements is one of the most interesting unknowns of modern cosmology. Sample variance forms a key component of the error on the local measurements, and will dominate the error budget in the future as more SNe Ia are observed. Many methods have been proposed to estimate sample variance in many contexts, and we compared results from a number of approximate methods to estimates from N-body simulations in a previous paper, confirming that sample variance for the Pantheon SNe Ia sample does not solve the Hubble tension. We now extend this analysis to include the more accurate analytic method based on calculating correlations between the radial peculiar velocities of SNe Ia, comparing this technique with results from numerical simulations. We consider the dependence of these errors on the linear power spectrum and how non-linear velocities contribute to the error. Using this technique, and matching sample variance errors from more approximate methods, we can define an effective volume for SNe Ia samples, finding that the Pantheon sample is equivalent to a top-hat sphere of radius $\sim220~h^{-1}$Mpc. We use this link between sample-variance errors to compute $\Delta H_{0}$ for idealised surveys with particular angular distributions of SNe Ia. For example, a half-sky survey at the Pantheon depth has the potential to suppress the sample variance of $H_{0}$ to $\sim0.1$ km s$^{-1}$Mpc$^{-1}$, a significant improvement compared with the current result. Finally, we consider the strength of large-scale velocity power spectrum required to explain the Hubble tension using sample variance, finding it requires an extreme model well beyond that allowed by other observations., Comment: 12 pages, 6 figures. Analyses updated and extended, minor to moderate updates due to referees' comments, replaced to match the accepted version by PRD

Constraints on the linear growth rate, $f\sigma_8$, using small scale redshift space distortion measurements have a significant statistical advantage over those made on large scales. However, these measurements need to carefully disentangle the linear and non-linear information when interpreting redshift space distortions in terms of $f\sigma_8$. It is particularly important to do this given that some previous measurements found a significant deviation from the expectation based on the $\Lambda$CDM model constrained by Planck CMB data. We construct a new emulator-based model for small scale galaxy clustering with scaling parameters for both the linear and non-linear velocities of galaxies, allowing us to isolate the linear growth rate. We train the emulator using simulations from the AbacusCosmos suite, and apply it to data from the extended Baryon Oscillation Spectroscopic Survey (eBOSS) luminous red galaxy sample. We obtain a value of $f\sigma_8(z=0.737)=0.368\pm0.041$, in 2.3-$\sigma$ tension with the Planck 2018 $\Lambda$CDM expectation, and find less dependence on the minimum measurement scale than previous analyses., Comment: 14 pages, 9 figures, submitted to MNRAS

The High Latitude Spectroscopic Survey (HLSS) is the reference baseline spectroscopic survey for NASA's Nancy Grace Roman space telescope, measuring redshifts of $\sim 10$M H$\alpha$ emission line galaxies over a $2000$ deg$^2$ footprint at $z=1-2$. In this work, we use a realistic Roman galaxy mock catalogue to explore optimal phenomenological modeling of the measured power spectrum. We consider two methods for modeling the redshift-space distortions (Kaiser squashing and another with a window function on $\beta$ that selects out the coherent radial infall pairwise velocities, $M_A$ and $M_B$, respectively), two models for the nonlinear impact of baryons that smears the BAO signal (a fixed ratio between the smearing scales in the perpendicular and parallel dimensions and another where these smearing scales are kept as a free parameters, P$_{dw}(k|k_*)$ and P$_{dw}(k|\Sigma_\perp,\Sigma_\parallel)$, respectively), and two analytical emulations of nonlinear growth (one employing the halo model and another formulated from simulated galaxy clustering of a semi-analytical model, $F_{HM}$ and $F_{SAM}$, respectively). We find that the best model combination employing $F_{HM}$ is $P_{dw}(k|k_*)*F_{HM}*M_B$, while the best combination employing $F_{SAM}$ is $P_{dw}(k|k_*)*F_{SAM}*M_B$, which leads to unbiased measurements of cosmological parameters. We compare these to the Effective Field Theory of Large-Scale Structure perturbation theory model $P_{EFT}(k|\Theta)$, and find that our simple phenomenological models are comparable across the entire redshift range for $k_{max}=0.25$ and $0.3$ $h$/Mpc. We expect the tools that we have developed to be useful in probing dark energy and testing gravity using Roman in an accurate and robust manner., Comment: 17 pages, 3 figures, 4 tables, Accepted to MNRAS 23 May 2023

There is untapped cosmological information in galaxy redshift surveys in the non-linear regime. In this work, we use the AEMULUS suite of cosmological $N$-body simulations to construct Gaussian process emulators of galaxy clustering statistics at small scales ($0.1-50 \: h^{-1}\,\mathrm{Mpc}$) in order to constrain cosmological and galaxy bias parameters. In addition to standard statistics -- the projected correlation function $w_\mathrm{p}(r_\mathrm{p})$, the redshift-space monopole of the correlation function $\xi_0(s)$, and the quadrupole $\xi_2(s)$ -- we emulate statistics that include information about the local environment, namely the underdensity probability function $P_\mathrm{U}(s)$ and the density-marked correlation function $M(s)$. This extends the model of AEMULUS III for redshift-space distortions by including new statistics sensitive to galaxy assembly bias. In recovery tests, we find that the beyond-standard statistics significantly increase the constraining power on cosmological parameters of interest: including $P_\mathrm{U}(s)$ and $M(s)$ improves the precision of our constraints on $\Omega_m$ by 27%, $\sigma_8$ by 19%, and the growth of structure parameter, $f \sigma_8$, by 12% compared to standard statistics. We additionally find that scales below $\sim6 \: h^{-1}\,\mathrm{Mpc}$ contain as much information as larger scales. The density-sensitive statistics also contribute to constraining halo occupation distribution parameters and a flexible environment-dependent assembly bias model, which is important for extracting the small-scale cosmological information as well as understanding the galaxy-halo connection. This analysis demonstrates the potential of emulating beyond-standard clustering statistics at small scales to constrain the growth of structure as a test of cosmic acceleration., Comment: Published in the Astrophysical Journal; updated to match journal version

Recent local measurements of the Hubble constant made using supernovae have delivered a value that differs by $\sim$5$\sigma$ (statistical error) from predictions using the Cosmic Microwave Background (CMB), or using Baryon Acoustic Oscillations (BAO) and Big-Bang Nucleosynthesis (BBN) constraints, which are themselves consistent. The effective volume covered by the supernovae is small compared to the other probes, and it is therefore interesting to consider whether sample variance (often also called cosmic variance) is a significant contributor to the offset. We consider four ways of calculating the sample variance: (i) perturbation theory applied to the luminosity distance, which is the most common method considered in the literature; (ii) perturbation of cosmological parameters, as is commonly used to alleviate super-sample covariance in sets of N-body simulations; (iii) a new method based on the variance between perturbed spherical top-hat regions; (iv) using numerical N-body simulations. All give consistent results showing that, for the Pantheon supernova sample, sample variance can only lead to fluctuations in $H_0$ of order $\pm1$ km s$^{-1}$Mpc$^{-1}$ or less. While this is not in itself a new result, the agreement between the methods used adds to its robustness. Furthermore, it is instructive to see how the different methods fit together. We also investigate the internal variance of the $H_{0}$ measurement using SH0ES and Pantheon data. By searching for an offset between measurements in opposite hemispheres, we find that the direction coincident with the CMB dipole has a higher $H_{0}$ measurement than the opposite hemisphere by roughly 4 km s$^{-1}$Mpc$^{-1}$. We compare this with a large number of simulations and find that the size of this asymmetry is statistically likely, but the preference of direction may indicate that further calibration is needed., Comment: 13 pages, 5 figures; comments welcome

We analyze clustering measurements of BOSS galaxies using a simulation-based emulator of two-point statistics. We focus on the monopole and quadrupole of the redshift-space correlation function, and the projected correlation function, at scales of $0.1\sim60~h^{-1}$Mpc. Although our simulations are based on $w$CDM with general relativity (GR), we include a scaling parameter of the halo velocity field, $\gamma_f$, defined as the amplitude of the halo velocity field relative to the GR prediction. We divide the BOSS data into three redshift bins. After marginalizing over other cosmological parameters, galaxy bias parameters, and the velocity scaling parameter, we find $f\sigma_{8}(z=0.25) = 0.413\pm0.031$, $f\sigma_{8}(z=0.4) = 0.470\pm0.026$ and $f\sigma_{8}(z=0.55) = 0.396\pm0.022$. Compared with Planck observations using a flat $\Lambda$CDM model, our results are lower by $1.9\sigma$, $0.3\sigma$ and $3.4\sigma$ respectively. These results are consistent with other recent simulation-based results at non-linear scales, including weak lensing measurements of BOSS LOWZ galaxies, two-point clustering of eBOSS LRGs, and an independent clustering analysis of BOSS LOWZ. All these results are generally consistent with a combination of $\gamma_f^{1/2}\sigma_8\approx 0.75$. We note, however, that the BOSS data is well fit assuming GR, i.e. $\gamma_f=1$. We cannot rule out an unknown systematic error in the galaxy bias model at non-linear scales, but near-future data and modeling will enhance our understanding of the galaxy--halo connection, and provide a strong test of new physics beyond the standard model., Comment: 30 pages, 23 figures and 4 tables; ApJ accepted, updated to the match the accepted version

ISCEA (Infrared Satellite for Cosmic Evolution Astrophysics) is a small astrophysics mission whose Science Goal is to discover how galaxies evolved in the cosmic web of dark matter at cosmic noon. Its Science Objective is to determine the history of star formation and its quenching in galaxies as a function of local density and stellar mass when the Universe was 3-5 Gyrs old (1.2 10 Mpc in each of 50 protocluster (cluster and cosmic web) fields at 1.2 < z < 2.1. ISCEA measures the star formation quenching factor in those fields, and galaxy kinematics with a precision < 50 km/s to deduce the 3D spatial distribution in each field. ISCEA will transform our understanding of galaxy evolution at cosmic noon. ISCEA is a small satellite observatory with a 30cm equivalent diameter aperture telescope with a FoV of 0.32 deg^2, and a multi-object spectrograph with a digital micro-mirror device (DMD) as its programmable slit mask. ISCEA will obtain spectra of 1000 galaxies simultaneously at an effective resolving power of R=1000, with 2.8"x2.8" slits, over the NIR wavelength range of 1.1 to 2.0 microns, a regime not accessible from the ground without large gaps in coverage. ISCEA will achieve a pointing accuracy of <= 2" FWHM over 200s. ISCEA will be launched into a Low Earth Orbit, with a prime mission of 2.5 years. ISCEA's space-qualification of DMDs opens a new window for spectroscopy from space, enabling revolutionary advances in astrophysics., Comment: 40 pages, 31 figures

The Nancy Grace Roman Space Telescope will conduct a High Latitude Spectroscopic Survey (HLSS) over a large volume at high redshift, using the near-IR grism (1.0-1.93 $\mu$m, $R=435-865$) and the 0.28 deg$^2$ wide field camera. We present a reference HLSS which maps 2000 deg$^2$ and achieves an emission line flux limit of 10$^{-16}$ erg/s/cm$^2$ at 6.5$\sigma$, requiring $\sim$0.6 yrs of observing time. We summarize the flowdown of the Roman science objectives to the science and technical requirements of the HLSS. We construct a mock redshift survey over the full HLSS volume by applying a semi-analytic galaxy formation model to a cosmological N-body simulation, and use this mock survey to create pixel-level simulations of 4 deg$^2$ of HLSS grism spectroscopy. We find that the reference HLSS would measure $\sim$ 10 million H$\alpha$ galaxy redshifts that densely map large scale structure at $z=1-2$ and 2 million [OIII] galaxy redshifts that sparsely map structures at $z=2-3$. We forecast the performance of this survey for measurements of the cosmic expansion history with baryon acoustic oscillations and the growth of large scale structure with redshift space distortions. We also study possible deviations from the reference design, and find that a deep HLSS at $f_{\rm line}>7\times10^{-17}$erg/s/cm$^2$ over 4000 deg$^2$ (requiring $\sim$1.5 yrs of observing time) provides the most compelling stand-alone constraints on dark energy from Roman alone. This provides a useful reference for future optimizations. The reference survey, simulated data sets, and forecasts presented here will inform community decisions on the final scope and design of the Roman HLSS., Comment: 29 pages, 8 figures, ApJ submitted

We simulate the emission line galaxy properties for the Nancy Grace Roman Space Telescope by creating a 4 deg$^{2}$ galaxy mock catalog using Galacticus, a semi-analytic model (SAM) for galaxy formation and evolution. The simulated galaxy properties include emission line luminosity, equivalent width (EW), broad band photometry, and spectral energy distribution (SED). We compare this mock catalog with the latest observational data from WISP and MOSDEF. We find that our Galacticus model makes predictions consistent with observational data, over a wide redshift range for multiple emission lines. We use this validated galaxy mock catalog to forecast the photometric completeness of H$\alpha$ and [OIII] emission lines as a function of line flux cut and H band magnitdue cut, for both Roman and Euclid-like surveys. Our prediction for the photometric completeness of a Euclid-like survey is in good agreement with previous work based on WISP data. Our predictions for the photometric completeness of possible implementations of the Roman High Latitude Wide Area Spectroscopic Survey (HLWASS) provides key input for the planning of such a survey. We find that at H = 24, a Euclid-like survey to the line flux limit of 2$\times 10^{-16}\,$erg/s/cm$^2$ is 97\% complete, while a Roman HLWASS to the line flux limit of $10^{-16}\,$erg/s/cm$^2$ is only 94.6\% complete (it becomes 98\% complete at H = 25)., Comment: 12 pages, 12 figures; comments welcome

We measure the small-scale clustering of the Data Release 16 extended Baryon Oscillation Spectroscopic Survey Luminous Red Galaxy sample, corrected for fibre-collisions using Pairwise Inverse Probability weights, which give unbiased clustering measurements on all scales. We fit to the monopole and quadrupole moments and to the projected correlation function over the separation range $7-60\,h^{-1}$Mpc with a model based on the Aemulus cosmological emulator to measure the growth rate of cosmic structure, parameterized by $f\sigma_8$. We obtain a measurement of $f\sigma_8(z=0.737)=0.408\pm0.038$, which is $1.4\sigma$ lower than the value expected from 2018 Planck data for a flat $\Lambda$CDM model, and is more consistent with recent weak-lensing measurements. The level of precision achieved is 1.7 times better than more standard measurements made using only the large-scale modes of the same sample. We also fit to the data using the full range of scales $0.1-60\,h^{-1}$Mpc modelled by the Aemulus cosmological emulator and find a $4.5\sigma$ tension in the amplitude of the halo velocity field with the Planck+$\Lambda$CDM model, driven by a mismatch on the non-linear scales. This may not be cosmological in origin, and could be due to a breakdown in the Halo Occupation Distribution model used in the emulator. Finally, we perform a robust analysis of possible sources of systematics, including the effects of redshift uncertainty and incompleteness due to target selection that were not included in previous analyses fitting to clustering measurements on small scales., Comment: 19 pages, 16 figures, published in MNRAS

We present measurements of the linear galaxy bias of H$\alpha$ and [OIII] emission line galaxies (ELGs) for the High Latitude Spectroscopic Survey (HLSS) of Nancy Grace Roman Space Telescope, using galaxy mocks constructed using semi-analytical model for galaxy formation, {\it Galacticus}, with a large cosmic volume and redshift coverage. We compute the two-point statistics of galaxies in configuration space and measure linear bias within scales of $10\sim50 h^{-1}$Mpc. We adopt different selection algorithms to investigate the impact of the Roman line flux cut, as well as the effect of dust model used to calibrate Galacticus, on the bias measurements. We consider galaxies with H$\alpha$ and [OIII] emissions over the redshift range $1

We present a realistic 2000 deg$^{2}$ H$\alpha$ galaxy mock catalog with $1

The weak lensing magnification of Type Ia Supernovae (SNe Ia) is sensitive to the clustering of matter, and provides an independent cosmological probe complementary to SN Ia distance measurements. The Nancy Grace Roman Space Telescope is uniquely sensitive to this measurement as it can discover high redshift SNe Ia and measure them with high precision. We present a methodology for reconstructing the probability distribution of the weak lensing magnification $\mu$ of SNe Ia, $p(\mu)$, from observational data, and using it to constrain cosmological parameters. We find that the reconstructed $p(\mu)$ can be fitted accurately by a stretched Gaussian distribution, and used to measure the variance of $\mu$, $\xi_\mu$, which can be compared to theoretical predictions in a likelihood analysis. Applying our methodology to a set of realistically simulated SNe Ia expected from the Roman Space Telescope, we find that using the weak lensing magnification of the SNe Ia constrains a combination of matter density $\Omega_m$ and matter clustering amplitude $\sigma_8$. SN Ia distances alone lead to a better than 1\% measurement of $\Omega_m$. The combination of SN Ia weak lensing magnification and distance measurements result in a $\sim$ 10\% measurement on $\sigma_8$. The SNe Ia from Roman will be powerful in constraining the cosmological model., Comment: 9 pages, 8 figures; comments welcome

As a physical and sufficient compression of the full CMB data, the CMB distance priors, or shift parameters, have been widely used and provide a convenient way to include CMB data when obtaining cosmological constraints. In this paper, we revisit this data vector and examine its stability under different cosmological models. We find that the CMB distance priors are an accurate substitute for the full CMB data when probing dark energy dynamics. This is true when the primordial power spectrum model is directly generalized from the power spectrum of the model used in the derivation of the distance priors from the CMB data. We discover a difference when a non-flat model with the untilted primordial inflation power spectrum is used to measure the distance priors. This power spectrum is a radical change from the more conventional tilted primordial power spectrum and violates fundamental assumptions for the reliability of the CMB shift parameters. We also investigate the performance of CMB distance priors when the sum of neutrino masses $\sum m_{\nu}$ and the effective number of relativistic species $N_{\text{eff}}$ are allowed to vary. Our findings are consistent with earlier results: the neutrino parameters can change the measurement of the sound horizon from CMB data, and thus the CMB distance priors. We find that when the neutrino model is allowed to vary, the cold dark matter density $\omega_{c}$ and $N_{\text{eff}}$ need to be included in the set of parameters that summarize CMB data, in order to reproduce the constraints from the full CMB data. We present an updated and expanded set of CMB distance priors which can reproduce constraints from the full CMB data within $1\sigma$, and are applicable to models with massive neutrinos, as well as non-standard cosmologies., Comment: 11 pages, 8 figures, 1 table; updated to match the version accepted for publication in JCAP

Models of the spatial distribution of dark matter halos must achieve new levels of precision and accuracy in order to satisfy the requirements of upcoming experiments. In this work, we present a halo bias emulator for modeling the clustering of halos on large scales. It incorporates the cosmological dependence of the bias beyond the mapping of halo mass to peak height. The emulator makes substantial improvements in accuracy compared to the widely used Tinker et al. (2010) model. Halos in this work are defined using an overdensity criteria of 200 relative to the mean background density. Halo catalogs are produced for 40 N-body simulations as part of the Aemulus project at snapshots from z=3 to z=0. The emulator is trained over the mass range $6\times10^{12}-7\times10^{15}\ h^{-1}M_{\odot}$. Using an additional suite of 35 simulations, we determine that the precision of the emulator is redshift dependent, achieving sub-percent levels for a majority of the redshift range. Two additional simulation suites are used to test the ability of the emulator to extrapolate to higher and lower masses. Our high-resolution simulation suite is used to develop an extrapolation scheme in which the emulator asymptotes to the Tinker et al. (2010) model at low mass, achieving ~3% accuracy down to $10^{11}\ h^{-1}M_{\odot}$. Finally, we present a method to propagate emulator modeling uncertainty into an error budget. Our emulator is made publicly available at \url{https://github.com/AemulusProject/bias_emulator}., Comment: 14 pages, 7 figures, Submitted to ApJ

We perform a simulation with Galacticus, a semi-analytical galaxy formation model, to predict the number counts of H$\alpha$ and [OIII] emitting galaxies. With a state-of-the-art N-body simulation, UNIT, we first calibrate Galacticus with the current observation of H$\alpha$ luminosity function. The resulting model coupled with a dust attenuation model, can reproduce the current observations, including the H$\alpha$ luminosity function from HiZELS and number density from WISP. We extrapolate the model prediction to higher redshift and the result is found to be consistent with previous investigations. We then use the same galaxy formation model to predict the number counts for [OIII] emitting galaxies. The result provides further validation of our galaxy formation model and dust model. We present number counts of H$\alpha$ and [OIII] emission line galaxies for three different line flux limits: $5\times10^{-17}$erg/s/cm$^{2}$, $1\times10^{-16}$erg/s/cm$^{2}$ (6.5$\sigma$ nominal depth for WFIRST GRS), and $2\times10^{-16}$erg/s/cm$^{2}$ (3.5$\sigma$ depth of Euclid GRS). At redshift $2

Weak lensing of Type Ia supernovae (SNe Ia) is a systematic uncertainty in the use of SNe Ia as standard candles, as well as an independent cosmological probe, if the corresponding magnification distribution can be extracted from data. We study the peak brightness distribution of SNe Ia in the Pantheon sample, and find that the high $z$ sub-sample shows distinct weak lensing signatures compared to the low $z$ subsample: a long tail at the bright end due to high magnifications and a shift of the peak brightness toward the faint end, consistent with findings from earlier work. We have developed a technique to reconstruct the weak lensing magnification distribution of SNe Ia, $p(\mu)$, from the measured SN Ia flux distribution, and applied it to the Pantheon sample. We find that $p(\mu)$ can be reconstructed at a significance better than 2$\sigma$ for the subsample of SNe Ia at $z>0.7$ (124 SNe Ia), and at a lower significance for the SNe Ia at $z>0.9$ (49 SNe Ia), due to the small number of SNe Ia at high redshifts. The large number of $z>1$ SNe Ia from future surveys will enable the use of $p(\mu)$ reconstructed from SNe Ia as an independent cosmological probe., Comment: 7 pages, 5 figures, comments welcome

We present a systematic analysis of the cosmological constraints from the "Pantheon Sample" of 1048 Type Ia Supernovae (SNe Ia) in the redshift range $0.01

Using the N-body simulations of the AEMULUS Project, we construct an emulator for the non-linear clustering of galaxies in real and redshift space. We construct our model of galaxy bias using the halo occupation framework, accounting for possible velocity bias. The model includes 15 parameters, including both cosmological and galaxy bias parameters. We demonstrate that our emulator achieves ~ 1% precision at the scales of interest, 0.1

The rapidly growing statistical precision of galaxy surveys has lead to a need for ever-more precise predictions of the observables used to constrain cosmological and galaxy formation models. The primary avenue through which such predictions will be obtained is suites of numerical simulations. These simulations must span the relevant model parameter spaces, be large enough to obtain the precision demanded by upcoming data, and be thoroughly validated in order to ensure accuracy. In this paper we present one such suite of simulations, forming the basis for the AEMULUS Project, a collaboration devoted to precision emulation of galaxy survey observables. We have run a set of 75 (1.05 h^-1 Gpc)^3 simulations with mass resolution and force softening of 3.51\times 10^10 (Omega_m / 0.3) ~ h^-1 M_sun and 20 ~ h^-1 kpc respectively in 47 different wCDM cosmologies spanning the range of parameter space allowed by the combination of recent Cosmic Microwave Background, Baryon Acoustic Oscillation and Type Ia Supernovae results. We present convergence tests of several observables including spherical overdensity halo mass functions, galaxy projected correlation functions, galaxy clustering in redshift space, and matter and halo correlation functions and power spectra. We show that these statistics are converged to 1% (2%) for halos with more than 500 (200) particles respectively and scales of r>200 ~ h^-1 kpc in real space or k ~ 3 h Mpc^-1 in harmonic space for z\le 1. We find that the dominant source of uncertainty comes from varying the particle loading of the simulations. This leads to large systematic errors for statistics using halos with fewer than 200 particles and scales smaller than k ~ 4 h^-1 Mpc. We provide the halo catalogs and snapshots detailed in this work to the community at https://AemulusProject.github.io., Comment: 16 pages, 12 figures, 3 Tables Project website: https://aemulusproject.github.io/

Existing models for the dependence of the halo mass function on cosmological parameters will become a limiting source of systematic uncertainty for cluster cosmology in the near future. We present a halo mass function emulator and demonstrate improved accuracy relative to state-of-the-art analytic models. In this work, mass is defined using an overdensity criteria of 200 relative to the mean background density. Our emulator is constructed from the AEMULUS simulations, a suite of 40 N-body simulations with snapshots from z=3 to z=0. These simulations cover the flat wCDM parameter space allowed by recent Cosmic Microwave Background, Baryon Acoustic Oscillation and Type Ia Supernovae results, varying the parameters w, Omega_m, Omega_b, sigma_8, N_{eff}, n_s, and H_0. We validate our emulator using five realizations of seven different cosmologies, for a total of 35 test simulations. These test simulations were not used in constructing the emulator, and were run with fully independent initial conditions. We use our test simulations to characterize the modeling uncertainty of the emulator, and introduce a novel way of marginalizing over the associated systematic uncertainty. We confirm non-universality in our halo mass function emulator as a function of both cosmological parameters and redshift. Our emulator achieves better than 1% precision over much of the relevant parameter space, and we demonstrate that the systematic uncertainty in our emulator will remain a negligible source of error for cluster abundance studies through at least the LSST Year 1 data set., Comment: https://aemulusproject.github.io/

We investigate the possibility of detecting in redshift surveys a hemispherical power asymmetry similar to that first reported in CMB observations. We assume the hemispherical asymmetry arises from a linear gradient in comoving coordinates in the perturbation amplitude. We predict the resulting clustering of galaxy or galaxy cluster tracers using an excursion set approach; doing so accounts for the variation of both the underlying clustering and the tracer bias. Based on the predicted variation of the clustering of tracers, we perform a Fisher matrix forecast of the galaxy clustering amplitude and calculate the statistical significance for ideal surveys and planned surveys. The results indicate that the DESI galaxy survey would be able to detect this signal with higher than $3\sigma$ significance if the asymmetry does exist. We also investigate the amplitude and scale dependence of the above result. The DESI galaxy survey can probe the dipole amplitude higher than 0.04, which correspond to a $\pm4\%$ difference of the temperature fluctuation along and opposite the dipole direction, at least at the $2\sigma$ level. Additionally, we investigate a modulation of the power spectrum that exhibits asymmetry only for large scales. This modulation is potentially detectable. For Milky Way galaxy mass tracers, the scale-dependent modulation yields a larger change in the large scale power spectrum than does a scale-independent modulation, because the former does not alter the bias., Comment: 8 pages, 7 figures, updated to match the ApJ accepted version

We compare a large suite of theoretical cosmological models to observational data from the cosmic microwave background, baryon acoustic oscillation measurements of expansion, Type Ia SNe measurements of expansion, redshift space distortion measurements of the growth of structure, and the local Hubble constant. Our theoretical models include parametrizations of dark energy as well as physical models of dark energy and modified gravity. We determine the constraints on the model parameters, incorporating the redshift space distortion data directly in the analysis. To determine whether models can be ruled out, we evaluate the $p$ value (the probability under the model of obtaining data as bad or worse than the observed data). In our comparison, we find the well known tension of H$_0$ with the other data; no model resolves this tension successfully. Among the models we consider, the large scale growth of structure data does not affect the modified gravity models as a category particularly differently than dark energy models; it matters for some modified gravity models but not others, and the same is true for dark energy models. We compute predicted observables for each model under current observational constraints, and identify models for which future observational constraints will be particularly informative., Comment: 35 pages, 33 figures, updated to match the ApJ accepted version

We present the first scientific results from the luminous red galaxy sample (LRG) of the extended Baryon Oscillation Spectroscopic Survey (eBOSS). We measure the small and intermediate scale clustering from a sample of more than 61,000 galaxies in the redshift range $0.6 < z < 0.9$. We interpret these measurements in the framework of the Halo Occupation Distribution. The bias of eBOSS LRGs is $2.30 \pm 0.03$, with a satellite fraction of $13\pm3$\% and a mean halo mass of $2.5\times10^{13}h^{-1}M_{\odot}$. These results are consistent with expectations, demonstrating that eBOSS galaxies will be reliable tracers of large scale structure at $z\sim 0.7$. The eBOSS galaxy bias implies a scatter of luminosity at fixed halo mass, $\sigma_{\log L}$, of 0.19 dex. Using the clustering of massive galaxies from BOSS-CMASS, BOSS-LOWZ, and SDSS, we find that $\sigma_{\log L}=0.19$ is consistent with observations over the full redshift range that these samples cover. The addition of eBOSS to previous surveys allows investigation of the evolution of massive galaxies over the past $\sim 7$ Gyr., Comment: 15 pages, 16 figures; changes following referee's comments; matches the accepted version

We present cosmological results from the final galaxy clustering data set of the Baryon Oscillation Spectroscopic Survey, part of the Sloan Digital Sky Survey III. Our combined galaxy sample comprises 1.2 million massive galaxies over an effective area of 9329 deg^2 and volume of 18.7 Gpc^3, divided into three partially overlapping redshift slices centred at effective redshifts 0.38, 0.51, and 0.61. We measure the angular diameter distance DM and Hubble parameter H from the baryon acoustic oscillation (BAO) method after applying reconstruction to reduce non-linear effects on the BAO feature. Using the anisotropic clustering of the pre-reconstruction density field, we measure the product DM*H from the Alcock-Paczynski (AP) effect and the growth of structure, quantified by f{\sigma}8(z), from redshift-space distortions (RSD). We combine measurements presented in seven companion papers into a set of consensus values and likelihoods, obtaining constraints that are tighter and more robust than those from any one method. Combined with Planck 2015 cosmic microwave background measurements, our distance scale measurements simultaneously imply curvature {\Omega}_K =0.0003+/-0.0026 and a dark energy equation of state parameter w = -1.01+/-0.06, in strong affirmation of the spatially flat cold dark matter model with a cosmological constant ({\Lambda}CDM). Our RSD measurements of f{\sigma}_8, at 6 per cent precision, are similarly consistent with this model. When combined with supernova Ia data, we find H0 = 67.3+/-1.0 km/s/Mpc even for our most general dark energy model, in tension with some direct measurements. Adding extra relativistic species as a degree of freedom loosens the constraint only slightly, to H0 = 67.8+/-1.2 km/s/Mpc. Assuming flat {\Lambda}CDM we find {\Omega}_m = 0.310+/-0.005 and H0 = 67.6+/-0.5 km/s/Mpc, and we find a 95% upper limit of 0.16 eV/c^2 on the neutrino mass sum., Comment: 38 pages, 20 figures, 12 tables. Submitted to MNRAS. Clustering data and likelihoods are made available at https://sdss3.org/science/boss_publications.php

The Extended Baryon Oscillation Spectroscopic Survey (eBOSS) will conduct novel cosmological observations using the BOSS spectrograph at Apache Point Observatory. Observations will be simultaneous with the Time Domain Spectroscopic Survey (TDSS) designed for variability studies and the Spectroscopic Identification of eROSITA Sources (SPIDERS) program designed for studies of X-ray sources. eBOSS will use four different tracers to measure the distance-redshift relation with baryon acoustic oscillations (BAO). Using more than 250,000 new, spectroscopically confirmed luminous red galaxies at a median redshift z=0.72, we project that eBOSS will yield measurements of $d_A(z)$ to an accuracy of 1.2% and measurements of H(z) to 2.1% when combined with the z>0.6 sample of BOSS galaxies. With ~195,000 new emission line galaxy redshifts, we expect BAO measurements of $d_A(z)$ to an accuracy of 3.1% and H(z) to 4.7% at an effective redshift of z= 0.87. A sample of more than 500,000 spectroscopically-confirmed quasars will provide the first BAO distance measurements over the redshift range 0.92.1; these new data will enhance the precision of $d_A(z)$ and H(z) by a factor of 1.44 relative to BOSS. Furthermore, eBOSS will provide improved tests of General Relativity on cosmological scales through redshift-space distortion measurements, improved tests for non-Gaussianity in the primordial density field, and new constraints on the summed mass of all neutrino species. Here, we provide an overview of the cosmological goals, spectroscopic target sample, demonstration of spectral quality from early data, and projected cosmological constraints from eBOSS., Comment: 39 pages, 13 figures, accepted AJ, revised for consistency with accepted version