Arxiv:2108.02207V1 [Astro-Ph.CO] 4 Aug 2021 Without Many Localised FRB Sources [26]

Cosmology with the moving lens eﬀect

Selim C. Hotinli,1 Kendrick M. Smith,2 Mathew S. Madhavacheril,2 and Marc Kamionkowski1 1Department of Physics & Astronomy, Johns Hopkins University, Baltimore, MD 21218, USA 2Perimeter Institute for Theoretical Physics, 31 Caroline St N, Waterloo, ON N2L 2Y5, Canada (Dated: August 6, 2021) Velocity fields can be reconstructed at cosmological scales from their influence on the correlation between the cosmic microwave background and large-scale structure. Effects that induce such correlations include the kinetic Sunyaev Zel’dovich (kSZ) effect and the moving-lens effect, both of which will be measured to high precision with upcoming cosmology experiments. Galaxy measurements also provide a window into measuring velocities from the effect of redshift-space distortions (RSDs). The information that can be accessed from the kSZ or RSDs, however, is limited by astrophysical uncertainties and systematic effects, which may significantly reduce our ability to constrain cosmological parameters such as fσ8. In this paper, we show how the large-scale transverse-velocity field, which can be reconstructed from measurements of the moving-lens effect, can be used to measure fσ8 to high precision.

I. INTRODUCTION rameter, defined such that brsd = 1 in the absence of anisotropic selection effects.1 Anisotropic selection Next generation cosmic microwave background (CMB) effects could significantly degrade measurements of experiments such as the Simons Observatory (SO) [1,2] f σ8 from future galaxy surveys, due to degeneracy and CMB-S4 [3], and galaxy surveys such as DESI [4] with brsd. In this paper, we propose a complementary and the Vera Rubin Observatory (VRO) [5] will gener- method for measuring f σ8 in which all degeneracies ate a wealth of new data with unprecedented precision with astrophysical bias parameters are broken. on small scales. Correlations between CMB anisotropies We will show that the temperature fluctuations and the galaxy density carry valuable cosmological in- sourced by the transverse velocities of large-scale struc- formation about the largest scales, creating novel op- ture (the moving-lens effect [22]) can be used to mea- portunities for inference [e.g.6–15]. These correlations sure fσ8 and to break the kSZ optical-depth degener- that are induced by interactions of CMB photons with acy. These transverse velocities are inferred by cross- the intervening large-scale structure include the thermal correlation of the effect they induce in the CMB with and kinetic Sunyaev Zel’dovich effects [16–20] and the galaxy surveys [e.g. 32–34]. The moving lens effect integrated Sachs-Wolfe effects [21], which includes the sourced by individual objects is expected to be small moving-lens effect [22]. Among these, the kinetic Sun- and can be easily confused with other effects. However, yaev Zel’dovich (kSZ) and the moving-lens effects depend combining the signal from many objects allows a tomo- on the peculiar velocities of cosmological structures. graphic reconstruction of the transverse-velocity field on large scales, analogous and complementary to kSZ to- A key product of velocity measurements is the combi- mography, which provides the radial-velocity fields [e.g. nation of the linear-theory growth rate f and the ampli- 23, 25]. Note that the transverse-velocities measured tude σ of matter fluctuations on the scale of 8h−1Mpc. 8 from moving-lens tomography will be biased similar to In principle, fσ can be determined by measurements of 8 the case of kSZ tomography. The moving-lens bias b the kSZ effect [e.g. 11–13, 23–29], or by measurements of ⊥ can be measured through galaxy-galaxy or galaxy-CMB redshift-space distortions (RSDs) on galaxy surveys [30]. lensing, however, since b is proportional to the galaxy- The cosmological information available from the kSZ ef- ⊥ matter cross-correlation, in contrast to the kSZ bias fect, however, is limited by the degeneracy of this signal b which is proportional to the galaxy-electron cross- with the optical depth of galaxies (except in the case of k correlation, which is difficult to measure independently primordial non-Gaussianity [12]). While fast radio burst arXiv:2108.02207v1 [astro-ph.CO] 4 Aug 2021 without many localised FRB sources [26]. (FRB) sources may provide a way to break some of these biases by measuring the galaxy-electron correlation [26], Small-scale CMB anisotropies are expected to be dom- the prospect of localising enough FRBs in the near future inated by the kSZ effect due to Compton scattering of is challenging. the CMB photon off free electrons that have non-zero peculiar velocities with respect to the CMB rest frame. Furthermore, recent studies show evidence of This induces a shift in the CMB temperature in the di- anisotropic selection effects [e.g. 31] for galaxy RSD rection of the free-electron radial velocity. The moving- measurements. Such effects are degenerate with redshift- lens effect is due to time-evolving gravitational potentials space distortions, in the sense that the large-scale bias on the line of sight. This induces additional tempera- of the galaxy field becomes (b +b fµ2), rather than g rsd ture anisotropies in the CMB, via the non-linear ISW 2 ˆ simply (bg +fµ ), where µ = rˆ k. Here, rˆ is the unit · 1 vector in the radial direction, kˆ is the three-dimensional According to Ref. [31], (brsd − 1) can be of order 0.1 − 0.2 in subcatalogs of SDSS. Fourier (unit) wavevector and brsd is a new bias pa- 2

(or Rees-Sciama) eﬀect. This signal, which is a dipole 104 pattern oriented along the object’s transverse velocity, is expected to be smaller than kSZ. Even so, it will be 103 measured with high signal-to-noise with upcoming CMB CIB ν =150GHz experiments [32–34]. ) 2 tSZ π 10 }

In what follows, we demonstrate that the transverse (2 late-kSZ / ` velocities from the moving-lens effect, when combined C reio-kSZ with galaxy measurements, can provide competitive con- 101 + 1) ML ` straints on the fσ parameter, comparable to the sce- ( 8 ` S4 nario where the RSD bias is measured externally to high 100 SO precision. This parameter is useful for studying a large HD range of physics, including dark energy [35], modified 1 10− gravity [36] and the effects of neutrino mass [37]. Pre- 3 4 cision measurement of fσ8 also allows using the kSZ to 10 10 learn about astrophysics, such as the characteristics of `, multipoles the electron density profiles around halos, by breaking FIG. 1. CMB signal and foregrounds at the frequency the degeneracy between the electron-scattering optical ν = 150 GHz. Total (lensed) CMB signal is shown with the depth and the growth rate [25]. black solid line. We show the anticipated post-component sep- This paper is organised as follows: We begin with a aration noise (including foreground residuals) for three exper- description of kSZ and moving-lens tomography in Sec- iments; CMB-S4 (S4), Simons Observatory (SO) and CMB- HD (HD); along with the most significant foreground com- tionII. We show the dependence of the velocity field ponents (before foreground cleaning) in temperature maps on cosmology and discuss the anticipated RSD bias in (shown at 150 GHz). These are the cosmic infrared back- SectionIII. We introduce a formalism for the three- ground (CIB) the thermal Sunyaev Zel’dovich effect (tSZ), velocity in SectionIV. We introduce a Fisher analy- late-time kSZ effect (due to free energetic electrons in in- sis in SectionIIIB. We provide results from our fore- tergalactic mediums where z ∈ [0, 6]), the kSZ signal from casts using experimental specifications anticipated for patchy-reionization and the moving-lens (ML) effect. next generation CMB and galaxy surveys in SectionV. We conclude with a discussion in SectionVI. Through- out, we use the standard ΛCDM model with cosmolog- projected onto the radial direction. Most of the cos- 2 2 mological information is contained in the radial-velocity ical parameters Ωb h , Ωcdmh ,As, ns, τ set equal to 0.022, 0.12, 2.2 {10−9, 3.04, 0.965, 0.06 , respectively.} field vk while ne depends primarily on astrophysics and { × } non-linear large-scale structure; see Ref. [25] for a detailed discussion of this point. We demonstrate the II. VELOCITY RECONSTRUCTION anticipated late-time kSZ signal in Fig.1, along with other foregrounds and the forecasted noise levels for next- generation CMB experiments. The component sepa- A. kSZ tomography rated CMB noise due to foreground residuals, including the cosmic infrared background (CIB), radio sources, Kinetic Sunyaev Zel’dovich tomography aims to ex- the thermal Sunyaev Zel’dovich effect (tSZ), and the tract cosmological information from the kSZ effect by us- kSZ signal are shown in Fig.1 for the Simons Obser- ing measurements of the CMB and a tracer of the electron vatory, CMB-S4 and CMB-HD experimental specifica- density, such as a galaxy survey, to reconstruct the radial- tions (these were produced using the orphics2 code as velocity field. The temperature anisotropy induced by described in AppendixA). We use the electron model the kSZ effect from a large-scale structure (LSS) shell at described in Ref. [38]. redshift z = z∗ is The (inverse) noise on the velocity reconstruction is Z L given by [25] ΘkSZ(θ) = K(z∗) dr qk(x), (1) 0 2 Z 2 ! 1 K∗ ksdks Pge(ks) = 2 TT,obs , (3) where x χ?θ + rrˆ, χ? is the conformal distance to the N (k ) χ 2π obs k L ∗ Pgg (ks)C` LSS shell,≡rˆ is the unit vector in the radial direction, θ `=kχ∗ is the angular direction on the sky, Θ(θ) = ∆T (θ)/T¯ is where k is the three-dimensional Fourier wavevector and the fractional CMB temperature fluctuation, and the integral is over small-scale Fourier modes kS. We −τ(z) 2 represent large-scale modes with an ‘L’ subscript. Here, K(z) = σT ne,0xe(z)e (1 + z) , (2) − Pge(k) is the power-spectrum of the galaxy-electron cor- obs is the radial weight function with units of 1/Mpc. Here, relation, Pgg (k) is the observed galaxy power spectrum TT,obs τ(z) is the optical depth at redshift z, σT is the Thom- and C` is the observed CMB spectrum. We show son cross-section, ne is the free electron number-density 2 github.com/msyriac/orphics and qk(x) = δe(x)vk(x) is the electron-momentum field, 3 the reconstruction noise from kSZ tomography in Fig.2. satisfy faH vˆ = b µ δ (k) , (9) k k k m B. Moving-lens tomography p 2 faH vˆ⊥ = b⊥ 1 µ δm(k) , (10) − k Gravitational potentials that evolve in time induce a temperature modulation on the CMB known as the inte- where H is the Hubble parameter, a is the scale factor, grated Sachs-Wolfe (ISW) effect, which has the form δm(k) is the matter overdensity, f = d ln D(a)/d ln a is the growth rate and the bias parameters bk and b⊥ are Z L true fid due to the mismatch between Pge and Pge , and sim- ΘML(θ) = 2 dr Φ(˙ x) , (4) − 0 ilarly for Pgm. Here, D(a) is the linear-theory growth factor for the matter spectrum that parametrises the where Φ is the gravitational potential. The ISW effect time evolution of the matter power-spectra, Pmm(a) = can be sourced by peculiar velocities of potentials, lead- 2 D (a)Pmm(a = 1). The combination of the galaxy and ing to the moving-lens effect, which has the form the velocity satisfies

Z L 2 2 Pgg(k, µ) = (bg + brsdfµ ) Pmm(k) , (11) ΘML(θ) = 2 dr ∇⊥Φ(x) v⊥(x) , (5) − 0 · faH 2 Pgv(k, µ) = b (bg + brsdfµ )Pmm(k) , (12) where v⊥(x) is the peculiar (comoving) transverse veloc- k ity. Our focus is on the large-scale velocity field, where we faH 2 anticipate that the true velocity is linear and curl-free, P (k, µ) = b2 P (k) , (13) vv k mm such that we can define a transverse velocity potential Υ(x) as v⊥(x) = ∇⊥Υ(x). We demonstrate the antici- p 2 where b b µ, b⊥ 1 µ for the radial and the trans- pated moving lens (ML) signal in Fig.1. The ML signal is ∈ { k − } verse velocities, respectively. Here, bg is the linear galaxy calculated using the velocity-reconstruction pipeline de- bias, P is the galaxy-velocity cross power spectrum and scribed in Ref. [39] for the redshift range z 0.1, 5 . gv ∈ { } Pgg is the galaxy auto power spectrum. Above, we in- Moving-lens tomography can be used to measure the 2 cluded the RSD effect [30] fµ with a bias factor brsd, large-scale transverse-velocity field. In the limit k which we discuss next. These equations demonstrate the kS, `/χ∗, we find the noise on the velocity reconstruction dependence of observables on cosmological parameters. to be ! 1 2µ2 k2 Z k3 dk P (k )2 = ⊥ L S S ΦΦ S , Υˆ Υˆ 2 obs ΘΘ,obs B. Redshift space distortions N (kL) χ∗ 2π P (kS)C ΦΦ ` `=kS χ∗ (6) In the absence of measurement biases (and on large- p 2 ˆ obs scales) the effect of RSDs is given by the Kaiser formula, where µ⊥ = 1 µ , µ = k rˆ, PΦΦ (k, µ) = PΦΦ(k) + − · which states that the galaxy profile ug is modified as NΦΦ(k, µ), and 2 2 4 2 2 ug(k) 1 + βµ ug(k) , (14) NΦΦ(k, µ) = [(1 + z)ρm,0/2] /[k bg ng W (k, µ) ] , (7) → with the photo-z term, where β = f/bg. Equivalently, the galaxy-density field in Fourier space takes the form 2 2 2 2 W (k, µ) = exp µ k σz /2H∗ , (8) 2 {− } δg(k) = (bg + µ f)δm(k) + noise . (15) where σz is the photo-z error and H∗ is the Hubble parameter evaluated at redshift z . Here, ρ = 3Ω H2, Extracting cosmological information from fδm(k) as- ∗ m,0 M 0 sumes that one can calibrate a Fourier mode measured ΩM is the matter fraction and H0 is the Hubble constant. We show the shot noise on the gravitational potential in from the galaxy density field to the comoving k, subject Fig.2 for the galaxy surveys we consider, together with only to the knowledge of the distance scales, Da(z) and the reconstruction noise from moving-lens tomography. H(z), for kk and k⊥, respectively. This can be done by using the BAO signature and the Alcock-Paczynski [40] method, internally for a galaxy survey, for example. In reality, however, the galaxies are not randomly dis- III. COSMOLOGY FROM VELOCITIES tributed and their orientations show correlations with their local large-scale structure environment (tidal fields, A. Cross-correlation signal density, dust, etc.). The appearance of the galaxy depends on the line-of- Where linear theory is valid, the reconstructed velocity sight, expressed in the galaxy’s frame of reference. This fields are proportional to the cosmic growth rate. They introduces an anisotropic selection function when mea- 4

9 10− kSZ ML : VROY10 4 10 kSZ ML : DESI 10 10− ] 3 ] 3 11 10−

102 ) [Mpc k ) [Mpc ( k

( 12 ΦΦ 10− vv P P 4

k DESI 0 13 10 10− VROY1 µ = 0.1 µ = 0.9 µ = 0.4 VROY10 14 3 2 3 2 10− 3 2 1 0 10− 10− 10− 10− 10− 10− 10− 10 k [1/Mpc] k [1/Mpc] k [1/Mpc]

FIG. 2. (Left) Comparison of the velocity signal to reconstruction noise at redshift z = 1. The power spectrum of the norm |v| of the 3-velocity is shown with solid black curves. Coloured lines correspond to reconstruction-noise spectra, calculated from the anticipated moving-lens (ML) and kSZ signals. Two galaxy-survey configurations are shown; these are DESI and Rubin Observatory Year 10 (VROY10). The CMB-S4 experimental specifications are used for all noise spectra. The two subplots correspond to nearly-transverse (µ = 0.1) and nearly-radial (µ = 0.9) velocity reconstruction. (Right) The gravitational potential signal (solid black line) is shown together with the shot noise for the galaxy surveys we consider in this paper. In all plots, the rise in the noise on small-scales for Rubin is due to scatter from photometric redshifts. suring galaxy clustering, which depends on the intrinsic where alignment of galaxies with its local environment [41]. In P fid(k ) this work, we represent such effects with an RSD bias as ge S 1 F (kS) = kS obs TT,obs , (19) Pgg (kS) C 2 `=kS χ δg(k) = (bg + brsdµ f)δm(k) + noise . (16) which depends on the difference between the true galaxy- Evidence for an RSD bias has recently been found at true fid electron correlation, Pge , and the fiducial model, Pge . 5σ in Ref. [31] where authors find brsd variations on the order of 10 20% within the BOSS CMASS and LOWZ samples [42−]. Next, we demonstrate the improvement on A similar bias can also appear in the transverse- cosmological constraints from the moving-lens effect in velocity reconstruction, as: the presence of an RSD bias. true vˆ⊥(k) = b⊥v⊥ (k) (20) + noise, with power spectrum N⊥(k, µ) ,

true where v⊥ is the true transverse-velocity amplitude. For IV. FISHER ANALYSIS the moving-lens reconstruction, this bias is due to the uncertain cross-correlation between the gravitational potential and the density tracer. For a galaxy field tracer, Following Ref. [25], we use a simplified ‘snapshot’ ge- for example, this bias can be written as: ometry throughout this paper to represent the observed R true Universe, which we take to be a 3D box with comoving dkSF (kS)P (kS) 3 Φg side length L and volume V = L , at a constant redshift b⊥ = R fid , (21) dkSF (kS)PΦg (kS) z∗ = 1 corresponding to a distance χ∗ 3400 Mpc. ' The optical-depth degeneracy can be represented on where F (k ) = k P fid(k )/[P obs(k )CTT,obs] . Here, S S Φg S gg S `=kS χ large scales with a bias on the radial-velocity measure- true PΨg (k) is the true galaxy-potential correlation and ment, in the form, fid PΨg(k) is the fiducial model. true vˆk(k) = bkvk (k) (17) + noise with power spectrum Nk(k) , We take the true velocity field to be curl-free, i.e. where vtrue is the true radial-velocity amplitude. The true ˆ k vj (k) = ikjs(k) , (22) radial velocity bias bk depends only on the small-scale modelling of the electron-galaxy correlation. It can be where s(k) is a scalar field with power spectrum Pvv(k) = 2 calculated as [25]: (faH/k) Pmm(k). To avoid covariance matrices which are not full rank, we combinev ˆk andv ˆ⊥ into a “unified” R dk F (k )P true(k ) S S ge S reconstructionv î (a three-vector). Then, the two-point bk = R fid , (18) dkSF (kS)Pge (kS) function ofv î can be written as the sum of signal and 5 noise contributions from the velocity, as: DESI, σ(brsd)/brsd = 5% 0 ∗ 3 3 0 4% vî(k)v ˆj(k ) = Sij(k) + Nij(k) (2π) δ (k k ) . − (23) 3% The average (shown with brackets) is over realisations of 8 the three-velocity. 2% /fσ ) Given a wavenumber ki and line-of-sight directionr î, 8 fσ one natural choice of basis for the 3-by-3 matrices Sij and ( σ 1% Nij is the “radial/tangential” basis, wherer î is a basis vector: 2 1.4 arcmin beam 10− 30 arcsec beam  1   kµ  p 2 DESI, b = 1 rî =  0  k =  k 1 µ  . (24) rsd 0 − 1 0 1 0 10− 10 10 ∆ [µK arcmin] T − In what follows, we add the galaxy field to Eq. (23), pro- FIG. 3. Constraints on fσ8 from moving-lens tomogra- moting Sij(k) and Nij(k) to 4-by-4 matrices. phy and galaxy clustering. Coloured lines correspond to constraints from the reconstructed transverse-velocity fields, In the radial/tangential basis, the noise covariance and their cross-correlation with the galaxy field on large scales [including the galaxy field, g(k)] can be read off from (with no RSD bias prior). A range of CMB experiment noise Eqs. (18) and (21): levels (∆T ) are shown and two beam sizes, appropriate for CMB-S4 and CMB-HD. We scale the noise levels for all of Nij(k) (25) the frequency bands of these experiments. The ILC calcula- 2 = diag Nk(k, µ),N⊥(k, µ),N⊥(k, µ), 1/[W (k, µ)ng] , tion is discussed in AppendixA. Blue (orange) markers show { } anticipated noise levels of CMB-S4 (CMB-HD). We use Ru- where ng is the galaxy number-count density in a given bin Observatory Year 10 (VROY10) for moving-lens tomog- box. To compute the signal covariance, we combine raphy. Note that here the radial-velocity reconstruction only Eqs. (18) and (21) to get marginally improves the constraints on fσ8: the constraints are dominated by the correlation of the transverse-velocity true field and galaxy clustering. We consider a 1% prior on the vî(k) = Bij(k)vj (k) + (noise) , (26) transverse-velocity bias [defined in Eq. (21)], from an exter- where the 4-by-4 bias matrix Bij(k) is given by nal measurement of the matter-galaxy correlation on small scales, such as galaxy and lensing correlations, for example. Bij(k) = diag bk, b⊥, b⊥, bg . (27) The 1% constraints on the transverse-velocity bias falls near { } the scenario where b⊥ is fixed and not marginalised. Finally, true By Eq. (22), the signal covariance of vi is we show constraints on the fσ8 parameter from measurements of the RSD effect, using the spectroscopic DESI galaxy sur- true true 0 ∗ ˆ ˆ 3 3 0 vey. The dotted gray horizontal lines corresponds to 5%−1% vi (k)vj (k ) = kikjPv(k) (2π) δ (k k ) , (28) − priors on the RSD bias from top to bottom. The solid gray and the signal covariance ofv î is horizontal line correspond to fixing brsd to 1. Recent studies show evidence for the case brsd 6= 1 [31]. Sij = cicjPvv(k) , (29) where 2D integral,

Z ∞ 2 Z 1 ci = (ck, c⊥, 0, cg) , (30) V k dk Fab = 2 dµ (35) c b µ , (31) 2 0 4π −1 k ≡ k p 2 ∂S(k, µ) ∂S(k, µ) c⊥ b⊥ 1 µ , (32) Tr C−1(k, µ) C−1(k, µ) . ≡ − 2 × ∂πa ∂πb cg (bg + brsdfµ ) k/(faH) . (33) ≡ In our forecasts we consider the following 6 parameters: Given parameters π , , π , the N-by-N Fisher ma- 1 N fσ8, bgσ8, bk, b⊥, brsd . trix is ··· { } Z 3 V d k −1 ∂S(k) −1 ∂S(k) Fab = 3 Tr C (k) C (k) , V. RESULTS 2 (2π) ∂πa ∂πb (34) where C = S + N is the total covariance. We can sim- We consider a range of CMB experiments with varying plify a little by noting that the signal and noise covari- noise, including SO and CMB-S4, as well as the futur- ance matrices in Eqs. (26) and (29) only depend on k istic CMB-HD survey. The galaxy-survey (shot) noise through (k, µ). This lets us reduce the 3D integral to a depends on the number-density function, ng(z). We con- 6 sider three LSS experiments: DESI, Vera Rubin Obser- In TableI, we show the detection signal-to-noise ra- vatory Year 1 (VROY1) and Year 10 (VROY10). The tio (SNR) of the reconstructed velocity and galaxy den- parameters for VRO are based on [43, 44] and those for sity cross-correlation Pvgˆ (k) from our experimental con- DESI are based on [45]. For DESI, we take a 116 Gpc3 figurations. The velocities are reconstructed either form −4 3 box with galaxy survey density ng = 1.4 10 /Mpc the moving-lens (ML) effect or the kSZ tomography. × 3 and σz = 0. For VROY1, we take a 113.4 Gpc box, We find similar SNR forecasts compared to previous re- −3 3 with galaxy survey density ng = 6.9 10 /Mpc . For sults [32, 34]. Nevertheless, note that our calculation of VROY10, we take a 180Gpc3 box, with× galaxy survey the transverse velocity noise differs from earlier work by −2 3 3 density ng = 1.2 10 /Mpc . For both VROY1 and the use of the simplified box formalism here. × VROY10, we consider photo-z errors satisfying σz = 0.06 We show forecasts on fσ8 in Fig.3. We find that at z = 1. For DESI, VROY1, VROY10 , we consider a moving-lens tomography, together with the galaxy den- { } galaxy bias of 1.5, 1.7, 1.6 , respectively. We take all the sity, can constrain fσ8 to high precision, comparable to boxes to be centered{ at z}= 1. For all forecasts, we set the scenario where the RSD bias is known up to few per- the maximum wavenumber we include in our analysis to cent accuracy. For SO and VROY10, improvements on kmax = 0.1/Mpc. fσ8 relative to RSD can be achieved if the RSD bias The total CMB power spectrum gets contributions is uncertain at the level of σ(brsd) > 5%. Using S4 from weak gravitational lensing, the kSZ effect (both and VROY10 further improves the benefit from using from reionization and late times) and the experimental the transverse velocities, while constraints from CMB- noise, which we take to satisfy HD and VROY10 improves upon the scenario where the 2 RSD bias is known up to around 1 2% accuracy. 2 `(` + 1)θFWHM − N` = ∆T exp . (36) Note that in our forecasts we considered a 1% prior on 8 ln 2 the transverse velocity field. This can potentially be pro- Note that when forecasting the CMB noise for these ex- vided from external measurements of the galaxy-matter periments, we also include the contribution to the power cross-correlation on small scales (from cross-correlation spectra from residual foregrounds after foreground clean- of galaxy-lensing measurements, for example) as defined < ing we discuss in AppendixA. Fig.1 demonstrates how in Eq. (21). For upcoming experiments, we find 1% the CMB noise and different residual foregrounds com- prior on the transverse-velocity bias recovers the cosmic-∼ pare with the CMB signal. variance limit with the transverse-velocity bias fixed as b⊥ = 1. Fig.4 shows the constraints on bk from a combination of moving-lens tomography, kSZ tomography, and galaxy DESI, σ(brsd)/brsd = 5% clustering. Moving-lens measurement of fσ8 removes the 4% 3% TABLE I. The detection signal-to-noise (SNR) of the (re- k 2% constructed) velocity and galaxy-density cross-correlation: /b

) P (k). Velocities are reconstructed from the kSZ and

k vgˆ b

( moving-lens (ML) tomography for various CMB and LSS ex- σ 1% periments. VROY1 and VROY10 refer to Vera Rubin Obser- 10 2 − vatory Year 1 and 10 respectively. 1.4 arcmin beam kSZ SNR CMB ML SNR CMB 30 arcsec beam LSS SO S4 HD LSS SO S4 HD DESI, brsd = 1 DESI 231 414 1170 DESI 7 8 62 1 0 1 VROY1 116 210 669 VROY1 16 28 85 10− 10 10 ∆ [µK arcmin] VROY10 123 228 775 VROY10 20 37 112 T − FIG. 4. Constraints on bk from combination of moving- lens tomography, kSZ tomography and galaxy clus- 3 tering. Coloured lines correspond to constraints from the Our forecasts for the moving-lens SNR differs from Refs. [32, 34], reconstructed transverse- and radial-velocity fields, and their in various ways: First, Refs [32, 34] included a map-based reconstruction on the two-sphere and with 11 redshift bins in the cross-correlation with the galaxy field on large scales (with range z ∈ [0, 3], unlike our single-box formalism which is simpler no RSD bias prior); similar to Figure3, with 1% prior on the but potentially more optimistic than the former. Furthermore, transverse-velocity bias. We use Rubin Observatory Year 10 forecasts in Refs [32, 34] did not include the adverse effects of (VROY10) survey for kSZ and moving-lens tomography. Gray residual foregrounds after foreground cleaning (except the ef- lines correspond to constraints on bk from measurements of fects of kSZ and weak lensing), while including delensed CMB the RSD effect and kSZ tomography, using the spectroscopic spectra with ideal noise. Delensing noticeably improves the de- DESI galaxy survey and the CMB-S4 experimental specifica- tection prospects of reconstruction at the low-noise limit, hence tions. The dotted gray lines corresponds to 5% − 1% priors important for experiments like CMB-S4 and CMB-HD, for ex- on the RSD bias from top to bottom. The solid grey line ample. Despite these differences, we find matching SNR forecasts for CMB-S4 and SO. correspond to fixing brsd to 1. 7

degeneracy can be broken by relatively robust measure- ) 8 HD: VROY10 VROY1 ments such as correlations of the galaxy density and fσ

( (galaxy or CMB) weak lensing. Overcoming the former

σ S4: VROY10 VROY1 “optical depth” degeneracy, however, is likely more diffi- SO: VROY10 VROY1 cult. Galaxy clustering is also sensitive to cosmic growth 101 due to the RSD effect [30], but future measurements may be biased by selection effects. Using moving-lens recon- structions that combine the Vera Rubin Observatory and CMB surveys on the other hand provides strong constraints on cosmological parameters such as fσ8 at the level of < 5%, 3% and 2% for the Simons Observatory, CMB-S4∼ and CMB-HD experiments respectively.

fractional improvement in 0 We focused on systematics in the RSD effect, but mov- 10 1 0 10− 10 ing forward, a more detailed analysis involving the mod- σ(brsd)/brsd (external prior) elling of all relevant effects will give further insight into the benefit of moving-lens tomography. These effects can FIG. 5. The fractional improvement in σ(fσ8) from the moving-lens effect, relative to the constraint from be separated into two categories. First, the theory pre- biased galaxy RSDs. The ratio of constraints from DESI dictions from our models are often gauge dependent and RSDs (with varying RSD bias prior) and the combined mea- defined in comoving coordinates, rather than in redshift surement of biased galaxy RSDs and transverse velocities is space. Mapping from the theory to gauge-independent shown for various experimental configurations. observables (such as the galaxy power spectrum or CMB temperature) requires a multitude of correction terms on top of the density fluctuations, which take into ac- degeneracy between the radial-velocity bias bk and fσ8, suffered by the radial velocities reconstructed with kSZ count the photon geodesics in an inhomogeneous Uni- verse. These include general relativistic effects such as tomography. This allows kSZ to measure directly bk, al- lowing the possibility of astrophysical inference. Doppler (magnification) terms, potential (Sachs-Wolfe, integrated Sachs-Wolfe, time delay) terms and weak lens- In Fig.5, we show the ratio of the errors on fσ8 ob- tained from the DESI RSDs and from forecasts com- ing, as well as anisotropies in the mapping from real space bining VRO galaxies and the moving-lens tomography. to redshift space induced by peculiar velocities (RSDs). The RSDs are taken to be biased for both calculations Second, galaxy clustering measurements are subject to selection effects, where detecting particular galaxies with with a varying prior assumption on σ(brsd) (shown on the x axis). Our results suggest that moving-lens tomog- various properties may be more or less likely than on average. These effects include tidal alignments [41], lensing raphy improves the fσ8 constraints for σ(brsd)/brsd > few%. Note that in the absence of prior knowledge on∼ magnification [46], and redshift-evolution biases [e.g. 47]. These biases can complicate the mapping between theory the RSD bias [or for σ(brsd) 5% or worse], moving- lens tomography may play an essential∼ role in our ability and observed quantities and should be explored in future to constrain cosmology from velocity reconstruction. So work. far, the only other observable that can constrain fσ8 (to- The upcoming decade will host a wealth of new and gether with the kSZ tomography) is the FRB dispersion high-quality data. These developments will open many measurements, which allow the breaking of the optical- novel opportunities for cosmological inference. In par- depth degeneracy of the kSZ signal [26]. Our forecasts ticular, using the CMB as a cosmological backlight, suggest the moving-lens effect may provide a constraining secondary fluctuations induced by the interaction be- power comparable to measurement of 105 FRBs, de- tween CMB photons and large-scale structure allow new pending on the uncertainty on the FRB∼ dispersion mea- methods like transverse-velocity reconstruction from the surement [26]. moving-lens effect, as described in this paper. These new observables provide the different tracers that can potentially improve the constraints on cosmological parame- VI. DISCUSSION ters, beyond the limits imposed by selection effects, astrophysical uncertainties, and cosmic variance. In this work, we have explored one of the cosmological applications of reconstructing the transverse-velocity field from the moving-lens effect. Velocity fields are di- VII. ACKNOWLEDGEMENTS rect tracers of the combination of cosmic growth rate and the amplitude of matter fluctuations and contain valu- We thank J. Colin Hill for comments and suggestions able cosmological information. While velocity measure- on the ILC noise calculation. SCH is supported by the ments from the kSZ and the moving-lens effects are bi- Horizon Fellowship from Johns Hopkins University. SCH ased due to uncertain electron-galaxy and matter-galaxy was also supported by a postdoctoral grant from Impe- cross-correlation on small-scales, respectively, the latter rial College London and a Perimeter Visiting Fellowship. 8

SCH thanks Matthew Johnson and Juan Cayuso for use- ment of Canada through the Department of Innovation, ful conversations. MK was supported by NSF Grant Science and Industry Canada and by the Province of On- No. 1818899 and the Simons Foundation. Research at tario through the Ministry of Colleges and Universities. Perimeter Institute is supported in part by the Govern-

[1] Peter Ade, James Aguirre, Zeeshan Ahmed, Simone [14] Zhen Pan and Matthew C. Johnson, “Forecasted con- Aiola, Aamir Ali, David Alonso, Marcelo A. Alvarez, straints on modified gravity from Sunyaev-Zel’dovich Kam Arnold, Peter Ashton, Jason Austermann, and tomography,” Phys. Rev. D 100, 083522 (2019), et al., “The Simons Observatory: science goals and fore- arXiv:1906.04208 [astro-ph.CO]. casts,” J. Cosmology Astropart. Phys. 2019, 056 (2019), [15] Selim C. Hotinli, James B. Mertens, Matthew C. arXiv:1808.07445 [astro-ph.CO]. Johnson, and Marc Kamionkowski, “Probing cor- [2] Adrian Lee, Maximilian H. Abitbol, Shunsuke Adachi, related compensated isocurvature perturbations using Peter Ade, James Aguirre, Zeeshan Ahmed, Simone scale-dependent galaxy bias,” Phys. Rev. D100, 103528 Aiola, Aamir Ali, David Alonso, Marcelo A. Alvarez, (2019), arXiv:1908.08953 [astro-ph.CO]. and et al., “The Simons Observatory,” in Bulletin of the [16] Ya. B. Zeldovich and R. A. Sunyaev, “The Interaction of American Astronomical Society, Vol. 51 (2019) p. 147, Matter and Radiation in a Hot-Model Universe,” Ap&SS arXiv:1907.08284 [astro-ph.IM]. 4, 301–316 (1969). [3] Kevork N. Abazajian et al., “CMB-S4 Science Book, [17] Y. B. Zel’Dovich, “Gravitational instability: an approxi- First Edition,” arXiv e-prints , arXiv:1610.02743 (2016), mate theory for large density perturbations.” A&A 500, arXiv:1610.02743 [astro-ph.CO]. 13–18 (1970). [4] Amir Aghamousa et al. (DESI), “The DESI Experiment [18] R. A. Sunyaev and I. B. Zeldovich, “Microwave back- Part I: Science,Targeting, and Survey Design,” (2016), ground radiation as a probe of the contemporary struc- arXiv:1611.00036 [astro-ph.IM]. ture and history of the universe,” ARA&A 18, 537–560 [5] LSST Science Collaboration et al., “LSST Science Book, (1980). Version 2.0,” ArXiv e-prints (2009), arXiv:0912.0201 [19] R. A. Sunyaev and Ya. B. Zeldovich, “The Observations [astro-ph.IM]. of Relic Radiation as a Test of the Nature of X-Ray Ra- [6] Uroˇs Seljak, “Extracting Primordial Non-Gaussianity diation from the Clusters of Galaxies,” Comments on As- without Cosmic Variance,” Phys. Rev. Lett. 102, 021302 trophysics and Space Physics 4, 173 (1972). (2009), arXiv:0807.1770 [astro-ph]. [20] S. Y. Sazonov and R. A. Sunyaev, “Microwave polar- [7] Pengjie Zhang and Matthew C. Johnson, “Testing eter- ization in the direction of galaxy clusters induced by nal inflation with the kinetic Sunyaev Zel’dovich ef- the CMB quadrupole anisotropy,” MNRAS 310, 765– fect,” J. Cosmology Astropart. Phys. 2015, 046 (2015), 772 (1999), arXiv:astro-ph/9903287 [astro-ph]. arXiv:1501.00511 [astro-ph.CO]. [21] R. K. Sachs and A. M. Wolfe, “Perturbations of a Cosmo- [8] Arka Banerjee, Bhuvnesh Jain, Neal Dalal, and Jessie logical Model and Angular Variations of the Microwave Shelton, “Tests of neutrino and dark radiation models Background,” ApJ 147, 73 (1967). from galaxy and CMB surveys,” J. Cosmology Astropart. [22] M. Birkinshaw and S. F. Gull, “A test for transverse mo- Phys. 2018, 022 (2018), arXiv:1612.07126 [astro-ph.CO]. tions of clusters of galaxies,” Nature 302, 315–317 (1983). [9] Marcel Schmittfull and Uroˇs Seljak, “Parameter con- [23] Anne-Sylvie Deutsch, Matthew C. Johnson, Moritz straints from cross-correlation of CMB lensing with Münchmeyer, and Alexandra Terrana, “Polarized Sun- galaxy clustering,” Phys. Rev. D 97, 123540 (2018), yaev Zel’dovich tomography,” J. Cosmology Astropart. arXiv:1710.09465 [astro-ph.CO]. Phys. 2018, 034 (2018), arXiv:1705.08907 [astro-ph.CO]. [10] Chirag Modi, Martin White, and Zvonimir Vlah, “Mod- [24] Anne-Sylvie Deutsch, Emanuela Dimastrogiovanni, Mat- eling CMB lensing cross correlations with CLEFT,” teo Fasiello, Matthew C. Johnson, and Moritz J. Cosmology Astropart. Phys. 2017, 009 (2017), Münchmeyer, “Primordial gravitational wave phe- arXiv:1706.03173 [astro-ph.CO]. nomenology with polarized Sunyaev Zel’dovich tomogra- [11] Anne-Sylvie Deutsch, Emanuela Dimastrogiovanni, phy,” Phys. Rev. D 100, 083538 (2019), arXiv:1810.09463 Matthew C. Johnson, Moritz Münchmeyer, and Alexan- [astro-ph.CO]. dra Terrana, “Reconstruction of the remote dipole and [25] Kendrick M. Smith, Mathew S. Madhavacheril, Moritz quadrupole fields from the kinetic Sunyaev Zel’dovich Münchmeyer, Simone Ferraro, Utkarsh Giri, and and polarized Sunyaev Zel’dovich effects,” Phys. Rev. D Matthew C. Johnson, “KSZ tomography and the bis- 98, 123501 (2018), arXiv:1707.08129 [astro-ph.CO]. pectrum,” arXiv e-prints , arXiv:1810.13423 (2018), [12] Moritz Münchmeyer, Mathew S. Madhavacheril, Si- arXiv:1810.13423 [astro-ph.CO]. mone Ferraro, Matthew C. Johnson, and Kendrick M. [26] Mathew S. Madhavacheril, Nicholas Battaglia, Smith, “Constraining local non-Gaussianities with ki- Kendrick M. Smith, and Jonathan L. Sievers, “Cosmol- netic Sunyaev-Zel’dovich tomography,” Phys. Rev. D ogy with kSZ: breaking the optical depth degeneracy with 100, 083508 (2019), arXiv:1810.13424 [astro-ph.CO]. Fast Radio Bursts,” arXiv e-prints , arXiv:1901.02418 [13] Juan I. Cayuso and Matthew C. Johnson, “Towards (2019), arXiv:1901.02418 [astro-ph.CO]. testing CMB anomalies using the kinetic and polarized [27] Dagoberto Contreras, Matthew C. Johnson, and Sunyaev-Zel’dovich effects,” Phys. Rev. D 101, 123508 James B. Mertens, “Towards detection of relativistic (2020), arXiv:1904.10981 [astro-ph.CO]. effects in galaxy number counts using kSZ tomogra- 9

phy,” J. Cosmology Astropart. Phys. 2019, 024 (2019), lende Prieto, James Annis, Stephen Bailey, Christophe arXiv:1904.10033 [astro-ph.CO]. Balland, and et al., “The DESI Experiment Part I: [28] Selim C. Hotinli and Matthew C. Johnson, “Re- Science,Targeting, and Survey Design,” arXiv e-prints constructing large scales at cosmic dawn,” (2020), , arXiv:1611.00036 (2016), arXiv:1611.00036 [astro- arXiv:2012.09851 [astro-ph.CO]. ph.IM]. [29] Utkarsh Giri and Kendrick M. Smith, “Exploring KSZ [46] Lam Hui, Enrique Gaztañaga, and Marilena Loverde, velocity reconstruction with N-body simulations and the “Anisotropic magnification distortion of the 3D galaxy halo model,” arXiv e-prints , arXiv:2010.07193 (2020), correlation. I. Real space,” Phys. Rev. D 76, 103502 arXiv:2010.07193 [astro-ph.CO]. (2007), arXiv:0706.1071 [astro-ph]. [30] Nick Kaiser, “Clustering in real space and in redshift [47] David Alonso, Philip Bull, Pedro G. Ferreira, Roy space,” MNRAS 227, 1–21 (1987). Maartens, and Mário G. Santos, “Ultra-large-scale [31] Andrej Obuljen, Will J. Percival, and Neal Dalal, “De- Cosmology in Next-generation Experiments with Single tection of anisotropic galaxy assembly bias in BOSS Tracers,” ApJ 814, 145 (2015), arXiv:1505.07596 [astro- DR12,” J. Cosmology Astropart. Phys. 2020, 058 (2020), ph.CO]. arXiv:2004.07240 [astro-ph.CO]. [48] Kevork Abazajian et al., “CMB-S4 Science Case, Ref- [32] Selim C. Hotinli, Joel Meyers, Neal Dalal, Andrew H. erence Design, and Project Plan,” arXiv e-prints Jaffe, Matthew C. Johnson, James B. Mertens, Moritz , arXiv:1907.04473 (2019), arXiv:1907.04473 [astro- Münchmeyer, Kendrick M. Smith, and Alexander ph.IM]. van Engelen, “Transverse Velocities with the Moving [49] Mathew S. Madhavacheril, Nicholas Battaglia, and Lens Effect,” Phys. Rev. Lett. 123, 061301 (2019), Hironao Miyatake, “Fundamental physics from fu- arXiv:1812.03167 [astro-ph.CO]. ture weak-lensing calibrated Sunyaev-Zel’dovich galaxy [33] Siavash Yasini, Nareg Mirzatuny, and Elena Pier- cluster counts,” Phys. Rev. D 96, 103525 (2017), paoli, “Pairwise Transverse Velocity Measurement with arXiv:1708.07502 [astro-ph.CO]. the Rees-Sciama Effect,” ApJ 873, L23 (2019), [50] J. Dunkley, E. Calabrese, J. Sievers, et al., “The At- arXiv:1812.04241 [astro-ph.CO]. acama Cosmology Telescope: likelihood for small-scale [34] Selim C. Hotinli, Matthew C. Johnson, and Joel Meyers, CMB data,” J. Cosmology Astropart. Phys. 2013, 025 “Optimal filters for the moving lens effect,” Phys. Rev. D (2013), arXiv:1301.0776 [astro-ph.CO]. 103, 043536 (2021), arXiv:2006.03060 [astro-ph.CO]. [51] G. Lagache, M. Béthermin,L. Montier, P. Serra, and [35] Eric V. Linder, “Cosmic growth history and expansion M. Tucci, “Impact of polarised extragalactic sources on history,” Phys. Rev. D 72, 043529 (2005), arXiv:astro- the measurement of CMB B-mode anisotropies,” A&A ph/0507263 [astro-ph]. 642, A232 (2020), arXiv:1911.09466 [astro-ph.CO]. [36] Eric V. Linder and Robert N. Cahn, “Parameterized [52] Antony Lewis, Anthony Challinor, and Anthony beyond-Einstein growth,” Astroparticle Physics 28, 481– Lasenby, “Efficient Computation of Cosmic Mi- 488 (2007), arXiv:astro-ph/0701317 [astro-ph]. crowave Background Anisotropies in Closed Friedmann- [37] J. R. Bond and A. S. Szalay, “The collisionless damping Robertson-Walker Models,” ApJ 538, 473–476 (2000), of density fluctuations in an expanding universe,” ApJ arXiv:astro-ph/9911177 [astro-ph]. 274, 443–468 (1983). [53] Hyunbae Park, Paul R. Shapiro, Eiichiro Komatsu, Il- [38] N. Battaglia, “The tau of galaxy clusters,” J. Cosmol- ian T. Iliev, Kyungjin Ahn, and Garrelt Mellema, ogy Astropart. Phys. 2016, 058 (2016), arXiv:1607.02442 “The Kinetic Sunyaev-Zel’dovich Effect as a Probe [astro-ph.CO]. of the Physics of Cosmic Reionization: The Effect [39] Juan Cayuso, Richard Bloch, Selim C. Hotinli, and of Self-regulated Reionization,” ApJ 769, 93 (2013), Matthew Johnson, “Upcoming,” . arXiv:1301.3607 [astro-ph.CO]. [40] C. Alcock and B. Paczynski, “An evolution free test for non-zero cosmological constant,” Nature 281, 358 (1979). [41] Christopher M. Hirata, “Tidal alignments as a contami- Appendix A: Internal Linear Combination noise nant of redshift space distortions,” MNRAS 399, 1074– curves 1087 (2009), arXiv:0903.4929 [astro-ph.CO]. [42] Franco D. Albareti et al., “The 13th Data Release of the Sloan Digital Sky Survey: First Spectroscopic Data We describe here the procedure used to analytically from the SDSS-IV Survey Mapping Nearby Galaxies calculate the effective total noise in CMB maps for the Si- at Apache Point Observatory,” ApJS 233, 25 (2017), mons Observatory (SO), CMB-S4 and CMB-HD surveys, arXiv:1608.02013 [astro-ph.GA]. guided by the Internal Linear Combination (ILC) algo- [43] A. P. Hearin, A. R. Zentner, and Z. Ma, “General re- rithm that would be used in a realistic analysis. Our ap- quirements on matter power spectrum predictions for proach here is approximate and the resulting noise curves cosmology with weak lensing tomography,” J. Cosmol- are only meant to be roughly representative of the fore- ogy Astropart. Phys. 4, 034 (2012), arXiv:1111.0052. cast performance of these surveys; we do however find [44] The LSST Dark Energy Science Collaboration, Rachel 10-20% agreement of the SO and CMB-S4 results with Mandelbaum, Tim Eifler, RenéeHloˇzek,Thomas Col- lett, Eric Gawiser, and others Scolnic, “The LSST Dark the official noise curves in [1] and [48], respectively. In Energy Science Collaboration (DESC) Science Require- particular, for SO, we find better than 10% agreement ments Document,” arXiv e-prints , arXiv:1809.01669 with the official forecasts for almost all relevant multi- (2018), arXiv:1809.01669 [astro-ph.CO]. poles. We employ this analytic approach (as opposed to [45] DESI Collaboration, Amir Aghamousa, Jessica Aguilar, the simulation-based work in [1] and [48]) so as to be able Steve Ahlen, Shadab Alam, Lori E. Allen, Carlos Al- to easily propagate the results of varying the white noise 10

Beam FWHM Noise RMS include contributions from tSZ, clustered CIB and con- (µK-arcmin) servative (large) levels of Poisson CIB foregrounds at SO S4 HD SO S4 HD these frequencies as well as the black-body late-time kSZ 39 GHz 5.10 5.10 36.300 36 12.4 3.4 following the approach in [49], which is based on fits to 93 GHz 2.20 2.20 15.300 8 2.0 0.6 ACT data from [50], but we do not include the tSZ-CIB 145 GHz 1.40 1.40 10.000 10 2.0 0.6 correlation. We include radio sources in the 39, 93 and 225 GHz 1.00 1.00 6.600 22 6.9 1.9 145 GHz channels using the flux-limit-dependent radio 280 GHz 0.90 0.90 5.400 54 16.7 4.6 source power model from [51], where for both SO and S4 we assume flux limits in those channels of 10, 7 and 10 TABLE II. Inputs to ILC noise: The beam and noise RMS mJy respectively. For HD, we assume lower flux limits parameters we assume for survey configurations roughly cor- of 2, 1 and 1 mJy respectively. In addition, we include responding to Simons Observatory (baseline), CMB-S4 and the black-body lensed CMB contribution calculated us- CMB-HD. We do not include the 30 GHz channel and do not include atmospheric noise since the kSZ and moving lens ing CAMB [52] as well as the reionization kSZ signal from information is primarily in scales ` > 2000. [53]. We do not include the contribution from the moving lens signal itself since it is sub-dominant to the sum of the above. RMS in the CMB experiments. We assume that each of these surveys measures the Given the above sky model involving signal compo- millimeter sky at 39, 93, 145, 225 and 280 GHz with nents correlated across frequencies indexed by i (and un- beam FWHM and white noise RMS as shown in Table correlated beam-deconvolved noise), the total covariance II. We do not include large-scale atmospheric noise but between two frequency channels at each multipole ` is do not expect this to make a difference for our forecasts ij C` . The final minimum-variance (standard ILC) noise given the availability of Planck data at 44, 100, 143 and for the black-body CMB+kSZ signal that includes all 217 GHz as well as the fact that most of the kSZ and contributions (including the sample variance in the CMB moving lens information is derived from small scales with −1 hP −1iji ` > 2000. and kSZ) is then given by N` = ij C ` . These For the signal contributions in the millimeter sky, we curves are shown in Figure1.