Momentum power spectrum of SDSS galaxies by massE cosmic ruler: 2.1x improvement in measure of growth rate

Peculiar motion of galaxies probes the structure growth in the Universe. In this study we employ the galaxy stellar mass-binding energy (massE) relation with only two nuisance parameters to build the largest peculiar-velocity (PV) catalog to date, consisting of 229,890 ellipticals from the main galaxy sample (MGS) of the Sloan Digital Sky Survey (SDSS). We quantify the distribution of the massE-based distances in individual narrow redshift bins (dz=0.005), and then estimate the PV of each galaxy based on its offset from the Gaussian mean of the distribution. As demonstrated with the Uchuu-SDSS mock data, the derived PV and momentum power spectra are insensitive to accurate calibration of the massE relation itself, enabling measurements out to a redshift of 0.2, well beyond the current limit of z=0.1 using other galaxy scaling laws. We then measure the momentum power spectrum and demonstrate that it remains almost unchanged if varying significantly the redshift bin size within which the distance is measured, as well as the intercept and slope of the massE relation, respectively. By fitting the spectra using the perturbation theory model with four free parameters, f{\sigma}8 is constrained to f{\sigma}8 =0.459+0.068-0.069 over {\Delta}z=0.02-0.2, 0.416+0.074-0.076 over {\Delta}z=0.02-0.1 and 0.526+0.133-0.143 over {\Delta}z=0.1-0.2. The error of f{\sigma}8 is 2.1 times smaller than that by the redshift space distortion (RSD) of the same sample. A Fisher-matrix forecast illustrates that the constraint on f{\sigma}8 from the massE-based PV can potentially exceed that from the stage-IV RSD in late universe (z<0.5).


INTRODUCTION
The large-scale structure of galaxies, as a consequence of the structure growth of matter in the expanding Universe, is a powerful way to understand the cosmos (Peebles 1980).The average excess number of galaxies surrounding a galaxy as a function of separation, i.e., the correlation function of the galaxy number density or its Fourier analog-power spectrum, has been widely used to characterize the large scale structure (Kaiser 1987;Eisenstein et al. 2005;Beutler et al. 2012;Alam et al. 2017Alam et al. , 2021)).This includes studying the baryonic acoustic oscillation around a separation of 100 Mpc/h, as well as examining the redshift space distortion and the full broadband correlation function.
In addition to the number density, galaxies experience peculiar motion under the gravitational potential of matter.Unlike the galaxy number density that is a biased tracer of ⋆ E-mail: yong@nju.edu.cnmatter, peculiar velocity (PV) is a more direct tracer at least on linear scales (Koda et al. 2014;Howlett 2019).When the density perturbation is small, the over-density evolves with cosmic time as: δ(r ,a)=ρ(r ,a)/ ρ(r ,a)-1 ∝ D(a) where r is a comoving coordinate, a is the scale factor and D(a) is the linear growth factor.The continuity equation then gives where the Hubble constant H(a) = da/dt/a and the growth rate f (a)= dlnD(a)/dlna.The equation indicates that on linear scales the PV divergence is directly related to the growth rate f (a) of matter.Under the General Relativity (GR), f (a) = Ωm(a) 0.55 (Lahav et al. 1991;Linder & Cahn 2007), and any deviation from it could imply the break down of GR.By taking the fact that the Fourier transform of a quantity is proportional to the transform of its divergence, F (∇•v ) =ik • F (v ).As a result, in addition to the density field, the PV field of a large scale structure offers a complementary way to constrain the structure growth or test GR by quantifying the velocity power spectrum.In practice, we only measure the PV along the line of sight to each observed galaxy, which thus gives a line-of-sight mass-weighted velocity power spectrum or momentum power spectrum.With a PV as denoted by vp, the observed redshift z obs follows where zcos is the cosmic redshift and c is the speed of light.Measurements of PV require measurements of observed redshift as well as cosmological-independent measurements of distances that infer zcos.The latter has been achieved through galaxy scaling law including the Tully-Fisher relation for spiral galaxies (Tully & Fisher 1977), the fundamental plane (FP) (Dressler et al. 1987) and surface brightness fluctuation for elliptical galaxies (Tonry & Schneider 1988;Tonry et al. 2001) as well as SN Ia (Riess et al. 1998).
Studies of momentum power spectrum has significant development both in theories and observations.Through distribution function and perturbation theory, non-linear effects on the momentum power spectrum has been characterized analytically (Vlah et al. 2012(Vlah et al. , 2013;;Okumura et al. 2014;Saito et al. 2014), which makes possible the comparison to observations over a full band of spatial scales.In observations, the technique of the density power spectrum (Feldman et al. 1994) has been developed (Koda et al. 2014;Hand et al. 2017Hand et al. , 2018;;Howlett 2019) and applied to several PV catalog to estimate the momentum power spectrum (e.g.Qin et al. 2019).These measurements, along with velocity correlation function or velocity-density cross correlations (Johnson et al. 2014;Adams & Blake 2017, 2020;Wang et al. 2021;Turner et al. 2023;Lai et al. 2023), have produced reasonable constraints on the growth rate.
Exiting PV surveys based on galaxy scaling law are limited to very low redshift, i.e., not beyond a redshift of 0.1, while those with SN Ia can go to higher redshift but its low number statistics make it not competitive.For example, the Two Micron All-Sky Survey (2MASS) Tully-Fisher survey covers a redshift range up to 0.033 with about 2000 objects (Masters et al. 2008).The Six-degree-Field Galaxy Survey peculiarvelocity (6dFGSv) is a FP-based PV survey up to a redshift of 0.053 with about 9000 objects (Magoulas et al. 2012;Springob et al. 2014).The SDSS FP-based peculiar velocity sample extends to a redshift of 0.1 with ∼ 3.4×10 4 objects (Howlett et al. 2022).Cosmicflows-IV data-set compiles a total number of ∼ 5.5×10 4 objects with all kinds of distance calibration (Tully et al. 2022), including those with surface brightness fluctuation and SN Ia etc.
Although the measurement error of PV increases with redshift, both the increasing volume at a larger redshift and the sampling on larger spatial scales can improve the constraint on the growth rate f (e.g.Koda et al. 2014).The Tully-Fisher relation requires spatially-resolved measurement of kinematic maps which is difficult beyond z=0.1.Although FP only requires single-fiber measurement of velocity dispersion, the FP relation contains three free parameters plus an additional nuisance parameter that quantifies the offset from the FP plane for galaxies with different properties (Magoulas et al. 2012;Howlett et al. 2022).
In Shi et al. (2021), the stellar mass of a galaxy has been found to tightly correlate with the galaxy binding energy within the effective radius, which is referred as the massE relation.This relation offers a new cosmic ruler with only two nuisance parameters (Shi et al. 2022), in contrast to the FP that has four nuisance parameters.The application of the massE cosmic ruler to ellipticals of the SDSS main galaxy sample (MGS) offer distance measurement with an accuracy as low as 0.35% for the redshift range of 0.05-0.2,which proves the existence of dark energy at 7-σ under flat-ΛCDM (Shi et al. 2022).In this study, we use the massE to measure PV of the SDSS MGS elliptical sample and obtain the momentum power spectrum out to z=0.2.Throughout the study, we adopt the Planck 2018 flat Λ-CDM cosmology (Planck Collaboration et al. 2020) with h=0.6736,Ωm=0.3153 and Ω b =0.0493 and σ8=0.811.

The massE cosmic ruler as the distance estimator
The distance measurement of the SDSS MGS elliptical sample with the massE cosmic ruler is detailed in Shi et al. (2022).Briefly, the massE ruler estimates the angular diameter distance of an elliptical galaxy with stellar mass well above 10 8 M⊙ as: where Re,as is the apparent effective radius in arcsec, σe is the velocity dispersion within Re,as, Mstar,1Mpc is the galaxy stellar mass placed at a distance of 1 Mpc and zcos is the cosmic redshift of the galaxy.Besides the above four observables, the two nuisance parameters in Equation 3 are β and SD 0 , where β is the slope of the massE relationship, and SD 0 is the intercept of the relation and is related to the absolute calibration of observables and is degenerate with the local Hubble constant H0.By fitting to the low-redshift sample compiled in Shi et al. (2021), Shi et al. (2022) estimated: and ln(D0/Mpc) = (26.06β− 5.54) ± 0.005. (5) The angular diameter distance shown in the above equation can be converted the luminosity distance through DL=DA(1+z) 2 .For a flat cosmology, the angular distance is related to the comoving distance through Dc=DA(1+z).
In practice, the observed redshift z obs is used to approximate zcos.Although the effect on the distance is negligible, the effect on the PV is noticeable as shown in Equation 10.By assuming zcos=z obs , the shift in the comoving distance can be derived as: in which we define κ=( 2β+0.25 2β−0.25 )1 c , and assume that vp is a small number as compared to the speed of light.
As shown later, the PV is not sensitive to the accurate calibration of the local massE relationship, i.e., β and SD 0 can vary significantly while the derived PV and momentum power spectrum remains almost unchanged.

The MGS-elliptical sample
The procedure to construct the PV catalog is shown in Table 1 that lists the number of objects for each step.We start with the SDSS Large-Scale-Structure (LSS) safe-0 catalog (Blanton et al. 2005) whose completeness and selection function have been quantified 1 .The redshift in the catalog is converted from the heliocentric frame to the CMB frame following Fixsen et al. (1996).
To define the elliptical sub-sample, we first cross-match with the MPA-JHU catalog with a search radius of 1.5 arcsec to query the stellar mass, Petrosion half-light radius in the r band (PETROR50_R), velocity dispersion (V_DISP) and the fraction of de Vaucouleurs component in the r band (FRACDEV_R) from the MPA-JHU catalog (Kauffmann et al. 2003;Brinchmann et al. 2004).The fraction of the successful match is 0.938, which will be included in the galaxy weight when calculating the momentum power spectrum.The SDSS MGS-elliptical sample is then defined as FRACDEV_R > 0.80, PETROR50_R > 0, along with V_DISP > 0 and 0.02 < z obs < 0.2.The lower redshift limit is adopted so that even in the lowest redshift bin of [0.02, 0.021], 95% of galaxies still have stellar masses larger than 3×10 8 M⊙ to ensure the validity of the distance calculation with Equation 3. The upper redshift limit is set to ensure a reasonable high galaxy number density.The footprint of the MGS-elliptical sample is shown in Figure 1

Measurements of peculiar velocities
To measure PV for the MGS-elliptical sample, we first derive quantities as needed by the distance calculator given by Equation 3. As detailed in Shi et al. (2022), the apparent size Re,as is derived from the Petrosian half-light radius PETROR50_R by adopting the cosmological model of the MPA-JHU catalog, and the flux-like Mstar,1Mpc is converted from the stellar mass in the MPA-JHU catalog using the same cosmological model.The velocity dispersion within the fiber aperture from that catalog is corrected to the effective radius to give σe.
In practice it is not recommended to estimate vp directly from Equation 2, because the log-normal distribution of the distance errors would make errors of vp non-Gaussian, thus complicating the error analysis of momentum power spectra and cosmological parameter constraints.In previous studies (Adams & Blake 2017), it has been demonstrated that vp can be estimated from the logarithmic ratio of the comoving  distance at z=z obs and that at z=zcos, i.e., η gal = ln( Dc(z obs ) Dc(zcos) ): where with requirements of vp/c ≪ 1 and αvp ≪ 1.The typical vp is about 300 km s −1 so that the first condition is satisfied.
The second condition requires z obs to be reasonably larger than 0, e.g.αvp ∼ 0.05 at z obs =0.02 and vp=300 km s −1 .
In the above equation, the true comoving distance of a galaxy, i.e., Dc(zcos), is measured directly from the massE cosmic ruler.In practice we replace zcos with z obs in the massE equation, which causes a systematic shift in ln(Dc(zcos)) as derived in Equation 6.To account for this shift, an updated equation to derive PV for the massE cosmic ruler is given as where κ (defined in Equation 6) dominates over α at increas- ing redshift and reduces PV errors accordingly compared to Equation 8.The remaining key is to obtain unbiased estimate of the comoving distance at the observed redshift of the galaxy, i.e., Dc(z obs ).Our overall strategy is to divide the interested redshift range into small bins and estimate Dc(z obs ) therein.The redshift bin size needs to balance two facts.On the one hand, galaxies in each redshift bin are composed of those located at slightly lower true distances with vp > 0 and those located at slightly larger distances with vp < 0. As a result, an increasing redshift bin size contains a larger number of galaxies to converge their mean true distance ⟨ln(Dc(zcos))⟩ to the distance at z obs , i.e., ln(Dc(z obs )).On the other hand, the redshift bin cannot be too large so that Dc(z obs ) is no longer a constant over the bin.
Our fiducial measurements adopt a redshift bin size of dz obs =5×10 −3 .In each redshift bin, as shown in Figure 2(a) for an example, the ln(Dc(zcos)) distribution is more or less normal.But because a small fraction of the sample has failed measurements in the velocity dispersion or galaxy sizes, there are outliers.We remove outliers in the following ways: for those within three times the standard deviation of the bestfit Gaussian, we include all of them; for those outside, we randomly include part of them so that the number of objects follows the Gaussian expectation at its ln(Dc,z obs ) bin.The new distribution is shown in Figure 2 (b).The exclusion of these outliers decreases the sample completeness and increases the error of the momentum power spectrum.The η bin is then defined as the minus of the offset of a galaxy from the Gaussian mean in that redshift bin.Because the redshift bin size is not infinitely small, we correct additional small offset between the observed redshift and the mean redshift of galaxies within the bin: where ⟨z⟩ is the mean redshift of galaxies within the given redshift bin.Besides the fiducial redshift bin size, we show that if varying dz obs from 1×10 −3 to 1×10 −2 , the result remains unchanged as discussed later.
Our final PV catalog, listed in Table 1 and illustrated in Figure 3, contains 229,890 galaxies, representing the largest number of objects in any PV catalog thus far.It is seven times larger than the previously largest one (Howlett et al. 2022) and even exceeds the expected number of PV objects generated by applying the FP plus Tully-Fisher methods to the DESI survey (Saulder et al. 2023).At z < 0.1 we include the whole SDSS MGS early-type galaxies with few selection, while the SDSS-FP sample has additional selections to have higher PV precision for individual galaxies, such as a brighter flux limit, galaxy inclination angle, Hα equivalent width etc.As a result, our sample at z <0.1 contains a factor of 2.5 more galaxies than the SDSS-FP one.Another important difference is because of our ability to apply the massE method up to a redshift of 0.2, which is well beyond the limit (∼0.1) of the Tully-Fisher and FP methods.The redshift limit of the latter two methods is not a rigid boundary based on the first principle, but rather derived from empirical usage.The Tully-Fisher method requires spatially-resolved kinematic measurements, limiting its redshift range.Although the FP can rely on single fiber measurements, its zero point is a function of galaxy properties (Magoulas et al. 2012;Howlett et al. 2022).While adjustments can be made to the FP zero point based on galaxy properties, even a small remaining systematic shift may lead to a large offset in the PV estimation at relatively high redshifts.On the hand, the massE contains solely two free parameters that correspond to the intercept and slope of the massE relationship, respectively.Consequently, we divide the redshift range into narrow individual bins, whereby PV is measured by the offset from the Gaussian mean of galaxies in narrow redshift bins, largely removing the necessity of accurately measuring the zero point throughout the entire redshift range.We will present the investigation using mock data-sets Three dotted lines mark the absolute luminosity of the magnitude cut (m r,limit =17.6) at three redshifts.(d), the distribution of the difference between the intrinsic distance and observed distance at z=0.2, which are derived from the intrinsic magnitude and observed magnitude, respectively.The filled histogram represents the one with intrinsic apparent magnitude below m r,limit , while open one for the observed apparent magnitude below m r,limit .(e), the fractional difference between intrinsic and observed distribution (the denominator is the observed distribution) as a function of the offset from the mean in terms of standard deviation of the distribution at z=0.2.The result contains ten simulations.The solid line is the best-fitted linear function with zero intercept.(f ), the slope of above linear function as a function of redshift, where the solid line is the best-fitted power law function.
to support the above statement in § 4.4.We do not see any apparent limit for massE to extend beyond a redshift of 0.2.However, given the rapid drop in the galaxy number density of the SDSS MGS toward higher redshift, in this study, we define our sample to be below a redshift of 0.2.For the momentum power spectrum measurement, we select those that have completeness fgotten ≥ 0.9 and lie within the contiguous region of the North Galactic cap.The corresponding footprint area is 6725 deg 2 as illustrated in Fig. 1(b).

The velocity zero point
In the above approach to calculate PV, we use the Gaussian mean of the distance distribution within a small redshift bin as the velocity zero point.This assumption is correct as along as the large scale structure within the redshift bin is not coherent so that no bulk motion exists.The best way to assure this is an all-sky survey that covers both northern and southern hemisphere.However the SDSS is not large enough especially for galaxies below z=0.1.Fortunately, the SDSS MGS have additional coverage along the equatorial plane besides a continuous area in the northern hemisphere.As a result, the Gaussian mean is measured only using galaxies in two great circles of the sky sphere, i.e., those along the equatorial plane plus a stripe along longitude as enclosed by the dashed boxes in Figure 1 (a).We demonstrate this with mock data in § 4.3.

The correction for the flux limit
The SDSS MGS is a flux-limited sample.Due to the luminosity function, sources that are intrinsically faint and below the flux limit have a higher likelihood of scattering to be above the limit, compared to sources that are intrinsically bright, which have a lower chance of falling below the limit.As a result, some fraction of sources that are above the flux limit should be intrinsically faint, thus overestimating its stellar mass and underestimating its distance.
The log distribution of r-band magnitude errors as a function of the magnitude is shown for the entire photometric catalog from mr=14 to 24 in Figure 5(a).We fit a double Gaussian function to the log distribution of magnitude errors in magnitude bins with ∆mr=0.4mag.An example of the fitting is illustrated in Figure 5(b) for mr=17.8-18.2.The SDSS MGS sample for the large scale structure has a flux limit around mr=17.6.Below the redshift of 0.1, this flux limit still probes those below the characteristic magnitude as shown in Figure 5(c).
We randomly draw 3×10 6 absolute magnitudes from the luminosity function for a given redshift.The distribution ranges from Mr=−24 to Mr=23.25-DM,where DM is the distance module, with 23.25 corresponding to the SDSS image 3-σ depth.Next, we assign the flux error to each absolute magnitude following Figure 5(a), and calculate the difference in the magnitude (∆Mr) between the intrinsic and observed ones.Since the mass-to-light ratio is a function of color but SDSS MGS is a purely flux-limited sample, this magnitude difference corresponds to a difference in distance such that ∆lnDc = ln(10) 2.5 We then obtain the ∆lnDc distribution for two separate samples: one with intrinsic apparent magnitudes below the SDSS MGS bright limit and the other with observed apparent magnitudes below the limit, as shown in Figure 5(d) for z=0.2.The difference between these two distributions arises from the fact that the latter is affected by the flux limit cut.
In Figure 5(e), we show the fractional difference (the denominator is the observed distribution) in individual bins of ∆lnDc/ std(lnDc).Here, std(lnDc) represents the standard deviation of ∆lnDc, and the offset is expressed in terms of standard deviations.This approach ensures that the above trend remains unchanged if lnDc experiences increasing Gaussian errors as in the real case.After conducting 10 simulations for a given redshift, we fit a linear function to the above trend with zero intercept.We obtain the slope of the trend for z=0.05, 0.1, 0.15 and 0.2, and then fit a power law to interpolate at any given redshift.To correct the bias introduced by the flux limit for each redshift bin in real case, we select a portion of galaxies according to its offset from the Gaussian mean and the linear function at given redshift, and reverse its sign of the offset.

The methodology
The line-of-sight mass-weighted momentum field is described as where δ(r ) is the galaxy over-density at the location vector r .For a catalog with discrete galaxies, it can be written as (Feldman et al. 1994;Hand et al. 2017;Howlett 2019): where w(r ) is the galaxy total weight that is ≥ 1, ng(r ) is the observed galaxy number density, and A is the normalization given in Equation 18 below.The momentum power spectrum P P (k )=⟨F P (k )F P (k ′ )⟩, the Fourier analog of the two-point correlation function 2 , can expand through Legendre polynomials: where µ = k • r h and r h = (r + r ′ )/2.The estimator for the P P ℓ (k) is given by (Feldman et al. 1994;Hand et al. 2017): where Ω k is the solid angle in Fourier space.By assuming the local plane-parallel approximation of ( k the integral over F P (r ) and F P (r ′ ) in the above equation can be separated into independent integrals (Yamamoto et al. 2006).By further decomposing the k •r into a product of spherical harmonics, the calculation of the above equation can be done with a small number of Fast Fourier Transforms (FFTs), which significantly saves the computing time (Hand et al. 2017).In this study, we are only interested in ℓ=0, i.e. the angle-averaged momentum power spectrum P P ℓ=0 (k), so that the local plane-parallel is exactly valid.

The implementation of the methodology
The python code nbodykit (Hand et al. 2018) implements the above calculation for the density power spectrum.We make some minor modification of their python codes FKPCatalog and ConvolvedFFTPower to produce FKPMomentumCatalog  and MomentumConvolvedFFTPower for the momentum power spectrum calculation.Here we list the key inputs in order to run the code: (1) w mpa−jhu -the completeness weight due to the failure match of the LSS catalog with the MPA-JHU catalog: the reciprocal of 0.938 as derived in Table 1 is taken as w mpa−jhu for all galaxies.
(2) wv p -the completeness weight for the vp measurement that is due to the failure in measurements of σe: we divide the RA, DEC and z-range of our sample into 10 bins, respectively.In each 3-D cell, the completeness is calculated as the ratio of the number of objects with available vp and the total number.The completeness is found to be almost a constant with the median value of 0.977 and a standard deviation of 0.029.As a result we assign wv p as a reciprocal of the median completeness for all galaxies.
(3) w sel -the weight for the selection function of the SDSS main galaxy LSS sample.This number is available in the catalog, which is 1/fgotten.
(4) The observed number density n(z) as a function of redshift: this function is shown in Figure 4(a).We interpolate the function to get the number density at the redshift of each galaxy.
(5) The ⟨v 2 p ⟩z as a function of redshift: this function is shown in Fig. 4(b).Similarly we interpolate it to obtain the velocity variance at the redshift of each galaxy.
(6) The FKP weight where P P is set to be 10 9 h −3 Mpc 3 km 2 s −2 , which is roughly the momentum power spectrum at our median spatial scale of 0.1 h/Mpc.(7) The range of spatial scales as bracketed by kmin and kmax: kmin is crucial for the measurement of structure growth that is sensitive to large spatial scales.We set it to be 0.02 h/Mpc.Our test with mock data shows that the momentum power spectrum is recovered well above kmin.For kmax, we test several values and found that the growth rate is insensitive to it.
(8) With the combined completeness weight w complete =w mpa−jhu wv p w sel and FKP weight wFKP, we ran FKPMomentumCatalog and MomentumConvolvedFFTPower implemented in nbodykit to get the momentum power spectrum which is further subtracted by the shot-noise contribution: where A is the normalization as given by Given wtot=w complete wFKP, the effective redshift of the sample is given by where i and j are members of all galaxy pairs in the sample.
For our whole sample, z eff is found to be 0.0959.

Mock galaxy catalogs and true momentum power spectrum
We use 32 Uchuu-SDSS mock galaxy catalogs that are constructed from the Uchuu 2.1 trillion N-body simulation (Dong-Páez et al. 2022).The catalogs match the depth, footprint, selection function etc of the SDSS MGS in the North Cap region.To construct the mock catalog, the r-band absolute magnitude Mr was assigned to a halo through subhalo abundance matching (SHAM) techniques.The g-r color was further assigned randomly through the probability distribution of the color-magnitude diagram.With Mr, g-r color and observed redshift, k-d tree was used to find the closest real galaxy for each mock galaxy, after which properties of real galaxies were attached to mock galaxies.For our use, we need one additional parameter -FRACDEV, which is obtained with k-d tree match also based on Mr, g-r color and observed redshift.
We define elliptical galaxies of each mock catalog following the same criteria as the observed one: (1) FRACDEV ≥ 0.8 and (2) 0.02 < z obs < 0.2.With true peculiar velocities obtained through Equation 2, the true momentum power spectrum is then calculated using nbodykit.The median of intrinsic momentum power spectra of mocks is presented as solid lines in Figure 8.

Mimicking PV measurements
To produce the PV for the mock galaxy catalog, we carry out the following simulation as outlined in Figure 6: (1) The apparent half-light radius in arcsec (R mock e,as ) for each mock galaxy: we first plot the observed physical halflight radius as a function of Mr of the SDSS MGS-elliptical sample in Figure 7(a).Here we include ellipticals with z obs =0.01-0.2 to extend the dynamic range.As shown in the figure, we fit the median and standard deviations of half-light radii at given Mr with second-order polynomial functions, respectively.We then assign the intrinsic physical size to each mock galaxy based on the above median trend along with its mock Mr. The associated errors are assigned based on standard deviations of trends (dashed lines in Figure 6 (a)).The mock physical size is then further converted to the angular size R mock e,as based on the mock true angular diameter distance.
(2) The stellar mass normalized at 1 Mpc (M mock star,1Mpc ) for each mock galaxy: we first plot the observed mass-to-light ratio in terms of logMstar+Mr/2.5 as a function of g-r color based on the SDSS MGS in Figure 7(b).Here we include all types of galaxies in order to extend the dynamic range.Again, second-order polynomial functions are fitted to the median and standard deviations of the relationship, respectively.We then assign an intrinsic stellar mass to each mock galaxy based on its mock g-r color and Mr using the median trend, while the observed errors are assigned based on the fitting to the standard deviations of the relationship (dashed lines in Figure 6 (b)).The stellar mass is then converted to that at 1 Mpc (M mock star,1Mpc ) based on the mock true luminosity distance.
(3) The velocity dispersion (σ mock e ) for each mock galaxy: we first use Equation 3 to calculate an intrinsic velocity dispersion from the intrinsic M mock star,1Mpc and R mock e,as along with the mock true distance, for which the slope β of the massE relationship is fixed to the best-fit value in Equation 4.About 10 4 objects of the PV sample have duplicate measurements of the velocity dispersion in the MPA-JHU catalog.As shown in Figure 6 (c), the distribution of the velocity dispersion errors can be described by two Gaussian profiles.We thus assign each galaxies following the error distribution in the figure .(4) The mock distance for each mock galaxy: based on the mock R mock e,as , M mock star,1Mpc and σ mock e with errors, we calculate the mock distance through Equation 3.Here we perturb the β and the corresponding D0 with their errors, respectively.
(5) The mock PV for each mock galaxy: we then follow the same procedure as the observed one to obtain the PV for each mock galaxy.For each of 31 mocks from mock-0 to mock-30, we run 60 simulations so that in total we have 1860 simulations for each observed set of PV measurement.Note that the last mock (mock-31) is used to test the zero velocity point.Each simulation includes random perturbation of the intrinsic R mock e,as , M mock star,1Mpc and σ mock e with their error distributions, and one perturbation of β and D0, respectively.
By construction, the mock catalog does not have failure measurements, so we do not remove any outliers.We confirm that the velocity errors produced by mocks are similar to the observed ones at different redshift.

The velocity zero point
In this section, we examine the impact of the velocity zero point on the momentum power spectrum.In our method to measure the PV, the zero velocity represents the Gaussian mean of the distance distribution within a narrow redshift bin.We first present the case where each mock uses its own distance distribution to calculate the Gaussian mean.As a result, the zero point is derived from galaxies in an about 7000 deg 2 area in the northern hemisphere.As shown in the first row of Figure 8, in the low redshift range of ∆z=0.02-0.10 the derived momentum power spectrum is systematically lower than the intrinsic one, while the one in ∆z=0.10-0.20 overall recovers the intrinsic one very well except for small k that is below 0.02 h/Mpc.This demonstrates that because of  .The first row: for each mock (mock 0 to mock 30), the zero velocity point (Gaussian mean) is determined from galaxies in that particular mock; the second row: the zero velocity point is determined from galaxies in that particular mock plus all galaxies in mock 31; the bottom row: for each mock, the zero velocity is determined from part of galaxies in mock 31 and part of galaxies in that particular mock to mimick galaxies in yellow dashed boxes in Figure 1.For details, see § 4.3.From left to right, the momentum power spectra are presented for three redshift ranges.The solid line is median of the intrinsic spectra of 31 mocks (mock 0 to mock 30), and symbols are the median of 31 mocks calculated with the massE-based PVs.
the limited volume at low redshift a small bulk motion still exists and affects our methodology.
The Uchuu-SDSS provides 32 independent mocks to simulate the SDSS MGS in the northern hemisphere.In the second case, we use the last mock (mock-31) to emulate a SDSS-MGS-like survey in the southern hemisphere.For each of the remaining mocks, we combine their galaxies with those from mock 31 to measure the Gaussian mean.This mean is then used to derive PV for the galaxies in that particular mock.Note that galaxies from mock 31 are solely used for determining the Gaussian mean and are not used in the calculation of the momentum power spectra.As shown in the middle panel of Figure 8, the intrinsic momentum power spectra for three different redshift ranges are well recovered.This indicates, even at redshifts below z=0.1, the effect of the velocity zero point on the momentum power spectrum is negligible if large enough survey is conducted.
Unfortunately, the SDSS MGS lacks a comparable survey in the southern hemisphere.But it does contain galaxies along the equatorial plane.In the third case, we thus combine part of galaxies from mock 31 with part of galaxies from each of the remaining mocks to calculate the Gaussian mean for that mock.For the part from mock 31, we select galaxies with a declination below 5 degrees, which roughly yields a similar sky area to the part of the SDSS MGS that lies in the equatorial plane but not in the continuous northern regions.This sky area is represented by the two orange boxes along the equatorial plane in Figure 1.For each of the remaining mocks, we select their galaxies along the longitude plane that is similar to the stripe in Figure 1.The above two parts of galaxies lie on two perpendicular great circles of the sky sphere, and thus effectively reduce the coherent bulk motion.As shown in the lower panel of Figure 8, in both the high and low redshift bins, the intrinsic momentum power spectra are recovered reasonably well.We thus use this strategy for both the observation data and mock data.

Comparisons between massE-based measurements and true ones in mock data-sets
To assess the validity of our methodology, we conduct two types of comparison between the massE-based measurements and the true ones in mock data-sets.The first one is the PV of individual galaxies, while the second focuses on the momentum power spectrum, including both the monopole and quadrupole.Figure 9 presents a comparison between the derived PV and the true velocity in three redshift bins for galaxies, where the derived PV is the mean value after 2000 simulations.The two velocities in all three redshift bins exhibit an almost one-to-one correlation.This suggests that, on a statistical level, the massE effectively recovers the true velocity.A small offset is observed in the zero velocity, which represents a random error on the zero velocity point.This offset adds to the shot noise (Howlett 2019), resulting in a slight increase in the momentum power spectrum above 0.1 h/Mpc, as seen in Figure 10.However, this should not affect the constraint on the growth rate, as it is not sensitive to the high k regime.
To further evaluate the velocity field, Figure 10 presents a comparison of the momentum power spectra between the massE-based measurements and true ones.The solid line indicates the mean value from 31 mocks, while the symbols represents the mean values of massE-based measurements for 31 mocks with each mock consisting of 60 simulations.The left panel of the figure displays the monopole of the momentum power spectrum for three different redshift ranges.The massE-based measurement well recovers the true momentum power spectrum except for a slight increment above k ∼ 0.1 h/Mpc as stated above.Although the survey volume of mock data is not large enough to measure the quadrupole of the momentum power spectrum, by examining the average of simulations it can give insights into whether the velocity field is recovered by the massE.Note that the quadrupole power spectrum is affected by the wide angle effect but for the SDSS MGS survey volume it is minimal compared to errors in the measurement (e.g.Castorina & White 2018).Additionally, when comparing the true quadrupole to the one derived from massE, both measurements should be affected by the wide-angle effect to the same extent.As shown in the right panel of the figure, across the whole redshift range (∆z=0.02-0.20), the massE-based measurement follows the true quadrupole momentum with no systematic shift.For two redshift sub-ranges, only at high k values exceeding 0.2 h/Mpc, the massE-based measurements appear to be higher than the true values, but note that the error is large too.
We perform additional tests, taking into account that the intrinsic massE (used to derive the velocity dispersion) does not show any redshift evolution but we mistakenly incorporate redshift evolution in the massE when calculating the PV (see Fig. 6 for two relations).The first case assumes the intercept in the second massE has a strong redshift evolution with a form of SD 0 (z)=SD 0 (1+z) 5 .As shown in Figure 11, the best fit to PVs of individual galaxies in three redshift bins still show good one-to-one correlations with slightly increase in errors.The upper panel of Figure 13 shows that both monopole and quadrupole of the momentum power spectrum over ∆z=0.02-0.2 is well recovered too.The second case assumes the slope in the second massE has a redshift evolution in a form of β(z)=β(1+z) 0.4 .At z=0.2 the evolution results in an approximately 20-σ offset from our best estimate of β as listed in Equation 4. Figure 12 shows that best fits to PVs of individual galaxies in three redshift bins also exhibit oneto-one correlations, but with larger errors.The lower panel of Figure 13 indicates the mean values of both the monopole and quadrupole momentum power spectra derived from massE still closely follow the intrinsic one.to the true median of 31 Uchuu-SDSS mocks, our observed momentum power spectrum is slightly higher at low k but the deviation becomes larger at k ≳ 0.07 h/Mpc.The difference could be the result that on small spatial scales nonlinear effects affect the power spectrum, while the Uchuu-SDSS simulation is purely N-body (Ishiyama et al. 2021).
In the figure, we overlay the result of TNG-300 hydrodynamics simulation that include baryons (Nelson et al. 2019).For this simulation, we use the snapshot at z=0 and include all sub-halos with stellar masses larger than 10 10 M⊙.The calculation of the momentum power spectrum is obtained through FFTPower implemented in nbodykit.Given a specific simulation box, we can have three lines of sight, each giving three independent measurements, as shown in the figure.The disparity between TNG and our observation above 0.07 h/Mpc is now much smaller.It remains unclear what causes the difference between our observations and Uchuu-SDSS mock data.The momentum power spectrum at high k is dominated by the product of the density and velocity (e.g., see Fig. 13 in Howlett ( 2019)).The Uchuu-SDSS mock is designed to have the same galaxy bias as the observed one, a fact we confirm by comparing the density power spectra of our observations and the Uchuu-SDSS mock.Accordingly the difference at high-k should not be attributed to different galaxy biases between the two.It is probable that the N-body simulation, such as Uchuu, does not account adequately for non-linear effects as compared with TNG hydrodynamic simulation.Fortunately, the constraint on the growth rate is not particularly sensitive to the maximum wavenumber kmax.
Figure 14(b) and (c) presents two additional momentum power spectra for two sub-redshift ranges of ∆z obs =0.10-0.20 and 0.02-0.10,respectively.Similar to the case of the full redshift range, the momentum power spectra within two subranges exhibit a slightly elevation when compared to the Uchuu-SDSS result at low k.However, they show considerably larger disparities above k=0.07h/Mpc, but tend to match more or less the result of the TNG hydrodynamic simulation at high k. Figure 13.The same as Figure 10 but only for the entire redshift range and the case that the intrinsic massE in the mock has no redshift evolution while the one used to derive PV has a strong redshift evolution in the intercept and slope, respectively.

Robustness of the measurement
In this section, we discuss the result if varying some parameters that we use to derive the momentum power spectrum including the massE relation itself and the redshift bin dz obs .In § 4.4, we confirm with mock data that the massE-based PV is insensitive to the accurate calibration of the massE relationship.Here we further demonstrate this with the observational data.In Figure 16(a), we overlay two momentum power spectra with the fiducial one: one is the case that the normalization SD 0 of the massE relation evolves strongly with redshift as SD 0 ∝ (1+z) 5 , and another is the one that the slope β of the relationship evolves strongly with the redshift as β ∝ (1+z) 0.4 .For the β, such a redshift evolution gives a value at z=0.2 that deviates from the observed one by 20-σ.In the first case, any alteration in the value of SD 0 causes all distances to shift uniformly with the same amount, thereby maintaining the offset of a distance from the Gaussian mean.As a result, the momentum power spectrum remains the same regardless of the magnitude of the SD 0 shift.In the second case, the momentum power spectrum changes a little bit but well within the error bar, which still benefits from our method.For galaxies within a narrow redshift bin, their physical properties including sizes, velocity dispersions and stellar masses exhibit a narrow range.If there is a change in the slope β, it only introduces a small amount of scatter in the distribution of distances.Therefore Figure 16(a) demonstrates that the momentum power spectrum with the massE-based PV is not sensitive to the accurate calibration of the massE relation.Figure 16(b) presents the momentum power spectra whose PV is measured with redshift bin sizes of dz obs =0.001 and dz obs =0.01, as compared to the fiducial one with dz obs =0.005.As shown in the figure, the three spectra are almost the same, indicating that the PV measurements is not sensitive to dz obs which can vary by a factor of 10.

The effective survey volume
Before detailed modeling, we first roughly evaluate the statistical power of our sample in constraining a cosmological parameter, given that the signal to noise of the parameter constrained at a given k is more or less proportional to the square root of the survey effective volume at that k (e.g.Koda et al. 2014).Similar to the definition of the density power spectrum (Feldman et al. 1994;Tegmark 1997), an effective volume for the momentum power spectrum can be written as where n(z) and ⟨v 2 p ⟩z are the galaxy number density and PV variance as a function of the redshift, respectively.The above equation indicates that for high galaxy number densities n(z) ≫ ⟨v 2 p ⟩z/P P (k), the effective volume is just the geometrical volume and thus increases with the cube of the redshift.On the other hand, for low galaxy number densities, because ⟨v 2 p ⟩z roughly scales with the square of the redshift at z < 0.2, the effective volume only scales linearly with the redshift.
To calculate the effective volume for our sample, we adopt the median of 31 Uchuu-SDSS mocks for P P (k).The n(z) and ⟨v 2 p ⟩z are shown in Figure 4 for our sample.We derived V eff (k) for the whole redshift range as well as two redshift sub-ranges as shown in Figure 17, along with 6dFGSv and SDSS-FP surveys.The 6dFGSv data has a spatial coverage of 17 000 deg 2 and a redshift range of 0.0001 to 0.0534 (Scrimgeour et al. 2016).We use (cz/(1 + z)ση) 2 to get its velocity variance by setting ση=0.324(Scrimgeour et al. 2016).The galaxy number density is almost a constant around 1.2×10 −3 h 3 Mpc −3 .We use the same P P (k) as for our sample.The SDSS FP-based PV sample covers a redshift range of 0.0033 to 0.1 and an area of 7016 deg 2 (Howlett et al. 2022).We use the above equation for the velocity variance with ση=0.23 (Howlett et al. 2022), measure n(z) from their catalog and also adopt the median of Uchuu-SDSS mocks for P P (k).
The integral quantity N k,modes = 0.07 0.02 k 2 V eff dk gives the total number of k modes on linear scales where the growth rate is sensitive to.As compared to the 6dFGSv and SDSS FP-based PV catalogs, our sample offers a factor of about 2.4 and 6.7 increase in the above quantity, respectively.As shown in the figure, even for z<0.1, our effective volume is still larger than the SDSS-FP catalog, with about a factor of 1.7 larger in terms of N k,modes .This is because of our larger sample even at z < 0.1 as presented in § 2.3.Note that, although our sample extends to z obs =0.2, the galaxy number density above z obs =0.1 of the SDSS MGS drops rapidly.

MCMC fitting
We use the model that is based on the distribution function and Eulerian perturbation theory as developed by Vlah et al. (2012Vlah et al. ( , 2013)); Okumura et al. (2014) and Saito et al. (2014).The analytic formula of the model is compiled in the appendix of Howlett (2019).Here we give the basic one:  where The term P11, P12, P13 and P22 can be calculated from the linear power spectrum while containing free parameters such as the growth rate f , the linear-scale bias b1, the non-linear bias b2, the non-linear velocity dispersion σv etc.Through these parameters the model accounts for the effect of smallscale non-linear motion on the power spectrum.Before running the MCMC fitting, we first quantify the effect of the window function for the momentum power spectrum measurement.Due to the survey volume as well as the sample selection as a function of both spatial positions and redshift, the power spectrum suffers from the effect of the window function.We evaluate this effect using the python code pypower.CatalogSmoothWindow3 with following steps: (1) a random elliptical catalog: we start with a random catalog available for the SDSS LSS sample4 (Blanton et al. 2005).We first follow the observed fraction of ellipticals in the LSS sample to randomly pick objects from the above catalog to define the random elliptical catalog.By fitting the observed redshift distribution of the MGS-elliptical sample with multiple Gaussian profiles, we assign a random redshift to those in the random elliptical catalog.After that, we remain those in the North Cap region and limit the redshift range to 0.02-0.2.
(2) the window function: based on the above random catalog,       0.02 and 0.1 h/Mpc.We thus carry out the MCMC fitting by incorporating correction factors into the theoretical power spectrum.The correction factors for two redshift sub-ranges of ∆z=0.02-0.10 and ∆z=0.10-0.20 are also produced following the same procedure, and their differences from the above one are small.
As discussed in Qin et al. (2019), the distribution of the momentum power spectrum at given k is not always normal as shown in Fig. 19(a), and it is thus suggested to first Box-Cox transform (Box & Cox 1964) the momentum power spectrum to Z for which the MCMC fitting is carried out: Here, because of large errors in P (k) some P(k) is negative so that we add a P shift which is set to be the minimum negative value of the distribution.The code scipy.stats.boxcox is used to find the best λ for each transformation.We set the background cosmological parameter to the Planck 2018 result with σ8(z=0)=0.811(Planck Collaboration et al. 2020), based on which the linear power spectrum at z=0 is calculated as an input of the theoretical power spectrum through nbodykit using the CLASS transfer function.In this case, if we multiply the derived f with σ8(z=0), it gives f σ8 at z = z eff .Following Qin et al. (2019), we set the following four free parameters: f , b1, b2 and σv.The MCMC fitting has been carried out through Python code emcee (Foreman-Mackey et al. 2013).Priors of four parameters are listed in Table 2.
We set kmin = 0.02 h/Mpc and conduct the MCMC fitting for kmax values of 0.15, 0.20 and 0.25 h/Mpc, respectively, for galaxies within ∆z=0.02-0.2.The corner plots of three fitting results are shown in Figure 20.The growth rate f σ8 remains unaffected by the non-linear parameters b2 and σv.While there are some correlations between f σ8 and b1 at a given kmax, the best-fit values of f σ8 exhibit no sensitivity to kmax.This is consistent with the expected advantage of the momentum power spectrum that is affected little by the bias and non-linear effects.Figure 22 offers a further comparison in f σ8 between observations and Uchuu-SDSS mocks, revealing that our observations and the mock data conform to the predictions of General Relativity (GR) and display no trends with respect to kmax.The errors from observation are slightly larger than those from mocks.This is reasonable because the median momentum power spectrum of mocks at each k almost precisely lies on the theoretical prediction.As shown in Figure 23, the best fit is slightly below the median of data points, which is caused by the fact that the covariance matrix contains significant non-diagonal parts.This is also the reason that the best-fitted f σ8 of mocks are not exactly but slightly smaller than the GR's prediction.
As shown in Figure 20 and Figure 21, there are noticeable discrepancies in the best-fit b1 between observations and mocks.This discrepancy is indicative of a disparity in the momentum power spectrum at high k regime between the two as shown in Figure 14.While both observations and mocks have a similar galaxy bias (∼ 1.5) as seen by their density power spectrum, the difference in the above result is likely due to the velocity field on small scales.As stated above, the comparison between Uchuu-SDSS and TNG indicates that Uchuu-SDSS as a N-body simulation may not fully account for non-linear effects.In addition, the observations exhibit higher values than TNG above k= 0.1 h/Mpc by about 20-30%.There may be two reasons for this: one is that there is some random zero-point offsets in velocity similar to those seen in mock simulations (see Figure 9), which adds a constant to the shot noise (Howlett 2019).This addition predominantly affects the high-k regime but becomes negligible at low-k; another is that since the mocks underestimate the power at high-k, the corresponding error at high-k is also underestimated.Nevertheless, these factors should not affect our conclusion regarding f σ8, as it primarily relies on the low k regime.
We also carry out fit to the momentum power spectra of two sub-redshift ranges of ∆z obs =0.02-0.1 and 0.1-0.2, which gives f σ8 of 0.416 +0.074 −0.076 and 0.526 +0.133 −0.148 , respectively, for kmax=0.20 h/Mpc.As shown in Figure 23, they are consistent with the result of the entire redshift range, which yields a f σ8 value of 0.459 +0.068 −0.069 , well within the 1-σ confidence interval.The inverse variance of f σ8 for the whole redshift range is almost the same as the sum of those of two redshift sub-ranges too.These results offer additional evidence for the robustness of our results.
Table 2 lists χ 2 /d.o.f.values for three redshift ranges, which are between 2 and 3. We suspect that this is due to the fact that the mock data under-estimate the momentum power spectra and associated errors at high k values as stated in § 5.1.To test this for ∆z=0.02-0.2,we increase the error by a factor of 1.5 for k above 0.1 h/Mpc and corresponding nondiagonal parts in the covariance matrix.By carrying out the MCMC fitting again, it is found that the best-fit f σ8 remains unchanged but now χ 2 /d.o.f. is 36/32.This further demonstrates that the growth rate measured from momentum power spectrum is insensitive to momentum power spectra on small scales.density power spectrum and momentum power spectrum of two data-sets together.This constraint is similar to that by Turner et al. (2023) using the combined velocity and density data of the 6dFGSv survey.If adopting the momentum power spectrum of the 6dFGSv survey alone, Qin et al. (2019) estimated 68% confidence range to be 0.226 (see their Table 2).Our measurement from ∆z=0.02 − 0.20 offers a factor of 1.7 higher accuracy.This improvement in accuracy is somewhat less than that expected by the square root of the increase in the effective volume (×6.7), but are consistent within 30%.This is reasonable given different mocks or fitting techniques (Eisenstein et al. 2005).Lai et al. (2023) measured the growth rate using the maximum-likelihood fields method through both density and PV data of the SDSS-FP survey (Howlett et al. 2022).The derived growth rate is 0.405 +0.076 −0.071 (stat)±0.009(sys), whose confidence range is consistent with our low-z one (0.416 +0.074 −0.076 for ∆z=0.02-0.1).They did not report the result when only using PV data.However, as mentioned earlier, although both our catalog and the SDSS-FP catalog are based on the same parent sample, our massE-based PV catalog includes 2.5 more objects than the SDSS-FP catalog at z < 0.1.This could account for our PV-only result exhibiting a similar level of error as theirs of combined PV and density data.Saulder et al. (2023) have provided the predicted number density and PV error for the DESI survey, from which they estimated a 68% confidence range for f σ8, which is 0.126 and 0.092 for kmax values of 0.1 h/Mpc and 0.2 h/Mpc, respectively.We measure their effective volume and find that their number of k-modes between 0.02 and 0.07 h/Mpc is 1.7 times more than our entire redshift range.This corresponds to a 1.3-fold enhancement in their precision, consistent with our constraint of f σ8 at a value of 0.459 +0.068 −0.069 at kmax=0.2 h/Mpc.Howlett et al. (2015) measured the redshift space distortion of the SDSS MGS sample through two-point correlation function, resulting in a derived f σ8 of 0.49 +0.15  −0.14 .Our massE-based PV measurement of the same sample thus offers a factor of 2.1 improvement in the error of f σ8.In Figure 24, we overlay our measurement on most recent collection of f σ8 that has been measured with redshift space distortion (Alam et al. 2021).

Fisher-matrix forecast for the massE-based PV
Our work has demonstrated that the massE can serve as a new cosmic ruler with only two nuisance parameters, allowing for the probing of the PV beyond the current redshift limit of 0.1.For future perspective, we envision that there is no apparent limit for the PV measurement with the massE to higher redshift, as it is insensitive to the accurate calibration of the massE relation itself.In addition, because of κ term in Equation 10 that is used to convert the distance measurement to PV, the PV variance scales with the redshift slower than z/(1 + z) by a few times, especially above z=0.2,as shown in Figure 25(a).As a result, although measuring the momentum power spectrum still requires a high galaxy number density above z=0.2, it has a potential to provide strong constraints on the growth rate in the late universe where dark energy dominates.To illustrate this quantitatively, we use the Fisher-matrix forecast code5 (Howlett et al. 2017;da Cunha et al. 2017) with a survey area of 14 k square degree, a distance error of ση=20%, a galaxy number density of 8×10 −3 h 3 /Mpc 3 and Equation 10 for the PV error.As shown in Figure 25(a), PV-only data yields stronger constraints compared to the RSD only for this specific data-set where kmax is set to be 0.2 h/Mpc.On the other hand, the momentum power spectrum mainly relies on linear scales to extract cosmological information, complementing the RSD.The combination of these two offers significant improvement in the growth rate, potentially exceeding the DESI's constraints by a factor of 1.7 to 4.7 in terms of errors. (a).

Figure 1 .
Figure 1.(a), the footprint of our massE-based SDSS MGS-PV catalog.Galaxies in three dotted boxes are used to measure the velocity zero point.(b), the north cap portion of the PV catalog for the momentum power spectrum measurement.

Figure 2 .
Figure 2. (a), the distribution of the massE-based comoving distance in one redshift bin.The solid red curve is the best-fit Gaussian distribution and two dotted lines mark the three standard deviations.(b), the new distribution after removing outliers as detailed in § 2.3.

Figure 4 .
Figure 4. (a), the number density of galaxies in our SDSS massEbased PV catalog, where symbols are the data and curves are the interpolated one.(b), the PV variance as a function of redshift, where symbols are the data and the curve is the interpolated one.

Figure 5 .
Figure5.(a), the r-band magnitude errors versus the apparent magnitude for the entire SDSS photometric catalog.(b), the distribution of r-band magnitude errors in one magnitude bin, along with the best-fitted double Gaussian profiles.(c), the r-band luminosity function with the characteristic brightness labeled by a star symbol.Three dotted lines mark the absolute luminosity of the magnitude cut (m r,limit =17.6) at three redshifts.(d), the distribution of the difference between the intrinsic distance and observed distance at z=0.2, which are derived from the intrinsic magnitude and observed magnitude, respectively.The filled histogram represents the one with intrinsic apparent magnitude below m r,limit , while open one for the observed apparent magnitude below m r,limit .(e), the fractional difference between intrinsic and observed distribution (the denominator is the observed distribution) as a function of the offset from the mean in terms of standard deviation of the distribution at z=0.2.The result contains ten simulations.The solid line is the best-fitted linear function with zero intercept.(f ), the slope of above linear function as a function of redshift, where the solid line is the best-fitted power law function.
Obs. R e vs. M r Fig.5-a Mock M r Obs.M/L vs. g-r Fig.5-b Mock g-r Mock M r massE relation 1 Mock s e Obs.log(s e ) Fig.7-c Mock M star, 1Mpc with error Mock s e with error Mock R e, as with error Mock v p Mock R e, as Mock M star, 1Mpc massE relation 2

Figure 6 .Figure 7
Figure 6.The flow chart to mimic the observed PV for Uchuu-SDSS mock galaxies.The first massE relation is the intrinsic one, and the second relation is used to measure PV.We test the cases if they are different in § 4.4.

Figure 8
Figure8.The first row: for each mock (mock 0 to mock 30), the zero velocity point (Gaussian mean) is determined from galaxies in that particular mock; the second row: the zero velocity point is determined from galaxies in that particular mock plus all galaxies in mock 31; the bottom row: for each mock, the zero velocity is determined from part of galaxies in mock 31 and part of galaxies in that particular mock to mimick galaxies in yellow dashed boxes in Figure1.For details, see § 4.3.From left to right, the momentum power spectra are presented for three redshift ranges.The solid line is median of the intrinsic spectra of 31 mocks (mock 0 to mock 30), and symbols are the median of 31 mocks calculated with the massE-based PVs.

Figure 9 .
Figure 9.The mean of the massE-based PV after 2000 simulations versus the true PVs in mocks.The lower panel for each redshift shows the residual between two velocities.

Figure 10 .
Figure10.The median of the intrinsic momentum power spectra of Uchuu-SDSS mock catalogs (solid line), as compared to its 1860 simulations of PV measurements that mimic observations (symbols are the median and error bars are the square root of the diagonal covariance).

Figure 11 .
Figure11.The same as Figure9but now the intrinsic massE in the mock has no redshift evolution, while the one used to derive PV has a strong redshift evolution in the intercept with S D 0 ∝ (1 + z) 5 .

Figure 12 .
Figure12.The same as Figure9but now the intrinsic massE in the mock has no redshift evolution, while the one used to derive PV has a strong redshift evolution in the slope with β ∝ (1 + z) 0.4 .

Figure 14 .
Figure 14.(a), the momentum power spectrum based on our SDSS massE-based PV catalog for the entire redshift range ∆z=0.02-0.2.Symbols are the observed data points.Error bars indicate the square root of the diagonal elements of the covariance.Results of N-body (Uchuu-SDSS) and hydrodynamic (TNG) simulations are also shown.(b), the same as (a) but for the redshift range ∆z=0.02-0.1.(c), the same as (a) but for the redshift range ∆z=0.02-0.1.

Figure 16 .
Figure 16.(a), the momentum power spectra obtained with different massE relations, including the best-fit one, a strong redshift evolution in the slope of the relation as well a strong redshift evolution in the intercept of the relation.(b), the momentum power spectrum obtained with three different sizes of redshift bins.
Figure 17.The effective volume of our PV catalog as compared to the literature PV catalogs.The calculation follows Equation 20.

Figure 18
Figure18.(a,) the window function of our PV catalog in the configuration space for our sample within ∆z=0.02-0.20.(b), the ratio of the theoretical spectrum convolved with effects of both wide angle and window function to the pure theoretical spectrum.The theoretical spectrum at high k resolution is calculated following Equation 21 at k from 0.001 to 10 with a bin size of 0.005.

Figure 20 .k
Figure 20.The corner plot of four free parameters used to fit the observed momentum power spectrum.Three results corresponding to three kmax values are overlaid.

Figure 22 .Figure 23 .
Figure 22.The best-fit f σ 8 values of observations and mocks within ∆z=0.02-0.2 for three different kmax values, where k min is fixed at 0.02 h/Mpc.
Figure 25.(a), the error variance of PV as a function of redshift for different formulae.The massE-based PV follows σ 2 η /(α + κ) 2 , where ση is the distance error of ln(Dc), α and κ are defined in Equation 10. (b), the Fisher matrix forecast for the constraint on f σ 8 for a survey with an area of 14 k square degree and average galaxy number density of 8×10 −3 h 3 /Mpc 3 .The DESI's forecast is from DESI Collaboration et al. (2016).

Table 1 .
The steps to select the sample.
Figure 15.The covariance matrix of the derived momentum power spectrum in Figure 14.