Exact CMB B-mode power spectrum from anisotropic cosmic birefringence

Toshiya Namikawa Center for Data-Driven Discovery, Kavli IPMU (WPI), UTIAS, The University of Tokyo, Kashiwa, 277-8583, Japan
(May 1, 2024)
Abstract

We calculate the cosmic microwave background (CMB) B𝐵Bitalic_B-mode power spectrum resulting from anisotropic cosmic birefringence, without relying on the thin approximation of the last scattering surface. Specifically, we consider the influence of anisotropic cosmic birefringence arising from massless axion-like particles. Comparing our results to those obtained using the thin approximation, we observe a suppression in the amplitude of the B𝐵Bitalic_B-mode power spectrum by approximately an order of magnitude at large angular scales (10less-than-or-similar-to10\ell\lesssim 10roman_ℓ ≲ 10) and by a factor of two at small angular scales (100greater-than-or-equivalent-to100\ell\gtrsim 100roman_ℓ ≳ 100) when not employing the thin approximation. We also constrain the amplitude of the angular power spectrum of the scale-invariant anisotropic cosmic birefringence using the SPTpol B𝐵Bitalic_B-mode power spectrum. We find that the amplitude is constrained as ACB×104=1.030.97+0.91(2σ)subscript𝐴CBsuperscript104subscriptsuperscript1.030.910.972𝜎A_{\rm CB}\times 10^{4}=1.03^{+0.91}_{-0.97}\,(2\,\sigma)italic_A start_POSTSUBSCRIPT roman_CB end_POSTSUBSCRIPT × 10 start_POSTSUPERSCRIPT 4 end_POSTSUPERSCRIPT = 1.03 start_POSTSUPERSCRIPT + 0.91 end_POSTSUPERSCRIPT start_POSTSUBSCRIPT - 0.97 end_POSTSUBSCRIPT ( 2 italic_σ ). The numerical code is publicly available here.

I Introduction

Cosmic birefringence — a rotation of the linear polarization plane of the cosmic microwave background (CMB) as they travel through space [1] — is a key observational effect on CMB as it provides a way to search for parity-violating physics in cosmology [2]. Recent analyses of the cross-correlation between the even-parity E𝐸Eitalic_E-modes and odd-parity B𝐵Bitalic_B-modes in the Planck polarization map suggest a tantalizing hint of cosmic birefringence [3, 4, 5, 6, 7]. Additionally, a polarization measurement of the Crab Nebula by POLARBEAR suggests a nonzero cosmic birefringence [8]. The phenomenon can be induced by a pseudoscalar field, such as axionlike particles (ALPs), coupled with electromagnetic fields through the Chern-Simons term, gϕγϕFμνF~μν/4subscript𝑔italic-ϕ𝛾italic-ϕsuperscript𝐹𝜇𝜈~𝐹𝜇𝜈4\mathcal{L}\supset-g_{\phi\gamma}\phi F^{\mu\nu}\tilde{F}{\mu\nu}/4caligraphic_L ⊃ - italic_g start_POSTSUBSCRIPT italic_ϕ italic_γ end_POSTSUBSCRIPT italic_ϕ italic_F start_POSTSUPERSCRIPT italic_μ italic_ν end_POSTSUPERSCRIPT over~ start_ARG italic_F end_ARG italic_μ italic_ν / 4, where gϕγ𝑔italic-ϕ𝛾g{\phi\gamma}italic_g italic_ϕ italic_γ represents the coupling constant, ϕitalic-ϕ\phiitalic_ϕ denotes an ALP field, Fμνsuperscript𝐹𝜇𝜈F^{\mu\nu}italic_F start_POSTSUPERSCRIPT italic_μ italic_ν end_POSTSUPERSCRIPT stands for the electromagnetic field tensor, and F~μνsubscript~𝐹𝜇𝜈\tilde{F}_{\mu\nu}over~ start_ARG italic_F end_ARG start_POSTSUBSCRIPT italic_μ italic_ν end_POSTSUBSCRIPT is its dual.

Multiple works have explored cosmic birefringence induced by the ALP field of dark energy [9, 10, 11, 12, 13, 14, 15, 16], early dark energy [12, 17, 18], dark matter [19, 20, 21], and by topological defects [22, 23, 24, 25], as well as by potential signatures of quantum gravity [26, 27]. Anticipated advancements in CMB experiments, including BICEP [28, 29], Simons Array [30], Simons Observatory [31], CMB-S4 [32], and LiteBIRD [33], are projected to substantially reduce polarization noise, thus enhancing sensitivity to cosmic birefringence signals.

Several recent studies have shown that the time evolution of the pseudoscalar fields during recombination and reionization significantly changes the EB𝐸𝐵EBitalic_E italic_B power spectrum [19, 34, 35, 36, 17, 18, 37, 38, 39]. Measuring the spectral shape of the EB𝐸𝐵EBitalic_E italic_B power spectrum will provide tomographic information on the pseudoscalar fields. We can also constrain the late-time evolution by observing the polarized Sunyaev-Zel’dovich effect [40, 41] and galaxies [42, 43]. This tomographic method can avoid the degeneracies with the instrumental miscalibration angle [44, 45, 46, 47] and half-wave plate nonidealities [48].

If the ALP fields fluctuate, δϕ𝛿italic-ϕ\delta\phiitalic_δ italic_ϕ, which has a spatial variation, the polarization rotation will be anisotropic (i.e., have different values in different directions in the sky); indeed, anisotropies in the cosmic birefringence are produced by many of the types of beyond-the-Standard-Model physics (see, e.g., [9, 49, 50, 51, 52, 53, 54]). In this paper, we investigate the impact of the time evolution of δϕ𝛿italic-ϕ\delta\phiitalic_δ italic_ϕ on the calculation of anisotropic cosmic birefringence. Specifically, we focus on the CMB B𝐵Bitalic_B-mode power spectrum which is most sensitive to anisotropic cosmic birefringence among the CMB two-point correlations.

The B𝐵Bitalic_B-mode power spectrum arising from anisotropic cosmic birefringence has often been computed using the approximation of the thin last-scattering surface (LSS). This approximation ignores the time evolution of δϕ𝛿italic-ϕ\delta\phiitalic_δ italic_ϕ during the recombination (and reionization) epoch to simplify the B𝐵Bitalic_B-mode power spectrum [55, 56, 57, 51, 58, 20, 59, 60, 61]. Our main objective is to demonstrate the significance of considering the thickness of the LSS when computing the B𝐵Bitalic_B-mode power spectrum. It’s worth noting that only a few studies have incorporated the time evolution of δϕ𝛿italic-ϕ\delta\phiitalic_δ italic_ϕ to derive the B𝐵Bitalic_B-mode power spectrum without the thin LSS approximation. Pospelov et al. (2009) [62] (hereafter, P09) derived the expression for the B𝐵Bitalic_B-mode power spectrum from time-varying anisotropic cosmic birefringence and numerically computed the power spectrum. Greco et al. (2024) [63] also derived the expression for the B𝐵Bitalic_B-mode power spectrum but did not numerically compute it due to numerical challenges. We further compare our results with the previous numerical computations by P09.

This paper is organized as follows. In section II, we review cosmic birefringence and CMB B𝐵Bitalic_B-mode power spectrum. In section III, we derive the B𝐵Bitalic_B-mode power spectrum induced by anisotropic cosmic birefringence without the thin LSS approximation and evaluate the B𝐵Bitalic_B-mode power spectrum numerically. Section IV is devoted to summary and conclusion.

II Cosmic birefringence

II.1 Cosmic birefringence from massless scalar fields

We first briefly review cosmic birefringence. The presence of the pseudoscalar fields rotates the polarization plane of CMB photons. This rotation angle in general depends on the line-of-sight direction, 𝒏^^𝒏\hat{\bm{n}}over^ start_ARG bold_italic_n end_ARG. The rotation angle for photons emitted at conformal time η𝜂\etaitalic_η in a line-of-sight direction 𝒏^^𝒏\hat{\bm{n}}over^ start_ARG bold_italic_n end_ARG is given by [64, 65, 66] 111We here ignore the gravitational lensing effect which changes the trajectory of the CMB photons. Even if we include the lensing effect, the isotropic cosmic birefringence remains unchanged [67, 68]. On the other hand, the anisotropic cosmic birefringence is modified since the rotation angle depends on δϕ𝛿italic-ϕ\delta\phiitalic_δ italic_ϕ at a slightly different position from χ𝒏^𝜒^𝒏\chi\hat{\bm{n}}italic_χ over^ start_ARG bold_italic_n end_ARG. This remapping effect would be, however, negligible in the anisotropic cosmic birefringence if its angular power spectrum is red; the lensing remapping effect changes any red spectrum only slightly [69]. Since the model for the anisotropic cosmic birefringence in this paper has a steep red spectrum, we ignore the lensing effect on the birefringence angle throughout the paper.

β(η,𝒏^)=g2[ϕ¯(η0)ϕ¯(η)δϕ(η,χ𝒏^)].𝛽𝜂^𝒏𝑔2delimited-[]¯italic-ϕsubscript𝜂0¯italic-ϕ𝜂𝛿italic-ϕ𝜂𝜒^𝒏\displaystyle\beta(\eta,\hat{\bm{n}})=\frac{g}{2}[\bar{\phi}(\eta_{0})-\bar{% \phi}(\eta)-\delta\phi(\eta,\chi\hat{\bm{n}})]\,.italic_β ( italic_η , over^ start_ARG bold_italic_n end_ARG ) = divide start_ARG italic_g end_ARG start_ARG 2 end_ARG [ over¯ start_ARG italic_ϕ end_ARG ( italic_η start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT ) - over¯ start_ARG italic_ϕ end_ARG ( italic_η ) - italic_δ italic_ϕ ( italic_η , italic_χ over^ start_ARG bold_italic_n end_ARG ) ] . (1)

Here, χ=η0η𝜒subscript𝜂0𝜂\chi=\eta_{0}-\etaitalic_χ = italic_η start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT - italic_η is the comoving distance of the source. The rotation angle is decomposed into the isotropic and anisotropic components. Anisotropies of the rotation angle are given by

α(η,𝒏^)gϕ2δϕ(η,χ𝒏^).𝛼𝜂^𝒏subscript𝑔italic-ϕ2𝛿italic-ϕ𝜂𝜒^𝒏\displaystyle\alpha(\eta,\hat{\bm{n}})\equiv-\frac{g_{\phi}}{2}\delta\phi(\eta% ,\chi\hat{\bm{n}})\,.italic_α ( italic_η , over^ start_ARG bold_italic_n end_ARG ) ≡ - divide start_ARG italic_g start_POSTSUBSCRIPT italic_ϕ end_POSTSUBSCRIPT end_ARG start_ARG 2 end_ARG italic_δ italic_ϕ ( italic_η , italic_χ over^ start_ARG bold_italic_n end_ARG ) . (2)

Throughout this paper, we consider the massless pseudoscalar field as the source of the anisotropic cosmic birefringence since the time evolution of the pseudoscalar field is well known and the model can produce a large signal for the anisotropic cosmic birefringence. In this case, the isotropic components become zero and the rotation angle is given by the anisotropic component, β(η,𝒏^)=α(η,𝒏^)𝛽𝜂^𝒏𝛼𝜂^𝒏\beta(\eta,\hat{\bm{n}})=\alpha(\eta,\hat{\bm{n}})italic_β ( italic_η , over^ start_ARG bold_italic_n end_ARG ) = italic_α ( italic_η , over^ start_ARG bold_italic_n end_ARG ) [50]. The evolution equation for the massless pseudoscalar field fluctuations is given by [50]:

d2δϕdη2+2dδϕdη+k2δϕ=0,superscriptd2𝛿italic-ϕdsuperscript𝜂22d𝛿italic-ϕd𝜂superscript𝑘2𝛿italic-ϕ0\displaystyle\dfrac{{\rm d}^{2}\delta\phi}{{\rm d}\eta^{2}}+2\mathcal{H}\dfrac% {{\rm d}\delta\phi}{{\rm d}\eta}+k^{2}\delta\phi=0\,,divide start_ARG roman_d start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT italic_δ italic_ϕ end_ARG start_ARG roman_d italic_η start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT end_ARG + 2 caligraphic_H divide start_ARG roman_d italic_δ italic_ϕ end_ARG start_ARG roman_d italic_η end_ARG + italic_k start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT italic_δ italic_ϕ = 0 , (3)

where =(da/dη)/ad𝑎d𝜂𝑎\mathcal{H}=({\rm d}a/{\rm d}\eta)/acaligraphic_H = ( roman_d italic_a / roman_d italic_η ) / italic_a with a𝑎aitalic_a being the scale factor. We introduce the transfer function to characterize the evolution of the fluctuations:

δϕ(η,𝒌)=T(k,η)δϕini(𝒌).𝛿italic-ϕ𝜂𝒌𝑇𝑘𝜂𝛿subscriptitalic-ϕini𝒌\displaystyle\delta\phi(\eta,\bm{k})=T(k,\eta)\delta\phi_{\rm ini}(\bm{k})\,.italic_δ italic_ϕ ( italic_η , bold_italic_k ) = italic_T ( italic_k , italic_η ) italic_δ italic_ϕ start_POSTSUBSCRIPT roman_ini end_POSTSUBSCRIPT ( bold_italic_k ) . (4)

For the massless pseudoscalar case, during the matter domination epoch, we have

T(k,η)𝑇𝑘𝜂\displaystyle T(k,\eta)italic_T ( italic_k , italic_η ) =3j1(kη)kη.absent3subscript𝑗1𝑘𝜂𝑘𝜂\displaystyle=3\frac{j_{1}(k\eta)}{k\eta}\,.= 3 divide start_ARG italic_j start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT ( italic_k italic_η ) end_ARG start_ARG italic_k italic_η end_ARG . (5)

The transfer function satisfies T(k,η)1𝑇𝑘𝜂1T(k,\eta)\to 1italic_T ( italic_k , italic_η ) → 1 at the super-horizon scale. We define the correlation of the pseudoscalar fluctuations as [70]

δϕini(𝒌)δϕini(𝒌)=2π2k3𝒫ϕ(k)(2π)3δ(3)(𝒌𝒌).delimited-⟨⟩𝛿subscriptsuperscriptitalic-ϕini𝒌𝛿subscriptitalic-ϕinisuperscript𝒌2superscript𝜋2superscript𝑘3subscript𝒫italic-ϕ𝑘superscript2𝜋3superscript𝛿3𝒌superscript𝒌\displaystyle\langle\delta\phi^{*}_{\rm ini}(\bm{k})\delta\phi_{\rm ini}(\bm{k% }^{\prime})\rangle=\frac{2\pi^{2}}{k^{3}}{\cal P}_{\phi}(k)(2\pi)^{3}\delta^{(% 3)}(\bm{k}-\bm{k}^{\prime})\,.⟨ italic_δ italic_ϕ start_POSTSUPERSCRIPT ∗ end_POSTSUPERSCRIPT start_POSTSUBSCRIPT roman_ini end_POSTSUBSCRIPT ( bold_italic_k ) italic_δ italic_ϕ start_POSTSUBSCRIPT roman_ini end_POSTSUBSCRIPT ( bold_italic_k start_POSTSUPERSCRIPT ′ end_POSTSUPERSCRIPT ) ⟩ = divide start_ARG 2 italic_π start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT end_ARG start_ARG italic_k start_POSTSUPERSCRIPT 3 end_POSTSUPERSCRIPT end_ARG caligraphic_P start_POSTSUBSCRIPT italic_ϕ end_POSTSUBSCRIPT ( italic_k ) ( 2 italic_π ) start_POSTSUPERSCRIPT 3 end_POSTSUPERSCRIPT italic_δ start_POSTSUPERSCRIPT ( 3 ) end_POSTSUPERSCRIPT ( bold_italic_k - bold_italic_k start_POSTSUPERSCRIPT ′ end_POSTSUPERSCRIPT ) . (6)

The dimensionless primordial power spectrum for the pseudoscalar fields originated from the initial vacuum fluctuations is given by [50]:

𝒫ϕ(k)subscript𝒫italic-ϕ𝑘\displaystyle{\cal P}_{\phi}(k)caligraphic_P start_POSTSUBSCRIPT italic_ϕ end_POSTSUBSCRIPT ( italic_k ) =(HI2π)2,absentsuperscriptsubscript𝐻𝐼2𝜋2\displaystyle=\left(\frac{H_{I}}{2\pi}\right)^{2}\,,= ( divide start_ARG italic_H start_POSTSUBSCRIPT italic_I end_POSTSUBSCRIPT end_ARG start_ARG 2 italic_π end_ARG ) start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT , (7)

where HIsubscript𝐻𝐼H_{I}italic_H start_POSTSUBSCRIPT italic_I end_POSTSUBSCRIPT is the expansion rate during inflation.

II.2 B-mode power spectrum with the thin LSS approximation

Here, we introduce the B𝐵Bitalic_B-mode power spectrum with the thin LSS approximation. The polarization plane of the CMB photon from the LSS is rotated during propagation, and the Stokes Q𝑄Qitalic_Q and U𝑈Uitalic_U polarization map is distorted as

[Q±iU](𝒏^)=[Q¯±iU¯](𝒏^)e±2iα(η,𝒏^),delimited-[]plus-or-minus𝑄i𝑈^𝒏delimited-[]plus-or-minus¯𝑄i¯𝑈^𝒏superscripteplus-or-minus2i𝛼subscript𝜂^𝒏\displaystyle[Q\pm{\rm i}\hskip 0.50003ptU](\hat{\bm{n}})=[\bar{Q}\pm{\rm i}% \hskip 0.50003pt\bar{U}](\hat{\bm{n}})\hskip 1.00006pt{\rm e}^{\pm 2{\rm i}% \hskip 0.35002pt\alpha(\eta_{*},\hat{\bm{n}})}\,,[ italic_Q ± roman_i italic_U ] ( over^ start_ARG bold_italic_n end_ARG ) = [ over¯ start_ARG italic_Q end_ARG ± roman_i over¯ start_ARG italic_U end_ARG ] ( over^ start_ARG bold_italic_n end_ARG ) roman_e start_POSTSUPERSCRIPT ± 2 roman_i italic_α ( italic_η start_POSTSUBSCRIPT ∗ end_POSTSUBSCRIPT , over^ start_ARG bold_italic_n end_ARG ) end_POSTSUPERSCRIPT , (8)

where the bar denotes the CMB anisotropies without the rotation. The E𝐸Eitalic_E- and B𝐵Bitalic_B-mode polarization are then obtained by [71]

Em±iBm=d2𝒏^Ym±2(𝒏^)[Q±iU](𝒏^),plus-or-minussubscript𝐸𝑚isubscript𝐵𝑚superscriptd2^𝒏subscriptsuperscriptsubscript𝑌𝑚plus-or-minus2^𝒏delimited-[]plus-or-minus𝑄i𝑈^𝒏\displaystyle E_{\ell m}\pm{\rm i}\hskip 0.50003ptB_{\ell m}=-\int\!\!\,{\rm d% }^{2}\hat{\bm{n}}\,\,{}_{\pm 2}Y_{\ell m}^{*}(\hat{\bm{n}})[Q\pm{\rm i}\hskip 0% .50003ptU](\hat{\bm{n}})\,,italic_E start_POSTSUBSCRIPT roman_ℓ italic_m end_POSTSUBSCRIPT ± roman_i italic_B start_POSTSUBSCRIPT roman_ℓ italic_m end_POSTSUBSCRIPT = - ∫ roman_d start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT over^ start_ARG bold_italic_n end_ARG start_FLOATSUBSCRIPT ± 2 end_FLOATSUBSCRIPT italic_Y start_POSTSUBSCRIPT roman_ℓ italic_m end_POSTSUBSCRIPT start_POSTSUPERSCRIPT ∗ end_POSTSUPERSCRIPT ( over^ start_ARG bold_italic_n end_ARG ) [ italic_Q ± roman_i italic_U ] ( over^ start_ARG bold_italic_n end_ARG ) , (9)

where Ym±2(𝒏^)subscriptsubscript𝑌𝑚plus-or-minus2^𝒏{}_{\pm 2}Y_{\ell m}(\hat{\bm{n}})start_FLOATSUBSCRIPT ± 2 end_FLOATSUBSCRIPT italic_Y start_POSTSUBSCRIPT roman_ℓ italic_m end_POSTSUBSCRIPT ( over^ start_ARG bold_italic_n end_ARG ) is the spin-2222 spherical harmonics. Similarly with the spin-00 (scalar) spherical harmonics, Ym(𝒏^)subscript𝑌𝑚^𝒏Y_{\ell m}(\hat{\bm{n}})italic_Y start_POSTSUBSCRIPT roman_ℓ italic_m end_POSTSUBSCRIPT ( over^ start_ARG bold_italic_n end_ARG ), the harmonic coefficients of the rotation angle is defined as

αLM(η)subscript𝛼𝐿𝑀subscript𝜂\displaystyle\alpha_{LM}(\eta_{*})italic_α start_POSTSUBSCRIPT italic_L italic_M end_POSTSUBSCRIPT ( italic_η start_POSTSUBSCRIPT ∗ end_POSTSUBSCRIPT ) =d2𝒏^YLM(𝒏^)α(η,𝒏^).absentsuperscriptd2^𝒏superscriptsubscript𝑌𝐿𝑀^𝒏𝛼subscript𝜂^𝒏\displaystyle=\int\!\!\,{\rm d}^{2}\hat{\bm{n}}\,\,Y_{LM}^{*}(\hat{\bm{n}})% \alpha(\eta_{*},\hat{\bm{n}})\,.= ∫ roman_d start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT over^ start_ARG bold_italic_n end_ARG italic_Y start_POSTSUBSCRIPT italic_L italic_M end_POSTSUBSCRIPT start_POSTSUPERSCRIPT ∗ end_POSTSUPERSCRIPT ( over^ start_ARG bold_italic_n end_ARG ) italic_α ( italic_η start_POSTSUBSCRIPT ∗ end_POSTSUBSCRIPT , over^ start_ARG bold_italic_n end_ARG ) . (10)

Expanding the right-hand side of Eq. 8 to first order in the rotation angle and transforming the Stokes Q/U𝑄𝑈Q/Uitalic_Q / italic_U parameters to E/B𝐸𝐵E/Bitalic_E / italic_B modes using Eq. 9, we obtain the distorted B𝐵Bitalic_B modes by the anisotropic cosmic birefringence as

Bmsubscript𝐵𝑚\displaystyle B_{\ell m}italic_B start_POSTSUBSCRIPT roman_ℓ italic_m end_POSTSUBSCRIPT =B¯m+LMm(1)m(LmMm)absentsubscript¯𝐵𝑚subscript𝐿𝑀superscriptsuperscript𝑚superscript1𝑚matrix𝐿superscript𝑚𝑀superscript𝑚\displaystyle=\bar{B}_{\ell m}+\sum_{LM\ell^{\prime}m^{\prime}}(-1)^{m}\begin{% pmatrix}\ell&L&\ell^{\prime}\\ -m&M&m^{\prime}\end{pmatrix}= over¯ start_ARG italic_B end_ARG start_POSTSUBSCRIPT roman_ℓ italic_m end_POSTSUBSCRIPT + ∑ start_POSTSUBSCRIPT italic_L italic_M roman_ℓ start_POSTSUPERSCRIPT ′ end_POSTSUPERSCRIPT italic_m start_POSTSUPERSCRIPT ′ end_POSTSUPERSCRIPT end_POSTSUBSCRIPT ( - 1 ) start_POSTSUPERSCRIPT italic_m end_POSTSUPERSCRIPT ( start_ARG start_ROW start_CELL roman_ℓ end_CELL start_CELL italic_L end_CELL start_CELL roman_ℓ start_POSTSUPERSCRIPT ′ end_POSTSUPERSCRIPT end_CELL end_ROW start_ROW start_CELL - italic_m end_CELL start_CELL italic_M end_CELL start_CELL italic_m start_POSTSUPERSCRIPT ′ end_POSTSUPERSCRIPT end_CELL end_ROW end_ARG )
×αLM(η)(WLE¯m+WL+B¯m).absentsubscript𝛼𝐿𝑀subscript𝜂subscriptsuperscript𝑊𝐿superscriptsubscript¯𝐸superscriptsuperscript𝑚subscriptsuperscript𝑊𝐿superscriptsubscript¯𝐵superscriptsuperscript𝑚\displaystyle\times\alpha_{LM}(\eta_{*})\left(-W^{-}_{\ell L\ell^{\prime}}\bar% {E}_{\ell^{\prime}m^{\prime}}+W^{+}_{\ell L\ell^{\prime}}\bar{B}_{\ell^{\prime% }m^{\prime}}\right)\,.× italic_α start_POSTSUBSCRIPT italic_L italic_M end_POSTSUBSCRIPT ( italic_η start_POSTSUBSCRIPT ∗ end_POSTSUBSCRIPT ) ( - italic_W start_POSTSUPERSCRIPT - end_POSTSUPERSCRIPT start_POSTSUBSCRIPT roman_ℓ italic_L roman_ℓ start_POSTSUPERSCRIPT ′ end_POSTSUPERSCRIPT end_POSTSUBSCRIPT over¯ start_ARG italic_E end_ARG start_POSTSUBSCRIPT roman_ℓ start_POSTSUPERSCRIPT ′ end_POSTSUPERSCRIPT italic_m start_POSTSUPERSCRIPT ′ end_POSTSUPERSCRIPT end_POSTSUBSCRIPT + italic_W start_POSTSUPERSCRIPT + end_POSTSUPERSCRIPT start_POSTSUBSCRIPT roman_ℓ italic_L roman_ℓ start_POSTSUPERSCRIPT ′ end_POSTSUPERSCRIPT end_POSTSUBSCRIPT over¯ start_ARG italic_B end_ARG start_POSTSUBSCRIPT roman_ℓ start_POSTSUPERSCRIPT ′ end_POSTSUPERSCRIPT italic_m start_POSTSUPERSCRIPT ′ end_POSTSUPERSCRIPT end_POSTSUBSCRIPT ) . (11)

Here, we define

WL±subscriptsuperscript𝑊plus-or-minus𝐿superscript\displaystyle W^{\pm}_{\ell L\ell^{\prime}}italic_W start_POSTSUPERSCRIPT ± end_POSTSUPERSCRIPT start_POSTSUBSCRIPT roman_ℓ italic_L roman_ℓ start_POSTSUPERSCRIPT ′ end_POSTSUPERSCRIPT end_POSTSUBSCRIPT =±2iζ±pL(2+1)(2L+1)(2+1)4π(L202),absentplus-or-minus2isuperscript𝜁plus-or-minussubscriptsuperscript𝑝minus-or-plus𝐿superscript212𝐿12superscript14𝜋matrix𝐿superscript202\displaystyle=\pm 2{\rm i}\hskip 0.50003pt\zeta^{\pm}p^{\mp}_{\ell L\ell^{% \prime}}\sqrt{\frac{(2\ell+1)(2L+1)(2\ell^{\prime}+1)}{4\pi}}\begin{pmatrix}% \ell&L&\ell^{\prime}\\ 2&0&-2\end{pmatrix}\,,= ± 2 roman_i italic_ζ start_POSTSUPERSCRIPT ± end_POSTSUPERSCRIPT italic_p start_POSTSUPERSCRIPT ∓ end_POSTSUPERSCRIPT start_POSTSUBSCRIPT roman_ℓ italic_L roman_ℓ start_POSTSUPERSCRIPT ′ end_POSTSUPERSCRIPT end_POSTSUBSCRIPT square-root start_ARG divide start_ARG ( 2 roman_ℓ + 1 ) ( 2 italic_L + 1 ) ( 2 roman_ℓ start_POSTSUPERSCRIPT ′ end_POSTSUPERSCRIPT + 1 ) end_ARG start_ARG 4 italic_π end_ARG end_ARG ( start_ARG start_ROW start_CELL roman_ℓ end_CELL start_CELL italic_L end_CELL start_CELL roman_ℓ start_POSTSUPERSCRIPT ′ end_POSTSUPERSCRIPT end_CELL end_ROW start_ROW start_CELL 2 end_CELL start_CELL 0 end_CELL start_CELL - 2 end_CELL end_ROW end_ARG ) , (12)

with ζ+=1superscript𝜁1\zeta^{+}=1italic_ζ start_POSTSUPERSCRIPT + end_POSTSUPERSCRIPT = 1, ζ=isuperscript𝜁i\zeta^{-}={\rm i}\hskip 0.50003ptitalic_ζ start_POSTSUPERSCRIPT - end_POSTSUPERSCRIPT = roman_i, and pL=[1(1)+L+]/2subscriptsuperscript𝑝minus-or-plus𝐿superscriptdelimited-[]minus-or-plus1superscript1𝐿superscript2p^{\mp}_{\ell L\ell^{\prime}}=[1\mp(-1)^{\ell+L+\ell^{\prime}}]/2italic_p start_POSTSUPERSCRIPT ∓ end_POSTSUPERSCRIPT start_POSTSUBSCRIPT roman_ℓ italic_L roman_ℓ start_POSTSUPERSCRIPT ′ end_POSTSUPERSCRIPT end_POSTSUBSCRIPT = [ 1 ∓ ( - 1 ) start_POSTSUPERSCRIPT roman_ℓ + italic_L + roman_ℓ start_POSTSUPERSCRIPT ′ end_POSTSUPERSCRIPT end_POSTSUPERSCRIPT ] / 2 is the result of the parity symmetry. The B𝐵Bitalic_B-mode power spectrum with the thin LSS approximation is given by (e.g. [55])

CBBsuperscriptsubscript𝐶𝐵𝐵\displaystyle C_{\ell}^{BB}italic_C start_POSTSUBSCRIPT roman_ℓ end_POSTSUBSCRIPT start_POSTSUPERSCRIPT italic_B italic_B end_POSTSUPERSCRIPT =|Bm|2=L12+1CLαα(η,η)absentdelimited-⟨⟩superscriptsubscript𝐵𝑚2subscriptsuperscript𝐿121superscriptsubscript𝐶𝐿𝛼𝛼subscript𝜂subscript𝜂\displaystyle=\langle|B_{\ell m}|^{2}\rangle=\sum_{\ell^{\prime}L}\frac{1}{2% \ell+1}C_{L}^{\alpha\alpha}(\eta_{*},\eta_{*})= ⟨ | italic_B start_POSTSUBSCRIPT roman_ℓ italic_m end_POSTSUBSCRIPT | start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT ⟩ = ∑ start_POSTSUBSCRIPT roman_ℓ start_POSTSUPERSCRIPT ′ end_POSTSUPERSCRIPT italic_L end_POSTSUBSCRIPT divide start_ARG 1 end_ARG start_ARG 2 roman_ℓ + 1 end_ARG italic_C start_POSTSUBSCRIPT italic_L end_POSTSUBSCRIPT start_POSTSUPERSCRIPT italic_α italic_α end_POSTSUPERSCRIPT ( italic_η start_POSTSUBSCRIPT ∗ end_POSTSUBSCRIPT , italic_η start_POSTSUBSCRIPT ∗ end_POSTSUBSCRIPT )
×[|WL|2CE¯E¯+|WL+|2CB¯B¯]absentdelimited-[]superscriptsuperscriptsubscript𝑊𝐿superscript2superscriptsubscript𝐶superscript¯𝐸¯𝐸superscriptsuperscriptsubscript𝑊𝐿superscript2superscriptsubscript𝐶superscript¯𝐵¯𝐵\displaystyle\times[|W_{\ell L\ell^{\prime}}^{-}|^{2}C_{\ell^{\prime}}^{\bar{E% }\bar{E}}+|W_{\ell L\ell^{\prime}}^{+}|^{2}C_{\ell^{\prime}}^{\bar{B}\bar{B}}]× [ | italic_W start_POSTSUBSCRIPT roman_ℓ italic_L roman_ℓ start_POSTSUPERSCRIPT ′ end_POSTSUPERSCRIPT end_POSTSUBSCRIPT start_POSTSUPERSCRIPT - end_POSTSUPERSCRIPT | start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT italic_C start_POSTSUBSCRIPT roman_ℓ start_POSTSUPERSCRIPT ′ end_POSTSUPERSCRIPT end_POSTSUBSCRIPT start_POSTSUPERSCRIPT over¯ start_ARG italic_E end_ARG over¯ start_ARG italic_E end_ARG end_POSTSUPERSCRIPT + | italic_W start_POSTSUBSCRIPT roman_ℓ italic_L roman_ℓ start_POSTSUPERSCRIPT ′ end_POSTSUPERSCRIPT end_POSTSUBSCRIPT start_POSTSUPERSCRIPT + end_POSTSUPERSCRIPT | start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT italic_C start_POSTSUBSCRIPT roman_ℓ start_POSTSUPERSCRIPT ′ end_POSTSUPERSCRIPT end_POSTSUBSCRIPT start_POSTSUPERSCRIPT over¯ start_ARG italic_B end_ARG over¯ start_ARG italic_B end_ARG end_POSTSUPERSCRIPT ]
=4L(2+1)(2L+1)4π(L202)2absent4subscriptsuperscript𝐿2superscript12𝐿14𝜋superscriptmatrix𝐿superscript2022\displaystyle=4\sum_{\ell^{\prime}L}\frac{(2\ell^{\prime}+1)(2L+1)}{4\pi}% \begin{pmatrix}\ell&L&\ell^{\prime}\\ 2&0&-2\end{pmatrix}^{2}= 4 ∑ start_POSTSUBSCRIPT roman_ℓ start_POSTSUPERSCRIPT ′ end_POSTSUPERSCRIPT italic_L end_POSTSUBSCRIPT divide start_ARG ( 2 roman_ℓ start_POSTSUPERSCRIPT ′ end_POSTSUPERSCRIPT + 1 ) ( 2 italic_L + 1 ) end_ARG start_ARG 4 italic_π end_ARG ( start_ARG start_ROW start_CELL roman_ℓ end_CELL start_CELL italic_L end_CELL start_CELL roman_ℓ start_POSTSUPERSCRIPT ′ end_POSTSUPERSCRIPT end_CELL end_ROW start_ROW start_CELL 2 end_CELL start_CELL 0 end_CELL start_CELL - 2 end_CELL end_ROW end_ARG ) start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT
×CLαα(η,η)[pL+CE¯E¯+pLCB¯B¯].absentsuperscriptsubscript𝐶𝐿𝛼𝛼subscript𝜂subscript𝜂delimited-[]subscriptsuperscript𝑝𝐿superscriptsuperscriptsubscript𝐶superscript¯𝐸¯𝐸subscriptsuperscript𝑝𝐿superscriptsuperscriptsubscript𝐶superscript¯𝐵¯𝐵\displaystyle\times C_{L}^{\alpha\alpha}(\eta_{*},\eta_{*})[p^{+}_{\ell L\ell^% {\prime}}C_{\ell^{\prime}}^{\bar{E}\bar{E}}+p^{-}_{\ell L\ell^{\prime}}C_{\ell% ^{\prime}}^{\bar{B}\bar{B}}]\,.× italic_C start_POSTSUBSCRIPT italic_L end_POSTSUBSCRIPT start_POSTSUPERSCRIPT italic_α italic_α end_POSTSUPERSCRIPT ( italic_η start_POSTSUBSCRIPT ∗ end_POSTSUBSCRIPT , italic_η start_POSTSUBSCRIPT ∗ end_POSTSUBSCRIPT ) [ italic_p start_POSTSUPERSCRIPT + end_POSTSUPERSCRIPT start_POSTSUBSCRIPT roman_ℓ italic_L roman_ℓ start_POSTSUPERSCRIPT ′ end_POSTSUPERSCRIPT end_POSTSUBSCRIPT italic_C start_POSTSUBSCRIPT roman_ℓ start_POSTSUPERSCRIPT ′ end_POSTSUPERSCRIPT end_POSTSUBSCRIPT start_POSTSUPERSCRIPT over¯ start_ARG italic_E end_ARG over¯ start_ARG italic_E end_ARG end_POSTSUPERSCRIPT + italic_p start_POSTSUPERSCRIPT - end_POSTSUPERSCRIPT start_POSTSUBSCRIPT roman_ℓ italic_L roman_ℓ start_POSTSUPERSCRIPT ′ end_POSTSUPERSCRIPT end_POSTSUBSCRIPT italic_C start_POSTSUBSCRIPT roman_ℓ start_POSTSUPERSCRIPT ′ end_POSTSUPERSCRIPT end_POSTSUBSCRIPT start_POSTSUPERSCRIPT over¯ start_ARG italic_B end_ARG over¯ start_ARG italic_B end_ARG end_POSTSUPERSCRIPT ] . (13)

Here, we define the angular power spectrum of the anisotropic rotation fields as [70]

CLαα(η,η)superscriptsubscript𝐶𝐿𝛼𝛼𝜂superscript𝜂\displaystyle C_{L}^{\alpha\alpha}(\eta,\eta^{\prime})italic_C start_POSTSUBSCRIPT italic_L end_POSTSUBSCRIPT start_POSTSUPERSCRIPT italic_α italic_α end_POSTSUPERSCRIPT ( italic_η , italic_η start_POSTSUPERSCRIPT ′ end_POSTSUPERSCRIPT ) =αLM(η)αLM(η)absentdelimited-⟨⟩subscript𝛼𝐿𝑀𝜂subscriptsuperscript𝛼𝐿𝑀superscript𝜂\displaystyle=\langle\alpha_{LM}(\eta)\alpha^{*}_{LM}(\eta^{\prime})\rangle= ⟨ italic_α start_POSTSUBSCRIPT italic_L italic_M end_POSTSUBSCRIPT ( italic_η ) italic_α start_POSTSUPERSCRIPT ∗ end_POSTSUPERSCRIPT start_POSTSUBSCRIPT italic_L italic_M end_POSTSUBSCRIPT ( italic_η start_POSTSUPERSCRIPT ′ end_POSTSUPERSCRIPT ) ⟩
=4πdlnk𝒫ϕ(k)uL(k,η)uL(k,η),absent4𝜋d𝑘subscript𝒫italic-ϕ𝑘subscript𝑢𝐿𝑘𝜂subscript𝑢𝐿𝑘superscript𝜂\displaystyle=4\pi\int\!\!\,{\rm d}\ln k\,\,{\cal P}_{\phi}(k)u_{L}(k,\eta)u_{% L}(k,\eta^{\prime})\,,= 4 italic_π ∫ roman_d roman_ln italic_k caligraphic_P start_POSTSUBSCRIPT italic_ϕ end_POSTSUBSCRIPT ( italic_k ) italic_u start_POSTSUBSCRIPT italic_L end_POSTSUBSCRIPT ( italic_k , italic_η ) italic_u start_POSTSUBSCRIPT italic_L end_POSTSUBSCRIPT ( italic_k , italic_η start_POSTSUPERSCRIPT ′ end_POSTSUPERSCRIPT ) , (14)

where we define

uL(k,η)=gϕ2jL(k(η0η))T(k,η).subscript𝑢𝐿𝑘𝜂subscript𝑔italic-ϕ2subscript𝑗𝐿𝑘subscript𝜂0𝜂𝑇𝑘𝜂\displaystyle u_{L}(k,\eta)=\frac{g_{\phi}}{2}j_{L}(k(\eta_{0}-\eta))T(k,\eta)\,.italic_u start_POSTSUBSCRIPT italic_L end_POSTSUBSCRIPT ( italic_k , italic_η ) = divide start_ARG italic_g start_POSTSUBSCRIPT italic_ϕ end_POSTSUBSCRIPT end_ARG start_ARG 2 end_ARG italic_j start_POSTSUBSCRIPT italic_L end_POSTSUBSCRIPT ( italic_k ( italic_η start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT - italic_η ) ) italic_T ( italic_k , italic_η ) . (15)

At large-angular scales, L100less-than-or-similar-to𝐿100L\lesssim 100italic_L ≲ 100, the transfer function is approximately unity, and the angular power spectrum of the anisotropic birefringence angle is given by [50]:

CLαα2πL(L+1)(gϕ2)2(HI2π)22πL(L+1)ACB,similar-to-or-equalssuperscriptsubscript𝐶𝐿𝛼𝛼2𝜋𝐿𝐿1superscriptsubscript𝑔italic-ϕ22superscriptsubscript𝐻𝐼2𝜋22𝜋𝐿𝐿1subscript𝐴CB\displaystyle C_{L}^{\alpha\alpha}\simeq\frac{2\pi}{L(L+1)}\left(\frac{g_{\phi% }}{2}\right)^{2}\left(\frac{H_{I}}{2\pi}\right)^{2}\equiv\frac{2\pi}{L(L+1)}A_% {\rm CB}\,,italic_C start_POSTSUBSCRIPT italic_L end_POSTSUBSCRIPT start_POSTSUPERSCRIPT italic_α italic_α end_POSTSUPERSCRIPT ≃ divide start_ARG 2 italic_π end_ARG start_ARG italic_L ( italic_L + 1 ) end_ARG ( divide start_ARG italic_g start_POSTSUBSCRIPT italic_ϕ end_POSTSUBSCRIPT end_ARG start_ARG 2 end_ARG ) start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT ( divide start_ARG italic_H start_POSTSUBSCRIPT italic_I end_POSTSUBSCRIPT end_ARG start_ARG 2 italic_π end_ARG ) start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT ≡ divide start_ARG 2 italic_π end_ARG start_ARG italic_L ( italic_L + 1 ) end_ARG italic_A start_POSTSUBSCRIPT roman_CB end_POSTSUBSCRIPT , (16)

where we introduce a constant, ACBsubscript𝐴CBA_{\rm CB}italic_A start_POSTSUBSCRIPT roman_CB end_POSTSUBSCRIPT, and use

0dlnxjL2(x)=12L(L+1).superscriptsubscript0d𝑥subscriptsuperscript𝑗2𝐿𝑥12𝐿𝐿1\displaystyle\int_{0}^{\infty}\!\!\,{\rm d}\ln x\,\,j^{2}_{L}(x)=\frac{1}{2L(L% +1)}\,.∫ start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT start_POSTSUPERSCRIPT ∞ end_POSTSUPERSCRIPT roman_d roman_ln italic_x italic_j start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT start_POSTSUBSCRIPT italic_L end_POSTSUBSCRIPT ( italic_x ) = divide start_ARG 1 end_ARG start_ARG 2 italic_L ( italic_L + 1 ) end_ARG . (17)

II.3 Other contributions

Refer to caption
Figure 1: Ratio of the B𝐵Bitalic_B-mode power spectra with and without the higher-order terms, 𝒪[(CLαα)2]𝒪delimited-[]superscriptsubscriptsuperscript𝐶𝛼𝛼𝐿2\mathcal{O}[(C^{\alpha\alpha}_{L})^{2}]caligraphic_O [ ( italic_C start_POSTSUPERSCRIPT italic_α italic_α end_POSTSUPERSCRIPT start_POSTSUBSCRIPT italic_L end_POSTSUBSCRIPT ) start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT ] (“higher-order”) and with and without the lensing effect (“lensing”).

In the above calculation, we made some simplifications. For example, we ignored the contributions from the higher-order terms of 𝒪[(CLαα)2]𝒪delimited-[]superscriptsuperscriptsubscript𝐶𝐿𝛼𝛼2\mathcal{O}[(C_{L}^{\alpha\alpha})^{2}]caligraphic_O [ ( italic_C start_POSTSUBSCRIPT italic_L end_POSTSUBSCRIPT start_POSTSUPERSCRIPT italic_α italic_α end_POSTSUPERSCRIPT ) start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT ] to the B𝐵Bitalic_B-mode power spectrum. We also ignore the lensing contributions.

To see the impact of the higher-order terms and lensing, we compare the cases with and without the higher-order terms or the lensing effect using the public software package class_rot developed by Cai & Guan (2021) [60] with the thin LSS approximation. We find that these contributions are less than a few percent at 1000less-than-or-similar-to1000\ell\lesssim 1000roman_ℓ ≲ 1000 and we can approximate B𝐵Bitalic_B-mode power spectrum with the perturbed spectrum without lensing. On the other hand, at very small scales, 3000greater-than-or-equivalent-to3000\ell\gtrsim 3000roman_ℓ ≳ 3000, both contributions become important.

III B-mode power spectrum from anisotropic cosmic birefringence

Here, we derive the expression for the B𝐵Bitalic_B-mode power spectra from time-varying anisotropic cosmic birefringence without the thin LSS approximation.

III.1 B-mode power spectrum without the thin LSS approximation

We follow Ref. [72] where they derive the B𝐵Bitalic_B-mode power spectrum from the Faraday rotation due to the primordial magnetic fields without the thin LSS approximation. The derivation is based on the total angular momentum method of Ref. [73].

The evolution equation for the 𝒒𝒒\bm{q}bold_italic_q Fourier mode of P±=Q±iUsubscript𝑃plus-or-minusplus-or-minus𝑄i𝑈P_{\pm}=Q\pm{\rm i}\hskip 0.50003ptUitalic_P start_POSTSUBSCRIPT ± end_POSTSUBSCRIPT = italic_Q ± roman_i italic_U is given by [72, 36]

dP±dη(η,𝒒,𝒏^)+iqμP±(η,𝒒,𝒏^)dsubscript𝑃plus-or-minusd𝜂𝜂𝒒^𝒏i𝑞𝜇subscript𝑃plus-or-minus𝜂𝒒^𝒏\displaystyle\dfrac{{\rm d}P_{\pm}}{{\rm d}\eta}(\eta,\bm{q},\hat{\bm{n}})+{% \rm i}\hskip 0.50003ptq\mu P_{\pm}(\eta,\bm{q},\hat{\bm{n}})divide start_ARG roman_d italic_P start_POSTSUBSCRIPT ± end_POSTSUBSCRIPT end_ARG start_ARG roman_d italic_η end_ARG ( italic_η , bold_italic_q , over^ start_ARG bold_italic_n end_ARG ) + roman_i italic_q italic_μ italic_P start_POSTSUBSCRIPT ± end_POSTSUBSCRIPT ( italic_η , bold_italic_q , over^ start_ARG bold_italic_n end_ARG )
=τ˙[P±(η,𝒒,𝒏^)+6P(0)(η,𝒒)4π5Y20±2(𝒏^)]absent˙𝜏delimited-[]subscript𝑃plus-or-minus𝜂𝒒^𝒏6superscript𝑃0𝜂𝒒4𝜋5subscriptsubscript𝑌20plus-or-minus2^𝒏\displaystyle=\dot{\tau}\left[-P_{\pm}(\eta,\bm{q},\hat{\bm{n}})+\sqrt{6}P^{(0% )}(\eta,\bm{q})\sqrt{\frac{4\pi}{5}}{}_{\pm 2}Y_{20}(\hat{\bm{n}})\right]= over˙ start_ARG italic_τ end_ARG [ - italic_P start_POSTSUBSCRIPT ± end_POSTSUBSCRIPT ( italic_η , bold_italic_q , over^ start_ARG bold_italic_n end_ARG ) + square-root start_ARG 6 end_ARG italic_P start_POSTSUPERSCRIPT ( 0 ) end_POSTSUPERSCRIPT ( italic_η , bold_italic_q ) square-root start_ARG divide start_ARG 4 italic_π end_ARG start_ARG 5 end_ARG end_ARG start_FLOATSUBSCRIPT ± 2 end_FLOATSUBSCRIPT italic_Y start_POSTSUBSCRIPT 20 end_POSTSUBSCRIPT ( over^ start_ARG bold_italic_n end_ARG ) ]
idαdη(η,𝒏^)P±(η,𝒒,𝒏^).minus-or-plusid𝛼d𝜂𝜂^𝒏subscript𝑃plus-or-minus𝜂𝒒^𝒏\displaystyle\mp{\rm i}\hskip 0.50003pt\dfrac{{\rm d}\alpha}{{\rm d}\eta}(\eta% ,\hat{\bm{n}})P_{\pm}(\eta,\bm{q},\hat{\bm{n}})\,.∓ roman_i divide start_ARG roman_d italic_α end_ARG start_ARG roman_d italic_η end_ARG ( italic_η , over^ start_ARG bold_italic_n end_ARG ) italic_P start_POSTSUBSCRIPT ± end_POSTSUBSCRIPT ( italic_η , bold_italic_q , over^ start_ARG bold_italic_n end_ARG ) . (18)

Here, μ=𝒏^𝒒^𝜇^𝒏^𝒒\mu=\hat{\bm{n}}\cdot\hat{\bm{q}}italic_μ = over^ start_ARG bold_italic_n end_ARG ⋅ over^ start_ARG bold_italic_q end_ARG, and we define the differential CMB optical depth, τ˙aneσT˙𝜏𝑎subscript𝑛esubscript𝜎T\dot{\tau}\equiv an_{\rm e}\sigma_{\rm T}over˙ start_ARG italic_τ end_ARG ≡ italic_a italic_n start_POSTSUBSCRIPT roman_e end_POSTSUBSCRIPT italic_σ start_POSTSUBSCRIPT roman_T end_POSTSUBSCRIPT, with nesubscript𝑛en_{\rm e}italic_n start_POSTSUBSCRIPT roman_e end_POSTSUBSCRIPT and σTsubscript𝜎T\sigma_{\rm T}italic_σ start_POSTSUBSCRIPT roman_T end_POSTSUBSCRIPT being the electron number density and the Thomson-scattering cross-section, respectively. The quantity, P(0)=(Θ26E2)/10superscript𝑃0subscriptΘ26subscript𝐸210P^{(0)}=(\Theta_{2}-\sqrt{6}E_{2})/10italic_P start_POSTSUPERSCRIPT ( 0 ) end_POSTSUPERSCRIPT = ( roman_Θ start_POSTSUBSCRIPT 2 end_POSTSUBSCRIPT - square-root start_ARG 6 end_ARG italic_E start_POSTSUBSCRIPT 2 end_POSTSUBSCRIPT ) / 10, is the polarization source with Θ2subscriptΘ2\Theta_{2}roman_Θ start_POSTSUBSCRIPT 2 end_POSTSUBSCRIPT and E2subscript𝐸2E_{2}italic_E start_POSTSUBSCRIPT 2 end_POSTSUBSCRIPT being the quadrupole of temperature and E𝐸Eitalic_E-mode polarization from the scalar perturbations [73]. The above equation has the following solution:

P±(𝒒,𝒏^)subscript𝑃plus-or-minus𝒒^𝒏\displaystyle P_{\pm}(\bm{q},\hat{\bm{n}})italic_P start_POSTSUBSCRIPT ± end_POSTSUBSCRIPT ( bold_italic_q , over^ start_ARG bold_italic_n end_ARG ) =0η0dηgv(η)6P(0)(η,𝒒)[1±2iα(η,𝒏^)]absentsuperscriptsubscript0subscript𝜂0differential-d𝜂subscript𝑔v𝜂6superscript𝑃0𝜂𝒒delimited-[]plus-or-minus12i𝛼𝜂^𝒏\displaystyle=-\int_{0}^{\eta_{0}}\!\!\,{\rm d}\eta\,\,g_{\rm v}(\eta)\sqrt{6}% P^{(0)}(\eta,\bm{q})\left[1\pm 2{\rm i}\hskip 0.50003pt\alpha(\eta,\hat{\bm{n}% })\right]= - ∫ start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT start_POSTSUPERSCRIPT italic_η start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT end_POSTSUPERSCRIPT roman_d italic_η italic_g start_POSTSUBSCRIPT roman_v end_POSTSUBSCRIPT ( italic_η ) square-root start_ARG 6 end_ARG italic_P start_POSTSUPERSCRIPT ( 0 ) end_POSTSUPERSCRIPT ( italic_η , bold_italic_q ) [ 1 ± 2 roman_i italic_α ( italic_η , over^ start_ARG bold_italic_n end_ARG ) ]
×(i)4π(2+1)ϵ(0)(q(η0η))Y0±2(𝒏^)\displaystyle\times\sum_{\ell}(-{\rm i}\hskip 0.50003pt)^{\ell}\sqrt{4\pi(2% \ell+1)}\epsilon^{(0)}_{\ell}(q(\eta_{0}-\eta)){}_{\pm 2}Y_{\ell 0}(\hat{\bm{n% }})× ∑ start_POSTSUBSCRIPT roman_ℓ end_POSTSUBSCRIPT ( - roman_i ) start_POSTSUPERSCRIPT roman_ℓ end_POSTSUPERSCRIPT square-root start_ARG 4 italic_π ( 2 roman_ℓ + 1 ) end_ARG italic_ϵ start_POSTSUPERSCRIPT ( 0 ) end_POSTSUPERSCRIPT start_POSTSUBSCRIPT roman_ℓ end_POSTSUBSCRIPT ( italic_q ( italic_η start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT - italic_η ) ) start_FLOATSUBSCRIPT ± 2 end_FLOATSUBSCRIPT italic_Y start_POSTSUBSCRIPT roman_ℓ 0 end_POSTSUBSCRIPT ( over^ start_ARG bold_italic_n end_ARG )
+𝒪(α2),𝒪superscript𝛼2\displaystyle+\mathcal{O}(\alpha^{2})\,,+ caligraphic_O ( italic_α start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT ) , (19)

where gv(η)subscript𝑔v𝜂g_{\rm v}(\eta)italic_g start_POSTSUBSCRIPT roman_v end_POSTSUBSCRIPT ( italic_η ) is the visibility function and

ϵ(0)(x)=38(+2)!(2)!j(x)x2.superscriptsubscriptitalic-ϵ0𝑥3822subscript𝑗𝑥superscript𝑥2\displaystyle\epsilon_{\ell}^{(0)}(x)=\sqrt{\frac{3}{8}\frac{(\ell+2)!}{(\ell-% 2)!}}\frac{j_{\ell}(x)}{x^{2}}\,.italic_ϵ start_POSTSUBSCRIPT roman_ℓ end_POSTSUBSCRIPT start_POSTSUPERSCRIPT ( 0 ) end_POSTSUPERSCRIPT ( italic_x ) = square-root start_ARG divide start_ARG 3 end_ARG start_ARG 8 end_ARG divide start_ARG ( roman_ℓ + 2 ) ! end_ARG start_ARG ( roman_ℓ - 2 ) ! end_ARG end_ARG divide start_ARG italic_j start_POSTSUBSCRIPT roman_ℓ end_POSTSUBSCRIPT ( italic_x ) end_ARG start_ARG italic_x start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT end_ARG . (20)

Next, we use the total angular momentum method to derive the expression for the B𝐵Bitalic_B-mode power spectrum. We decompose the polarization fields as

P±(𝒒,𝒏^)subscript𝑃plus-or-minus𝒒^𝒏\displaystyle P_{\pm}(\bm{q},\hat{\bm{n}})italic_P start_POSTSUBSCRIPT ± end_POSTSUBSCRIPT ( bold_italic_q , over^ start_ARG bold_italic_n end_ARG ) =(i)4π2+1absentsubscriptsuperscripti4𝜋21\displaystyle=\sum_{\ell}(-{\rm i}\hskip 0.50003pt)^{\ell}\sqrt{\frac{4\pi}{2% \ell+1}}= ∑ start_POSTSUBSCRIPT roman_ℓ end_POSTSUBSCRIPT ( - roman_i ) start_POSTSUPERSCRIPT roman_ℓ end_POSTSUPERSCRIPT square-root start_ARG divide start_ARG 4 italic_π end_ARG start_ARG 2 roman_ℓ + 1 end_ARG end_ARG
×m=(E(m)(𝒒)±iB(m)(𝒒))Ym±2(𝒏^),\displaystyle\times\sum_{m=-\ell}^{\ell}(E_{\ell}^{(m)}(\bm{q})\pm{\rm i}% \hskip 0.50003ptB_{\ell}^{(m)}(\bm{q})){}_{\pm 2}Y_{\ell m}(\hat{\bm{n}})\,,× ∑ start_POSTSUBSCRIPT italic_m = - roman_ℓ end_POSTSUBSCRIPT start_POSTSUPERSCRIPT roman_ℓ end_POSTSUPERSCRIPT ( italic_E start_POSTSUBSCRIPT roman_ℓ end_POSTSUBSCRIPT start_POSTSUPERSCRIPT ( italic_m ) end_POSTSUPERSCRIPT ( bold_italic_q ) ± roman_i italic_B start_POSTSUBSCRIPT roman_ℓ end_POSTSUBSCRIPT start_POSTSUPERSCRIPT ( italic_m ) end_POSTSUPERSCRIPT ( bold_italic_q ) ) start_FLOATSUBSCRIPT ± 2 end_FLOATSUBSCRIPT italic_Y start_POSTSUBSCRIPT roman_ℓ italic_m end_POSTSUBSCRIPT ( over^ start_ARG bold_italic_n end_ARG ) , (21)

with 𝒒^=𝒛^^𝒒^𝒛\hat{\bm{q}}=\hat{\bm{z}}over^ start_ARG bold_italic_q end_ARG = over^ start_ARG bold_italic_z end_ARG. Note that B(m)(𝒒)superscriptsubscript𝐵𝑚𝒒B_{\ell}^{(m)}(\bm{q})italic_B start_POSTSUBSCRIPT roman_ℓ end_POSTSUBSCRIPT start_POSTSUPERSCRIPT ( italic_m ) end_POSTSUPERSCRIPT ( bold_italic_q ) vanishes in the absence of cosmic birefringence, and the scalar perturbations alone lead to m=0𝑚0m=0italic_m = 0. We then invert the above equation with the spin-2 spherical harmonics and obtain:

B(m)(𝒒)superscriptsubscript𝐵𝑚𝒒\displaystyle B_{\ell}^{(m)}(\bm{q})italic_B start_POSTSUBSCRIPT roman_ℓ end_POSTSUBSCRIPT start_POSTSUPERSCRIPT ( italic_m ) end_POSTSUPERSCRIPT ( bold_italic_q ) =i2+14π1(i)14π(21+1)ϵ1(0)absentsuperscripti214𝜋subscriptsubscript1superscriptisubscript14𝜋2subscript11subscriptsuperscriptitalic-ϵ0subscript1\displaystyle={\rm i}\hskip 0.50003pt^{\ell}\sqrt{\frac{2\ell+1}{4\pi}}\sum_{% \ell_{1}}(-{\rm i}\hskip 0.50003pt)^{\ell_{1}}\sqrt{4\pi(2\ell_{1}+1)}\epsilon% ^{(0)}_{\ell_{1}}= roman_i start_POSTSUPERSCRIPT roman_ℓ end_POSTSUPERSCRIPT square-root start_ARG divide start_ARG 2 roman_ℓ + 1 end_ARG start_ARG 4 italic_π end_ARG end_ARG ∑ start_POSTSUBSCRIPT roman_ℓ start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT end_POSTSUBSCRIPT ( - roman_i ) start_POSTSUPERSCRIPT roman_ℓ start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT end_POSTSUPERSCRIPT square-root start_ARG 4 italic_π ( 2 roman_ℓ start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT + 1 ) end_ARG italic_ϵ start_POSTSUPERSCRIPT ( 0 ) end_POSTSUPERSCRIPT start_POSTSUBSCRIPT roman_ℓ start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT end_POSTSUBSCRIPT
×d2𝒏^[Y102(𝒏^)Ym2(𝒏^)+Y102(𝒏^)Ym2(𝒏^)]\displaystyle\times\int\!\!\,{\rm d}^{2}\hat{\bm{n}}\,\,[{}_{2}Y_{\ell_{1}0}(% \hat{\bm{n}}){}_{2}Y_{\ell m}^{*}(\hat{\bm{n}})+{}_{-2}Y_{\ell_{1}0}(\hat{\bm{% n}}){}_{-2}Y_{\ell m}^{*}(\hat{\bm{n}})]× ∫ roman_d start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT over^ start_ARG bold_italic_n end_ARG [ start_FLOATSUBSCRIPT 2 end_FLOATSUBSCRIPT italic_Y start_POSTSUBSCRIPT roman_ℓ start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT ( over^ start_ARG bold_italic_n end_ARG ) start_FLOATSUBSCRIPT 2 end_FLOATSUBSCRIPT italic_Y start_POSTSUBSCRIPT roman_ℓ italic_m end_POSTSUBSCRIPT start_POSTSUPERSCRIPT ∗ end_POSTSUPERSCRIPT ( over^ start_ARG bold_italic_n end_ARG ) + start_FLOATSUBSCRIPT - 2 end_FLOATSUBSCRIPT italic_Y start_POSTSUBSCRIPT roman_ℓ start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT ( over^ start_ARG bold_italic_n end_ARG ) start_FLOATSUBSCRIPT - 2 end_FLOATSUBSCRIPT italic_Y start_POSTSUBSCRIPT roman_ℓ italic_m end_POSTSUBSCRIPT start_POSTSUPERSCRIPT ∗ end_POSTSUPERSCRIPT ( over^ start_ARG bold_italic_n end_ARG ) ]
×0η0dηgv(η)6P(0)(𝒒,η)α(η,𝒏^),\displaystyle\times\int_{0}^{\eta_{0}}\!\!\,{\rm d}\eta\,\,g_{\rm v}(\eta)% \sqrt{6}P^{(0)}(\bm{q},\eta)\alpha(\eta,\hat{\bm{n}})\,,× ∫ start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT start_POSTSUPERSCRIPT italic_η start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT end_POSTSUPERSCRIPT roman_d italic_η italic_g start_POSTSUBSCRIPT roman_v end_POSTSUBSCRIPT ( italic_η ) square-root start_ARG 6 end_ARG italic_P start_POSTSUPERSCRIPT ( 0 ) end_POSTSUPERSCRIPT ( bold_italic_q , italic_η ) italic_α ( italic_η , over^ start_ARG bold_italic_n end_ARG ) , (22)

where we use Eq. 19. The B𝐵Bitalic_B-mode power spectrum is defined as [73]

CBB=4π(2+1)2d3𝒒(2π)3m=B(m)(𝒒)B(m)(𝒒).superscriptsubscript𝐶𝐵𝐵4𝜋superscript212superscriptd3𝒒superscript2𝜋3superscriptsubscript𝑚delimited-⟨⟩superscriptsubscript𝐵𝑚𝒒superscriptsubscript𝐵𝑚𝒒\displaystyle C_{\ell}^{BB}=\frac{4\pi}{(2\ell+1)^{2}}\int\!\!\frac{\,{\rm d}^% {3}\bm{q}\,}{(2\pi)^{3}}\,\sum_{m=-\ell}^{\ell}\langle B_{\ell}^{(m)*}(\bm{q})% B_{\ell}^{(m)}(\bm{q})\rangle\,.italic_C start_POSTSUBSCRIPT roman_ℓ end_POSTSUBSCRIPT start_POSTSUPERSCRIPT italic_B italic_B end_POSTSUPERSCRIPT = divide start_ARG 4 italic_π end_ARG start_ARG ( 2 roman_ℓ + 1 ) start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT end_ARG ∫ divide start_ARG roman_d start_POSTSUPERSCRIPT 3 end_POSTSUPERSCRIPT bold_italic_q end_ARG start_ARG ( 2 italic_π ) start_POSTSUPERSCRIPT 3 end_POSTSUPERSCRIPT end_ARG ∑ start_POSTSUBSCRIPT italic_m = - roman_ℓ end_POSTSUBSCRIPT start_POSTSUPERSCRIPT roman_ℓ end_POSTSUPERSCRIPT ⟨ italic_B start_POSTSUBSCRIPT roman_ℓ end_POSTSUBSCRIPT start_POSTSUPERSCRIPT ( italic_m ) ∗ end_POSTSUPERSCRIPT ( bold_italic_q ) italic_B start_POSTSUBSCRIPT roman_ℓ end_POSTSUBSCRIPT start_POSTSUPERSCRIPT ( italic_m ) end_POSTSUPERSCRIPT ( bold_italic_q ) ⟩ . (23)

Substituting Eq. 22 into Eq. 23, the B𝐵Bitalic_B-mode power spectrum in the presence of the time-varying anisotropic cosmic birefringence is given by

CBBsubscriptsuperscript𝐶𝐵𝐵\displaystyle C^{BB}_{\ell}italic_C start_POSTSUPERSCRIPT italic_B italic_B end_POSTSUPERSCRIPT start_POSTSUBSCRIPT roman_ℓ end_POSTSUBSCRIPT =4π2+1m=dlnq𝒫(q)12i1+2(21+1)(22+1)absent4𝜋21superscriptsubscript𝑚d𝑞subscript𝒫𝑞subscriptsubscript1subscript2superscriptisubscript1subscript22subscript112subscript21\displaystyle=\frac{4\pi}{2\ell+1}\sum_{m=-\ell}^{\ell}\int\!\!\,{\rm d}\ln q% \,\,{\cal P}_{\cal R}(q)\sum_{\ell_{1}\ell_{2}}{\rm i}\hskip 0.50003pt^{-\ell_% {1}+\ell_{2}}\sqrt{(2\ell_{1}+1)(2\ell_{2}+1)}= divide start_ARG 4 italic_π end_ARG start_ARG 2 roman_ℓ + 1 end_ARG ∑ start_POSTSUBSCRIPT italic_m = - roman_ℓ end_POSTSUBSCRIPT start_POSTSUPERSCRIPT roman_ℓ end_POSTSUPERSCRIPT ∫ roman_d roman_ln italic_q caligraphic_P start_POSTSUBSCRIPT caligraphic_R end_POSTSUBSCRIPT ( italic_q ) ∑ start_POSTSUBSCRIPT roman_ℓ start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT roman_ℓ start_POSTSUBSCRIPT 2 end_POSTSUBSCRIPT end_POSTSUBSCRIPT roman_i start_POSTSUPERSCRIPT - roman_ℓ start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT + roman_ℓ start_POSTSUBSCRIPT 2 end_POSTSUBSCRIPT end_POSTSUPERSCRIPT square-root start_ARG ( 2 roman_ℓ start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT + 1 ) ( 2 roman_ℓ start_POSTSUBSCRIPT 2 end_POSTSUBSCRIPT + 1 ) end_ARG
×d2𝒏^[Y102(𝒏^)Ym2(𝒏^)+Y102(𝒏^)Ym2(𝒏^)]d2𝒏^[Y202(𝒏^)Ym2(𝒏^)+Y202(𝒏^)Ym2(𝒏^)]\displaystyle\times\int\!\!\,{\rm d}^{2}\hat{\bm{n}}\,\,[{}_{2}Y_{\ell_{1}0}^{% *}(\hat{\bm{n}}){}_{2}Y_{\ell m}(\hat{\bm{n}})+{}_{-2}Y_{\ell_{1}0}^{*}(\hat{% \bm{n}}){}_{-2}Y_{\ell m}(\hat{\bm{n}})]\int\!\!\,{\rm d}^{2}\hat{\bm{n}}^{% \prime}\,\,[{}_{2}Y_{\ell_{2}0}(\hat{\bm{n}}^{\prime}){}_{2}Y_{\ell m}^{*}(% \hat{\bm{n}}^{\prime})+{}_{-2}Y_{\ell_{2}0}(\hat{\bm{n}}^{\prime}){}_{-2}Y_{% \ell m}^{*}(\hat{\bm{n}}^{\prime})]× ∫ roman_d start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT over^ start_ARG bold_italic_n end_ARG [ start_FLOATSUBSCRIPT 2 end_FLOATSUBSCRIPT italic_Y start_POSTSUBSCRIPT roman_ℓ start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT start_POSTSUPERSCRIPT ∗ end_POSTSUPERSCRIPT ( over^ start_ARG bold_italic_n end_ARG ) start_FLOATSUBSCRIPT 2 end_FLOATSUBSCRIPT italic_Y start_POSTSUBSCRIPT roman_ℓ italic_m end_POSTSUBSCRIPT ( over^ start_ARG bold_italic_n end_ARG ) + start_FLOATSUBSCRIPT - 2 end_FLOATSUBSCRIPT italic_Y start_POSTSUBSCRIPT roman_ℓ start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT start_POSTSUPERSCRIPT ∗ end_POSTSUPERSCRIPT ( over^ start_ARG bold_italic_n end_ARG ) start_FLOATSUBSCRIPT - 2 end_FLOATSUBSCRIPT italic_Y start_POSTSUBSCRIPT roman_ℓ italic_m end_POSTSUBSCRIPT ( over^ start_ARG bold_italic_n end_ARG ) ] ∫ roman_d start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT over^ start_ARG bold_italic_n end_ARG start_POSTSUPERSCRIPT ′ end_POSTSUPERSCRIPT [ start_FLOATSUBSCRIPT 2 end_FLOATSUBSCRIPT italic_Y start_POSTSUBSCRIPT roman_ℓ start_POSTSUBSCRIPT 2 end_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT ( over^ start_ARG bold_italic_n end_ARG start_POSTSUPERSCRIPT ′ end_POSTSUPERSCRIPT ) start_FLOATSUBSCRIPT 2 end_FLOATSUBSCRIPT italic_Y start_POSTSUBSCRIPT roman_ℓ italic_m end_POSTSUBSCRIPT start_POSTSUPERSCRIPT ∗ end_POSTSUPERSCRIPT ( over^ start_ARG bold_italic_n end_ARG start_POSTSUPERSCRIPT ′ end_POSTSUPERSCRIPT ) + start_FLOATSUBSCRIPT - 2 end_FLOATSUBSCRIPT italic_Y start_POSTSUBSCRIPT roman_ℓ start_POSTSUBSCRIPT 2 end_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT ( over^ start_ARG bold_italic_n end_ARG start_POSTSUPERSCRIPT ′ end_POSTSUPERSCRIPT ) start_FLOATSUBSCRIPT - 2 end_FLOATSUBSCRIPT italic_Y start_POSTSUBSCRIPT roman_ℓ italic_m end_POSTSUBSCRIPT start_POSTSUPERSCRIPT ∗ end_POSTSUPERSCRIPT ( over^ start_ARG bold_italic_n end_ARG start_POSTSUPERSCRIPT ′ end_POSTSUPERSCRIPT ) ]
×0η0dη0η0dηs1(q,η)s2(q,η)α(η,𝒏^)α(η,𝒏^),\displaystyle\times\int_{0}^{\eta_{0}}\!\!\,{\rm d}\eta\,\,\int_{0}^{\eta_{0}}% \!\!\,{\rm d}\eta^{\prime}\,\,s_{\ell_{1}}(q,\eta)s_{\ell_{2}}(q,\eta^{\prime}% )\langle\alpha(\eta,\hat{\bm{n}})\alpha(\eta^{\prime},\hat{\bm{n}}^{\prime})% \rangle\,,× ∫ start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT start_POSTSUPERSCRIPT italic_η start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT end_POSTSUPERSCRIPT roman_d italic_η ∫ start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT start_POSTSUPERSCRIPT italic_η start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT end_POSTSUPERSCRIPT roman_d italic_η start_POSTSUPERSCRIPT ′ end_POSTSUPERSCRIPT italic_s start_POSTSUBSCRIPT roman_ℓ start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT end_POSTSUBSCRIPT ( italic_q , italic_η ) italic_s start_POSTSUBSCRIPT roman_ℓ start_POSTSUBSCRIPT 2 end_POSTSUBSCRIPT end_POSTSUBSCRIPT ( italic_q , italic_η start_POSTSUPERSCRIPT ′ end_POSTSUPERSCRIPT ) ⟨ italic_α ( italic_η , over^ start_ARG bold_italic_n end_ARG ) italic_α ( italic_η start_POSTSUPERSCRIPT ′ end_POSTSUPERSCRIPT , over^ start_ARG bold_italic_n end_ARG start_POSTSUPERSCRIPT ′ end_POSTSUPERSCRIPT ) ⟩ , (24)

where we relate the correlation of the polarization source to the primordial curvature power spectrum as

[P(0)(𝒒,η)]P(0)(𝒒,η)=2π2q3𝒫(q)S(q,η)S(q,η),delimited-⟨⟩superscriptdelimited-[]superscript𝑃0𝒒𝜂superscript𝑃0𝒒superscript𝜂2superscript𝜋2superscript𝑞3subscript𝒫𝑞𝑆𝑞𝜂𝑆𝑞superscript𝜂\displaystyle\langle[P^{(0)}(\bm{q},\eta)]^{*}P^{(0)}(\bm{q},\eta^{\prime})% \rangle=\frac{2\pi^{2}}{q^{3}}{\cal P}_{\cal R}(q)S(q,\eta)S(q,\eta^{\prime})\,,⟨ [ italic_P start_POSTSUPERSCRIPT ( 0 ) end_POSTSUPERSCRIPT ( bold_italic_q , italic_η ) ] start_POSTSUPERSCRIPT ∗ end_POSTSUPERSCRIPT italic_P start_POSTSUPERSCRIPT ( 0 ) end_POSTSUPERSCRIPT ( bold_italic_q , italic_η start_POSTSUPERSCRIPT ′ end_POSTSUPERSCRIPT ) ⟩ = divide start_ARG 2 italic_π start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT end_ARG start_ARG italic_q start_POSTSUPERSCRIPT 3 end_POSTSUPERSCRIPT end_ARG caligraphic_P start_POSTSUBSCRIPT caligraphic_R end_POSTSUBSCRIPT ( italic_q ) italic_S ( italic_q , italic_η ) italic_S ( italic_q , italic_η start_POSTSUPERSCRIPT ′ end_POSTSUPERSCRIPT ) , (25)

and define a projected polarization source function

s(q,η)gv(η)6S(η,q)ϵ(0)(q(η0η)).subscript𝑠𝑞𝜂subscript𝑔v𝜂6𝑆𝜂𝑞subscriptsuperscriptitalic-ϵ0𝑞subscript𝜂0𝜂\displaystyle s_{\ell}(q,\eta)\equiv g_{\rm v}(\eta)\sqrt{6}S(\eta,q)\epsilon^% {(0)}_{\ell}(q(\eta_{0}-\eta))\,.italic_s start_POSTSUBSCRIPT roman_ℓ end_POSTSUBSCRIPT ( italic_q , italic_η ) ≡ italic_g start_POSTSUBSCRIPT roman_v end_POSTSUBSCRIPT ( italic_η ) square-root start_ARG 6 end_ARG italic_S ( italic_η , italic_q ) italic_ϵ start_POSTSUPERSCRIPT ( 0 ) end_POSTSUPERSCRIPT start_POSTSUBSCRIPT roman_ℓ end_POSTSUBSCRIPT ( italic_q ( italic_η start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT - italic_η ) ) . (26)

Note that in deriving Eq. (24) we assume that the source of the anisotropic cosmic birefringence does not correlate with the curvature perturbations. This assumption is valid for the massless pseudoscalar field model considered in this paper but does not hold in general.

We use the following expression of the two-point correlation function of the anisotropic cosmic birefringence:

α(η,𝒏^)α(η,𝒏^)=LMYLM(𝒏^)YLM(𝒏^)CLαα(η,η).delimited-⟨⟩𝛼𝜂^𝒏𝛼superscript𝜂superscript^𝒏subscript𝐿𝑀subscriptsuperscript𝑌𝐿𝑀^𝒏subscript𝑌𝐿𝑀superscript^𝒏superscriptsubscript𝐶𝐿𝛼𝛼𝜂superscript𝜂\displaystyle\langle\alpha(\eta,\hat{\bm{n}})\alpha(\eta^{\prime},\hat{\bm{n}}% ^{\prime})\rangle=\sum_{LM}Y^{*}_{LM}(\hat{\bm{n}})Y_{LM}(\hat{\bm{n}}^{\prime% })C_{L}^{\alpha\alpha}(\eta,\eta^{\prime})\,.⟨ italic_α ( italic_η , over^ start_ARG bold_italic_n end_ARG ) italic_α ( italic_η start_POSTSUPERSCRIPT ′ end_POSTSUPERSCRIPT , over^ start_ARG bold_italic_n end_ARG start_POSTSUPERSCRIPT ′ end_POSTSUPERSCRIPT ) ⟩ = ∑ start_POSTSUBSCRIPT italic_L italic_M end_POSTSUBSCRIPT italic_Y start_POSTSUPERSCRIPT ∗ end_POSTSUPERSCRIPT start_POSTSUBSCRIPT italic_L italic_M end_POSTSUBSCRIPT ( over^ start_ARG bold_italic_n end_ARG ) italic_Y start_POSTSUBSCRIPT italic_L italic_M end_POSTSUBSCRIPT ( over^ start_ARG bold_italic_n end_ARG start_POSTSUPERSCRIPT ′ end_POSTSUPERSCRIPT ) italic_C start_POSTSUBSCRIPT italic_L end_POSTSUBSCRIPT start_POSTSUPERSCRIPT italic_α italic_α end_POSTSUPERSCRIPT ( italic_η , italic_η start_POSTSUPERSCRIPT ′ end_POSTSUPERSCRIPT ) . (27)

Substituting Eq. 27 into Eq. 24, and employing the formulas for the Wigner 3j3𝑗3j3 italic_j symbols in Appendix A, we finally obtain the B𝐵Bitalic_B-mode power spectrum as

CBB=4LpL+(2+1)(2L+1)4π(L202)2𝒞,LEE,superscriptsubscript𝐶𝐵𝐵4subscriptsuperscript𝐿subscriptsuperscript𝑝𝐿superscript2superscript12𝐿14𝜋superscriptmatrix𝐿superscript2022subscriptsuperscript𝒞𝐸𝐸superscript𝐿\displaystyle C_{\ell}^{BB}=4\sum_{\ell^{\prime}L}p^{+}_{\ell L\ell^{\prime}}% \frac{(2\ell^{\prime}+1)(2L+1)}{4\pi}\begin{pmatrix}\ell&L&\ell^{\prime}\\ 2&0&-2\end{pmatrix}^{2}{\cal C}^{EE}_{\ell^{\prime},L}\,,italic_C start_POSTSUBSCRIPT roman_ℓ end_POSTSUBSCRIPT start_POSTSUPERSCRIPT italic_B italic_B end_POSTSUPERSCRIPT = 4 ∑ start_POSTSUBSCRIPT roman_ℓ start_POSTSUPERSCRIPT ′ end_POSTSUPERSCRIPT italic_L end_POSTSUBSCRIPT italic_p start_POSTSUPERSCRIPT + end_POSTSUPERSCRIPT start_POSTSUBSCRIPT roman_ℓ italic_L roman_ℓ start_POSTSUPERSCRIPT ′ end_POSTSUPERSCRIPT end_POSTSUBSCRIPT divide start_ARG ( 2 roman_ℓ start_POSTSUPERSCRIPT ′ end_POSTSUPERSCRIPT + 1 ) ( 2 italic_L + 1 ) end_ARG start_ARG 4 italic_π end_ARG ( start_ARG start_ROW start_CELL roman_ℓ end_CELL start_CELL italic_L end_CELL start_CELL roman_ℓ start_POSTSUPERSCRIPT ′ end_POSTSUPERSCRIPT end_CELL end_ROW start_ROW start_CELL 2 end_CELL start_CELL 0 end_CELL start_CELL - 2 end_CELL end_ROW end_ARG ) start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT caligraphic_C start_POSTSUPERSCRIPT italic_E italic_E end_POSTSUPERSCRIPT start_POSTSUBSCRIPT roman_ℓ start_POSTSUPERSCRIPT ′ end_POSTSUPERSCRIPT , italic_L end_POSTSUBSCRIPT , (28)

where the distorted E𝐸Eitalic_E-mode power spectrum is defined as

𝒞,LEE4πdlnk𝒫ϕ(k)𝒞,LEE(k).subscriptsuperscript𝒞𝐸𝐸superscript𝐿4𝜋d𝑘subscript𝒫italic-ϕ𝑘subscriptsuperscript𝒞𝐸𝐸superscript𝐿𝑘\displaystyle{\cal C}^{EE}_{\ell^{\prime},L}\equiv 4\pi\int\!\!\,{\rm d}\ln k% \,\,{\cal P}_{\phi}(k){\cal C}^{EE}_{\ell^{\prime},L}(k)\,.caligraphic_C start_POSTSUPERSCRIPT italic_E italic_E end_POSTSUPERSCRIPT start_POSTSUBSCRIPT roman_ℓ start_POSTSUPERSCRIPT ′ end_POSTSUPERSCRIPT , italic_L end_POSTSUBSCRIPT ≡ 4 italic_π ∫ roman_d roman_ln italic_k caligraphic_P start_POSTSUBSCRIPT italic_ϕ end_POSTSUBSCRIPT ( italic_k ) caligraphic_C start_POSTSUPERSCRIPT italic_E italic_E end_POSTSUPERSCRIPT start_POSTSUBSCRIPT roman_ℓ start_POSTSUPERSCRIPT ′ end_POSTSUPERSCRIPT , italic_L end_POSTSUBSCRIPT ( italic_k ) . (29)

Here, we define the E𝐸Eitalic_E-mode power spectrum for each k𝑘kitalic_k as

𝒞,LEE(k)4πdlnq𝒫(q)[Δ,L(q,k)]2,subscriptsuperscript𝒞𝐸𝐸superscript𝐿𝑘4𝜋d𝑞subscript𝒫𝑞superscriptdelimited-[]subscriptΔsuperscript𝐿𝑞𝑘2\displaystyle{\cal C}^{EE}_{\ell^{\prime},L}(k)\equiv 4\pi\int\!\!\,{\rm d}\ln q% \,\,{\cal P}_{\cal R}(q)\left[\Delta_{\ell^{\prime},L}(q,k)\right]^{2}\,,caligraphic_C start_POSTSUPERSCRIPT italic_E italic_E end_POSTSUPERSCRIPT start_POSTSUBSCRIPT roman_ℓ start_POSTSUPERSCRIPT ′ end_POSTSUPERSCRIPT , italic_L end_POSTSUBSCRIPT ( italic_k ) ≡ 4 italic_π ∫ roman_d roman_ln italic_q caligraphic_P start_POSTSUBSCRIPT caligraphic_R end_POSTSUBSCRIPT ( italic_q ) [ roman_Δ start_POSTSUBSCRIPT roman_ℓ start_POSTSUPERSCRIPT ′ end_POSTSUPERSCRIPT , italic_L end_POSTSUBSCRIPT ( italic_q , italic_k ) ] start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT , (30)

with

Δ,L(q,k)=0η0dηs(q,η)uL(k,η).subscriptΔsuperscript𝐿𝑞𝑘superscriptsubscript0subscript𝜂0differential-d𝜂subscript𝑠superscript𝑞𝜂subscript𝑢𝐿𝑘𝜂\displaystyle\Delta_{\ell^{\prime},L}(q,k)=\int_{0}^{\eta_{0}}\!\!\,{\rm d}% \eta\,\,s_{\ell^{\prime}}(q,\eta)u_{L}(k,\eta)\,.roman_Δ start_POSTSUBSCRIPT roman_ℓ start_POSTSUPERSCRIPT ′ end_POSTSUPERSCRIPT , italic_L end_POSTSUBSCRIPT ( italic_q , italic_k ) = ∫ start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT start_POSTSUPERSCRIPT italic_η start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT end_POSTSUPERSCRIPT roman_d italic_η italic_s start_POSTSUBSCRIPT roman_ℓ start_POSTSUPERSCRIPT ′ end_POSTSUPERSCRIPT end_POSTSUBSCRIPT ( italic_q , italic_η ) italic_u start_POSTSUBSCRIPT italic_L end_POSTSUBSCRIPT ( italic_k , italic_η ) . (31)

Note that the thin LSS approximation leads to

Δ,L(q,k)uL(k,η)0η0dηs(q,η),similar-to-or-equalssubscriptΔsuperscript𝐿𝑞𝑘subscript𝑢𝐿𝑘subscript𝜂superscriptsubscript0subscript𝜂0differential-d𝜂subscript𝑠superscript𝑞𝜂\displaystyle\Delta_{\ell^{\prime},L}(q,k)\simeq u_{L}(k,\eta_{*})\int_{0}^{% \eta_{0}}\!\!\,{\rm d}\eta\,\,s_{\ell^{\prime}}(q,\eta)\,,roman_Δ start_POSTSUBSCRIPT roman_ℓ start_POSTSUPERSCRIPT ′ end_POSTSUPERSCRIPT , italic_L end_POSTSUBSCRIPT ( italic_q , italic_k ) ≃ italic_u start_POSTSUBSCRIPT italic_L end_POSTSUBSCRIPT ( italic_k , italic_η start_POSTSUBSCRIPT ∗ end_POSTSUBSCRIPT ) ∫ start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT start_POSTSUPERSCRIPT italic_η start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT end_POSTSUPERSCRIPT roman_d italic_η italic_s start_POSTSUBSCRIPT roman_ℓ start_POSTSUPERSCRIPT ′ end_POSTSUPERSCRIPT end_POSTSUBSCRIPT ( italic_q , italic_η ) , (32)

and we obtain

𝒞,LEEsubscriptsuperscript𝒞𝐸𝐸superscript𝐿\displaystyle{\cal C}^{EE}_{\ell^{\prime},L}caligraphic_C start_POSTSUPERSCRIPT italic_E italic_E end_POSTSUPERSCRIPT start_POSTSUBSCRIPT roman_ℓ start_POSTSUPERSCRIPT ′ end_POSTSUPERSCRIPT , italic_L end_POSTSUBSCRIPT 4πdlnk𝒫ϕ(k)[uL(k,η)]2similar-to-or-equalsabsent4𝜋d𝑘subscript𝒫italic-ϕ𝑘superscriptdelimited-[]subscript𝑢𝐿𝑘subscript𝜂2\displaystyle\simeq 4\pi\int\!\!\,{\rm d}\ln k\,\,{\cal P}_{\phi}(k)[u_{L}(k,% \eta_{*})]^{2}≃ 4 italic_π ∫ roman_d roman_ln italic_k caligraphic_P start_POSTSUBSCRIPT italic_ϕ end_POSTSUBSCRIPT ( italic_k ) [ italic_u start_POSTSUBSCRIPT italic_L end_POSTSUBSCRIPT ( italic_k , italic_η start_POSTSUBSCRIPT ∗ end_POSTSUBSCRIPT ) ] start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT
×4πdlnq𝒫(q)[0η0dηs(q,η)]2absent4𝜋d𝑞subscript𝒫𝑞superscriptdelimited-[]superscriptsubscript0subscript𝜂0differential-d𝜂subscript𝑠superscript𝑞𝜂2\displaystyle\times 4\pi\int\!\!\,{\rm d}\ln q\,\,{\cal P}_{\cal R}(q)\left[% \int_{0}^{\eta_{0}}\!\!\,{\rm d}\eta\,\,s_{\ell^{\prime}}(q,\eta)\right]^{2}× 4 italic_π ∫ roman_d roman_ln italic_q caligraphic_P start_POSTSUBSCRIPT caligraphic_R end_POSTSUBSCRIPT ( italic_q ) [ ∫ start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT start_POSTSUPERSCRIPT italic_η start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT end_POSTSUPERSCRIPT roman_d italic_η italic_s start_POSTSUBSCRIPT roman_ℓ start_POSTSUPERSCRIPT ′ end_POSTSUPERSCRIPT end_POSTSUBSCRIPT ( italic_q , italic_η ) ] start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT
=CLαα(η,η)CE¯E¯,absentsuperscriptsubscript𝐶𝐿𝛼𝛼subscript𝜂subscript𝜂superscriptsubscript𝐶superscript¯𝐸¯𝐸\displaystyle=C_{L}^{\alpha\alpha}(\eta_{*},\eta_{*})C_{\ell^{\prime}}^{\bar{E% }\bar{E}}\,,= italic_C start_POSTSUBSCRIPT italic_L end_POSTSUBSCRIPT start_POSTSUPERSCRIPT italic_α italic_α end_POSTSUPERSCRIPT ( italic_η start_POSTSUBSCRIPT ∗ end_POSTSUBSCRIPT , italic_η start_POSTSUBSCRIPT ∗ end_POSTSUBSCRIPT ) italic_C start_POSTSUBSCRIPT roman_ℓ start_POSTSUPERSCRIPT ′ end_POSTSUPERSCRIPT end_POSTSUBSCRIPT start_POSTSUPERSCRIPT over¯ start_ARG italic_E end_ARG over¯ start_ARG italic_E end_ARG end_POSTSUPERSCRIPT , (33)

where we use Eq. 14 and the expression for the primordial E𝐸Eitalic_E-mode power spectrum:

CE¯E¯=4πdlnq𝒫(q)[0η0dηs(q,η)]2.superscriptsubscript𝐶superscript¯𝐸¯𝐸4𝜋d𝑞subscript𝒫𝑞superscriptdelimited-[]superscriptsubscript0subscript𝜂0differential-d𝜂subscript𝑠superscript𝑞𝜂2\displaystyle C_{\ell^{\prime}}^{\bar{E}\bar{E}}=4\pi\int\!\!\,{\rm d}\ln q\,% \,{\cal P}_{\cal R}(q)\left[\int_{0}^{\eta_{0}}\!\!\,{\rm d}\eta\,\,s_{\ell^{% \prime}}(q,\eta)\right]^{2}\,.italic_C start_POSTSUBSCRIPT roman_ℓ start_POSTSUPERSCRIPT ′ end_POSTSUPERSCRIPT end_POSTSUBSCRIPT start_POSTSUPERSCRIPT over¯ start_ARG italic_E end_ARG over¯ start_ARG italic_E end_ARG end_POSTSUPERSCRIPT = 4 italic_π ∫ roman_d roman_ln italic_q caligraphic_P start_POSTSUBSCRIPT caligraphic_R end_POSTSUBSCRIPT ( italic_q ) [ ∫ start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT start_POSTSUPERSCRIPT italic_η start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT end_POSTSUPERSCRIPT roman_d italic_η italic_s start_POSTSUBSCRIPT roman_ℓ start_POSTSUPERSCRIPT ′ end_POSTSUPERSCRIPT end_POSTSUBSCRIPT ( italic_q , italic_η ) ] start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT . (34)

Substituting Eq. 33 into Eq. 28 coincides with the expression for the birefringence-induced B𝐵Bitalic_B-mode power spectrum with the thin LSS approximation of Eq. 13 but without the intrinsic B𝐵Bitalic_B-mode power spectrum.

III.2 Impact of the thin LSS approximation

Refer to captionRefer to caption
Figure 2: Comparison of the B𝐵Bitalic_B-mode power spectra with and without the thin approximation denoted as “exact” and “thin”, respectively (Left), and their ratio (Right). For the thin case, we also show the two-source approximation where the B𝐵Bitalic_B-mode power spectrum is the sum of the contributions from the recombination and reionization epochs. In Left, as a reference, we plot the E𝐸Eitalic_E-mode power spectrum in a black solid line and the lensing-induced B𝐵Bitalic_B-mode power spectrum in a black dashed line.
Refer to captionRefer to caption
Figure 3: The distorted E𝐸Eitalic_E-mode power spectrum defined in Eq. 29 divided by the anisotropic birefringence power spectrum. Since 𝒞,LEE/CLααCEEsubscriptsuperscript𝒞𝐸𝐸superscript𝐿superscriptsubscript𝐶𝐿𝛼𝛼superscriptsubscript𝐶superscript𝐸𝐸{\cal C}^{EE}_{\ell^{\prime},L}/C_{L}^{\alpha\alpha}\to C_{\ell^{\prime}}^{EE}caligraphic_C start_POSTSUPERSCRIPT italic_E italic_E end_POSTSUPERSCRIPT start_POSTSUBSCRIPT roman_ℓ start_POSTSUPERSCRIPT ′ end_POSTSUPERSCRIPT , italic_L end_POSTSUBSCRIPT / italic_C start_POSTSUBSCRIPT italic_L end_POSTSUBSCRIPT start_POSTSUPERSCRIPT italic_α italic_α end_POSTSUPERSCRIPT → italic_C start_POSTSUBSCRIPT roman_ℓ start_POSTSUPERSCRIPT ′ end_POSTSUPERSCRIPT end_POSTSUBSCRIPT start_POSTSUPERSCRIPT italic_E italic_E end_POSTSUPERSCRIPT in the thin approximation, we also plot the E𝐸Eitalic_E-mode power spectrum for comparison.

Here, we compare the B𝐵Bitalic_B-mode power spectra with and without the thin LSS approximation. To compute the exact B𝐵Bitalic_B-mode power spectrum, we first compute the distorted E𝐸Eitalic_E-mode power spectrum per wavenumber given in Eq. (30) by modifying the Boltzmann code CLASS [74]. We then numerically integrate 𝒞,LEE(k)subscriptsuperscript𝒞𝐸𝐸superscript𝐿𝑘{\cal C}^{EE}_{\ell^{\prime},L}(k)caligraphic_C start_POSTSUPERSCRIPT italic_E italic_E end_POSTSUPERSCRIPT start_POSTSUBSCRIPT roman_ℓ start_POSTSUPERSCRIPT ′ end_POSTSUPERSCRIPT , italic_L end_POSTSUBSCRIPT ( italic_k ) to obtain 𝒞,LEEsubscriptsuperscript𝒞𝐸𝐸superscript𝐿{\cal C}^{EE}_{\ell^{\prime},L}caligraphic_C start_POSTSUPERSCRIPT italic_E italic_E end_POSTSUPERSCRIPT start_POSTSUBSCRIPT roman_ℓ start_POSTSUPERSCRIPT ′ end_POSTSUPERSCRIPT , italic_L end_POSTSUBSCRIPT from Eq. (29). We finally obtain the B𝐵Bitalic_B-mode power spectrum from Eq. (28) where we employ WignerFamilies to compute the Wigner 3j symbols. We choose the amplitude of 𝒫ϕ(k)subscript𝒫italic-ϕ𝑘{\cal P}_{\phi}(k)caligraphic_P start_POSTSUBSCRIPT italic_ϕ end_POSTSUBSCRIPT ( italic_k ) so that ACB=3×105subscript𝐴CB3superscript105A_{\rm CB}=3\times 10^{-5}italic_A start_POSTSUBSCRIPT roman_CB end_POSTSUBSCRIPT = 3 × 10 start_POSTSUPERSCRIPT - 5 end_POSTSUPERSCRIPT which corresponds to the current bound from observations [75, 76, 77].

Figure 2 shows the B𝐵Bitalic_B-mode power spectra with and without the thin LSS approximation. We ignore the primordial B𝐵Bitalic_B-mode power spectrum in Eq. 13. The amplitude of the B𝐵Bitalic_B-mode power spectrum is suppressed by an order of magnitude at <1010\ell<10roman_ℓ < 10, and by approximately a factor of two at >100100\ell>100roman_ℓ > 100. We also show the B𝐵Bitalic_B-mode power spectrum with the thin LSS approximation, (13), but separating the contribution from the recombination and reionization, i.e.,

CBB,rec+CBB,rei=4LpL+(2+1)(2L+1)4πsuperscriptsubscript𝐶𝐵𝐵recsuperscriptsubscript𝐶𝐵𝐵rei4subscriptsuperscript𝐿subscriptsuperscript𝑝𝐿superscript2superscript12𝐿14𝜋\displaystyle C_{\ell}^{BB,\rm rec}+C_{\ell}^{BB,\rm rei}=4\sum_{\ell^{\prime}% L}p^{+}_{\ell L\ell^{\prime}}\frac{(2\ell^{\prime}+1)(2L+1)}{4\pi}italic_C start_POSTSUBSCRIPT roman_ℓ end_POSTSUBSCRIPT start_POSTSUPERSCRIPT italic_B italic_B , roman_rec end_POSTSUPERSCRIPT + italic_C start_POSTSUBSCRIPT roman_ℓ end_POSTSUBSCRIPT start_POSTSUPERSCRIPT italic_B italic_B , roman_rei end_POSTSUPERSCRIPT = 4 ∑ start_POSTSUBSCRIPT roman_ℓ start_POSTSUPERSCRIPT ′ end_POSTSUPERSCRIPT italic_L end_POSTSUBSCRIPT italic_p start_POSTSUPERSCRIPT + end_POSTSUPERSCRIPT start_POSTSUBSCRIPT roman_ℓ italic_L roman_ℓ start_POSTSUPERSCRIPT ′ end_POSTSUPERSCRIPT end_POSTSUBSCRIPT divide start_ARG ( 2 roman_ℓ start_POSTSUPERSCRIPT ′ end_POSTSUPERSCRIPT + 1 ) ( 2 italic_L + 1 ) end_ARG start_ARG 4 italic_π end_ARG
×(L202)2x=rei,recCLαα(ηx,ηx)CE¯E¯,x,absentsuperscriptmatrix𝐿superscript2022subscript𝑥reirecsuperscriptsubscript𝐶𝐿𝛼𝛼subscript𝜂𝑥subscript𝜂𝑥superscriptsubscript𝐶superscript¯𝐸¯𝐸𝑥\displaystyle\qquad\times\begin{pmatrix}\ell&L&\ell^{\prime}\\ 2&0&-2\end{pmatrix}^{2}\sum_{x={\rm rei},{\rm rec}}C_{L}^{\alpha\alpha}(\eta_{% x},\eta_{x})C_{\ell^{\prime}}^{\bar{E}\bar{E},x}\,,× ( start_ARG start_ROW start_CELL roman_ℓ end_CELL start_CELL italic_L end_CELL start_CELL roman_ℓ start_POSTSUPERSCRIPT ′ end_POSTSUPERSCRIPT end_CELL end_ROW start_ROW start_CELL 2 end_CELL start_CELL 0 end_CELL start_CELL - 2 end_CELL end_ROW end_ARG ) start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT ∑ start_POSTSUBSCRIPT italic_x = roman_rei , roman_rec end_POSTSUBSCRIPT italic_C start_POSTSUBSCRIPT italic_L end_POSTSUBSCRIPT start_POSTSUPERSCRIPT italic_α italic_α end_POSTSUPERSCRIPT ( italic_η start_POSTSUBSCRIPT italic_x end_POSTSUBSCRIPT , italic_η start_POSTSUBSCRIPT italic_x end_POSTSUBSCRIPT ) italic_C start_POSTSUBSCRIPT roman_ℓ start_POSTSUPERSCRIPT ′ end_POSTSUPERSCRIPT end_POSTSUBSCRIPT start_POSTSUPERSCRIPT over¯ start_ARG italic_E end_ARG over¯ start_ARG italic_E end_ARG , italic_x end_POSTSUPERSCRIPT , (35)

where we approximate that CE¯E¯,reisuperscriptsubscript𝐶superscript¯𝐸¯𝐸reiC_{\ell^{\prime}}^{\bar{E}\bar{E},{\rm rei}}italic_C start_POSTSUBSCRIPT roman_ℓ start_POSTSUPERSCRIPT ′ end_POSTSUPERSCRIPT end_POSTSUBSCRIPT start_POSTSUPERSCRIPT over¯ start_ARG italic_E end_ARG over¯ start_ARG italic_E end_ARG , roman_rei end_POSTSUPERSCRIPT and CE¯E¯,recsuperscriptsubscript𝐶superscript¯𝐸¯𝐸recC_{\ell^{\prime}}^{\bar{E}\bar{E},{\rm rec}}italic_C start_POSTSUBSCRIPT roman_ℓ start_POSTSUPERSCRIPT ′ end_POSTSUPERSCRIPT end_POSTSUBSCRIPT start_POSTSUPERSCRIPT over¯ start_ARG italic_E end_ARG over¯ start_ARG italic_E end_ARG , roman_rec end_POSTSUPERSCRIPT are obtained from the E𝐸Eitalic_E-mode power spectrum at 10superscript10\ell^{\prime}\leq 10roman_ℓ start_POSTSUPERSCRIPT ′ end_POSTSUPERSCRIPT ≤ 10 and 10<10superscript10<\ell^{\prime}10 < roman_ℓ start_POSTSUPERSCRIPT ′ end_POSTSUPERSCRIPT, respectively. The exact B𝐵Bitalic_B-mode power spectrum is still different from this two-source approximation.

Figure 3 shows the ratio, 𝒞,LEE/CLααsubscriptsuperscript𝒞𝐸𝐸superscript𝐿superscriptsubscript𝐶𝐿𝛼𝛼{\cal C}^{EE}_{\ell^{\prime},L}/C_{L}^{\alpha\alpha}caligraphic_C start_POSTSUPERSCRIPT italic_E italic_E end_POSTSUPERSCRIPT start_POSTSUBSCRIPT roman_ℓ start_POSTSUPERSCRIPT ′ end_POSTSUPERSCRIPT , italic_L end_POSTSUBSCRIPT / italic_C start_POSTSUBSCRIPT italic_L end_POSTSUBSCRIPT start_POSTSUPERSCRIPT italic_α italic_α end_POSTSUPERSCRIPT for different L𝐿Litalic_L. In the thin LSS limit, this ratio becomes equivalent to the E𝐸Eitalic_E-mode power spectrum. The ratio is substantially suppressed at low superscript\ell^{\prime}roman_ℓ start_POSTSUPERSCRIPT ′ end_POSTSUPERSCRIPT. If the B𝐵Bitalic_B-mode multipole is 10less-than-or-similar-to10\ell\lesssim 10roman_ℓ ≲ 10 in Eq. 28, the contribution to the B𝐵Bitalic_B-mode power spectrum comes from Lsimilar-to-or-equalssuperscript𝐿\ell^{\prime}\simeq Lroman_ℓ start_POSTSUPERSCRIPT ′ end_POSTSUPERSCRIPT ≃ italic_L due to the property of the Wigner 3j3𝑗3j3 italic_j symbol. Since the rotation field is significant at low-L𝐿Litalic_L, the significant reduction of the amplitude in the low-\ellroman_ℓ B𝐵Bitalic_B-mode power spectrum in Fig. 2 would be the result of the substantial suppression appeared in 𝒞,LEE/CLααsubscriptsuperscript𝒞𝐸𝐸superscript𝐿superscriptsubscript𝐶𝐿𝛼𝛼{\cal C}^{EE}_{\ell^{\prime},L}/C_{L}^{\alpha\alpha}caligraphic_C start_POSTSUPERSCRIPT italic_E italic_E end_POSTSUPERSCRIPT start_POSTSUBSCRIPT roman_ℓ start_POSTSUPERSCRIPT ′ end_POSTSUPERSCRIPT , italic_L end_POSTSUBSCRIPT / italic_C start_POSTSUBSCRIPT italic_L end_POSTSUBSCRIPT start_POSTSUPERSCRIPT italic_α italic_α end_POSTSUPERSCRIPT at low-superscript\ell^{\prime}roman_ℓ start_POSTSUPERSCRIPT ′ end_POSTSUPERSCRIPT. On the other hand, the B𝐵Bitalic_B-modes at high-\ellroman_ℓ are also produced by high-superscript\ell^{\prime}roman_ℓ start_POSTSUPERSCRIPT ′ end_POSTSUPERSCRIPT E𝐸Eitalic_E-modes, and the suppression in the B𝐵Bitalic_B-mode power spectrum becomes mild.

III.3 Comparison with previous works

Our results show that the B𝐵Bitalic_B-mode power spectrum is significantly suppressed if we appropriately include the time evolution of the pseudoscalar fields. On the other hand, P09 [62] shows that the B𝐵Bitalic_B-mode power spectrum at low-\ellroman_ℓ is rather enhanced significantly. We here discuss possible sources of the difference between our results and the previous attempt to compute the B𝐵Bitalic_B-mode power spectrum with time-varying anisotropic cosmic birefringence by P09.

P09 derives the B𝐵Bitalic_B-mode power spectrum by extending the formula of Ref. [71]. Choosing a Fourier wavenumber of a single scalar perturbation mode being parallel to the z𝑧zitalic_z axis, cosmic birefringence induces the Stokes U𝑈Uitalic_U polarization as U(𝒏^,η)=2α(𝒏^,η)Q¯(𝒏^,η)𝑈^𝒏𝜂2𝛼^𝒏𝜂¯𝑄^𝒏𝜂U(\hat{\bm{n}},\eta)=2\alpha(\hat{\bm{n}},\eta)\bar{Q}(\hat{\bm{n}},\eta)italic_U ( over^ start_ARG bold_italic_n end_ARG , italic_η ) = 2 italic_α ( over^ start_ARG bold_italic_n end_ARG , italic_η ) over¯ start_ARG italic_Q end_ARG ( over^ start_ARG bold_italic_n end_ARG , italic_η ) at conformal time η𝜂\etaitalic_η at the lowest order of α𝛼\alphaitalic_α. We then obtain the B𝐵Bitalic_B-mode as

Bmsubscript𝐵𝑚\displaystyle B_{\ell m}italic_B start_POSTSUBSCRIPT roman_ℓ italic_m end_POSTSUBSCRIPT =(2)!(+2)!d2𝒏^Ym(𝒏^)ð2+ð¯22U(𝒏^)absent22superscriptd2^𝒏superscriptsubscript𝑌𝑚^𝒏superscriptitalic-ð2superscript¯italic-ð22𝑈^𝒏\displaystyle=-\sqrt{\frac{(\ell-2)!}{(\ell+2)!}}\int\!\!\,{\rm d}^{2}\hat{\bm% {n}}\,\,Y_{\ell m}^{*}(\hat{\bm{n}})\frac{\eth^{2}+\bar{\eth}^{2}}{2}U(\hat{% \bm{n}})= - square-root start_ARG divide start_ARG ( roman_ℓ - 2 ) ! end_ARG start_ARG ( roman_ℓ + 2 ) ! end_ARG end_ARG ∫ roman_d start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT over^ start_ARG bold_italic_n end_ARG italic_Y start_POSTSUBSCRIPT roman_ℓ italic_m end_POSTSUBSCRIPT start_POSTSUPERSCRIPT ∗ end_POSTSUPERSCRIPT ( over^ start_ARG bold_italic_n end_ARG ) divide start_ARG italic_ð start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT + over¯ start_ARG italic_ð end_ARG start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT end_ARG start_ARG 2 end_ARG italic_U ( over^ start_ARG bold_italic_n end_ARG )
=32(2)!(+2)!d2𝒏^Ym(𝒏^)d3𝒒absent3222superscriptd2^𝒏subscriptsuperscript𝑌𝑚^𝒏superscriptd3𝒒\displaystyle=\frac{3}{2}\sqrt{\frac{(\ell-2)!}{(\ell+2)!}}\int\!\!\,{\rm d}^{% 2}\hat{\bm{n}}\,\,Y^{*}_{\ell m}(\hat{\bm{n}})\int\!\!\,{\rm d}^{3}\bm{q}\,\,= divide start_ARG 3 end_ARG start_ARG 2 end_ARG square-root start_ARG divide start_ARG ( roman_ℓ - 2 ) ! end_ARG start_ARG ( roman_ℓ + 2 ) ! end_ARG end_ARG ∫ roman_d start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT over^ start_ARG bold_italic_n end_ARG italic_Y start_POSTSUPERSCRIPT ∗ end_POSTSUPERSCRIPT start_POSTSUBSCRIPT roman_ℓ italic_m end_POSTSUBSCRIPT ( over^ start_ARG bold_italic_n end_ARG ) ∫ roman_d start_POSTSUPERSCRIPT 3 end_POSTSUPERSCRIPT bold_italic_q
×0η0dηð2+ð¯22eixμα(𝒏^,η)gv(η)Π(q,η)ξ(𝒒),\displaystyle\times\int_{0}^{\eta_{0}}\!\!\,{\rm d}\eta\,\,\frac{\eth^{2}+\bar% {\eth}^{2}}{2}\hskip 1.00006pt{\rm e}^{{\rm i}\hskip 0.35002ptx\mu}\alpha(\hat% {\bm{n}},\eta)g_{\rm v}(\eta)\Pi(q,\eta)\xi(\bm{q})\,,× ∫ start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT start_POSTSUPERSCRIPT italic_η start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT end_POSTSUPERSCRIPT roman_d italic_η divide start_ARG italic_ð start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT + over¯ start_ARG italic_ð end_ARG start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT end_ARG start_ARG 2 end_ARG roman_e start_POSTSUPERSCRIPT roman_i italic_x italic_μ end_POSTSUPERSCRIPT italic_α ( over^ start_ARG bold_italic_n end_ARG , italic_η ) italic_g start_POSTSUBSCRIPT roman_v end_POSTSUBSCRIPT ( italic_η ) roman_Π ( italic_q , italic_η ) italic_ξ ( bold_italic_q ) , (36)

where ðitalic-ð\ethitalic_ð and ð¯¯italic-ð\bar{\eth}over¯ start_ARG italic_ð end_ARG are the spin rising and lowering operators, respectively, x=q(η0η)𝑥𝑞subscript𝜂0𝜂x=q(\eta_{0}-\eta)italic_x = italic_q ( italic_η start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT - italic_η ), Π(q,η)Π𝑞𝜂\Pi(q,\eta)roman_Π ( italic_q , italic_η ) is the polarization source [71], and ξ(𝒒)𝜉𝒒\xi(\bm{q})italic_ξ ( bold_italic_q ) is the primordial fluctuations. P09 replaces the spin operators as (ð2+ð¯2)/2=m2(1+x2)2x2superscriptitalic-ð2superscript¯italic-ð22superscript𝑚2superscript1superscriptsubscript𝑥22superscript𝑥2(\eth^{2}+\bar{\eth}^{2})/2=m^{2}-(1+\partial_{x}^{2})^{2}x^{2}( italic_ð start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT + over¯ start_ARG italic_ð end_ARG start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT ) / 2 = italic_m start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT - ( 1 + ∂ start_POSTSUBSCRIPT italic_x end_POSTSUBSCRIPT start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT ) start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT italic_x start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT and obtains the B𝐵Bitalic_B-modes induced by cosmic birefringence as

Bmsubscript𝐵𝑚\displaystyle B_{\ell m}italic_B start_POSTSUBSCRIPT roman_ℓ italic_m end_POSTSUBSCRIPT =32(2)!(+2)!d2𝒏^d3𝒒d3𝒌Ym(𝒏^)absent3222superscriptd2^𝒏superscriptd3𝒒superscriptd3𝒌subscriptsuperscript𝑌𝑚^𝒏\displaystyle=\frac{3}{2}\sqrt{\frac{(\ell-2)!}{(\ell+2)!}}\int\!\!\,{\rm d}^{% 2}\hat{\bm{n}}\,\,\int\!\!\,{\rm d}^{3}\bm{q}\,\,\int\!\!\,{\rm d}^{3}\bm{k}\,% \,Y^{*}_{\ell m}(\hat{\bm{n}})= divide start_ARG 3 end_ARG start_ARG 2 end_ARG square-root start_ARG divide start_ARG ( roman_ℓ - 2 ) ! end_ARG start_ARG ( roman_ℓ + 2 ) ! end_ARG end_ARG ∫ roman_d start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT over^ start_ARG bold_italic_n end_ARG ∫ roman_d start_POSTSUPERSCRIPT 3 end_POSTSUPERSCRIPT bold_italic_q ∫ roman_d start_POSTSUPERSCRIPT 3 end_POSTSUPERSCRIPT bold_italic_k italic_Y start_POSTSUPERSCRIPT ∗ end_POSTSUPERSCRIPT start_POSTSUBSCRIPT roman_ℓ italic_m end_POSTSUBSCRIPT ( over^ start_ARG bold_italic_n end_ARG )
×0η0dη[m2(1+x2)2x2]eixμ+iyν\displaystyle\times\int_{0}^{\eta_{0}}\!\!\,{\rm d}\eta\,\,[m^{2}-(1+\partial_% {x}^{2})^{2}x^{2}]\hskip 1.00006pt{\rm e}^{{\rm i}\hskip 0.35002ptx\mu+{\rm i}% \hskip 0.35002pty\nu}× ∫ start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT start_POSTSUPERSCRIPT italic_η start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT end_POSTSUPERSCRIPT roman_d italic_η [ italic_m start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT - ( 1 + ∂ start_POSTSUBSCRIPT italic_x end_POSTSUBSCRIPT start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT ) start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT italic_x start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT ] roman_e start_POSTSUPERSCRIPT roman_i italic_x italic_μ + roman_i italic_y italic_ν end_POSTSUPERSCRIPT
×gv(η)Π(q,η)T(k,η)ξ(𝒒)gϕ2δϕ(𝒌),absentsubscript𝑔v𝜂Π𝑞𝜂𝑇𝑘𝜂𝜉𝒒subscript𝑔italic-ϕ2𝛿italic-ϕ𝒌\displaystyle\times g_{\rm v}(\eta)\Pi(q,\eta)T(k,\eta)\xi(\bm{q})\frac{g_{% \phi}}{2}\delta\phi(\bm{k})\,,× italic_g start_POSTSUBSCRIPT roman_v end_POSTSUBSCRIPT ( italic_η ) roman_Π ( italic_q , italic_η ) italic_T ( italic_k , italic_η ) italic_ξ ( bold_italic_q ) divide start_ARG italic_g start_POSTSUBSCRIPT italic_ϕ end_POSTSUBSCRIPT end_ARG start_ARG 2 end_ARG italic_δ italic_ϕ ( bold_italic_k ) , (37)

where y=k(η0η)𝑦𝑘subscript𝜂0𝜂y=k(\eta_{0}-\eta)italic_y = italic_k ( italic_η start_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT - italic_η ) and ν=𝒏^𝒌^𝜈^𝒏^𝒌\nu=\hat{\bm{n}}\cdot\hat{\bm{k}}italic_ν = over^ start_ARG bold_italic_n end_ARG ⋅ over^ start_ARG bold_italic_k end_ARG. The above expression, however, ignores the operation of the spin operator to α(𝒏^,η)𝛼^𝒏𝜂\alpha(\hat{\bm{n}},\eta)italic_α ( over^ start_ARG bold_italic_n end_ARG , italic_η ). Thus, the expression for the B𝐵Bitalic_B-mode power spectrum does not reduce to the B𝐵Bitalic_B-mode power spectrum derived in the literature even in the thin LSS approximation. This discrepancy would be a source of the difference in the behavior of the B𝐵Bitalic_B-mode power spectrum at large angular scales.

III.4 Constraints on anisotropic cosmic birefringence with B-mode power spectrum

Refer to captionRefer to caption
Figure 4: Left: The angular power spectrum from the anisotropic cosmic birefringence and lensing with parameters corresponding to the best-fit values (orange solid), the lensing-only power spectrum (black dashed), and the bandpower of the SPTpol B𝐵Bitalic_B-mode power spectrum. For reference, we also show the B𝐵Bitalic_B-mode power spectrum from the isotropic polarization rotation of 0.630.630.63\,0.63deg with lensing (green solid) where the angle is estimated from Ref. [76]. Right: The posterior distribution of ACBsubscript𝐴CBA_{\rm CB}italic_A start_POSTSUBSCRIPT roman_CB end_POSTSUBSCRIPT with 1111, 2222, and 3333σ𝜎\sigmaitalic_σ regions.

We constrain ACBsubscript𝐴CBA_{\rm CB}italic_A start_POSTSUBSCRIPT roman_CB end_POSTSUBSCRIPT with the SPTpol B𝐵Bitalic_B-mode power spectrum [78] since that has the lowest measurement errors at small angular scales. We use the public data that contains the auto and cross spectra of 150150150150 and 90909090 GHz, the bandpower window function, bandpower covariance matrix, and the seven eigenvectors of the beam covariance matrix. We follow [78] to analyze the data; we model the theoretical B𝐵Bitalic_B-mode power spectrum that includes contributions from the anisotropic cosmic birefringence, lensing, Galactic dust, and Poisson sources. We then multiply the overall calibration parameters at each frequency and correct the beam error. In total, we have 15151515 parameters; ACBsubscript𝐴CBA_{\rm CB}italic_A start_POSTSUBSCRIPT roman_CB end_POSTSUBSCRIPT, the overall lensing amplitude, Galactic dust amplitude, Poisson source amplitude for the three auto and cross spectra, calibration parameters for the two frequencies, and the seven nuisance parameters for the beam error. We add the same Gaussian priors for the dust amplitude, calibration parameters, and beam errors as in [78]. We further add a Gaussian prior to the lensing amplitude with the mean of 1111 and the standard deviation of 0.0250.0250.0250.025 which corresponds to the current best uncertainty on the lensing amplitude by Planck [79]. We use the Markov Chain Monte Carlo (MCMC) package of emcee to fit the bandpowers to the above model.

Figure 4 shows the results of our constraint on ACBsubscript𝐴CBA_{\rm CB}italic_A start_POSTSUBSCRIPT roman_CB end_POSTSUBSCRIPT. We also show the angular power spectrum from the cosmic birefringence and lensing with parameters corresponding to the best-fit values. We obtain the 2σ2𝜎2\,\sigma2 italic_σ constraint as ACB×104=1.030.97+0.91subscript𝐴CBsuperscript104subscriptsuperscript1.030.910.97A_{\rm CB}\times 10^{4}=1.03^{+0.91}_{-0.97}italic_A start_POSTSUBSCRIPT roman_CB end_POSTSUBSCRIPT × 10 start_POSTSUPERSCRIPT 4 end_POSTSUPERSCRIPT = 1.03 start_POSTSUPERSCRIPT + 0.91 end_POSTSUPERSCRIPT start_POSTSUBSCRIPT - 0.97 end_POSTSUBSCRIPT. If we increase the standard deviation of the lensing amplitude in the Gaussian prior to 0.10.10.10.1, the constraint becomes ACB×104=1.131.10+1.06subscript𝐴CBsuperscript104subscriptsuperscript1.131.061.10A_{\rm CB}\times 10^{4}=1.13^{+1.06}_{-1.10}italic_A start_POSTSUBSCRIPT roman_CB end_POSTSUBSCRIPT × 10 start_POSTSUPERSCRIPT 4 end_POSTSUPERSCRIPT = 1.13 start_POSTSUPERSCRIPT + 1.06 end_POSTSUPERSCRIPT start_POSTSUBSCRIPT - 1.10 end_POSTSUBSCRIPT at 2σ2𝜎2\sigma2 italic_σ. Similar to Ref. [78], even if we decrease or increase the prior on the dust amplitude by 50%percent5050\%50 %, the constraint is almost unchanged and is mostly driven by the 150150150150 GHz auto power spectrum.

IV Summary and Discussion

We derived the B𝐵Bitalic_B-mode power spectrum induced by anisotropic cosmic birefringence without the thin LSS approximation. We numerically computed the B𝐵Bitalic_B-mode power spectrum and found that the power spectrum is suppressed by an order of magnitude at 10less-than-or-similar-to10\ell\lesssim 10roman_ℓ ≲ 10 and a factor of two at 100greater-than-or-equivalent-to100\ell\gtrsim 100roman_ℓ ≳ 100 compared to that with the thin LSS approximation. We also compared the results of P09 and the discrepancy in the behavior of the B𝐵Bitalic_B-mode power spectrum at low-\ellroman_ℓ would be the ignorance of the spin-operator to the rotation angle which also depends on the line-of-sight direction. We constrained the anisotropic birefringence with the SPTpol B𝐵Bitalic_B-mode power spectrum and obtained ACB×104=1.030.97+0.91(2σ)subscript𝐴CBsuperscript104subscriptsuperscript1.030.910.972𝜎A_{\rm CB}\times 10^{4}=1.03^{+0.91}_{-0.97}\,(2\sigma)italic_A start_POSTSUBSCRIPT roman_CB end_POSTSUBSCRIPT × 10 start_POSTSUPERSCRIPT 4 end_POSTSUPERSCRIPT = 1.03 start_POSTSUPERSCRIPT + 0.91 end_POSTSUPERSCRIPT start_POSTSUBSCRIPT - 0.97 end_POSTSUBSCRIPT ( 2 italic_σ ).

The SPTpol B𝐵Bitalic_B-mode power spectrum slightly prefers a nonzero amplitude of the anisotropic birefringence at more than 2σ2𝜎2\sigma2 italic_σ statistical significance. The results indicate that the SPTpol B𝐵Bitalic_B-mode power spectrum contains unknown systematics that was not characterized in Ref. [78], and/or suggests new physics.

For the former possibility, any E𝐸Eitalic_E-to-B𝐵Bitalic_B leakage could bias the B𝐵Bitalic_B-mode power spectrum. In Ref. [78], E𝐸Eitalic_E-to-B𝐵Bitalic_B leakages are mitigated with the measured EB𝐸𝐵EBitalic_E italic_B cross-power spectrum. Noise fluctuations in the EB𝐸𝐵EBitalic_E italic_B power spectrum could, however, lead to residuals of the leakage in the B𝐵Bitalic_B-mode power spectrum. One of the possible sources of the E𝐸Eitalic_E-to-B𝐵Bitalic_B leakage that can mimic the effect of anisotropic cosmic birefringence is the miscalibration angle. In Fig. 4, we also show the B𝐵Bitalic_B-mode power spectrum induced by the isotropic rotation of 0.630.630.63\,0.63deg which corresponds to the value estimated in Ref. [76]. The E𝐸Eitalic_E-to-B𝐵Bitalic_B leakage from this rotation effect is simultaneously mitigated by the method described in Ref. [78], but some residuals might still exist in the B𝐵Bitalic_B-mode power spectrum. In the parameter search, if we replace the B𝐵Bitalic_B-mode power spectrum from the anisotropic cosmic birefringence with that from the isotropic rotation, we obtain 0.840.840.84\,0.84deg for the best-fit value of the isotropic rotation angle. Thus, to explain the slight preference of the nonzero ACBsubscript𝐴CBA_{\rm CB}italic_A start_POSTSUBSCRIPT roman_CB end_POSTSUBSCRIPT by the residual B𝐵Bitalic_B-modes from isotropic polarization rotation, we need miscalibration-induced B𝐵Bitalic_B-modes comparable to that estimated in Ref. [76].

For the latter possibility, we note that ACBsubscript𝐴CBA_{\rm CB}italic_A start_POSTSUBSCRIPT roman_CB end_POSTSUBSCRIPT has been constrained by directly reconstructing the anisotropic birefringence angle from the correlations between CMB E𝐸Eitalic_E- and B𝐵Bitalic_B-modes at different angular scales [80, 81, 82]. The best constraint to date on ACBsubscript𝐴CBA_{\rm CB}italic_A start_POSTSUBSCRIPT roman_CB end_POSTSUBSCRIPT from the reconstruction is an order of magnitude smaller than that obtained in this work [83, 84, 85, 86, 75, 87, 88, 89, 77]. The constraints from the reconstruction, however, have ignored the effect of the time evolution of the anisotropic cosmic birefringence and are currently not possible to directly compare with our result. We leave a further investigation for the consistency of the constraint on anisotropic cosmic birefringence from the different observables for our future work.

Acknowledgements.
We thank Eiichiro Komatsu and Thomas Crawford for the helpful comments. This work is supported in part by JSPS KAKENHI Grant No. JP20H05859 and No. JP22K03682. Part of this work uses resources of the National Energy Research Scientific Computing Center (NERSC). The Kavli IPMU is supported by World Premier International Research Center Initiative (WPI Initiative), MEXT, Japan.

Appendix A Useful formula

The spin-weighted spherical harmonics are related to the Wigner 3j3𝑗3j3 italic_j symbols as [90]

d2𝒏^Y1m1s1(𝒏^)Y2m2s2(𝒏^)Y3m3s3(𝒏^)=(21+1)(22+1)(23+1)4π(123s1s2s3)(123m1m2m3).superscriptd2^𝒏subscriptsubscript𝑌subscript1subscript𝑚1subscript𝑠1^𝒏subscriptsubscript𝑌subscript2subscript𝑚2subscript𝑠2^𝒏subscriptsubscript𝑌subscript3subscript𝑚3subscript𝑠3^𝒏2subscript112subscript212subscript314𝜋matrixsubscript1subscript2subscript3subscript𝑠1subscript𝑠2subscript𝑠3matrixsubscript1subscript2subscript3subscript𝑚1subscript𝑚2subscript𝑚3\displaystyle\int\!\!\,{\rm d}^{2}\hat{\bm{n}}\,\,{}_{s_{1}}Y_{\ell_{1}m_{1}}(% \hat{\bm{n}}){}_{s_{2}}Y_{\ell_{2}m_{2}}(\hat{\bm{n}}){}_{s_{3}}Y_{\ell_{3}m_{% 3}}(\hat{\bm{n}})=\sqrt{\frac{(2\ell_{1}+1)(2\ell_{2}+1)(2\ell_{3}+1)}{4\pi}}% \begin{pmatrix}\ell_{1}&\ell_{2}&\ell_{3}\\ -s_{1}&-s_{2}&-s_{3}\end{pmatrix}\begin{pmatrix}\ell_{1}&\ell_{2}&\ell_{3}\\ m_{1}&m_{2}&m_{3}\end{pmatrix}\,.∫ roman_d start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT over^ start_ARG bold_italic_n end_ARG start_FLOATSUBSCRIPT italic_s start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT end_FLOATSUBSCRIPT italic_Y start_POSTSUBSCRIPT roman_ℓ start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT italic_m start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT end_POSTSUBSCRIPT ( over^ start_ARG bold_italic_n end_ARG ) start_FLOATSUBSCRIPT italic_s start_POSTSUBSCRIPT 2 end_POSTSUBSCRIPT end_FLOATSUBSCRIPT italic_Y start_POSTSUBSCRIPT roman_ℓ start_POSTSUBSCRIPT 2 end_POSTSUBSCRIPT italic_m start_POSTSUBSCRIPT 2 end_POSTSUBSCRIPT end_POSTSUBSCRIPT ( over^ start_ARG bold_italic_n end_ARG ) start_FLOATSUBSCRIPT italic_s start_POSTSUBSCRIPT 3 end_POSTSUBSCRIPT end_FLOATSUBSCRIPT italic_Y start_POSTSUBSCRIPT roman_ℓ start_POSTSUBSCRIPT 3 end_POSTSUBSCRIPT italic_m start_POSTSUBSCRIPT 3 end_POSTSUBSCRIPT end_POSTSUBSCRIPT ( over^ start_ARG bold_italic_n end_ARG ) = square-root start_ARG divide start_ARG ( 2 roman_ℓ start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT + 1 ) ( 2 roman_ℓ start_POSTSUBSCRIPT 2 end_POSTSUBSCRIPT + 1 ) ( 2 roman_ℓ start_POSTSUBSCRIPT 3 end_POSTSUBSCRIPT + 1 ) end_ARG start_ARG 4 italic_π end_ARG end_ARG ( start_ARG start_ROW start_CELL roman_ℓ start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT end_CELL start_CELL roman_ℓ start_POSTSUBSCRIPT 2 end_POSTSUBSCRIPT end_CELL start_CELL roman_ℓ start_POSTSUBSCRIPT 3 end_POSTSUBSCRIPT end_CELL end_ROW start_ROW start_CELL - italic_s start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT end_CELL start_CELL - italic_s start_POSTSUBSCRIPT 2 end_POSTSUBSCRIPT end_CELL start_CELL - italic_s start_POSTSUBSCRIPT 3 end_POSTSUBSCRIPT end_CELL end_ROW end_ARG ) ( start_ARG start_ROW start_CELL roman_ℓ start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT end_CELL start_CELL roman_ℓ start_POSTSUBSCRIPT 2 end_POSTSUBSCRIPT end_CELL start_CELL roman_ℓ start_POSTSUBSCRIPT 3 end_POSTSUBSCRIPT end_CELL end_ROW start_ROW start_CELL italic_m start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT end_CELL start_CELL italic_m start_POSTSUBSCRIPT 2 end_POSTSUBSCRIPT end_CELL start_CELL italic_m start_POSTSUBSCRIPT 3 end_POSTSUBSCRIPT end_CELL end_ROW end_ARG ) . (38)

In particular, we obtain

ImML1subscriptsuperscript𝐼𝐿subscript1𝑚𝑀\displaystyle I^{\ell L\ell_{1}}_{mM}italic_I start_POSTSUPERSCRIPT roman_ℓ italic_L roman_ℓ start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT end_POSTSUPERSCRIPT start_POSTSUBSCRIPT italic_m italic_M end_POSTSUBSCRIPT d2𝒏^[Y102(𝒏^)Ym2(𝒏^)+Y102(𝒏^)Ym2(𝒏^)]YLM(𝒏^)absentsuperscriptd2^𝒏delimited-[]subscriptsubscriptsuperscript𝑌subscript102^𝒏subscriptsubscript𝑌𝑚2^𝒏subscriptsubscriptsuperscript𝑌subscript102^𝒏subscriptsubscript𝑌𝑚2^𝒏subscriptsuperscript𝑌𝐿𝑀^𝒏\displaystyle\equiv\int\!\!\,{\rm d}^{2}\hat{\bm{n}}\,\,[{}_{2}Y^{*}_{\ell_{1}% 0}(\hat{\bm{n}}){}_{2}Y_{\ell m}(\hat{\bm{n}})+{}_{-2}Y^{*}_{\ell_{1}0}(\hat{% \bm{n}}){}_{-2}Y_{\ell m}(\hat{\bm{n}})]Y^{*}_{LM}(\hat{\bm{n}})≡ ∫ roman_d start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT over^ start_ARG bold_italic_n end_ARG [ start_FLOATSUBSCRIPT 2 end_FLOATSUBSCRIPT italic_Y start_POSTSUPERSCRIPT ∗ end_POSTSUPERSCRIPT start_POSTSUBSCRIPT roman_ℓ start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT ( over^ start_ARG bold_italic_n end_ARG ) start_FLOATSUBSCRIPT 2 end_FLOATSUBSCRIPT italic_Y start_POSTSUBSCRIPT roman_ℓ italic_m end_POSTSUBSCRIPT ( over^ start_ARG bold_italic_n end_ARG ) + start_FLOATSUBSCRIPT - 2 end_FLOATSUBSCRIPT italic_Y start_POSTSUPERSCRIPT ∗ end_POSTSUPERSCRIPT start_POSTSUBSCRIPT roman_ℓ start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT 0 end_POSTSUBSCRIPT ( over^ start_ARG bold_italic_n end_ARG ) start_FLOATSUBSCRIPT - 2 end_FLOATSUBSCRIPT italic_Y start_POSTSUBSCRIPT roman_ℓ italic_m end_POSTSUBSCRIPT ( over^ start_ARG bold_italic_n end_ARG ) ] italic_Y start_POSTSUPERSCRIPT ∗ end_POSTSUPERSCRIPT start_POSTSUBSCRIPT italic_L italic_M end_POSTSUBSCRIPT ( over^ start_ARG bold_italic_n end_ARG )
=2pL1+(1)M(21+1)(2+1)(2L+1)4π(L1202)(L1mM0).absent2subscriptsuperscript𝑝𝐿subscript1superscript1𝑀2subscript11212𝐿14𝜋matrix𝐿subscript1202matrix𝐿subscript1𝑚𝑀0\displaystyle=2p^{+}_{\ell L\ell_{1}}(-1)^{M}\sqrt{\frac{(2\ell_{1}+1)(2\ell+1% )(2L+1)}{4\pi}}\begin{pmatrix}\ell&L&\ell_{1}\\ -2&0&2\end{pmatrix}\begin{pmatrix}\ell&L&\ell_{1}\\ m&-M&0\end{pmatrix}\,.= 2 italic_p start_POSTSUPERSCRIPT + end_POSTSUPERSCRIPT start_POSTSUBSCRIPT roman_ℓ italic_L roman_ℓ start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT end_POSTSUBSCRIPT ( - 1 ) start_POSTSUPERSCRIPT italic_M end_POSTSUPERSCRIPT square-root start_ARG divide start_ARG ( 2 roman_ℓ start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT + 1 ) ( 2 roman_ℓ + 1 ) ( 2 italic_L + 1 ) end_ARG start_ARG 4 italic_π end_ARG end_ARG ( start_ARG start_ROW start_CELL roman_ℓ end_CELL start_CELL italic_L end_CELL start_CELL roman_ℓ start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT end_CELL end_ROW start_ROW start_CELL - 2 end_CELL start_CELL 0 end_CELL start_CELL 2 end_CELL end_ROW end_ARG ) ( start_ARG start_ROW start_CELL roman_ℓ end_CELL start_CELL italic_L end_CELL start_CELL roman_ℓ start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT end_CELL end_ROW start_ROW start_CELL italic_m end_CELL start_CELL - italic_M end_CELL start_CELL 0 end_CELL end_ROW end_ARG ) . (39)

Using

mM(L1mM0)(L2mM0)=121+1δ12,subscript𝑚𝑀matrix𝐿subscript1𝑚𝑀0matrix𝐿subscript2𝑚𝑀012subscript11subscript𝛿subscript1subscript2\displaystyle\sum_{mM}\begin{pmatrix}\ell&L&\ell_{1}\\ m&-M&0\end{pmatrix}\begin{pmatrix}\ell&L&\ell_{2}\\ m&-M&0\end{pmatrix}=\frac{1}{2\ell_{1}+1}\delta_{\ell_{1}\ell_{2}}\,,∑ start_POSTSUBSCRIPT italic_m italic_M end_POSTSUBSCRIPT ( start_ARG start_ROW start_CELL roman_ℓ end_CELL start_CELL italic_L end_CELL start_CELL roman_ℓ start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT end_CELL end_ROW start_ROW start_CELL italic_m end_CELL start_CELL - italic_M end_CELL start_CELL 0 end_CELL end_ROW end_ARG ) ( start_ARG start_ROW start_CELL roman_ℓ end_CELL start_CELL italic_L end_CELL start_CELL roman_ℓ start_POSTSUBSCRIPT 2 end_POSTSUBSCRIPT end_CELL end_ROW start_ROW start_CELL italic_m end_CELL start_CELL - italic_M end_CELL start_CELL 0 end_CELL end_ROW end_ARG ) = divide start_ARG 1 end_ARG start_ARG 2 roman_ℓ start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT + 1 end_ARG italic_δ start_POSTSUBSCRIPT roman_ℓ start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT roman_ℓ start_POSTSUBSCRIPT 2 end_POSTSUBSCRIPT end_POSTSUBSCRIPT , (40)

we obtain

mMImML1ImML2=4pL1+(2+1)(2L+1)4π(L1202)2δ12.subscript𝑚𝑀subscriptsuperscript𝐼𝐿subscript1𝑚𝑀subscriptsuperscript𝐼𝐿subscript2𝑚𝑀4subscriptsuperscript𝑝𝐿subscript1212𝐿14𝜋superscriptmatrix𝐿subscript12022subscript𝛿subscript1subscript2\displaystyle\sum_{mM}I^{\ell L\ell_{1}}_{mM}I^{\ell L\ell_{2}}_{mM}=4p^{+}_{% \ell L\ell_{1}}\frac{(2\ell+1)(2L+1)}{4\pi}\begin{pmatrix}\ell&L&\ell_{1}\\ -2&0&2\end{pmatrix}^{2}\delta_{\ell_{1}\ell_{2}}\,.∑ start_POSTSUBSCRIPT italic_m italic_M end_POSTSUBSCRIPT italic_I start_POSTSUPERSCRIPT roman_ℓ italic_L roman_ℓ start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT end_POSTSUPERSCRIPT start_POSTSUBSCRIPT italic_m italic_M end_POSTSUBSCRIPT italic_I start_POSTSUPERSCRIPT roman_ℓ italic_L roman_ℓ start_POSTSUBSCRIPT 2 end_POSTSUBSCRIPT end_POSTSUPERSCRIPT start_POSTSUBSCRIPT italic_m italic_M end_POSTSUBSCRIPT = 4 italic_p start_POSTSUPERSCRIPT + end_POSTSUPERSCRIPT start_POSTSUBSCRIPT roman_ℓ italic_L roman_ℓ start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT end_POSTSUBSCRIPT divide start_ARG ( 2 roman_ℓ + 1 ) ( 2 italic_L + 1 ) end_ARG start_ARG 4 italic_π end_ARG ( start_ARG start_ROW start_CELL roman_ℓ end_CELL start_CELL italic_L end_CELL start_CELL roman_ℓ start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT end_CELL end_ROW start_ROW start_CELL - 2 end_CELL start_CELL 0 end_CELL start_CELL 2 end_CELL end_ROW end_ARG ) start_POSTSUPERSCRIPT 2 end_POSTSUPERSCRIPT italic_δ start_POSTSUBSCRIPT roman_ℓ start_POSTSUBSCRIPT 1 end_POSTSUBSCRIPT roman_ℓ start_POSTSUBSCRIPT 2 end_POSTSUBSCRIPT end_POSTSUBSCRIPT . (41)

References

  • [1] W.-T. Ni Prog. Theor. Phys. Suppl. 172 (2008) 49–60, arXiv:0712.4082.
  • [2] E. Komatsu Nature Rev. Phys. 4 (2022), no. 7 452–469, arXiv:2202.13919.
  • [3] Y. Minami and E. Komatsu Phys. Rev. Lett.  125 (2020), no. 22 221301, arXiv:2011.11254.
  • [4] P. Diego-Palazuelos et al. Phys. Rev. Lett.  128 (2022), no. 9 091302, arXiv:2201.07682.
  • [5] J. R. Eskilt Astron. Astrophys.  662 (2022) A10, arXiv:2201.13347.
  • [6] J. R. Eskilt and E. Komatsu Phys. Rev. D 106 (2022), no. 6 063503, arXiv:2205.13962.
  • [7] J. R. Eskilt et al. arXiv:2305.02268.
  • [8] POLARBEAR Collaboration , S. Adachi et al. arXiv:2403.02096.
  • [9] S. M. Carroll Phys. Rev. Lett.  81 (1998) 3067–3070, astro-ph/9806099.
  • [10] S. Panda, Y. Sumitomo, and S. P. Trivedi Phys. Rev. D 83 (2011) 083506, arXiv:1011.5877.
  • [11] T. Fujita, Y. Minami, K. Murai, and H. Nakatsuka Phys. Rev. D 103 (2021), no. 6 063508, arXiv:2008.02473.
  • [12] T. Fujita, K. Murai, H. Nakatsuka, and S. Tsujikawa Phys. Rev. D 103 (2021), no. 4 043509, arXiv:2011.11894.
  • [13] G. Choi, W. Lin, L. Visinelli, and T. T. Yanagida Phys. Rev. D 104 (2021), no. 10 L101302, arXiv:2106.12602.
  • [14] I. Obata J. Cosmol. Astropart. Phys.  09 (2022) 062, arXiv:2108.02150.
  • [15] S. Gasparotto and I. Obata J. Cosmol. Astropart. Phys.  08 (2022), no. 08 025, arXiv:2203.09409.
  • [16] M. Galaverni, F. Finelli, and D. Paoletti Phys. Rev. D 107 (2023), no. 8 083529, arXiv:2301.07971.
  • [17] K. Murai, F. Naokawa, T. Namikawa, and E. Komatsu Phys. Rev. D 107 (2023) L041302, arXiv:2209.07804.
  • [18] J. R. Eskilt, L. Herold, E. Komatsu, K. Murai, T. Namikawa, and F. Naokawa Phys. Rev. Lett.  131 (2023), no. 12 121001, arXiv:2303.15369.
  • [19] F. Finelli and M. Galaverni Phys. Rev. D 79 (2009) 063002, arXiv:0802.4210.
  • [20] G.-C. Liu and K.-W. Ng Phys. Dark Univ. 16 (2017) 22–25, arXiv:1612.02104.
  • [21] M. A. Fedderke, P. W. Graham, and S. Rajendran Phys. Rev. D 100 (2019) 015040, arXiv:1903.02666.
  • [22] F. Takahashi and W. Yin J. Cosmol. Astropart. Phys.  04 (2021) 007, arXiv:2012.11576.
  • [23] N. Kitajima, F. Kozai, F. Takahashi, and W. Yin J. Cosmol. Astropart. Phys.  10 (2022) 043, arXiv:2205.05083.
  • [24] M. Jain, R. Hagimoto, A. J. Long, and M. A. Amin J. Cosmol. Astropart. Phys.  10 (2022) 090, arXiv:2208.08391.
  • [25] D. Gonzalez, N. Kitajima, F. Takahashi, and W. Yin Phys. Lett. B 843 (11, 2023) 137990, arXiv:2211.06849.
  • [26] R. C. Myers and M. Pospelov Phys. Rev. Lett.  90 (2003) 211601, hep-ph/0301124.
  • [27] A. Arvanitaki, S. Dimopoulos, S. Dubovsky, N. Kaloper, and J. March-Russell Phys. Rev. D 81 (2010) 123530, arXiv:0905.4720.
  • [28] J. Cornelison et al. Proc. SPIE Int. Soc. Opt. Eng. 12190 (2022) 121901X, arXiv:2207.14796.
  • [29] L. Moncelsi, P. A. R. Ade, Z. Ahmed, M. Amiri, D. Barkats, R. Basu Thakur, C. A. Bischoff, J. J. Bock, V. Buza, J. R. Cheshire, et al. Proc. SPIE Int. Soc. Opt. Eng. 11453 (2020) 1145314, arXiv:2012.04047.
  • [30] A. Suzuki et al. J. Low Temp. Phys. 184 (2016), no. 3-4 805–810, arXiv:1512.07299.
  • [31] The Simons Observatory Collaboration J. Cosmol. Astropart. Phys.  02 (2019) 056, arXiv:1808.07445.
  • [32] CMB-S4 Collaboration Astrophys. J.  926 (2022), no. 1 54, arXiv:2008.12619.
  • [33] LiteBIRD Collaboration Prog. Theor. Exp. Phys.  (2, 2022) arXiv:2202.02773.
  • [34] G. Sigl and P. Trivedi (11, 2018) arXiv:1811.07873.
  • [35] B. D. Sherwin and T. Namikawa Mon. Not. R. Astron. Soc.  520 (2021) 3298–3304, arXiv:2108.09287.
  • [36] H. Nakatsuka, T. Namikawa, and E. Komatsu Phys. Rev. D 105 (2022) 123509, arXiv:2203.08560.
  • [37] G.-C. Liu, S. Lee, and K.-W. Ng Phys. Rev. Lett.  97 (2006) 161303, astro-ph/0606248.
  • [38] S. Lee, G.-C. Liu, and K.-W. Ng Phys. Rev. D 89 (2014), no. 6 063010, arXiv:1307.6298.
  • [39] G. Gubitosi, M. Martinelli, and L. Pagano J. Cosmol. Astropart. Phys.  12 (2014), no. 2014 020, arXiv:1410.1799.
  • [40] N. Lee, S. C. Hotinli, and M. Kamionkowski Phys. Rev. D 106 (2022), no. 8 083518, arXiv:2207.05687.
  • [41] T. Namikawa and I. Obata Phys. Rev. D 108 (2023), no. 8 083510, arXiv:2306.08875.
  • [42] S. M. Carroll and G. B. Field Phys. Rev. Lett.  79 (1997) 2394–2397, astro-ph/9704263.
  • [43] W. W. Yin, L. Dai, J. Huang, L. Ji, and S. Ferraro, “A New Probe of Cosmic Birefringence Using Galaxy Polarization and Shapes”, 2024.
  • [44] E. Y. S. Wu et al. Phys. Rev. Lett.  102 (2009) 161302, arXiv:0811.0618.
  • [45] E. Komatsu et al. Astrophys. J.  192 (2011) 18, arXiv:1001.4538.
  • [46] J. P. Kaufman, N. J. Miller, M. Shimon, D. Barkats, C. Bischoff, I. Buder, B. G. Keating, J. M. Kovac, et al. Phys. Rev. D 89 (2014) 062006, arXiv:1312.7877.
  • [47] Planck Collaboration Astron. Astrophys.  596 (2016) A13, arXiv:1605.08633.
  • [48] M. Monelli, E. Komatsu, A. E. Adler, M. Billi, P. Campeti, N. Dachlythra, A. J. Duivenvoorden, J. E. Gudmundsson, and M. Reinecke arXiv:2211.05685.
  • [49] A. Lue, L. Wang, and M. Kamionkowski Phys. Rev. Lett.  83 (1999) 1506–1509, astro-ph/9812088.
  • [50] R. R. Caldwell, V. Gluscevic, and M. Kamionkowski Phys. Rev. D 84 (2011) 043504, arXiv:1104.1634.
  • [51] S. Lee, G.-C. Liu, and K.-W. Ng Phys. Lett. B 746 (2015) 406–409, arXiv:1403.5585.
  • [52] D. Leon, J. Kaufman, B. Keating, and M. Mewes Mod. Phys. Lett. A 32 (2017) 1730002, arXiv:1611.00418.
  • [53] L. Yin, J. Kochappan, T. Ghosh, and B.-H. Lee J. Cosmol. Astropart. Phys.  10 (2023) 007, arXiv:2305.07937.
  • [54] R. Z. Ferreira, S. Gasparotto, T. Hiramatsu, I. Obata, and O. Pujolas arXiv:2312.14104.
  • [55] M. Li and X. Zhang Phys. Rev. D 78 (2008) 103516, arXiv:0810.0403.
  • [56] M. Li and B. Yu J. Cosmol. Astropart. Phys.  06 (2013) 016, arXiv:1303.1881.
  • [57] W. Zhao and M. Li Phys. Rev. D 89 (2014), no. 10 103518, arXiv:1403.3997.
  • [58] S. Lee, G.-C. Liu, and K.-W. Ng The Universe 4 (2016), no. 4 29–44, arXiv:1912.12903.
  • [59] L. Pogosian, M. Shimon, M. Mewes, and B. Keating Phys. Rev. D 100 (2019), no. 2 023507, arXiv:1904.07855.
  • [60] H. Cai and Y. Guan Phys. Rev. D 105 (2022), no. 6 063536, arXiv:2111.14199.
  • [61] T. Namikawa, A. Naruko, R. Saito, A. Taruya, and D. Yamauchi J. Cosmol. Astropart. Phys.  10 (2021) 029, arXiv:2103.10639.
  • [62] M. Pospelov, A. Ritz, and C. Skordis Phys. Rev. Lett.  103 (2009) 051302, arXiv:0808.0673.
  • [63] A. Greco, N. Bartolo, and A. Gruppuso arXiv:2401.07079.
  • [64] D. Harari and P. Sikivie Phys. Rev. B 289 (1992) 67–72.
  • [65] S. M. Carroll, G. B. Field, and R. Jackiw Phys. Rev. D 41 (1990) 1231.
  • [66] S. M. Carroll and G. B. Field Phys. Rev. D 43 (Jun, 1991) 3789–3793.
  • [67] T. Namikawa Mon. Not. R. Astron. Soc.  506 (2021), no. 1 1250–1257, arXiv:2105.03367.
  • [68] F. Naokawa and T. Namikawa Phys. Rev. D 108 (2023), no. 6 063525, arXiv:2305.13976.
  • [69] A. Lewis and A. Challinor Phys. Rep. 429 (2006) 1–65, astro-ph/0601594.
  • [70] A. Greco, N. Bartolo, and A. Gruppuso J. Cosmol. Astropart. Phys.  05 (2023) 026, arXiv:2211.06380.
  • [71] M. Zaldarriaga and U. Seljak Phys. Rev. D 55 (1997) 1830–1840, astro-ph/9609170.
  • [72] L. Pogosian, A. P. S. Yadav, Y.-F. Ng, and T. Vachaspati Phys. Rev. D 84 (2011) 043530, arXiv:1106.1438. [Erratum: Phys.Rev.D 84, 089903 (2011)].
  • [73] W. Hu and M. J. White Phys. Rev. D 56 (1997) 596–615, astro-ph/9702170.
  • [74] D. Blas, J. Lesgourgues, and T. Tram J. Cosmol. Astropart. Phys.  2011 (2011), no. 7 034, arXiv:1104.2933.
  • [75] T. Namikawa et al. Phys. Rev. D 101 (2020), no. 8 083527, arXiv:2001.10465.
  • [76] F. Bianchini et al. Phys. Rev. D 102 (2020), no. 8 083504, arXiv:2006.08061.
  • [77] G. Zagatti, M. Bortolami, A. Gruppuso, P. Natoli, L. Pagano, and G. Fabbian arXiv:2401.11973.
  • [78] J. T. Sayre et al. Phys. Rev. D 101 (2020), no. 12 122003, arXiv:1910.05748.
  • [79] J. Carron, M. Mirmelstein, and A. Lewis J. Cosmol. Astropart. Phys.  09 (2022) 039, arXiv:2206.07773.
  • [80] M. Kamionkowski Phys. Rev. Lett.  102 (2009) 111302, arXiv:0810.1286.
  • [81] A. P. S. Yadav, M. Shimon, and B. G. Keating Phys. Rev. D 86 (2012), no. 8 083002, arXiv:1207.6640.
  • [82] T. Namikawa Phys. Rev. D 95 (2017), no. 4 043523, arXiv:1612.07855.
  • [83] V. Gluscevic, D. Hanson, M. Kamionkowski, and C. M. Hirata Phys. Rev. D 86 (2012) 103529, arXiv:1206.5546.
  • [84] POLARBEAR Collaboration Phys. Rev. D 92 (2015) 123509, arXiv:1509.02461.
  • [85] BICEP2 and Keck Array Collaborations Phys. Rev. D 96 (2017), no. 10 102003, arXiv:1705.02523.
  • [86] D. Contreras, P. Boubel, and D. Scott J. Cosmol. Astropart. Phys.  1712 (2017), no. 12 046, arXiv:1705.06387.
  • [87] A. Gruppuso, D. Molinari, P. Natoli, and L. Pagano J. Cosmol. Astropart. Phys.  11 (2020) 066, arXiv:2008.10334.
  • [88] M. Bortolami, M. Billi, A. Gruppuso, P. Natoli, and L. Pagano J. Cosmol. Astropart. Phys.  09 (2022) 075, arXiv:2206.01635.
  • [89] BICEP2 and Keck Array Collaborations Astrophys. J.  949 (2023), no. 2 43, arXiv:2210.08038.
  • [90] D. Varshalovich, A. Moskalev, and V. Kersonskii, Quantum Theory of Angular Momentum. World Scientific, 1989.