Is Lensing Amplitude Anomaly in the Cosmic Microwave Background the Evidence of Extremely Low Frequency Primordial Gravitational Wave?

Wenshuai Liu1 674602871@qq.com 1School of Physics, Henan Normal University, Xinxiang 453007, China
(June 19, 2024)
Abstract

Trajectories of photons of cosmic microwave background (CMB) from the surface of last scattering to us could be deflected by extremely low frequency primordial gravitational wave (PGW). With large scale structure (LSS) producing a smoothing of the acoustic peaks in the power spectrum of the CMB anisotropies through weak lensing, the presence of extremely low frequency PGW could enhance the effect of weak lensing on CMB due to the coupling of extremely low frequency PGW and LSS, thus, give rise to much more smoothing of the spectrum. This may be an natural explanation for the lensing amplitude anomaly observed by Planck, meaning that lensing amplitude anomaly may be the evidence of extremely low frequency PGW.

pacs:
98.80.-k,98.62.En

I Introduction

Inflation, which gives rise to not only a flat, homogeneous, and isotropic Universe but also seed perturbations growing to create large scale structure in the Universe (Grishchuk, 1976, 1977; Starobinsky, 1980; Linde, 1982), predicts the existence of primordial gravitational waves (PGWs) with a nearly scale invariant spectrum (Abbott & Wise, 1984; Starobinskii, 1985; Rubakov et al., 1982; Fabbri & Pollock, 1983; Starobinsky, 1979; Sahni, 1990; Allen, 1988). The detection of PGWs would confirm the inflationary scenario and determine the energy scale, thus, it is of great significance to find tools to detect PGWs. The traditional method to detect extremely low frequency PGWs with range of 1018superscript101810^{-18}10 start_POSTSUPERSCRIPT - 18 end_POSTSUPERSCRIPTHz-1016superscript101610^{-16}10 start_POSTSUPERSCRIPT - 16 end_POSTSUPERSCRIPTHz is the B-modes of polarization of cosmic microwave background (CMB) (Kamionkowski et al., 1997; Seljak & Zaldarriaga, 1997) induced by extremely low frequency PGWs through Thomson scattering. From observation, it shows that dust in our Milky Way can make foreground contamination, leading to challenges with the method of B-modes. Thus, finding an alternative method of detection is worthy of study. When the frequency of PGW is less than 1018superscript101810^{-18}10 start_POSTSUPERSCRIPT - 18 end_POSTSUPERSCRIPTHz, the wavelength can be larger than the horizon of the Universe. Here, PGWs with frequency less than 1018superscript101810^{-18}10 start_POSTSUPERSCRIPT - 18 end_POSTSUPERSCRIPTHz are also called as extremely low frequency PGWs.

A new method proposed by Liu Liu (2022) uses gravitational lens system with a non-aligned source-detector-observer configuration to detect extremely low frequency PGWs, which shows that the time delay in such gravitational lens system with perturbation from extremely low frequency PGWs could deviate from the one deduced from the theoretical model obviously. This means that such gravitational lens system could be used as a possible long-baseline detector of extremely low frequency PGWs. However, the deflector adopted in Liu (2022) is a point mass model which is unrealistic in strong lensing. Work of Liu Liu (2023) adopts gravitational lens system with singular isothermal sphere deflector. It shows from Liu Liu (2023) that, under the perturbation of extremely low frequency PGWs with arbitrary direction of propagation, the time delay from observation can strongly deviate from that deduced from the theory.

It shows from Liu Liu (2023) that extremely low frequency PGWs can affect gravitational lensing through perturbation on the trajectory of photon no matter what the gravitational lensing is strong lening or weak lensing, meaning that extremely low frequency PGWs could also have effect on the weak lensing of CMB through the coupling of extremely low frequency PGWs and LSS.

The CMB radiation undergoes weak lensing shown in Figure 1 by LSS along path since last scattering surface and the effect of weak lensing on CMB is a remapping of the CMB temperature anisotropies and the CMB polarization field, leading to a smoothing of the acoustic peaks in the power spectrum of the CMB anisotropies Seljak (1996). With the presence of extremely low frequency PGWs, the effect of weak lensing on CMB may be enhanced by the coupling of extremely low frequency PGW and LSS, thus, it may give rise to much more smoothing of the spectrum. Results from Planck Aghanim et al. (2020) show that the lensing-reconstruction power spectrum shows to be almost consistent with that expected for ΛΛ\Lambdaroman_ΛCDM models which fit the CMB spectra. However, the amount of smoothing observed in the CMB angular power spectrum is larger than the amount of the smoothing derived from the power spectrum of the reconstructed lensing potential, resulting the lensing parameter AL>1subscript𝐴𝐿1A_{L}>1italic_A start_POSTSUBSCRIPT italic_L end_POSTSUBSCRIPT > 1. The observed AL>1subscript𝐴𝐿1A_{L}>1italic_A start_POSTSUBSCRIPT italic_L end_POSTSUBSCRIPT > 1 or the observed smoothing which is larger than expected may be the evidence of extremely low frequency PGWs.

In addition to the lensing amplitude anomaly observed in CMB, other anomalies like Hubble tension, Ωk0subscriptΩ𝑘0\Omega_{k}\neq 0roman_Ω start_POSTSUBSCRIPT italic_k end_POSTSUBSCRIPT ≠ 0 and S8subscript𝑆8S_{8}italic_S start_POSTSUBSCRIPT 8 end_POSTSUBSCRIPT tension may be also the evidence of extremely low frequency PGWs since H0subscript𝐻0H_{0}italic_H start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT, ΩKsubscriptΩ𝐾\Omega_{K}roman_Ω start_POSTSUBSCRIPT italic_K end_POSTSUBSCRIPT and S8subscript𝑆8S_{8}italic_S start_POSTSUBSCRIPT 8 end_POSTSUBSCRIPT determined based on CMB can be affected by the weak lensing through the coupling of extremely low frequency PGWs and LSS.

II Weak lensing by the coupling of PGWs and LSS

After inflation, PGWs can be stretched to scale larger than the horizon, this makes PGWs with wavelength larger than the horizon not behave like waves but have an almost frozen amplitude, meaning that such PGWs like a constant anisotropy or shear. Such PGWs start oscillating inside the horizon after they re-enter the horizon, and the amplitude are damped due to expansion of the universe.

In order to only investigate the effect of extremely low frequency PGWs on the photon of CMB from the last scattering surface, we first give the metric of extremely low frequency PGWs as

hij=a0a[(uiujvivj)h++(uivj+viuj)h×]×cos(ωη𝐤𝐱)subscript𝑖𝑗subscript𝑎0𝑎delimited-[]subscript𝑢𝑖subscript𝑢𝑗subscript𝑣𝑖subscript𝑣𝑗subscriptsubscript𝑢𝑖subscript𝑣𝑗subscript𝑣𝑖subscript𝑢𝑗subscript𝜔𝜂𝐤𝐱h_{ij}=\frac{a_{0}}{a}[(u_{i}u_{j}-v_{i}v_{j})h_{+}+(u_{i}v_{j}+v_{i}u_{j})h_{% \times}]\times\cos(\omega\eta-\mathbf{k}\cdot\mathbf{x})italic_h start_POSTSUBSCRIPT italic_i italic_j end_POSTSUBSCRIPT = divide start_ARG italic_a start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT end_ARG start_ARG italic_a end_ARG [ ( italic_u start_POSTSUBSCRIPT italic_i end_POSTSUBSCRIPT italic_u start_POSTSUBSCRIPT italic_j end_POSTSUBSCRIPT - italic_v start_POSTSUBSCRIPT italic_i end_POSTSUBSCRIPT italic_v start_POSTSUBSCRIPT italic_j end_POSTSUBSCRIPT ) italic_h start_POSTSUBSCRIPT + end_POSTSUBSCRIPT + ( italic_u start_POSTSUBSCRIPT italic_i end_POSTSUBSCRIPT italic_v start_POSTSUBSCRIPT italic_j end_POSTSUBSCRIPT + italic_v start_POSTSUBSCRIPT italic_i end_POSTSUBSCRIPT italic_u start_POSTSUBSCRIPT italic_j end_POSTSUBSCRIPT ) italic_h start_POSTSUBSCRIPT × end_POSTSUBSCRIPT ] × roman_cos ( italic_ω italic_η - bold_k ⋅ bold_x ) (1)

where a𝑎aitalic_a is the scale factor and the present value is set to be a0=1subscript𝑎01a_{0}=1italic_a start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT = 1. The conformal time η=te+(z+L)𝜂subscript𝑡𝑒𝑧𝐿\eta=t_{e}+(z+L)italic_η = italic_t start_POSTSUBSCRIPT italic_e end_POSTSUBSCRIPT + ( italic_z + italic_L ) in order to approach the level of approximation. z𝑧zitalic_z is in form of conformal distance and conformal time. L=14000𝐿14000L=14000italic_L = 14000Gpc is the conformal distance between the last scattering surface and us. tesubscript𝑡𝑒t_{e}italic_t start_POSTSUBSCRIPT italic_e end_POSTSUBSCRIPT is the time when the photon was emitted at (x=0,y=0,z=L)formulae-sequence𝑥0formulae-sequence𝑦0𝑧𝐿(x=0,y=0,z=-L)( italic_x = 0 , italic_y = 0 , italic_z = - italic_L ) so that ωte𝜔subscript𝑡𝑒\omega t_{e}italic_ω italic_t start_POSTSUBSCRIPT italic_e end_POSTSUBSCRIPT acts as the initial phase, 𝐤=ω(sinθcosϕ,sinθsinϕ,cosθ)𝐤𝜔𝜃italic-ϕ𝜃italic-ϕ𝜃\mathbf{k}=\omega(\sin\theta\cos\phi,\sin\theta\sin\phi,\cos\theta)bold_k = italic_ω ( roman_sin italic_θ roman_cos italic_ϕ , roman_sin italic_θ roman_sin italic_ϕ , roman_cos italic_θ ) is the propagation vector, 𝐮=(sinϕ,cosϕ,0)𝐮italic-ϕitalic-ϕ0\mathbf{u}=(\sin\phi,-\cos\phi,0)bold_u = ( roman_sin italic_ϕ , - roman_cos italic_ϕ , 0 ), 𝐯=(cosθcosϕ,cosθsinϕ,sinθ)𝐯𝜃italic-ϕ𝜃italic-ϕ𝜃\mathbf{v}=(\cos\theta\cos\phi,\cos\theta\sin\phi,-\sin\theta)bold_v = ( roman_cos italic_θ roman_cos italic_ϕ , roman_cos italic_θ roman_sin italic_ϕ , - roman_sin italic_θ ), ω=2πf𝜔2𝜋𝑓\omega=2\pi fitalic_ω = 2 italic_π italic_f, f𝑓fitalic_f is the frequency of gravitational wave at present, h+subscripth_{+}italic_h start_POSTSUBSCRIPT + end_POSTSUBSCRIPT and h×subscripth_{\times}italic_h start_POSTSUBSCRIPT × end_POSTSUBSCRIPT are the amplitude of the two polarizations of the gravitational wave at present, respectively.

In Eq. (1), when the wavelength of PGW is larger than the horizon of the universe, it is frozen and does’t evolve, resulting η=te𝜂subscript𝑡𝑒\eta=t_{e}italic_η = italic_t start_POSTSUBSCRIPT italic_e end_POSTSUBSCRIPT in Eq. (1). Only when the wavelength of PGW is smaller than the horizon of the universe can η=te+(z+L)𝜂subscript𝑡𝑒𝑧𝐿\eta=t_{e}+(z+L)italic_η = italic_t start_POSTSUBSCRIPT italic_e end_POSTSUBSCRIPT + ( italic_z + italic_L ).

Refer to caption
Figure 1: The orange curve represents the last scattering surface and the black dot is the observer. The blue line shows the trajectory of CMB photon without perturbation of LSS and PGWs, and the black curve is the trajectory of CMB photon with perturbation of LSS and PGWs. The blue dashed line is the direction of CMB photon the observer sees.
Refer to caption
Figure 2: The orange line represents a small fraction of the last scattering surface, the big black dot is a cluster(In reality, there are many clusters between the the last scattering surface and the observer along line of sight. Here, we set one cluster in order to describe the effect of PGWs on lensing) and the small dot is the observer. The black line represents the trajectory of CMB photon with perturbation of the cluster but without perturbation of PGWs. The blue curve is the trajectory of CMB photon with perturbation of the cluster and PGWs. Due to the fact that PGW could deflect the trajectory of the photon, the closest distance between the photon and the cluster may be smaller with perturbation of PGWs compared with that without perturbation of PGWs, thus, the deflection by the cluster can be larger (α2>α1subscript𝛼2subscript𝛼1\alpha_{2}>\alpha_{1}italic_α start_POSTSUBSCRIPT 2 end_POSTSUBSCRIPT > italic_α start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT), resulting θ2>θ1subscript𝜃2subscript𝜃1\theta_{2}>\theta_{1}italic_θ start_POSTSUBSCRIPT 2 end_POSTSUBSCRIPT > italic_θ start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT.

To show the deflection of the photon of CMB due to extremely low frequency PGWs, the method to calculate the deflection angle is similar to the method of Liu Liu (2023) to calculate the angular position of the image of the deflector in the gravitational lens system. When we consider PGW with wavelength larger than (and almost equal to) the horizon of the universe, the frequencies we adopt are f=1019𝑓superscript1019f=10^{-19}italic_f = 10 start_POSTSUPERSCRIPT - 19 end_POSTSUPERSCRIPT, f=5×1019𝑓5superscript1019f=5\times 10^{-19}italic_f = 5 × 10 start_POSTSUPERSCRIPT - 19 end_POSTSUPERSCRIPT, and f=1018𝑓superscript1018f=10^{-18}italic_f = 10 start_POSTSUPERSCRIPT - 18 end_POSTSUPERSCRIPT, and we set the tensor to scalar ratio to be r=0.01𝑟0.01r=0.01italic_r = 0.01, r=0.005𝑟0.005r=0.005italic_r = 0.005, 0.0010.0010.0010.001 to calculate the root-mean-square of deflection angles. For r=0.01𝑟0.01r=0.01italic_r = 0.01, the root-mean-square of deflection angles are 6.93423×1076.93423superscript1076.93423\times 10^{-7}6.93423 × 10 start_POSTSUPERSCRIPT - 7 end_POSTSUPERSCRIPT, 3.38361×1063.38361superscript1063.38361\times 10^{-6}3.38361 × 10 start_POSTSUPERSCRIPT - 6 end_POSTSUPERSCRIPT and 6.15404×1066.15404superscript1066.15404\times 10^{-6}6.15404 × 10 start_POSTSUPERSCRIPT - 6 end_POSTSUPERSCRIPT when f=1019𝑓superscript1019f=10^{-19}italic_f = 10 start_POSTSUPERSCRIPT - 19 end_POSTSUPERSCRIPT, f=5×1019𝑓5superscript1019f=5\times 10^{-19}italic_f = 5 × 10 start_POSTSUPERSCRIPT - 19 end_POSTSUPERSCRIPT, and f=1018𝑓superscript1018f=10^{-18}italic_f = 10 start_POSTSUPERSCRIPT - 18 end_POSTSUPERSCRIPT. For r=0.005𝑟0.005r=0.005italic_r = 0.005, they are 4.90326×1074.90326superscript1074.90326\times 10^{-7}4.90326 × 10 start_POSTSUPERSCRIPT - 7 end_POSTSUPERSCRIPT, 2.39258×1062.39258superscript1062.39258\times 10^{-6}2.39258 × 10 start_POSTSUPERSCRIPT - 6 end_POSTSUPERSCRIPT and 4.35158×1064.35158superscript1064.35158\times 10^{-6}4.35158 × 10 start_POSTSUPERSCRIPT - 6 end_POSTSUPERSCRIPT when f=1019𝑓superscript1019f=10^{-19}italic_f = 10 start_POSTSUPERSCRIPT - 19 end_POSTSUPERSCRIPT, f=5×1019𝑓5superscript1019f=5\times 10^{-19}italic_f = 5 × 10 start_POSTSUPERSCRIPT - 19 end_POSTSUPERSCRIPT, and f=1018𝑓superscript1018f=10^{-18}italic_f = 10 start_POSTSUPERSCRIPT - 18 end_POSTSUPERSCRIPT. For 0.0010.0010.0010.001, they are 2.19273×1072.19273superscript1072.19273\times 10^{-7}2.19273 × 10 start_POSTSUPERSCRIPT - 7 end_POSTSUPERSCRIPT, 1.06996×1061.06996superscript1061.06996\times 10^{-6}1.06996 × 10 start_POSTSUPERSCRIPT - 6 end_POSTSUPERSCRIPT and 1.946023×1061.946023superscript1061.946023\times 10^{-6}1.946023 × 10 start_POSTSUPERSCRIPT - 6 end_POSTSUPERSCRIPT when f=1019𝑓superscript1019f=10^{-19}italic_f = 10 start_POSTSUPERSCRIPT - 19 end_POSTSUPERSCRIPT, f=5×1019𝑓5superscript1019f=5\times 10^{-19}italic_f = 5 × 10 start_POSTSUPERSCRIPT - 19 end_POSTSUPERSCRIPT, and f=1018𝑓superscript1018f=10^{-18}italic_f = 10 start_POSTSUPERSCRIPT - 18 end_POSTSUPERSCRIPT.

According to the above, it shows that extremely low frequency PGWs could strongly deflect the photon of CMB. Although the deflection of CMB due to LSS is about 2 arcminutes which is larger than the above, the coupling of extremely low frequency PGW and LSS may enhance the weak lensing shown in Figure 2 and the deflection induced by the coupling of extremely low frequency PGW and LSS could be larger than linear superposition of the deflection induced by LSS and that induced by extremely low frequency PGWs, resulting the observed smoothing of the spectrum may be larger than that derived from theory.

We only give a qualitative analysis of the effect of extremely low frequency PGWs on weak lensing, the detailed quantitative calculation should be conducted to confirm if extremely low frequency PGWs could explain the lensing amplitude anomaly. The detailed confirmation can be as follows.

Using the coupling of density perturbation derived from ΛΛ\Lambdaroman_ΛCDM and tensor perturbation determined by scalar to tensor ratio r𝑟ritalic_r to numerically generate the CMB map and calculate the CMB temperature power spectrum, we can infer the scalar to tensor ratio r𝑟ritalic_r by setting the calculated CMB temperature power spectrum to be consistent with the observed temperature power spectrum. If the inferred r𝑟ritalic_r is larger than the observed value, the lensing amplitude anomaly may be not due to PGWs. If the derived r𝑟ritalic_r is smaller than the observed value, the lensing potential should be reconstructed based on the generated CMB map in order to get the CMB temperature power spectrum and the amount of smoothing based on the power spectrum of the reconstructed lensing potential. If the resulting amount of smoothing is consistent with the amount of smoothing derived from the power spectrum of the reconstructed lensing potential in Planck Aghanim et al. (2020) and if the power spectrum of the reconstructed lensing potential based on the generated CMB map is consistent with that reconstructed in Planck Aghanim et al. (2020), it means that extremely low frequency PGWs may account for the lensing amplitude anomaly.

III DISCUSSION

The effect of extremely low frequency PGWs on weak lensing of CMB by LSS is discussed qualitatively. It shows that the effect of weak lensing on CMB may be enhanced by the coupling of extremely low frequency PGWs and LSS if extremely low frequency PGWs exist. Since weak lensing by LSS can lead to a smoothing of the acoustic peaks in the angular power spectrum of the CMB anisotropies, thus, weak lensing by the coupling of extremely low frequency PGWs and LSS may give rise to much more smoothing of the spectrum.

Observations by Planck Aghanim et al. (2020) show many anomalies such as lensing amplitude anomaly observed in CMB, Hubble tension, ΩK0subscriptΩ𝐾0\Omega_{K}\neq 0roman_Ω start_POSTSUBSCRIPT italic_K end_POSTSUBSCRIPT ≠ 0 and S8subscript𝑆8S_{8}italic_S start_POSTSUBSCRIPT 8 end_POSTSUBSCRIPT tension. From the point of view of weak lensing of CMB, the presence of extremely low frequency PGWs may account for the lensing amplitude anomaly since the coupling of extremely low frequency PGWs and LSS could enhance the effect of weak lensing of CMB. Due to the fact that H0subscript𝐻0H_{0}italic_H start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT, ΩKsubscriptΩ𝐾\Omega_{K}roman_Ω start_POSTSUBSCRIPT italic_K end_POSTSUBSCRIPT and S8subscript𝑆8S_{8}italic_S start_POSTSUBSCRIPT 8 end_POSTSUBSCRIPT determined based on CMB can be affected by the coupling of extremely low frequency PGWs and LSS, anomalies like Hubble tension, ΩK0subscriptΩ𝐾0\Omega_{K}\neq 0roman_Ω start_POSTSUBSCRIPT italic_K end_POSTSUBSCRIPT ≠ 0 and S8subscript𝑆8S_{8}italic_S start_POSTSUBSCRIPT 8 end_POSTSUBSCRIPT tension may also reflect the existence of extremely low frequency PGWs.

In this work, we give a qualitative investigation on how the coupling of extremely low frequency PGWs and LSS may enhance the effect of weak lensing on CMB and wonder if the lensing amplitude anomaly might be the evidence of extremely low frequency PGWs. A detailed quantitative verification on extremely low frequency PGWs accounting for lensing amplitude anomaly and other anomalies like Hubble tension, ΩK0subscriptΩ𝐾0\Omega_{K}\neq 0roman_Ω start_POSTSUBSCRIPT italic_K end_POSTSUBSCRIPT ≠ 0 and S8subscript𝑆8S_{8}italic_S start_POSTSUBSCRIPT 8 end_POSTSUBSCRIPT tension will be conducted in the near future.

References

  • Grishchuk (1976) Grishchuk, L. P. 1976, JETP Lett. 23, 293
  • Grishchuk (1977) Grishchuk, L. P. 1977, Sov. Phys. Usp. 20, 319
  • Starobinsky (1980) Starobinsky, A. A. 1980, PLB 91, 99
  • Linde (1982) Linde, A. D. 1982, PLB 108, 389
  • Abbott & Wise (1984) Abbott, L. F., & Wise, M. B. 1984, Nucl. Phys. B244, 541
  • Starobinskii (1985) Starobinskii, A. 1985, Sov. Astron. Lett. 11, 133
  • Rubakov et al. (1982) Rubakov, V. A., Sazhin, M.V., & Veryaskin, A.V. 1982, Phys. Lett. 115B, 189
  • Fabbri & Pollock (1983) Fabbri, R., & Pollock, M. D. 1983, Phys. Lett. 125B, 445
  • Starobinsky (1979) Starobinsky, A. A. 1979, JETP Lett. 30, 682
  • Sahni (1990) Sahni, V. 1990, Phys. Rev. D42, 453
  • Allen (1988) Allen, B. 1988, Phys. Rev. D37, 2078
  • Kamionkowski et al. (1997) Kamionkowski, M., Kosowsky, A., & Stebbins, A. 1997, Phys. Rev. Lett. 78, 2058
  • Seljak & Zaldarriaga (1997) Seljak, U., & Zaldarriaga, M. 1997, Phys. Rev. Lett. 78, 2054
  • Liu (2022) Liu, W., 2022, MNRAS, 517, 2769
  • Liu (2023) Liu, W., 2023, arXiv:2212.02221
  • Seljak (1996) Seljak, U., 1996, ApJ, 463, 1
  • Aghanim et al. (2020) Aghanim, N., Planck Collaboration VI. 2020, A&A, 641, A6