Vmode polarization of the cosmic microwave background is expected to be vanishingly small in the LambdaCDM model and hence usually ignored. Nonetheless several astrophysical effects as well as beyond standard model physics could produce it at a detectable level. A realistic halfwave plate  an optical element commonly used in CMB experiments to modulate the polarized signal  can provide sensitivity to V modes without significantly spoiling that to linear polarization. We assess this sensitivity for some newgeneration CMB experiments such as the LiteBIRD satellite the groundbased Simons Observatory and a CMBS4like experiment. We forecast the efficiency of these experiments to constrain the phenomenology of certain classes of BSM models inducing mixing of linear polarization states and generation of V modes in the CMB. We find that newgeneration experiments can improve current limits by 1to3 orders of magnitude depending on the data combination. The inclusion of Vmode information dramatically boosts the sensitivity to these BSM models.
['Astrophysics - Cosmology and Nongalactic Astrophysics']
Unveiling V Modes Enhancing CMB Sensitivity to BSM Physics with a NonIdeal HalfWave Plate
{'N. Raffuzzi, a,b M. Lembo, a,b S. Giardiello, c,a M. Gerbino, b M. Lattanzi, b P. Natoli, a,b and L. Pagano a,b,d': "- a Dipartimento di Fisica e Scienze della Terra, Università degli Studi di Ferrara, via Saragat 1, I-44122 Ferrara, Italy\n- b Istituto Nazionale di Fisica Nucleare, Sezione di Ferrara, via Saragat 1, I-44122 Ferrara, Italy\n- c School of Physics and Astronomy, Cardiff University, The Parade, Cardiff, Wales CF24 3AA, United Kingdom\n- d Institut d'Astrophysique Spatiale, CNRS, Univ. Paris-Sud, Université Paris-Saclay, Bât. 121, 91405 Orsay cedex, France  \nE-mail: nicolelia.raffuzzi@unife.it, margherita.lembo@unife.it, GiardielloS@cardiff.ac.uk  \nAbstract. V-mode polarization of the cosmic microwave background is expected to be vanishingly small in the Λ CDM model and, hence, usually ignored. Nonetheless, several astrophysical effects, as well as beyond standard model physics could produce it at a detectable level. A realistic half-wave plate - an optical element commonly used in CMB experiments to modulate the polarized signal - can provide sensitivity to V modes without significantly spoiling that to linear polarization. We assess this sensitivity for some new-generation CMB experiments, such as the LiteBIRD satellite, the ground-based Simons Observatory and a CMB-S4-like experiment. We forecast the efficiency of these experiments to constrain the phenomenology of certain classes of BSM models inducing mixing of linear polarization states and generation of V modes in the CMB. We find that new-generation experiments can improve current limits by 1-to-3 orders of magnitude, depending on the data combination. The inclusion of V-mode information dramatically boosts the sensitivity to these BSM models.", '1 Introduction': 'Observations of anisotropies in the cosmic microwave background (CMB) radiation proved a major observational channel for modern cosmology. The Planck satellite observed them in temperature and polarization with unprecedented precision, providing state-of-the-art constraints on cosmology and fundamental physics [1]. Ground-based experiments have complemented these observations, especially targeting the polarization signal [2-5]. CMB E-mode polarization remains a key focus for next-generation experiments, aiming to achieve a cosmicvariance-limited estimate of the optical depth to reionization τ [6], the least constrained parameter in the Λ CDM model. Also, observations of large-scale B-mode polarization are crucial for improving sensitivity to, or detecting, the tensor-to-scalar ratio r [6-8], ac proxy for primordial gravitational waves. Enhanced sensitivity to small-scale E-mode and B-mode anisotropies will further constrain parameters of Λ CDM and its extensions [7, 8]. Additionally, improved observations of polarized CMB will help explore BSM physics, such as deviations from standard electromagnetism [9], neutrino properties [6], dark sectors [10], and supersymmetries [11, 12].  \nThe circular polarization of the CMB, also known as V modes, is expected to be small in the standard cosmological model, as it is not produced by Thomson scattering during recombination and reionization. Several standard and non-standard physical mechanism can however source a degree of circular polarization in CMB photons, by converting linear polarization generated at the last scattering surface. Among the standard mechanisms, photonphoton scattering via Heisenberg-Euler interaction at recombination produces the strongest circular polarization [13, 14]. For a comprehensive overview of conventional mechanisms that can source V-mode polarization in the CMB radiation see e.g. Ref. [15]. In the realm of non-standard physics, extensions of QED suggest that V-modes arise in the presence of Lorentz-violating operators [9, 16, 17]. Additionally, Refs. [18, 19] show that a pseudo-scalar field or axion inflation leads to V-mode generation. Magneto-optic effects are another class of phenomena capable of generating CMB circular polarization[20, 21]. The detection of circular  \npolarization in the CMB has therefore the potential to offer further evidence for novel physics or more stringent constraints than those available from observations of linear polarization only.  \n̸  \nVast improvements are expected in the observation of the CMB polarized signal both from ground-based experiments, such as the currently running Simons Observatory (SO) [7] and the next-generation CMB experiment CMB-S4 [8], and from satellites such as the LiteBIRD mission [6]. Achieving the promised sensitivity to detect the faint cosmological signals hidden in polarization requires monitoring and minimizing systematic effects. Upcoming observations will benefit from improved scanning strategies, including the use of a rotating half-wave plate (HWP) as a polarization modulator for SO and LiteBIRD. Previous studies have already demonstrated the effectiveness of a rotating HWP in mitigating 1/f noise, [22], and reducing temperature-to-polarization leakage resulting from pair-differencing of orthogonal detectors [23, 24]. However, despite the addition of a HWP is remarkably useful, the inclusion of another optical element in the acquisition chain brings along additional systematic effects, which in turn may compromise the final scientific product. One example is a HWP with a non-ideal phase shift, i.e., a HWP inducing a phase difference = 180 · between the two components of the electromagnetic wave traversing it. Interestingly enough, this deviation, possibly degrading the sensitivity to linear polarization, can also cause a coupling between total intensity and circular polarization. This would allow to investigate the presence of V modes in the CMB radiation [25]. In this paper, we exploit this possibility showing how forthcoming CMB experiments, equipped with HWPs, can be sensitive to circular polarization and, in turn, provide valuable information on a specific class of BSM physical models. We investigate a simple mechanism named the Generalized Faraday Effect (GFE), i.e. the mixing of CMB polarization states, including a partial conversion of linear into circular polarization [26, 27].  \nThe paper is structured as follows: Section 2 introduces the mathematical formalism (Jones calculus [28]) describing the optical behaviour of a (realistic) HWP, the instrumental systematic effect induced by a non-ideal phase shift, as well as the expected sensitivity to V modes; Section 3 illustrates the specific mechanism generating circular polarization and its cosmological phenomenology, while Section 4 describes the analysis method adopted in this work; in Section 5, we present the results of the analysis, forecasting the performance of future CMB experiments. Finally, Section 6 provides conclusions. Appendix A includes a full set of plots with posterior probability distributions of the key cosmological parameters for this work and for each experiment configuration considered in this analysis.', '2 Sensitivity to V modes': 'In this section, we give an overview of the mathematical framework to model a HWP with non-ideal phase shift β (for a detailed description, see, e.g., [29, 30]). We also provide an estimate of the sensitivity to V modes for a generic CMB experiment employing a non-ideal HWP.', '2.1 The non-ideal HWP': 'An ideal half-wave plate (HWP) is an optical device that induces a π -phase shift between the two orthogonal components of the incident wave. In the Jones formalism, the matrix describing the behavior of an ideal HWP is:  \nJ Ideal HWP = ( 1 0 0 -1 ) .  \nA realistic HWP can be described by the more general Jones matrix [29]:  \nJ HWP = ( 1 + h 1 ζ 1 e iχ 1 ζ 2 e iχ 2 -(1 + h 2 ) e iβ ) , (2.1)  \nwhere all but one entries of the matrix (e.g., J 11 ) are complex. The parameters h 1 , 2 are real and negative-defined, and indicate deviations from unitary transmission of the two polarized components, E x,y , of the incident light due to absorption. The parameters ζ 1 , 2 and χ 1 , 2 correspond to the amplitudes and phases of the off-diagonal terms that are responsible for cross-polarization, i.e., mixing between the two orthogonal components of the wave. The parameter β represents the departure from the ideal phase shift of π . In this work, we focus on the effects of a non-vanishing β and how it can be used to measure circular polarization. In the following, we assume for simplicity that the effects due to absorption and cross-polarization are negligible and fix the relative parameters in Eq. 2.1 to zero. The Jones matrix of a HWP with non-ideal phase shift and spinning with angular velocity ω is:  \nJ HWP ( θ ) = R T ( θ ) J HWP R ( θ ) = ( J 11 ( θ ) J 12 ( θ ) J 21 ( θ ) J 22 ( θ ) ) , with R ( θ ) = ( cos θ sin θ -sin θ cos θ ) , (2.2)  \nwhere θ = ωt is the time-dependent angle between the HWP fast axis and the x-axis in the chosen reference frame, and:  \nJ 11 ( θ ) = cos 2 θ -e iβ sin 2 θ J 12 ( θ ) = J 21 ( θ ) = ( 1 + e iβ ) cos θ sin θ J 22 ( θ ) = -e iβ cos 2 θ +sin 2 θ. (2.3)  \nIn the case under study, the complete optical chain traversed by the incoming signal includes a rotating HWP, two orthogonal polarizers, and the detector. Eventually, the Jones matrix of the full optical chain 1 is:  \nJ x,y ( θ ) = J pol , ( x,y ) R T ( θ ) J HWP R ( θ ) , with J pol , x = ( 1 0 0 0 ) , J pol , y = ( 0 0 0 1 ) . (2.4)  \nSo far, we have considered a scenario where the incoming wave is fully polarized and has a quasi-monochromatic frequency. However, the CMB signal is only partially polarized, and a CMB experiment measures the time-averaged incident intensity. Therefore, it is more natural to express the signal in terms of the Stokes vector s = ( T , Q , U , V ) , and employ the Müller formalism [31], which also simplifies tracking the impact on each Stokes parameter.  \nGiven a Jones matrix, it is straightforward to compute the corresponding 4 × 4 Müller matrix. Explicitly, the Müller elements, for the x-oriented polarizer, are:  \nM x TT = 1 2 ( | J 11 | 2 + | J 12 | 2 ) = 1 2 M x TQ = 1 2 ( | J 11 | 2 -| J 12 | 2 ) = 1 2 Cos ( β 2 ) 2 Cos(4 θ ) M x TU = Re[(J 11 J ∗ 12 )] = 1 2 Cos ( β 2 ) 2 Sin(4 θ ) M x TV = Im[(J 11 J ∗ 12 )] = -1 2 Sin( β ) Sin(2 θ ) . (2.5)  \nThe y-oriented elements can be obtained by applying a 90-degree rotation ( 4 θ → 4 θ + π/ 2 ) to the above expressions. Finally, the total power collected by the detector is:  \nd obs = Σ i = x,y ( M i TT T +M i TQ Q +M i TU U +M i TV V ) . (2.6)  \nFrom Eqs. 2.5 - 2.6, it is clear that, when β = 0 , there is no contribution to the observed power from a non-vanishing V-mode signal. Conversely, the departure of the phase shift from the ideal value of π offers a handle to the detection of V modes.', '2.2 Sensitivity to Stokes parameters': "We now quantify the sensitivity to V modes of a CMB experiment equipped with a HWP with non-ideal phase shift. For the sake of simplicity, we consider Gaussian, stationary and uncorrelated instrumental noise of variance σ 2 pix in real space. This simplifies the expression of the noise covariance matrix in pixel space 2 .  \nThe square root of the diagonal elements of the inverse noise covariance in pixel space provides an estimate of the experimental noise per pixel to each Stokes parameter. In the case of a HWP with non-ideal phase shift, the noise - which we normalize to the total intensity T one - is a function of β :  \n( σ Q , U ( β ) σ T ) pix = √ 2 Cos 2 ( β/ 2) ( σ V ( β ) σ T ) pix = √ 2 Sin( β ) . (2.8)  \nIn Figure 1, we show the sensitivity to the Stokes parameters Q , U , V , normalized to the sensitivity to the Stokes T , as a function of the phase shift β . The expected noise ratio for ideal CMB experiments is ( σ Q , U /σ T ) pix = √ 2 . The non-ideal phase shift leads to a mixing of power among the Q , U , V Stokes parameters. This can be understood from the fact that the sensitivity to Q and U decreases (i.e., the noise ratio increases, in red in Fig. 1) while that of V (blue) increases for larger values of β . However, while the sensitivity to V rapidly changes with β , the degradation of the sensitivity to Q , U significantly deviates from the ideal case only for fairly large values of β ≳ 20 · (for comparison, typical values for β vary from a few degrees up to a few tenths [32-34]).  \nIt is also instructive to work out the sensitivity to V modes in harmonic space for a given experimental setup, as described in terms of the linear polarization noise level σ Q , angular resolution θ FWHM , observed sky fraction f sky and non-ideality HWP parameter β .  \nˆ m pix = ( A T N -1 A ) -1 pix ( A T N -1 d ) pix , (2.7)  \nwhere N = ⟨ nn T ⟩ is the noise covariance matrix in time domain and ( A T N -1 A ) -1 pix is the noise covariance matrix in pixel space. For the simple noise properties considered in this work, i.e. N = diag ( σ 2 ) , the noise covariance matrix in a pixel is simply given by σ 2 pix = σ 2 ( A T A ) -1 pix . This becomes σ 2 pix = 4 σ 2 N hit diag ( 1 , 2 / cos( β/ 2) 4 , 2 / cos( β/ 2) 4 , 2 / sin( β ) 2 ) , where N hit is the number of hits on that pixel, since we assume a very homogeneous and redundant scanning strategy. See, e.g. [30] and references therein for further details.  \nFigure 1 : Noise ratio of the Stokes parameters for linear polarization σ Q , U , pix ( β ) (red) and circular polarization σ V , pix ( β ) (blue), normalized to the total intensity noise, as a function of the HWP non-ideal phase shift β . Typical values for β range from a few degrees to a few tenths [32-34]. For ideal CMB experiments, ( σ Q,U /σ T ) pix = √ 2 for polarization-sensitive detectors. In presence of a non-ideal phase shift, the noise level for circular polarization σ V , pix ( β ) decreases rapidly with increasing β , while the noise level for linear polarization σ Q , U , pix ( β ) does not increase significantly until very large values of β are reached. As a result, a slightly non-ideal HWP (small value of β ) can help gain sensitivity to V modes without significantly degrading the sensitivity to linear polarization.  \n<!-- image -->  \nNote that in the rest of this section we will denote with σ P the nominal linear polarization noise, corresponding to an ideal HWP plate, i.e. σ Q = σ Q ( β = 0) = √ 2 σ T . In this way factors of sin β and cos β will appear explicitly in the equations.  \nThe observed power spectrum of the signal V modes ˆ C VV ℓ has an associated variance Σ 2 ℓ :  \nΣ 2 ℓ = 2 (2 ℓ +1) f sky [ C VV , fid ℓ + N ℓ B 2 ℓ 1 Sin 2 ( β ) ] , (2.9)  \nwhere C VV , fid ℓ is the true underlying power spectrum of the signal, N ℓ is the noise power spectrum of the linear polarization, assuming an ideal half-wave plate, and B ℓ is the harmonic equivalent of a Gaussian beam of width θ FWHM . We assume a white noise power spectrum N ℓ = σ 2 Q , and use a fiducial Gaussian approximation for the likelihood L ( C VV ℓ ) = Pr( ˆ C VV ℓ | C VV ℓ ) of a theoretical power spectrum C VV ℓ :  \n-2 ln L ( C VV ℓ ) = ( C VV ℓ -ˆ C VV ℓ ) 2 Σ 2 ℓ , (2.10)  \nFor the purpose of our analysis, we take ˆ C VV ℓ = ⟨ ˆ C VV ℓ ⟩ = C VV , fid ℓ . The inverse probability Pr( C VV ℓ | ˆ C VV ℓ ) is computed through the use of Bayes' theorem with a flat prior on Cℓ .  \nArmed with the above, we compute, for each multipole ℓ , the 95% upper limit on C VV ℓ in the assumption of vanishing V-mode signal, i.e., for C VV , fid ℓ = 0 , and for different experimental configurations. The results are shown in Fig. 2. Here, we show the results for a LiteBIRD-like satellite (black solid curve) and for a ground-based Simons Observatory/SO SAT 3 experiment (magenta dashed line), assuming β = 10 · in both cases. The other parameters describing each configuration can be read in Tabs. 1 (noise and angular resolution) and 2 (sky coverage). For comparison, the current 95% upper bounds from the balloon-borne experiment SPIDER [25] (blue and red) and the ground-based CLASS [35] (yellow) are also shown. A LiteBIRDlike experiment would improve current bounds by several orders of magnitude and extend the sensitivity over a broader range of angular scales. On the other hand, a ground-based experiment like SO will dominate at smaller scales.  \nIn Figure 3, we show iso-contours of constant integrated 95% sensitivity to a scaleinvariant VV spectrum, as a function of β and the instrumental noise level in linear polarization, σ Q . In the left (right) panel, the angular resolution and sky fraction are fixed to the values corresponding the LiteBIRD (SO) experiment. The two panels also correspond to different choices of the observed multipole range, again representatives of the aforementioned experiments: 2 ≤ ℓ ≤ 300 in the left panel and 50 ≤ ℓ ≤ 300 in the right panel. As expected, the sensitivity increases for smaller values of σ Q and large values of β . The stars in the plot indicate the expected performance of a LiteBIRD-like (left panel) and SO (right panel) experiment for β = 10 · .", '3 Production of V modes': 'In the previous section, we have quantified the sensitivity to V modes assuming a generic scaleinvariant spectrum. In the following, we work out the sensitivity of different combinations of CMB experiments to a V-mode signal generated in the phenomenological framework studied in [27]. We now briefly recall the basis of this formalism.  \nWe consider a simple class of models allowing for the in vacuo mixing of CMB polarization states during propagation, including the conversion of linear into circular polarization. The mixing of polarization states in a weakly anisotropic medium is dubbed Generalized Faraday Effect (see sec. 14.3 of [36] for a description of the general framework, and [26] for an application to extragalactic synchrotron sources), or GFE. As the name suggests, the usual Faraday rotation is a particular case of GFE. The radiative transfer equation for the linearly polarized radiation propagating in a weakly anisotropic non-absorbing medium reads [36]  \nd d s S = ρ × S (3.1)  \nwhere s is an affine parameter following the direction of wave propagation, and S = ( Q,U,V ) is the polarization Stokes vector, precessing about the direction of ρ = ( ρ Q , ρ U , ρ V ) , which is determined by the dielectric tensor of the medium and thus encodes its optical properties. If ˆ ρ ≡ ρ / | ρ | is directed along the V -axis in the polarization space ( ρ Q = ρ U = 0 ), then the Stokes vector S precesses about it. The linear polarization component, i.e. the projection of S on the ( Q, U ) plane, will thus rotate while mantaining a constant magnitude. This is the usual  \nV modes  -by-  sensitivity  \nFigure 2 : The figure illustrates the ℓ -byℓ (i.e. no binning) 95% C.L. sensitivity for current and future CMB experiments to a VV-spectrum. Assuming a non-ideal HWP with a phase shift of β = 10 · , the black curve represents the sensitivity for a LiteBIRD-like experiment and the dark magenta curve corresponds to the sensitivity for a SO SAT experiment. Details on the full experimental setup adopted to obtain these curves are in Tab. 1. For comparison, the current bounds on V modes at different frequencies and ℓ -ranges from the balloon-borne SPIDER [25] (blue and red) and the ground-based CLASS [35] (yellow) are also shown.  \n<!-- image -->  \nFaraday rotation, where no circular polarization is generated, and only the Q and U states mix among themselves. On the other hand, if ρ has non-vanishing Q and U components, a partial conversion of linear into circular polarization occurs, generating non-zero V modes.  \nIn Ref. [27], some of us have developed a formalism that allows to derive the modifications to the CMB spectra in temperature, linear and circular polarization, and their cross-correlations, given an effective dielectric tensor that describes the optical properties of the space traversed by CMB photons. This allows to be agnostic about the specific model beyond the mixing among CMB polarization components and just capture the basic phenomenology. Therefore, once the phenomenological parameters are constrained, it is possible to interpret the results in light of a specific physical model that could have induced GFE; see, e.g. [9] for a recent application.  \nHere we consider the case of an effective dieletric tensor that is homogeneous and independent of both the direction and magnitude of the radiation wave vector.  \nIn this specific case, parity-violating CMB cross spectra are vanishing, TB = EB = EV = BV = 0 (see Ref. [27] for details).  \nThe remaining CMB power spectra C XY ℓ (with X, Y = T, E, B, V ) can be expressed  \n4  \n4  \nFigure 3 : Contours of constant V-mode 95% sensitivity for a scale-invariant VV spectrum as function of non-ideal phase-shift β and nominal (i.e., assuming β = 0 ) linear polarization sensitivity σ Q . Left panel: the sensitivity is integrated over the multipole range 2 ≤ ℓ ≤ 300 for an angular resolution of 30 arcmin, similar to that achievable by LiteBIRD from space. For the LiteBIRD experiment, the expected noise level is σ Q ≃ 2 . 2 µ Karcmin , integrated over 22 frequency channels and 3-year observation time [6]. The star marks the combination of this noise level with a value of the nonideality parameter β = 10 · (chosen arbitrarily). Right panel: the sensitivity is integrated over 50 ≤ ℓ ≤ 300 for an angular resolution of 17 arcmin, similar to that achievable by the ground-based SO-SAT experiment. For the SO SAT, the noise level is σ Q ≃ 2 µ Karcmin , combining the 93 and 145 GHz bands [7]. The star marks the combination of this noise level with a value of the nonideality parameter β = 10 · (chosen arbitrarily).  \n<!-- image -->  \nin terms of the CMB spectra produced at the last-scattering surface ˜ C XY ℓ :  \nC TE ℓ = ( 1 -β 2 V + β 2 E 8 π ) ˜ C TE ℓ C EE ℓ = ( 1 -β 2 V + β 2 E 4 π ) ˜ C EE ℓ + β 2 V 4 π [ W (1) ℓ ˜ C EE ℓ + W (1) ℓ +1 ˜ C BB ℓ +1 + W (1) ℓ -1 ˜ C BB ℓ -1 ] C BB ℓ = ( 1 -β 2 V + β 2 E 4 π ) ˜ C BB ℓ + β 2 V 4 π [ W (1) ℓ ˜ C BB ℓ + W (1) ℓ +1 ˜ C EE ℓ +1 + W (1) ℓ -1 ˜ C EE ℓ -1 ] C VV ℓ = β 2 E π [ W (2) ℓ +2 ˜ C BB ℓ +2 + W (2) ℓ +1 ˜ C EE ℓ +1 + W (2) ℓ ˜ C BB ℓ + W (2) ℓ -1 ˜ C EE ℓ -1 + W (2) ℓ -2 ˜ C BB ℓ -2 ] , (3.2)  \nwhere W (1 , 2) are combinations of Wigner 3 j -symbols, and the phenomenological parameters β 2 E and β 2 V are, respectively, combinations of the diagonal and off-diagonal components of the effective susceptibility tensor [27].  \nFrom Eqs. 3.2 and Fig. 4, it is clear that the modification of the CMB polarization signal due to this particular realization of GFE has three main effects:  \n)  \n1  \n+  \n(  \n- · rescaling of the amplitude of TE , EE , BB spectra via the combination β 2 V + β 2 E ;\n- · mixing between E and B modes (cosmic birefringence), proportional to β 2 V ;\n- · sourcing of V modes from E and B modes, proportional to β 2 E .  \n1.5  \nFigure 4 : Impact of Generalized Faraday Effect (GFE) on polarization power spectra. The reference (blue) power spectra are computed using Planck best-fit values from TT, TE, EE, lowE + lensing [1] and vanishing GFE parameters, as reported in Tab. 3, while the modified spectra (orange) are calculated fixing β 2 V = 0 . 4 and β 2 E = 1 . 8 (chosen arbitrarily to make the GFE visible).  \n<!-- image -->', '4 Analisys method and datasets': "In this section, we present methodology and data employed in the likelihood analysis. We simulate the performance of different upcoming CMB experiments (satellite and ground-based) and forecast their sensitivity to the phenomenological GFE parameters β 2 E,V introduced in the previous section. To this scope, we perform a Monte Carlo Markov Chain (MCMC) analysis to obtain forecasted constraints on β 2 E and β 2 V jointly with the six Λ CDM parameters  \nand the tensor-to-scalar ratio r 4 . We use a customized version 5 of the CosmoMC [37] package to perform the analysis. We assume convergence of the MCMC chain by requiring that the Gelman-Rubin statistics R -1 ∼ 0 . 01 .  \nWe run forecasts for five combinations of data sets which simulate the performance of different CMB experiments:  \n- · a future LiteBIRD-like satellite experiment. We refer to this case as ' L ';\n- · LiteBIRD-like in combination with a simulated version of the Planck data, labeled ' LP ';\n- · a currently running Simons Observatory (SO) ground-based experiment, in combination with Planck, labeled ' PSO '. We include simulated performance of the observations with both the Large Aperture Telescope (LAT) which targets the smallest angular scales and the array of Small Aperture Telescopes (SATs) which target B modes at degree angular scales. The HWP is deployed on SATs only;\n- · Planck in combination with a next-generation ground-based experiment, i.e. CMBS4-like, labeled ' PS4 '. As for SO, we include information from both LAT and SATs. We note that the baseline design of CMB-S4, as described in [8, 38], does not include a HWP. Therefore, when referring to a CMB-S4-like experiment, we actually refer to the combination of expected sensitivity as determined by the sky coverage, angular resolution and noise properties reported in [8, 38], augmented with the use of a HWP in the telescopes targeting B modes to modulate the polarization signal;\n- · LiteBIRD-like in combination with CMB-S4-like, labeled ' LS4 '.  \nThe instrument characteristics of each simulated experiment are reported in Tab. 1. In the following, when referring to individual experiments, we will omit the 'like' suffix for brevity.  \nAs anticipated, for the sake of simplicity, we do not make use of the real data collected by Planck. We instead consider a simple white noise spectrum with 40 µ K-arcmin in temperature, and 80 µ K-arcmin in polarization [41], which reasonably reproduces the Planck performance. The noise spectrum for LiteBIRD is computed given the detector sensitivity and the angular resolution of LiteBIRD in the central frequency channels [6]. For SO [7], we make use of the publicly available noise curves 6 . The latter also include residual contributions from component separation, which dominates the cosmological signal both at the very small scales probed by the SO LAT, as well as at the intermediate scales probed by the SO SAT. CMB-S4 noise curves were computed referring to the central channels for both SAT and LAT [8].  \nWhen adding V-mode data in the analysis, the noise spectra of LiteBIRD, SO and CMB-S4 are further rescaled to account for the non-ideal phase shift, as detailed in Sec. 2.  \nTable 1 : Instrumental specifications for each experiment: sensitivity to total intensity in µ Karcmin (linear polarization sensitivity is obtained by multiplying this value by a √ 2 factor, except for Planck which is 2 as only half the detectors are polarized); beam size FWHM in arcmin; non-ideal phase shift in degrees (for those experiments that will deploy - or we simulate will deploy - a HWP [7, 38, 39]). The noise levels for SO SAT and LAT are reported for reference. They are coadded over the central 90 and 150 GHz channels, and the FWHM corresponds to the 145 GHz channel. For more details, see [7, 40]. In this work, we make use of the official SO noise curves.  \nIn practice, the rescaled noise curves become:  \nN EE ℓ B 2 ℓ → 1 B 2 ℓ N EE ℓ Cos 4 ( β/ 2) N BB ℓ B 2 ℓ → 1 B 2 ℓ N BB ℓ Cos 4 ( β/ 2) N VV ℓ B 2 ℓ → 1 B 2 ℓ N VV ℓ Sin 2 ( β ) . (4.1)  \nNote that Planck/CMB-S4 and LiteBIRD experiments observe largely overlapping sky fractions. While LiteBIRD dominates over Planck in B-mode sensitivity, the performance of the two experiments are different in T and E modes at different angular scales, depending on the respective noise level and angular resolution. A similar consideration applies to the combination ' LS4 ' as far as the T and E-fields ('TE' hereafter) are concerned. Therefore, for the 'TE' fields in the ' LP ' and ' LS4 ' cases, we use an effective inverse noise weighted curve:  \nN eff ℓ = [ ( N ℓ B 2 ℓ ) -1 X + ( N ℓ B 2 ℓ ) -1 LiteBIRD ] -1 , (4.2)  \nwhere X = Planck , CMB-S4. We can distinguish two regions where the two experiments perform differently: LiteBIRD dominates at larger scales, whereas Planck or CMB-S4 outperform LiteBIRD at smaller angular scales.  \nDifferently, when Planck is combined with SO/CMB-S4, information from the former is only included in the sky fraction and/or at the angular scales not observed by SO/CMB-S4(LAT). As a result, in the ' PSO / PS4 ' case, we make use of the individual noise curves for each experiment.  \nWe employ an exact likelihood in harmonic space, rescaled by a ℓ -independent sky fraction f sky to account for the experiment-specific sky coverage.  \n-2 ln L ℓ = f sky (2 ℓ +1) ( Tr [ ˆ S ℓ S -1 ℓ ] +lndet[ S ℓ ] ) + const. (4.3)  \nwhere ˆ S and S are matrices of the 'observed' (i.e., simulated, in our case) ˆ C ℓ and theoretical C ℓ power spectra, respectively (see, e.g., Sec. 2.2 of [42] for the definition of ˆ S and S ). When V modes are included in the likelihood analysis, the standard ˆ S and S matrices are extended to include the C VV ℓ . The sky fraction, range of angular scales and observed fields for each data combination 7 are reported in Tab. 2.  \nWe simulate the observed spectra ˆ C ℓ as the sum of a reference set of CMB spectra obtained with the Boltzmann code camb [43] from a fiducial set of cosmological parameters and the noise power spectrum in harmonic space: ˆ C ℓ = C fid ℓ +N ℓ / B 2 ℓ (see, e.g., Sec. 3.3 of [42] for the extension of the definition of ˆ S and S in Eq. 4.3 in presence of noise and beam smearing). We assume the fiducial model to be Λ CDM, i.e., vanishing V-mode signal and vanishing primordial tensor modes. The cosmological parameters corresponding to the fiducial model are reported in Tab. 3. A major source of variance in the observation of primordial B modes is the BB lensing contribution, which can be significantly reduced via appropriate de-lensing procedures. In order to quantify the effect of delensing in our results, we consider two cases: a fully lensed fiducial BB, i.e., no de-lensing, and a partly de-lensed fiducial BB. The two cases are obtained by properly rescaling the BB lensing contribution (in both the fiducial and the theory BB spectra) by a factor A L = 1 (fully lensed) and A L = 0 . 3 (70% delensing). The GFE-modified (tensorial) B modes are later summed to the A L -rescaled lensed B modes. The underlying assumption, also used in Ref. [27], posits that GFE exclusively affects primordial B modes, even though GFE and lensing should occur simultaneously as integrated effects along the line of sight. This assumption is valid for LiteBIRD and SO, see, e.g., [27]. However, for fourth-generation ground-based experiments such as CMB-S4, caution is advised regarding the interplay between lensing and GFE. In order to assess the validity of the assumption for CMB-S4, we compute the difference between the BB spectrum obtained as detailed above and the BB spectrum obtained by applying GFE to the sum of tensor and lensing BB spectra. We then compare this difference with the CMB-S4 noise curve. We find that, for reasonably small values of GFE parameters ( β 2 E,V ≲ 0 . 1 ), the difference in the predicted CMB spectra between the two approaches is well below the noise contribution. Furthermore, other studies in the literature, such as Ref. [44], have confirmed that lensed B modes have a negligible impact on cosmic birefringence.  \nTo better quantify the impact of a potential sensitivity to V modes, we compare results obtained with and without the inclusion of V modes in the datasets for each combination of experiments. We refer to the two cases as 'TEBV' and 'TEB', respectively. In the 'TEBV' case, we assume the non-ideal phase shift to be β = 10 · ; as already noted, this choice is in agreement with the expected values of the phase shift from simulations. In addition, as shown in Fig. 1, this value allows concurrently for a relatively acceptable noise level in V -Stokes parameter and negligible degradation of the noise level in Q / U -Stokes. In the 'TEB' case, we assume the phase shift to be β = 0 .  \nTable 2 : Data combinations considered in this work. For each experiment, we report the observed fields (temperature T, linear polarization E and B, circular polarization V) and the corresponding range of angular scales (multipoles), as well as the sky fraction observed ( f sky ). Throughout the main text, we refer to each combination with the corresponding label in the first column.  \nTable 3 : Fiducial values of cosmological parameters for the reference cosmological model assumed to simulate the observed CMB spectra in this analysis. They are taken from the Planck best-fit values from TT, TE, EE + lowE + lensing analysis [1].", '5 Results': "In this section, we present and discuss the results obtained for the various datasets and cases introduced in the previous section. These results are also compared with the current bounds  \nreported in literature [27, 45]. We focus on a subset of cosmological parameters: the optical depth ( τ ), the tensor-to-scalar ratio ( r ), and the GFE parameters ( β 2 E and β 2 V ) 8 . The first two parameters are linked to the primary science goals of experiments observing large and intermediate-scale polarization, such as LiteBIRD, SO, and CMB-S4, while the latter two are the main focus of this study. The 68% C.L. constraints on τ and the 95% C.L. upper bounds on the remaining three parameters are in Tab. 4 for both the fully lensed and partially delensed scenarios. Marginalized 1D and 2D posterior probabilities for r , τ , β 2 E and β 2 V can be found in the Figures in Appendix A. We compare the results obtained with and without the inclusion of V-mode data. In doing so, we expect to place more stringent limits on β 2 E and β 2 V . In addition, we test whether the partial degradation of sensitivity to linear polarization resulting from a non-ideal phase shift affects the constraints on τ and r .  \nFrom Tab. 4, we can draw the following considerations:  \n- · as known, the constraints on τ and r are driven by the ability of a given data combination to recover the large-scale E-mode signal and the large-to-intermediate B-mode signal respectively. In Tab. 4, compare, e.g., the improvement of the constraints on τ from the combinations including LiteBIRD with respect to those including Planck. Similarly, compare the lack of improvement of the constraints on r from the ' L ' dataset with respect to the ' LP ' combination.  \nIgnoring for a moment the impact of V-mode observations (but see the second bullet point below), the constraints on β 2 E are mostly driven by the ability to reconstruct the TE and EE spectra at intermediate and small scales, since β 2 E acts as an overall rescaling factor for the CMB spectra (see Eqs. 3.2). Indeed, by comparing the first rows of each Dataset box in Tab. 4, we note a steady improvement of the constraints on β 2 E when considering dataset which are more and more sensitive to small-scale E (compare, e.g., ' L ' and ' LP '; also, note the lack of improvement when comparing ' PS4 ' with ' LS4 ', as CMB-S4 dominates the sensitivity to E in both data combinations).  \nThe constraints on β 2 V are mostly driven by the sensitivity of each dataset to B modes, both via the overall rescaling of the amplitude (first term on the RHS in Eq. 3.2) and via the mixing with E modes (second term on the RHS) 9 .  \n- · the addition of V-mode data tightens the constraints on β 2 E by roughly an order of magnitude but has no effect on β 2 V (compare the two rows corresponding to each data label in Tab. 4; see also Appendix A). Indeed, C VV ℓ in Eq. 3.2 is solely affected by β 2 E that acts as an overall rescaling factor of a mixing of ˜ C EE ℓ and ˜ C BB ℓ . The inclusion of V modes, with the corresponding rescaling of the noise spectra as per Eq. 4.1 and Fig. 1, does not degrade the constraints on r and τ . This is one of the main findings of this work. This stands in contrast to the assertion made in [27] regarding the necessity of improving linear polarization for increased accuracy on GFE parameters. The inclusion of V-mode data proves to be an additional, significant source of constraining power,  \nbeyond what is provided by linear polarization alone. The sensitivity (i.e. the low VV noise) is enhanced by the non-ideal phase shift of the HWP (as shown in Fig. 1), where a relatively conservative value of β = 10 · was employed;  \n- · by reducing the contribution of signal variance in the observation of primordial B modes, the de-lensing procedure results in a factor of two improvement of the constraints on r and β 2 V (compare rows with the same data label in the top and bottom panels in Tab. 4). No improvement is seen for τ and β 2 E , as could be visually verified by comparing the left (no de-lensing) with the right (70% de-lensing) panel of triangular plots in Appendix A. This confirms that the constraining power on these parameters comes essentially from TE and EE spectra rather than from BB (in addition to VV for β 2 E ).\n- Improvements on r and β 2 V are expected, since they directly impact the C BB ℓ power spectrum as rescaling factors. In Eq. 3.2, we recognize an attenuating factor ( β 2 V + β 2 E ) / 4 π in front of the unmodified ˜ C BB ℓ . In the second term on the RHS, the amplitude of the combination of unmodified EE and BB spectra - with the ˜ C BB ℓ spectrum being subdominant compared to the unmodified ˜ C EE ℓ - is regulated by β 2 V . This requires β 2 V to suppress the power coming from the EE spectrum, thus restoring the expected amplitude of the observed tensorial B modes ( r = 0 ).\n- Also, de-lensing does not help improve the constraining power on the optical depth τ . In principle, information on τ can be also extracted from measurements of the reionization peak in the BB spectrum. However, the strong degeneracy between r and τ in that region as well as the residual lensing variance dominating also at lowℓ s for vanishingly small values of r severely limit the information content that can be extracted from BB to improve the constraints on τ . As a result, the constraining power for this parameter comes essentially from E modes only;\n- · when comparing the constraints on r from ' PSO ' and ' PS4 ' in the fully lensed case (top panel in Tab. 4), we note that the latter dataset provides - somehow unexpectedly - a looser bound on r than the former dataset (see also left panels of Figs. 7 and 8). This feature can be attributed to the fact that the two experiments driving the constraints on r in the two data combinations - namely, SO and CMB-S4 - are both cosmic variance limited in BB in the fully lensed scenario. With the instrumental noise being subdominant, the ability to jointly constrain the parameters affecting C BB ℓ is solely determined by volume effects in the parameter space explored in the MCMC analysis. More in detail, at first order, we can assume that the amplitude of C BB ℓ is determined by a combination of r , β 2 E and β 2 V (see Eqs. 3.2). When switching from ' PSO ' to ' PS4 ' data combination, the enhanced sensitivity of the latter to E modes leads to more stringent limits on β 2 E , thus allowing for a larger excursion of r in parameter space. For the same reason, the constraints on β 2 V from the two data combinations are virtually unchanged. In this case, the volume effect responsible for the degradation of the constraints on r is suppressed for β 2 V by the extra contribution of the E-B mixing term on the RHS in the equation for C BB ℓ .  \nConversely, for partial de-lensing (bottom panel of Tab. 4; see also right panels of Figs. 7 and 8), we observe the expected improvement of the constraint on r when moving from ' PSO ' to ' PS4 '. While ' PS4 ' remains cosmic variance limited, the impact of instrumental noise in ' PSO ' cannot be neglected and washes out the volume effect described above.  \nFinally, in Tab. 5 we quote again the 95% C.L. upper limits for the five analyzed cases combining TEBV modes, and compare them with current bounds. We note that current limits in the literature are obtained without a de-lensing procedure applied to data. Nevertheless, we report forecasts for both the fully lensed and the partially de-lensed scenarios, to also assess the benefit of de-lensing. The current upper limit on r is taken from the analysis of BICEP2, Keck Array, BICEP3, WMAP, and Planck data, marginalizing over a given model of Galactic dust and synchrotron contaminations [45]. As already mentioned, no de-lensing is applied to the data. Future experiments targeting r will see considerable improvement in the constraints by applying suitable de-lensing techniques (see, e.g., [8, 46, 47]). Our forecasts qualitatively confirm this expectation. Present bounds on β 2 E,V are taken from the analysis of [27] and were derived using temperature and linear polarization anisotropies from the Planck legacy release [48] and BICEP2/Keck 2015 [49] . The inclusion of current V-mode observations does not help improve the constraints on the GFE parameters [27]. Forecasted sensitivity to β 2 V is dominated by the sensitivity of upcoming experiments to B modes, with de-lensing playing a crucial role. As far as β 2 E is concerned, while de-lensing has minimal impact on the forecasted sensitivity, the inclusion of V modes from future experiments significantly improves the constraints. Our forecast shows that more stringent constraints ranging from 1 to 3 orders of magnitude can be achieved, depending on the considered parameter and dataset.", '6 Conclusions': "̸  \nIn this study, we investigated how the deployment of a realistic half-wave plate (HWP) inducing a non-ideal (i.e., = π ) phase shift, may affect the noise characteristics of CMB experiments, while providing sensitivity to the presence of a non-vanishing V-mode signal. We then forecasted the expected performance of new-generation CMB experiments, noting that our analysis did not consider foregrounds.  \nDeviations from the ideal phase shift modify the noise properties of a CMB experiment deploying a HWP by degrading the sensitivity to Q / U polarization. This could potentially spoil the ability to reach the primary science goals of future CMB experiments, such as a CVL determination of the optical depth τ and a tight constraint on/high-significance detection of the tensor-to-scalar ratio r . However, at the same time, deviations from the ideal performance of a HWP allow for the experiment to become sensitive to possibly non-vanishing V modes in the observed signal. This work showed how even a significant deviation from the ideal phase shift (of the order of 10 · ) can open a window to the investigation of beyond-the-standardmodel (BSM) cosmological scenarios while having no impact on the sensitivity to τ and r .  \nAs a test case, we focused on the phenomenology of a class of models leading to BSM electromagnetic effects on the CMB, the Generalized Faraday Effect (GFE) [27]. The GFE is responsible for E-B mixing as well as for sourcing circular polarization from the conversion of linear CMB polarization. These effects are modeled via two phenomenological parameters β 2 E and β 2 V - which modify the shape of the CMB spectra as they emerge from the last scattering surface (Eqs. 3.2). We forecasted the sensitivity to these parameters from currently running and future CMB experiments which will or might deploy a HWP. At the same time, we checked that the sensitivity to τ and r expected from these experiments is not degraded. We considered different combinations of simulated data from a LiteBIRD-like satellite experiment ( L ), from the ground-based Simons Observatory ( SO ) and a modified version (i.e., with the  \nTable 4 : Constraints on a subset of cosmological parameters sampled in this analysis: 68% C.L. constraints for the optical depth ( τ ), and 95% C.L. upper limits for the tensor-to-scalar ratio ( r ) and the GFE parameters ( β 2 E , β 2 V ). The results are shown assuming a lensing factor A L = 1 (no de-lensing, top table) and A L = 0 . 3 (70% de-lensing efficiency, bottom table), for different data combinations ( L , LP , PSO , PS4 , and LS4 ). For each data combination, we compare the results obtained with and without the inclusion of V-mode observations.  \ninclusion of a HWP) of a ground-based CMB-S4-like experiment ( S4 ) 10 . We also included information from Planck ( P ) observations where relevant.  \nBy comparing the results with current limits in the literature [27], we forecasted a 1-to-3 order-of-magnitude improvement of the constraints on β 2 V and β 2 E , depending on the dataset,  \nTable 5 : 95% C.L upper limits on the tensor-to-scalar ratio and GFE parameters ( β 2 E and β 2 V ) for the different data combinations considered in this analysis. These limits are obtained from the analysis of T, E, B and V modes, with a lensing factor A L = 1 (second, third and fourth columns) and A L = 0 . 3 (last three columns). In the last row, bounds on the same parameters from current data [27, 45] are also reported to ease the comparison, where no de-lensing procedures are included.  \nwith the most stringent limits achieved with the ' LS4 ' combination (see Tab. 5).  \nThe possibility to include observations of V modes, as enabled by the deployment of a realistic HWP, is key to obtain dramatic improvements in the constraints on β 2 E , ranging from a factor of a few for the combinations ' PSO ' and ' PS4 ' to an order of magnitude and more for the data combinations including LiteBIRD (see Tab. 4). Better sensitivity to intermediateto-small scale E modes is also an important driver of the improved constraints on β 2 E , as can be noted by comparing ' L ', ' LP ', ' PSO '/' PS4 '. Constraints on β 2 V are mostly insensitive to the inclusion of V modes in the analysis. Instead, they are driven by improved sensitivity to B modes. The constraints improve roughly by a factor of two when partial de-lensing (obtained as a 70% suppression of lensing power in the simulated BB data) is applied in the analysis (compare the top and bottom panels in Tab. 4). The stringent constraints on β 2 V and β 2 E will provide valuable insights on the fundamental physics governing the evolution of the Universe, particularly from the CMB last scattering to the present epoch. The GFE parameters can be mapped onto the elements of an effective 'cosmic susceptibility tensor', which describes the optical properties of the Universe seen as the medium through which CMB photons propagate. The mapping between β 2 V and β 2 E and the elements of the susceptibility tensor can be obtained once a specific physical model responsible for the GFE is specified. As such, constraints on β 2 V and β 2 E shed light on the nature of the GFE and its potential origins, such as magnetic fields [20, 21, 36, 50-52] or other BSM phenomena [9, 13, 14, 19, 53] that could have influenced the polarization state of the CMB. The methodology outlined in this work can be also applied to models capable of generating parity-breaking spectra (e.g. Ref. [9]). In conclusion, the results presented in this study indicate that the possibility of including observations of V modes in the cosmological analysis is a key tool to investigate more deeply specific aspects of fundamental physics which are otherwise less accessible, without spoiling the constraining power to other primary science targets. This emphasizes the necessity of fully exploit the potential of upcoming CMB experiments, including the extraction of V-mode data, which will provide unprecedented insights on the fundamental physics of the early Universe.", 'Acknowledgments': 'Wethank Alessandro Gruppuso for useful discussions during the preparation of the manuscript and for feedback on the final version of the paper. We acknowledge the financial support from the INFN InDark initiative and from the COSMOS network (www.cosmosnet.it) through the ASI (Italian Space Agency) Grants 2016-24-H.0 and 2016-24-H.1-2018, as well as 2020-9-HH.0 (participation in LiteBIRD, phase A). We acknowledge the use of camb [43], CosmoMC [37], GetDist [54], and the use of computing facilities at CINECA. MG is funded by the European Union (ERC, RELiCS, project number 101116027) and by the PRIN (Progetti di ricerca di Rilevante Interesse Nazionale) number 2022WJ9J33. SG acknowledges support from the Horizon 2020 ERC Starting Grant (Grant agreement No 849169) and from STFC and UKRI (grant numbers ST/W002892/1 and ST/X006360/1). This is not an official SO Collaboration paper.', 'References': "- [1] Planck collaboration, Planck 2018 results. VI. Cosmological parameters , Astron. Astrophys. 641 (2020) A6 [ 1807.06209 ].\n- [2] POLARBEAR collaboration, A Measurement of the Degree Scale CMB B -mode Angular Power Spectrum with POLARBEAR , Astrophys. J. 897 (2020) 55 [ 1910.02608 ].\n- [3] POLARBEAR collaboration, A Measurement of the CMB E -mode Angular Power Spectrum at Subdegree Scales from670 Square Degrees of POLARBEAR Data , Astrophys. J. 904 (2020) 65 [ 2005.06168 ].\n- [4] ACT collaboration, The Atacama Cosmology Telescope: DR4 Maps and Cosmological Parameters , JCAP 12 (2020) 047 [ 2007.07288 ].\n- [5] SPT-3G collaboration, Measurements of the E-mode polarization and temperature-E-mode correlation of the CMB from SPT-3G 2018 data , Phys. Rev. D 104 (2021) 022003 [ 2101.01684 ].\n- [6] LiteBIRD collaboration, Probing Cosmic Inflation with the LiteBIRD Cosmic Microwave Background Polarization Survey , PTEP 2023 (2023) 042F01 [ 2202.02773 ].\n- [7] Simons Observatory collaboration, The Simons Observatory: Science goals and forecasts , JCAP 02 (2019) 056 [ 1808.07445 ].\n- [8] CMB-S4 collaboration, CMB-S4: Forecasting Constraints on Primordial Gravitational Waves , Astrophys. J. 926 (2022) 54 [ 2008.12619 ].\n- [9] L. Caloni, S. Giardiello, M. Lembo, M. Gerbino, G. Gubitosi, M. Lattanzi et al., Probing Lorentz-violating electrodynamics with CMB polarization , JCAP 03 (2023) 018 [ 2212.04867 ].\n- [10] G. Choi, C.-T. Chiang and M. LoVerde, Probing Decoupling in Dark Sectors with the Cosmic Microwave Background , JCAP 06 (2018) 044 [ 1804.10180 ].\n- [11] I. Dalianis and Y. Watanabe, Probing the BSM physics with CMB precision cosmology: an application to supersymmetry , JHEP 02 (2018) 118 [ 1801.05736 ].\n- [12] PRISM collaboration, PRISM (Polarized Radiation Imaging and Spectroscopy Mission): An Extended White Paper , JCAP 02 (2014) 006 [ 1310.1554 ].\n- [13] I. Motie and S.-S. Xue, Euler-Heisenberg Lagrangian and Photon Circular Polarization , EPL 100 (2012) 17006 [ 1104.3555 ].\n- [14] R.F. Sawyer, Photon-photon interactions can be a source of CMB circular polarization , 1408.5434 .  \n- [15] P. Montero-Camacho and C.M. Hirata, Exploring circular polarization in the CMB due to conventional sources of cosmic birefringence , JCAP 08 (2018) 040 [ 1803.04505 ].\n- [16] S. Alexander, J. Ochoa and A. Kosowsky, Generation of Circular Polarization of the Cosmic Microwave Background , Phys. Rev. D 79 (2009) 063524 [ 0810.2355 ].\n- [17] D. Colladay and V.A. Kostelecky, Lorentz violating extension of the standard model , Phys. Rev. D 58 (1998) 116002 [ hep-ph/9809521 ].\n- [18] F. Finelli and M. Galaverni, Rotation of Linear Polarization Plane and Circular Polarization from Cosmological Pseudo-Scalar Fields , Phys. Rev. D 79 (2009) 063002 [ 0802.4210 ].\n- [19] S. Alexander, E. McDonough, A. Pullen and B. Shapiro, Physics Beyond The Standard Model with Circular Polarization in the CMB and CMB-21cm Cross-Correlation , JCAP 01 (2020) 032 [ 1911.01418 ].\n- [20] D. Ejlli, Magneto-optic effects of the cosmic microwave background , Nucl. Phys. B 935 (2018) 83 [ 1607.02094 ].\n- [21] D. Ejlli, On the CMB circular polarization: I. The Cotton-Mouton effect , Eur. Phys. J. C 79 (2019) 231 [ 1810.04947 ].\n- [22] B.R. Johnson et al., MAXIPOL: Cosmic Microwave Background Polarimetry Using a Rotating Half-Wave Plate , Astrophys. J. 665 (2007) 42 [ astro-ph/0611394 ].\n- [23] S.A. Bryan, T.E. Montroy and J.E. Ruhl, Modeling dielectric half-wave plates for cosmic microwave background polarimetry using a Mueller matrix formalism , Appl. Opt. 49 (2010) 6313 [ 1006.3359 ].\n- [24] ABS collaboration, Systematic effects from an ambient-temperature, continuously rotating half-wave plate , Rev. Sci. Instrum. 87 (2016) 094503 [ 1601.05901 ].\n- [25] SPIDER collaboration, A New Limit on CMB Circular Polarization from SPIDER , Astrophys. J. 844 (2017) 151 [ 1704.00215 ].\n- [26] M.P. Kennett and D.B. Melrose, Neutrino emission via the plasma process in a magnetized plasma , Phys. Rev. D 58 (1998) 093011 [ astro-ph/9901156 ].\n- [27] M. Lembo, M. Lattanzi, L. Pagano, A. Gruppuso, P. Natoli and F. Forastieri, Cosmic Microwave Background Polarization as a Tool to Constrain the Optical Properties of the Universe , Phys. Rev. Lett. 127 (2021) 011301 [ 2007.08486 ].\n- [28] R.C. Jones, A New Calculus for the Treatment of Optical SystemsI. Description and Discussion of the Calculus , J. Opt. Soc. Am. 31 (1941) 488.\n- [29] D. O'Dea, A. Challinor and B.R. Johnson, Systematic errors in cosmic microwave background polarization measurements , Mon. Not. Roy. Astron. Soc. 376 (2007) 1767 [ astro-ph/0610361 ].\n- [30] S. Giardiello et al., Detailed study of HWP non-idealities and their impact on future measurements of CMB polarization anisotropies from space , Astron. Astrophys. 658 (2022) A15 [ 2106.08031 ].\n- [31] W.C. Jones, T.E. Montroy, B.P. Crill, C.R. Contaldi, T.S. Kisner, A.E. Lange et al., Instrumental and Analytic Methods for Bolometric Polarimetry , Astron. Astrophys. 470 (2007) 771 [ astro-ph/0606606 ].\n- [32] G. Pisano et al., Development and application of metamaterial-based half-wave plates for the NIKA and NIKA2 polarimeters , Astron. Astrophys. 658 (2022) A24 [ 2006.12081 ].\n- [33] G. Pisano, G. Savini, P.A.R. Ade, V. Haynes and W.K. Gear, Achromatic half-wave plate for submillimeter instruments in cosmic microwave background astronomy: experimental characterization , Appl. Opt. 45 (2006) 6982.  \n- [34] G. Savini, G. Pisano and P.A.R. Ade, Achromatic half-wave plate for submillimeter instruments in cosmic microwave background astronomy: modeling and simulation , Appl. Opt. 45 (2006) 8907.\n- [35] I.L. Padilla et al., Two-year Cosmology Large Angular Scale Surveyor (CLASS) Observations: A Measurement of Circular Polarization at 40 GHz , 1911.00391 .\n- [36] D.B. Melrose, Collective plasma radiation processes , Ann. Rev. Astron. Astrophys. 29 (1991) 31.\n- [37] A. Lewis and S. Bridle, Cosmological parameters from cmb and other data: A monte carlo approach , Phys. Rev. D 66 (2002) 103511.\n- [38] K. Abazajian et al., CMB-S4 Science Case, Reference Design, and Project Plan , 1907.04473 .\n- [39] LiteBIRD collaboration, LiteBIRD: JAXA's new strategic L-class mission for all-sky surveys of cosmic microwave background polarization , Proc. SPIE Int. Soc. Opt. Eng. 11443 (2020) 114432F [ 2101.12449 ].\n- [40] Simons Observatory collaboration, The Simons Observatory: Astro2020 Decadal Project Whitepaper , Bull. Am. Astron. Soc. 51 (2019) 147 [ 1907.08284 ].\n- [41] Planck collaboration, Planck 2013 results. I. Overview of products and scientific results , Astron. Astrophys. 571 (2014) A1 [ 1303.5062 ].\n- [42] M. Gerbino, M. Lattanzi, M. Migliaccio, L. Pagano, L. Salvati, L. Colombo et al., Likelihood methods for CMB experiments , Front. in Phys. 8 (2020) 15 [ 1909.09375 ].\n- [43] A. Lewis and A. Challinor, 'CAMB: Code for Anisotropies in the Microwave Background.' Astrophysics Source Code Library, record ascl:1102.026, Feb., 2011.\n- [44] T. Namikawa, Exact CMB B-mode power spectrum from anisotropic cosmic birefringence , 2404.13771 .\n- [45] BICEP, Keck collaboration, Improved Constraints on Primordial Gravitational Waves using Planck, WMAP, and BICEP/Keck Observations through the 2018 Observing Season , Phys. Rev. Lett. 127 (2021) 151301 [ 2110.00483 ].\n- [46] CMB-S4 collaboration, CMB-S4: Iterative Internal Delensing and r Constraints , Astrophys. J. 964 (2024) 148 [ 2310.06729 ].\n- [47] K. Wolz et al., The Simons Observatory: pipeline comparison and validation for large-scale B-modes , 2302.04276 .\n- [48] Planck collaboration, Planck 2018 results. V. CMB power spectra and likelihoods , Astron. Astrophys. 641 (2020) A5 [ 1907.12875 ].\n- [49] BICEP2, Keck Array collaboration, BICEP2 / Keck Array x: Constraints on Primordial Gravitational Waves using Planck, WMAP, and New BICEP2/Keck Observations through the 2015 Season , Phys. Rev. Lett. 121 (2018) 221301 [ 1810.05216 ].\n- [50] S. De and H. Tashiro, Circular Polarization of the CMB: A probe of the First stars , Phys. Rev. D 92 (2015) 123506 [ 1401.1371 ].\n- [51] N. Lemarchand, J. Grain, G. Hurier, F. Lacasa and A. Ferté, Secondary CMB anisotropies from magnetized haloes - I. Power spectra of the Faraday rotation angle and conversion rate , Astron. Astrophys. 630 (2019) A149 [ 1810.09221 ].\n- [52] A. Cooray, A. Melchiorri and J. Silk, Is the cosmic microwave background circularly polarized? , Physics Letters B 554 (2003) 1.\n- [53] T.-C. Li, T. Zhu, W. Zhao and A. Wang, Power spectra and circular polarization of primordial gravitational waves with parity and Lorentz violations , 2403.05841 .\n- [54] A. Lewis, GetDist: a Python package for analysing Monte Carlo samples , 1910.13970 .", 'A Complete set of triangular plots': 'We show here the triangular plots for the five cases analyzed, displaying 1D and 2D posteriors for the parameters we consider relevant in our analysis, i.e. r , τ , β 2 E and β 2 V . The general behaviour of all plots looks similar, the inclusion of V-mode polarization substantially enhances our constraining power improving the bound on β 2 E . Furthermore, the addition of the de-lensing procedure is effective in reducing the upper limits of r and β 2 V (see the right panel of Figs. 5, 6, 7, 8, 9). In contrast, as expected, bounds on τ were found to be unaffected by both the de-lensing and the addition of V-mode data.  \nFigure 5 : One and two-dimensional posterior probability distributions for L dataset (see Tab. 2 for details). We report a subset of the total parameter space ( Λ CDM model + r + β 2 E + β 2 V ) showing TTTEEEBB in green and TTTEEEBBVV in magenta. The left panel omits de-lensing procedures ( A L = 1 ), while the right panel employs 70% de-lensing ( A L = 0 . 3 ).  \n<!-- image -->  \ne  \nv  \nr  \n2  \ne  \n2  \nv  \nTEB  \nTEBV  \n2  \n0.060  \n0.055  \n0.050  \n0.15  \ne  \n0.10  \n0.05  \n0.0010  \nv  \n0.0006  \n0.0002  \n2  \n2  \n4  \n×10  \n4  \n0.05  \n0.06  \n0.05  \n0.15  \n3  \n9  \n×10  \n4  \n2  \n0.060  \n0.055  \n0.050  \n0.10  \ne  \n0.05  \n0.0016  \nv  \n0.0010  \n0.0004  \n2  \n2  \n5  \n8  \nr  \nFigure 7 : One and two-dimensional posterior probability distributions for PSO dataset (see Tab. 2 for details). We report a subset of the total parameter space ( Λ CDM model + r + β 2 E + β 2 V ) showing TTTEEEBB in green and TTTEEEBBVV in magenta. The left panel omits de-lensing procedures ( A L = 1 ), while the right panel employs 70% de-lensing ( A L = 0 . 3 ).  \n<!-- image -->  \n0.05  \n0.06  \n0.03  \n0.10  \n2  \ne  \n0.001  \n2  \nv  \nFigure 6 : One and two-dimensional posterior probability distributions for LP dataset (see Tab. 2 for details). We report a subset of the total parameter space ( Λ CDM model + r + β 2 E + β 2 V ) showing TTTEEEBB in green and TTTEEEBBVV in magenta. The left panel omits de-lensing procedures ( A L = 1 ), while the right panel employs 70% de-lensing ( A L = 0 . 3 ).  \n<!-- image -->  \ne  \nv  \n0.002  \nr  \n2  \n0.048  \n0.064  \n0.05  \n0.10  \ne  \n0.003  \n2  \nv  \nTEB  \nTEBV  \n×10  \n4  \nTEB  \nTEBV  \n2  \n0.06  \n0.05  \n0.10  \ne  \n0.06  \n0.02  \n0.004  \n2  \nv  \n0.002  \nFigure 9 : One and two-dimensional posterior probability distributions for LS4 dataset (see Tab. 2 for details). We report a subset of the total parameter space ( Λ CDM model + r + β 2 E + β 2 V ) showing TTTEEEBB in green and TTTEEEBBVV in magenta. The left panel omits de-lensing procedures ( A L = 1 ), while the right panel employs 70% de-lensing ( A L = 0 . 3 ).  \n<!-- image -->  \ne  \nv  \ne  \nv  \nFigure 8 : One and two-dimensional posterior probability distributions for PS4 dataset (see Tab. 2 for details). We report a subset of the total parameter space ( Λ CDM model + r + β 2 E + β 2 V ) showing TTTEEEBB in green and TTTEEEBBVV in magenta. The left panel omits de-lensing procedures ( A L = 1 ), while the right panel employs 70% de-lensing ( A L = 0 . 3 ).  \n<!-- image -->  \ne  \nv  \ne  \nv'}