On the CMB circular polarization. I. The Cotton-Mouton effect

Generation of cosmic microwave background (CMB) elliptic polarization due to the Cotton-Mouton (CM) effect in a cosmic magnetic field is studied. We concentrate on the generation of CMB circular polarization and on the rotation angle of the CMB polarization plane from the decoupling time until at present. For the first time, a rather detailed analysis of the CM effect for an arbitrary direction of the cosmic magnetic field with respect to photon direction of propagation is done. Considering the CMB linearly polarized at the decoupling time, it is shown that the CM effect is one of the most substantial effects in generating circular polarization especially in the low part of the CMB spectrum. It is shown that in the frequency range $10^8$ Hz $\leq \nu_0\leq 10^9$ Hz, the degree of circular polarization of the CMB at present for perpendicular propagation with respect to the cosmic magnetic field is in the range $ 10^{-13}\lesssim P_C(t_0)\lesssim 7.65\times 10^{-7}$ or Stokes circular polarization parameter $2.7 \times 10^{-13}$ K $\lesssim |V(t_0)|\lesssim 2 \times 10^{-6}$ K for values of the cosmic magnetic field amplitude at present in the range $10^{-9}$ G $\lesssim B\lesssim 8\times 10^{-8}$ G. On the other hand, for not perpendicular propagation with respect to the cosmic magnetic field we find $10^{-15}\lesssim P_C(t_0)\lesssim 6\times 10^{-12}$ or $2.72 \times 10^{-15}$ K $\lesssim |V(t_0)| \lesssim 10^{-11}$ K, for the same values of the cosmic magnetic field amplitude and same frequency range. Estimates on the rotation angle of the CMB polarization plane $\delta\psi_0$ due to the CM effect and constraints on the cosmic magnetic field amplitude from current constraints on $\delta\psi_0$ due to a combination of the CM and Faraday effects are found.


Introduction
In the last two decades, there have been many established observational facts about the nature and properties of the CMB and their possible implications in cosmology. Among these, it has already been established the fact that the CMB has a linear polarization with a degree of polarization at present of the order P L (t 0 ) 10 −6 . This linear polarization is believed to have been generated at the decoupling time mostly due to the Thomson scattering of the CMB photons on electrons. In general, if the incident electromagnetic radiation has an isotropic intensity distribution, Thomson scattering does not generate a net linear polarization. In the specific case of the CMB the fact that linear polarization has been initially observed by DASI, WMAP and BOOMERANG collaborations [1] and then re-confirmed by other collaborations, implies that at the decoupling time the CMB intensity did not have an isotropic distribution, a fact which is widely confirmed from the observation of the CMB temperature anisotropy. Another important consequence of the Thomson scattering is that it does not generates circular polarization in the case when electrons are assumed to be unpolarized. Based on this fact, during these years it has been erroneously assumed, at least from the theoretical point of view, that the CMB does not have a circular polarization at all even though there have been initial studies that might support its existence [2] and also initial experimental efforts to detect it [3].
In the recent years there have been several other theoretical studies exploring the possibility of CMB circular polarization from standard and non-standard effects and also new experiments such as MIPOL [4] and SPIDER [5] aiming to detect it. The MIPOL [4] collaborations reported an upper limit on the degree of circular polarization at present of P C (t 0 ) 7 × 10 −5 − 5 × 10 −4 at the frequency 33 GHz and at angular scales between 8 • and 24 • . On the other hand, the SPIDER collaboration reported an upper limit on the CMB circular polarization power spectrum ( + 1)C V V /(2π) < 255(µK) 2 for multipole momenta 33 < < 307 at the CMB frequencies ν 0 = 95 GHz and ν 0 = 150 GHz. From the theoretical point of view, studies based on non-standard effects that generate circular polarization include; the interaction of the CMB with a vector field via a Chern-Simons term [6], non commutative geometry [7] and free photon-photon scattering due to the Euler-Heisenberg Lagrangian term [8]. On the other hand, some theoretical studies of standard effects include; the electron-positron scattering in magnetized plasma at the decoupling time [9], the propagation of the CMB photons in magnetic field of supernova remnants of the first stars [10], the scattering of the CMB photons with cosmic neutrino background [11] and also the alignment of the cosmological matter particles in the post-decoupling epoch which results in an anisotropic susceptibility matter tensor [12]. For a recent and not complete review of the CMB circular polarization see Ref. [13]. Apart from the circular polarization generation effects mentioned above, there is a class of effects called magneto-optic effects which generate CMB circular polarization as well. In Ref. [14] and Ref. [15], I studied the most important magneto-optic effects which can generate CMB circular polarization when the CMB interacts with large-scale cosmic magnetic fields. Among the effects which I studied one of them is a standard effect, namely the CM effect, and the other effects are non-standard and include the vacuum polarization in an external magnetic field due to one loop electron-positron, one loop millicharged fermion-antifermion and the photon-pseudoscalar mixing in a magnetic field. For all these effects to occur it is necessary the presence of a magnetic field which gives rise to birefringence effects due to the fact that each of the photon states acquires different indexes of refraction in the presence of the magnetized plasma.
While it is well known that it does exist a magnetic field in galaxies and galaxy clusters with an order of magnitude of few µG, it is still not known if such a field is present also in the intergalactic space. The only information that we have about intergalactic magnetic fields are only in forms of upper and lower limits on the field magnitude at the present epoch. The upper limits on the magnetic field amplitude are found from observations of the CMB temperature anisotropy and from the rotation angle of the CMB polarization plane due to the Faraday effect. The temperature anisotropy upper limit is usually stronger than the Faraday effect limit, as reported by the Planck collaboration [16], where the limit from CMB temperature anisotropy is B e0 3 nG while the limit from the Faraday effect is B e0 1380 nG. One important aspect of these limits is that they differ from each other roughly speaking by three orders of magnitude and most importantly these limits do not mutually exclude each other from the simple fact that they are model depended. For a general review on large-scale cosmic magnetic field see Ref. [17].
One key aspect which distinguishes the CMB linear polarization with the CMB circular polarization, is that the former being generated at the decoupling time due to the Thomson scattering does not depend on the CMB photon frequency because of the nature of Thomson scattering, while the latter in most cases strongly depends on the CMB frequency. Because of this frequency dependence of the circular polarization, there is in some sense a kind of uncertainty on how to use and interpret the current limits obtained by experiments such as MIPOL and SPIDER since their limits are usually derived by observing the CMB in a specific frequency and it is not known how much substantial could be the signal at other frequencies.
In order to study and detected CMB circular polarization, it is very important to first identify the circular polarization (possibly standard) effects that generate substantial CMB circular polarization and identify their frequency band where the signal is the strongest. So far, there has been a tendency in the literature to study the circular polarization in the high-frequency range, namely for frequencies above ten or few hundred GHz. This tendency has been partially influenced by the fact that most important CMB experiments such as WMAP and Planck operates at these frequencies and therefore their data in these frequencies might be useful in some way. In addition, there are some effects such as photon-photon scattering in a magnetic field [14] and the free photonphoton scattering [8], [12] which are linearly proportional to the CMB frequency and one might hope that the higher is the frequency the stronger is the circular polarization signal. Even though this is true, the signal for such effects is still very weak even at very high frequencies to be detected in the near future.
Based on the facts discussed above, it is rather logical to explore the CMB circular polarization at low frequencies and study the magnitude of the signal. In this work, I study such possibility and concentrate on the CM effect in a large-scale magnetic field. As we will see, the CM effect is proportional to the square of the magnetic field amplitude, B 2 , and inversely proportional to the third power of the CMB frequency, namely ν −3 . It is especially the scaling law with the frequency of ν −3 which makes the CM effect one of the most important effects in generating circular polarization of the CMB. I partially studied this effect in a previous work [14] where some estimates of the degree of circular polarization were made for a specific configuration of the magnetic field with the respect to the photon direction of propagation. In this work, I study the CM effect in details for an arbitrary configuration of the magnetic field direction and for arbitrary magnetic field amplitude profile. By generalizing the CM effect to an arbitrary direction of the magnetic field with respect to the observer's direction, the system of differential equations for the Stokes parameters has additional terms with respect to the case studied in Ref. [14]. In addition, I also study in details the impact that the CM effect has on the rotation angle of the CMB polarization plane and its interaction with the Faraday effect.
This paper is organized in the following way: in Sec. 2, I discuss in a concise way the propagation of the electromagnetic radiation in a magnetized plasma and derive the elements of the photon polarization tensor in a cold magnetized plasma. In Sec. 3, I derive the system of differential equations for the Stokes parameters in an expanding universe. In Sec. 4, I find perturbative solutions of the equations of motion in various regimes. In Sec. 5, I calculate in details the generation of the CMB circular polarization due to the CM effect at present. In Sec. 6, I study the rotation angle of the CMB polarization plane due to the CM effect alone and also due to a combination of the CM and Faraday effects. In Sec. 7, I conclude. In this work I use the metric with signature η µν = diag[1, −1, −1, −1] and work with the rationalized Lorentz-Heaviside natural units (k B = = c = ε 0 = µ 0 = 1) with e 2 = 4πα. In addition in this work we use the values of the cosmological parameters found by the Planck collaboration [18] with Ω Λ 0.68, Ω M 0.31, h 0 0.67 with zero spatial curvature with Ω κ = 0.

Propagation of the electromagnetic waves in a magnetized plasma
In this section we give a detailed description of the propagation of electromagnetic waves in a cold magnetized plasma. This description is useful because it would allow us to understand in details how electromagnetic waves propagate in a cold magnetized plasma and which are the most common effects which give rise to birefringence effects in the medium. In this section we use the same notation as in Ref. [19] where basics of propagation of the electromagnetic waves in a cold magnetized plasma are presented in the appendix.
When electromagnetic waves (photons) propagate in a medium several effects manifest which include dispersion, absorption and scattering of the electromagnetic radiation. In connection with the dispersion phenomena, the effects of the medium on the incident electromagnetic wave are usually described in terms of the photon polarization tensor Π ij (i, j = x, y, z) with components in a given cartesian coordinate system where the medium is at rest. Consequently, in a medium the free Maxwell equations in momentum space, in absence of external currents, get modified to for a plane electromagnetic wave travelling into the medium. Here ω is the photon energy and we used the expression k ij = ωn ij with k ij being the photon momentum tensor and n ij being the index of refraction tensor of photons in medium. We may see that the role of Π ij in (1) is to give to photons an "effective mass" in the medium.
In the case when the medium is isotropic, we have that n ij is a diagonal tensor with diagonal entries corresponding to the photon indexes of refraction in medium where n ii = 1. In the case when photons propagate in vacuum, we have that n ij = δ ij and we get the on-shell photon relation ω = k 2 where k is the photon wave-vector and Π ij = 0. The explicit expression of the photon polarization tensor Π ij depends on the induced currents that enter a given problem. In this work we are interested in a cold magnetized plasma which is quite common situation in astrophysics and cosmology. We assume that the magnetized plasma is with almost no collisions, globally neutral and homogeneous. In addition, there is not an external electric field, namely E e = 0 and the presence of the external magnetic field B e locally breaks the isotropy of the plasma since it singles out a preferred direction in a given region of space where the plasma is located.
In the cold magnetized plasma approximation, consider now an incident electromagnetic wave propagating along the observer's z axis which points to the East, in a magnetized plasma with external magnetic field vector B e = B en . Heren = [cos(Θ), sin(Θ) cos(Φ), sin(Θ) sin(Φ)] is a unit vector in the direction of the external magnetic field B e and Θ, Φ are, respectively, the polar and azimutal angles between the magnetic field B e and x and y axes. As shown in Ref. [19], the medium polarization vector P satisfies the equation of motion where E is the electric field of the incident electromagnetic wave, ω pl = 4παn e /m e is the plasma frequency, n e is the free electron number density and ω c = eB e /m e is the cyclotron frequency. In Eq. (2) the dot symbol (·) above P denotes the derivative with respect to the time.
Assume that the fields evolve in time harmonically at a given point x and then let us write where ω is the incident electromagnetic wave energy. By using the expressions in (3) into Eq. (2) and then solving for the components of P , after we get the following solution in terms of the incident electric field components E j , in the case when ω = 0 and ω = ω c where χ ij (ω) are the components of the electric susceptibility tensor , χ yx = χ * xy , The expressions for the components of χ ij in (5) are valid for an incident electromagnetic wave with an arbitrary direction of propagation with respect to B e . In addition, the components χ ij do not explicitly depend on x but only implicitly through B e (x, t) which enters in ω c . Another fact is that the expressions for χ ij in (5) are valid for arbitrary external magnetic profile B e (x, t). After these general comments about (5), let us find the components of the photon polarization tensor in a cold magnetized plasma. In order to do that we have to relate the components of χ ij with Π ij . It is well known that the components of the index of refraction tensor n ij are related to the relative permittivity tensor 1 ε ij through the relation n 2 ij = ε ij . On the other hand, the relative permittivity tensor ε ij is related to the electric susceptibility tensor χ ij , through the relation χ ij = ε ij − δ ij . By using these relations into the relation (1), we get Π ij = −ω 2 χ ij .
By using the expressions for χ ij in (5) into (6), we get Now by using the constraint (10) into the (i = x, y) components of the electric field in (9), we get the following equation for the transverse components of the electric field whereΠ ij for (i, j = x, y) is the effective photon polarization tensor of the transverse electromagnetic field in the cold magnetized plasmaΠ where the components of Π ij are given in (7). The effective expression for the polarization tensor in (12) takes into account the mixing of the longitudinal electromagnetic wave in plasma with the usual transverse electromagnetic waves. From expressions (12) and (7), we find the following expressions for the components ofΠ ij The expressions given in (13) are the most general form of the elements ofΠ ij in a cold magnatized plasma. As already mentioned above they take into account the mixing of the longitudinal electric field with the usual transverse electric field. We may note that this contribution in (13) is inversely proportional to ω 4 − ω 2 ω 2 pl − ω 2 ω 2 c + ω 2 c ω 2 pl sin 2 (Θ) sin 2 (Φ) = 0. The latter condition is satisfied as far as ω > 0 and where we must have ω 2 pl + ω 2 c ≥ 2ω c ω pl | sin(Φ) sin(Θ)| in order to have real and positive roots of the quadratic equation. Another important question to ask is for what minimum frequencies we have propagating transverse electromagnetic waves? This can be seen by requiring that all spatial derivatives on the left hand side in Eq. (9) are zero, namely a non propagating electric field in space. In that case we would have Π ij − ω 2 δ ij E j = 0 where nontrivial solution exist only if det(M ij ) = 0 with M ij ≡ Π ij − ω 2 δ ij . However, the solution of det(M ij ) = 0 in terms of ω would be quite complicated in the case when all components of Π ij = 0. For this reason it would be convenient to rotate the coordinate system in such a way that Φ = π/2 and Θ = π/2, namely B e is along the direction of propagation of the electromagnetic wave. Under a rotation of the coordinate system we have that M ij in the new coordinate system is related to the old M ij through M ij = R il R jm M lm and E j = R jk E k where R il is an orthogonal rotation matrix with unit determinant. In the rotated coordinate system the equation M ij E j = 0 becomes M ij E j = M ij (Φ = π/2, Θ = π/2)E j = 0. Consequently, the condition det(M ij ) = 0 is equivalent to det[M ij (Φ = π/2, Θ = π/2)] = 0. Now by requiring that det[M ij (Φ = π/2, Θ = π/2)] = 0 and after doing some algebra we find that the lower bound on the frequencies for propagation are ω > ω pl and ω > ±ω c /2 + ω 2 c + ω 2 pl /4 /2.

Solutions of the equations of motion of the Stokes parameters
In the previous section we derived the most general form of the elements of the photon polarization tensor in a cold magnetized plasma for arbitrary direction of propagation of the electromagnetic waves with respect to the external magnetic field B e . In this section, we focus on our attention on deriving the equations of motion of the Stokes parameters in an expanding universe and provide perturbative solutions of the equations of motion. As in the previous section, let us consider an electromagnetic wave propagating along the z direction in a cartesian reference system with wave vector k = (0, 0, k) in a cold magnetized plasma with arbitrary direction of the external magnetic field B e . The linearized equations of motion for the vector potential transverse components A x and A y in an unperturbed FRW metric for the CMB photons are given by [14] i∂ where A x and A y are respectively the transverse components of the vector potential A of the CMB photons with respect to the x and y axes, is the Hubble parameter, I is a 2 × 2 identity matrix and M is the mixing matrix which is given by where M x = −Π xx /(2ω), M y = −Π yy /(2ω) and M CF = −Π xy /(2ω). The term M CF = M C + iM F takes into account the combination of the CM and Faraday effects in a magnetized plasma. In order to describe the polarization of the light and more precisely in our case of the CMB photons, it is better to work with the Stokes parameters rather than the wave equation (15). The procedure in obtaining the equations of motion of the Stokes parameters has been presented in [14] and it consists on two steps; first write the equations of motion for the polarization density matrix ρ based on the wave equation (15) and second, express the polarization density matrix in terms of the Stokes parameters in order to get the equations of motion of the latter quantities. The equations of motion of the polarization density matrix in an unperturbed FRW metric are given by [14] ∂ρ ∂t where D = (3/2)H(t)I is the damping matrix which takes into account the damping of the electromagnetic waves in an expanding universe due to the Hubble friction. In our case the field mixing matrix M is Hermitian, namely M = M † since in our case we do not include any process which might change the number of photons due to decay or absorption in the medium 2 . Now by using the connection between the Stokes parameters and the polarization density matrix elements as shown in Ref. [20], see also the appendix of Ref. [14], we get the following equations of motion of the effective Stokes parameterṡ I(k,n, t) = −3H(t)I(k,n, t), where we have defined ∆M ≡ M y − M x with the dot sign above Stokes parameters indicating the time derivative with respect to the cosmological time t. For simplicity, in (18) we have dropped the symbols B e , Φ and Θ which do appear in the elements of M .
The system of linear differential equations (18) can be written in a more compact form asṠ(k,n, t) = A(k, t)S(k,n, t) where S = (I, Q, U, V ) T is the Stokes vector formed with the Stokes parameters and A(k, t) is the time dependent coefficient matrix which is given by In most cases is more convenient to express the quantities in A as a function of the photon temperature T rather than the cosmological time t, so, in this case one needs to express the time derivative in an expanding universe as ∂ t = −HT ∂ T in the equations of motion of the Stokes vector, namely S (k,n, T ) =Ã(k, T )S(k,n, T ). At this stage is more convenient to write the matrixÃ(k, T ) as the sum ofÃ(k, where in an expanding universe the wave-vector k = k(T ) is a function of the temperature T . We may note that with respect to the case when the direction of B e is in the xz plane as studied in Ref. [14], for arbitrary magnetic field direction, do appear the terms 2M C in the matrix B. The appearance of these terms which makes possible the mixing of the Q parameter with U and V parameters, complicate the situation with respect to the case when M C = 0.

Series solution of the polarization equations of motion
In the previous section, Sec. 3, we found the equations of motion of the Stokes parameters in an expanding universe for an arbitrary direction of the external magnetic field B e with respect to the electromagnetic wave direction of propagation. In this section we focus on our attention on perturbative solutions of the equations of motion in some limiting cases. Before aiming to find these solutions, it is very important to explicitly calculate each term which enters the matrix B(k, T ) since it will be very useful in what follows. Let us recall the definitions of M F ≡ −Im{Π xy }/(2ω), M C ≡ −Re{Π xy }/(2ω) and ∆M ≡ M y − M x = (Π xx −Π yy )/(2ω). Now by using the expressions of the photon polarization tensor given in (13) we get The expression for the elements of the matrix B given in (21), which are the most general ones for arbitrary magnetic field direction and magnitude, can be further simplified by making some reasonable assumptions on the parameters. Since in this work we concentrate on the CMB frequency spectrum we have that ω ω pl and ω ω c . In order to see this, let us calculate explicitly the numerical values of the parameters. The numerical value of the angular plasma frequency which enters the expressions in (21) can be written as ω pl = 5.64 × 10 4 n e /cm 3 (rad/s) or ν pl = ω pl /(2π) = 8976.33 n e /cm 3 (Hz) for the frequency. On the other hand the numerical value of the cyclotron angular frequency is given by ω c = 1.76 × 10 7 (B/G) (rad/s) or ν c = 2.8 × 10 6 (B/G) (Hz). However, in the case of CMB photons propagating in an expanding universe, we can express the time t in terms of the cosmological temperature T as t = t(T ) as we did in the previous section. Therefore, the conditions ω ω pl and ω ω c , in an expanding universe, are respectively satisfied when where we expressed ν(t) = ν 0 [a(t 0 )/a(t)] = ν 0 (T /T 0 ) with ν 0 being the frequency of the electromagnetic radiation at the present time t = t 0 at the temperature T = T 0 , a(t) being the universe expansion scale factor and B 0 = B(t 0 ) = B(T 0 ) is the magnetic field strength at the present time. Here we expressed the number density of free electrons as n e (t) = n e (T ) 0.76 n B (T 0 )X e (T )(T /T 0 ) 3 where n B (T 0 ) is the total baryon number density at the present time and X e (T ) is the ionization function of the free electrons. The factor of 0.76 takes into account the contribution of hydrogen atoms to the free electrons at the post decoupling time. By taking for example n B (T 0 ) 2.47 × 10 −7 cm −3 as given by the Planck collaboration [18], and expressing a(t 0 )/a(t) = T /T 0 , we can write the conditions (22) as Given the fact that the present day CMB photon frequencies are in the frequency part above ν 0 ≥ 10 8 Hz, the condition given in (23) is well satisfied for physically reasonable values of X e (T ) and B 0 . With these considerations in mind, we can simplify (21) for ω ω pl and ω ω c . Consequently, we can write the expressions in (21) as From the expressions (24) we may note that each expression within the square brackets is composed of a first term of trigonometric functions and a second term which is the product of trigonometric functions with terms ω 2 pl /ω 2 or ω 2 pl ω 2 c /ω 4 . However, since we are in the regime when ω ω pl and ω ω c , we also have ω 2 pl /ω 2 1 and ω 2 pl ω 2 c /ω 4 1. This fact tells us that in the case when the trigonometric functions in the first and second terms within the square brackets in (24) are different from zero, the second term is usually much smaller than the first term. In order to see this, let us consider the case when Θ = 0, namely when the magnetic field has components only along the x. In this case 2 1, so the contribution coming from the second term can be completely neglected. One can see that by making similar examples, the contribution of the second terms within the square brackets in (24), which arise due to the mixing of the longitudinal electromagnetic wave with the transverse waves, can be neglected with respect to the first terms. Consequently, in the regime studied in this work ω ω pl and ω ω c , we have that

Neumann series solutions
Here we present a Neumann series solutions of the equations of motion by making use of the perturbation theory. Let us concentrate on the full equation S (k,n, T ) = [B(k, T ) + (3/T )I 4×4 ] S(k,n, T ) and omit from now on the dependence of the Stokes vector onn and k and matrix B on k. From the equation of motion of the Stokes vector, the term 3/T is a term which takes into account the damping of the fields in an expanding universe.
In the case when there is not a magnetic field the solution of the equation It is worth to stress since now that the effective scaling of the Stokes vector in an expanding universe is not (T /T i ) 3 but (T /T i ) 2 as discussed in details in Ref. [14]. In the case when the magnetic field is present, namely when

T ). In this case, the equations of motion for S (T ) = [B(T ) + (3/T )I 4×4 ] S(T ) in components becomeS
The system of the first order of linear differential equations given in (26) cannot be solved exactly except in some particular cases. However, one of the main characteristic of a linear system of first order of differential equations is that its general solution is given byS(T ) =M (T )S(T i ), whereM (T ) is the solution matrix. Consequently if we put the general solutionS(T ) =M (T )S(T i ) in (26), we get that the solution matrixM satisfies the equatioñ with the initial conditions thatM lj (T i ) = I 4×4 . Therefore the solution of the system (26) is reduced to the solution of the differential equations for the matrixM lj in (27). The system of differential equations given in (27) can formally be solved as a convergent Neumann series in the case when the non zero elements of the matrix B lm (T ) satisfy In order to find the parameter space arising from the conditions T T i dT B lm (T ) < 1, we need to evaluate explicitly each element in the matrix B ij (T ). In each element in B ij (T ) enters the product H(T )T , where the Hubble parameter in the case of zero spatial curvature is given by where after the decoupling epoch the contribution of relativistic particles to the total energy density and consequently to the Hubble parameter can be safely neglected. In addition since the contribution of the cosmological parameter to the Hubble parameter is important only for low redshifts, we may approximate the Hubble parameter in our calculations as The conditions that The condition |M F (T )| < 1 is also satisfied by the following stronger condition where we used the fact that |sin(Θ) sin(Φ)| ≤ 1 in |M F (T )| < 1. So, the condition (31) is a stronger condition on the parameter space with respect to the case when the term |sin(Θ) sin(Φ)| is taken into account. In case when |sin(Θ) sin(Φ)| → 0, the condition |M F (T )| < 1 is in principle satisfied for any finite value of the parameters B e0 , ν 0 and T . On the other hand, the conditions |M C (T )| < 1 and |∆M(T )| < 1 are respectively satisfied by the much stronger conditions 6.05 × 10 31 Hz ν 0 where we used the fact | sin(2Θ) cos(Φ)| ≤ 1 in |M C (T )| < 1 and that sin 2 (Θ) cos 2 (Φ) − cos 2 (Θ) ≤ 1 in |∆M (T )| < 1. Again, in the cases when | sin(2Θ) cos(Φ)| → 0 and sin 2 (Θ) cos 2 (Φ) − cos 2 (Θ) → 0, the conditions |M C (T )| < 1 and |∆M(T )| < 1 are in principle satisfied for any finite values of the parameters B e0 , ν 0 and T .
In order to find the parameter space for the conditions T T i dT B lm (T ) < 1 it is necessary to know the expression for the free electron ionization function X e (T ). This function satisfies a complicated differential equations as shown in Ref. [21] and in general it is calculated by solving the differential equation numerically. In Fig. 1a the plots of X e (T ), X e (T )T 1/2 and X e (T )T 3/2 as a function of the CMB temperature T are shown. In the temperature interval 57.22 K≤ T ≤ 2970 K the curve of the ionization function X e (T ) is obtained by solving the differential equation for X e (T ) as given in Ref. [21], where the lower limit T = 57.22 K corresponds to the start of reionization epoch at redshift z ion ∼ 20 and the upper limit corresponds to the CMB decoupling temperature T i = 2970 K for redshift 1 + z 1090. The complete re-ionization is reached approximately at z ion 7. The evolution of X e (T ) in the temperature interval 21.8 K ≤ T ≤ 57.22 K has been obtained by a smooth interpolation of the curve X e (T ) in the interval 57.22 K ≤ T ≤ 2970 K with X e (T ) = 1 in the interval 2.725 K ≤ T ≤ 21.8 K. By using the numerical solutions found for X e (T ) as descibed above and plotted in Fig. 1a, we get the following values for 4.45 × 10 6 (K 5/2 ) and With these values of the integrals, the stronger conditions (31) and (32) are respectively satisfied when  Having found the parameter space where the condition for the convergence of the Neumann series certainly holds, we are at the position now to calculate the Stokes parameters. Since the Stokes vector is given byS(T ) = M (T )S(T i ), the only thing that we need to calculateS(T ) is to calculate the elements of the matrixM (T ) given in (28). Since we are in the regime where T T i dT B lm (T ) < 1, it will be sufficient for our purposes to truncate the Neumann series (28) at the second order. Now by looking at the structure of the matrix B(T ) in (20), we may note that B 1j = B j1 = B jj = 0 with the rest of the elements different from zero. Let us define for commodity Then the vanishing elements ofM areM 1j =M j1 = 0 while the non zero elements ofM are given bỹ (35) The matrix elements found in (35) allow us to find explicit expression for the Stokes vectorS(T ) up to second order in the perturbation theory. Consequently, the expression of the Stokes parameters are given bỹ It is very important to stress that the expressions found of the Stokes parameters in (36) are valid for arbitrary direction of the external magnetic field with respect to the photon propagation, namely for arbitrary Θ, Φ and for arbitrary profile of B e (T ) and n e (T ).

Power series solution for dominant Faraday effect
In the previous section we found Neumann series solutions in the case when the conditions |M F (T 0 )| < 1, |M C (T 0 )| < 1 and |∆M(T 0 )| < 1 are satisfied. However, we did not make any specific assumptions on the relative magnitude among |M F (T 0 )|, |M C (T 0 )| and |∆M(T 0 )|, namely we did not specify which of these terms is bigger than the others. As we can see from The conditions in (37) After by inserting the expansion (38) and B(T ) = B 1 (T ) + B 2 (T ) in the equation (27) and collecting the terms with the appropriate power in , we get the following matrix system of equations where for simplicity we suppressed the matrix elements indexes in B 1 , B 2 andM . The system of equations 39 has to be solved with the initial conditionsM In order to solve the system in (41) let us multiply the third equation with the imaginary unit (i) and after sum it with the second equation. Then we get (42) We may observe that (42) is a first order non homogeneous linear differential equation for the functionM Now by equating the real and imaginary parts of the left hand side of (43) with the those of the right hand side and directly integrating the first and the fourth equations in (43) The expressions in (44) allows us to find recursively the elements ofM ij (T ) to the order m + 1 in the case when the elements at the order m are known. Since we already know the elements ofM ij at the order m = 0 as given in (40), we can recursively calculate those at the order m + 1. Let us define for simplicity By using the definitions in (45) and the expressions in (44) we get the following expressions for m = 0 of the matrix elements of M (1) We can proceed in the same way to find the matrix elements of 2M (2) ij starting from the elements of M (1) ij given in (46). However, for our purposes it will be sufficient to consider only the elements ofM ij up to the first order in . By using the expressionS j (T ) =M jl (T )S l (T i ) M (0) jl (T ) + M (1) jl (T ) S l (T i ), we get for the elements of the Stokes vector the following expressionsĨ It is worth to remind that the expressions found in (47) are valid for M F (T ) = 0 or when | sin(Θ) sin(Φ)| = 0.

T ) S (T ). Then the equationS(T ) = [B 1 (T ) + B 2 (T )]S(T ) becomes
As we may note, so far we did not make any assumption on the matrix B 1 (T ) and in principle it can be also a null matrix depending on the situation. After doing some lengthy calculations we get We can solve Eq. (48) as a convergent Neumann series as we did in sec. 4.2 as far as T T i dT M ij (T ) < 1. By using the matrix expression (49) in (48), we get the following solution forS(T ) up to the first order We may note that the solution (51) exactly coincides with the solutions found in (47) which we found by using the regular perturbation theory. The solution (51) has been found without any restriction on the magnitude and sign of M F (T ) which is different from the result of Sec. 4.2 where we worked under the assumption that M F (T ) = 0, which for fixed and non zero values of B e0 , ν 0 and T , is equivalent to | sin(Θ) sin(Φ)| = 0. This fact tells us that the condition on the Faraday effect term M F (T ) = 0 is not necessary in order to find the solution (51) and that the condition M F (T ) = 0 comes out only in the regular perturbation theory. On the other hand, in this section in order to use the Neumann series expansion we required that T T i dT M ij (T ) < 1 or equivalently that |L

Degree of circular polarization
In the previous section we found perturbative solutions of the equations of motion of the Stokes parameters in two different regimes by using perturbation theory. In this section, we focus on our attention on generation of circular polarization, where in specific, we study the expected degree of circular polarization at present time and the expected rotation angle of the CMB polarization plane. We separate our analysis by first studying the solutions found in Sec. 4.1 and second study those found in Sec. 4.2. In what follows, we consider the evolution of the CMB polarization and rotation angle of the polarization plane starting from the decoupling epoch at the temperature T = T i until at the present time at the temperature T = T 0 . Moreover, we consider the CMB at the decoupling epoch partially polarized where it acquires only a linear polarization due to the Thomson scattering off the CMB photons on electrons with no initial circular polarization, namely Q i = 0, U i = 0 and V i = 0 as studied in Ref. [14].

Case when |M
Let us consider first the generation of the circular polarization and calculate its degree of polarization at present time where T = T 0 in the case when |M F (T 0 | < 1, |M C (T 0 | < 1 and |∆M (T 0 )| < 1. By using the expression for the Stokes parametersĨ(T ) andṼ (T ) found in (36), the degree of circular polarization of the CMB at present is given by It is quite convenient at this stage to normalize the CMB intensity at the decoupling time to unityĨ i = 1. In addition, we have thatĨ In what follows, we assume that V i = 0 if not specified otherwise. In order to calculate P C (T 0 ) we need to calculate explicitly the matrix elementsM 42 (T 0 ) and M 43 (T 0 ). For the first term entering inM 42 (T 0 ) we have while for the second term we have while the second term inM 43 (T 0 ) is given by Since all terms inM 42 andM 43 depend on the angles Θ and Φ and because some terms inM 42 andM 43 are equal to zero when averaged over the angles Θ and Φ, it is more convenient to calculate the root mean square of the degree of circular polarization instead of the mean value. By using the expressions (54)-(57) inM 42 andM 43 we get (by suppressing for the moment the units) One important thing about expression (58) is that the second terms proportional to Q 2 i and U 2 i must be in magnitude smaller that the first terms. The reason is because these terms correspond to second order terms in perturbation theory where their magnitudes must be smaller than the first order terms in order to have a convergent series. This fact implies that care must be used in order to choose the values of B e0 and ν 0 in order to evaluate P rms C (T 0 ). However, since we are in the regime where the constraints (33) must be satisfied, usually there is not reason to worry about since the values of the parameters ν 0 and B e0 that satisfy (33) automatically keep the magnitudes of the second order terms smaller than the first ones.
In Figs. 2 and 3 plots of the root mean square of the degree of circular polarization P rms C (T 0 ) as functions of the magnetic field amplitude B e0 and ν 0 are shown. In obtaining the plots we used the expression 58 where we expressed U i = rQ i with r being a parameter which can have either sign and which value is not a priori known. In addition, we have chosen those values of B e0 and ν 0 that satisfy the constraints (33). Usually if the stronger constraint on the Faraday effect term is satisfied, namely the first constraint on the left hand side in (33), the remaining two stronger constraints which arise from |M C (T 0 )| < 1 and |∆M(T 0 )| < 1 are also satisfied. We may observe from Figs.

Case when |M
In the case when |M F (T 0 )| = 0 and |M C (T 0 )| < 1, |∆M (T 0 )| < 1 the constraints on ν 0 and B e0 are much less stringent than in the previous section. In fact, for finite values of ν 0 and B e0 which interest us, the only possibility for the condition |M F (T 0 )| = 0 to hold is only when | sin(Θ) sin(Φ)| = 0 which occurs either when Θ = nπ or Φ = nπ with n ≥ 0. In both cases the direction of the magnetic field is perpendicular to the direction of photon propagation where M F (T 0 ) = 0 and M C (T 0 ) = 0. Consequently, the constraints in (33) reduce to only the constraint (Hz/ν 0 ) 3 (B e0 /G) 2 < 8.35 × 10 −39 which correspond to the stronger constraint on |∆M(T 0 )| < 1, namely the region within the black line in 1b.
In order to calculate the degree of circular polarization, let us use the results obtained in 4.3 that we found without any restriction on the magnitude of M F (T 0 ). For absent Faraday effect M F (T 0 ) = 0 and consequently a vanishing M C (T 0 ) = 0, the degree of circular polarization is given by (59) As we can see from (59) the degree of circular polarization for transverse magnetic field depend only on ∆M(T 0 ). The most important thing is that we do not have anymore the constraints on M (F, C) but only those on |∆M(T 0 )| < 1. In Figs. 4 and 5 plots of the degree of circular polarization for transverse magnetic field as a function of ν 0 , B e0 and |r| are shown. We may observe in Fig. 4 that for higher values of B e0 and lower values of ν 0 , the acquired degree of circular polarization of the CMB is quite substantial and be comparable with that of the linear polarization for some values of the parameters. For example, as we can see from 4, for B e0 = 8 × 10 −8 G, we get P C (T 0 ) 7.65 × 10 −7 for ν 0 = 10 8 Hz and P C (T 0 ) 7.65 × 10 −10 for ν 0 = 10 9 Hz. It is worth to point out that the expression (59) can also be obtained by using the perturbative approach used in the previous section.

Case when |M
In the case when |M F (T 0 )| ≥ 1 the situation is more complicated with respect to the previous cases. One aspect is that in (33) only the last two inequalities must be satisfied while the first inequality has not to be satisfied anymore. This fact tells us that the allowed region of parameters is that within the black line in Fig. 1b and that outside the region within the dotted line. In this case the degree of circular polarization can be calculated by using the results of Sec. 4.3 that we derived for arbitrary values of M F (T ) The expression (60) is valid for any value of M F (T ) and for |M C (T 0 )| < 1, |∆M(T 0 )| < 1 even though in this section we study the case when |M F (T )| ≥ 1. The main difficulty on calculating P C (T 0 ) analytically stands from the fact that in all terms M F , G C and ∆G enters the ionization function which does not have any known analytic expression. In order to find an analytic expression for P C , in this section we approximate X e (T ) X e whereX e Then we have that (62) where we used the identity cos(α − β) = cos(α) cos(β) + sin(α) sin(β). Let us define x ≡ AT 3/2 where dT = (2/3)A −2/3 x −1/3 dx and get  ν 0 (Hz) Be0=10 -9 G, |r|=0.1 Be0=10 -9 G, |r|=10

Rotation angle of the polarization plane
In the previous section we studied the generation of the CMB circular polarization by calculating explicitly P C (T 0 ) in various regimes. In this section, we focus on our attention on the rotation angle of the CMB polarization plane from the decoupling epoch until today. Apart from generating circular polarization, the CM effect also generates linear polarization with non zero Stokes parametersQ(T ) andŨ (T ). At a given cosmological temperature T the rotation angle of the polarization plane is given by where we must haveQ(T ) = 0. Let us write ψ(T ) = ψ(T i )+δψ(T ) where ψ(T i ) is the angle of the CMB polarization plane at the temperature T i = 2970 K corresponding to the decoupling time in the common reference frame used to study the CMB and ψ(T ) is the angle of the polarization plane at temperature T < T i . Here δψ(T ) is the amount of the rotation angle of the polarization plane from the decoupling time until at the time corresponding to the temperature T and it is the quantity which interests us. Since for the frequency range of interest in this work the magnitude of the effects which we study are in general small, namely |M F (T )| < 1, |∆M(T )| < 1, |M C (T )| < 1, and because experimentally δψ(T 0 ) is constrained to a small quantity (in radians), we expect that the rotation angle of the CMB polarization plane from decoupling epoch until at present to be a small quantity |δψ(T )| 1. In this case by using the trigonometric identity we can write  .
In order to calculate δψ(T 0 ) we need to calculate each matrix element in (70) at T = T 0 . Consequently, we have thatM where we numerically calculated T i T 0 dT X e (T )T 3/2 T i T dT X e (T )T 3/2 = 9.88 × 10 12 (K 5 ). We also get the following expression for On the other hand, we get the following expressions for where we used the numerically integrated value of The last expression to calculate is One important thing to observe about expressions (73) and (74) is that the second order of iteration terms must be smaller than unity because we are in the regime when 0 ≤ |M F (T )| < 1, |∆M(T )| < 1 and |M C (T )| < 1. Consequently one must choose the values of the parameters carefully in such way that such conditions are met. However, as far as the conditions in (33) are satisfied we do not have to worry about what we said above. On the other hand for the terms appearing in (72) and (71) the second order of iteration terms can be equal or even larger than the first order terms since we are in the situation where only 0 ≤ |M F (T )| < 1, |∆M(T )| < 1 and |M C (T )| < 1 must be met and no further conditions between the relative magnitude of these terms is needed at this stage. All told, let us consider first the case when M F (T 0 ) = 0 which occurs when Θ = 0, namely absent Faraday effect. In this case we also have B = 0. From expression (70) and dropping for simplicity the units we obtain where we neglected the sub-leading order term proportional to δψ 2 on the left hand side in expression (70). It is evident that since we are working under the condition |∆M(T 0 )| < 1 we have that the right hand side of (75) is less than one for any value of r. In Fig. 8 plots of the present epoch CMB rotation angle of the polarization plane δψ 0 = δψ(T 0 ) given by expression (75) are shown for various values of the parameters. We may note that substantial rotation of the polarization plane occurs only at low frequencies and for higher values of B e0 . On the other hand for higher values of the frequency and lower values of the magnetic field amplitude δψ 0 is extremely small.  ). There are several interesting facts about (76) that we need to clarify. First thing is that δψ(T 0 ) depends on the angles Θ and Φ through the temperature independent constants A, B and C. If we take for simplicity first the average over the angles Θ and Φ where A = 0, BC = 0 and then the absolute value we get | 2δψ(T 0 ) | 1 1 + r 2 r 1 − 1.6 × 10 6 (9 A 2 /(4X 2 e )) − 9.88 × 10 12 C 2 X −2 e 1 − 1.6 × 10 6 (9 A 2 /(4X 2 e )) − 9.88 × 10 12 B 2 X −2 where we used again geometric expansion and neglected the terms proportional to B 2 − C 2 A 2 and to B 2 − C 2 B 2 . Therefore in the case when we take the average value over the angles Θ and Φ, the average value of the rotation angle is still dominated by the CM effect and the contribution of the Faraday effect is completely negligible at the second order in perturbation theory. Second thing, is that in the case when we do not take the average value over the angles we have that δψ(T 0 ) is dominated by the Faraday effect term and we can derive from (76) where again we used geometric series expansion and used the fact in the case when |M F (T 0 )| < 1 we have that also that |∆M(T 0 )|, |M C (T 0 )| |M F (T 0 )|. Indeed, the latter condition can be easily understood because the values of |∆M(T 0 )| and |M C (T 0 )| for those values of the parameters that satisfy the condition |M F (T 0 )| < 1 are much smaller than the values of |M F (T 0 )| as far as the angles Θ and Φ are not zero.
In case when |M F (T )| > 1, |∆M(T )| < 1 and |M C (T )| < 1, we cannot use anymore the same expressions that we used above because of the fact that |M F (T )| > 1. In this case we can use the expressions for the Stokes parameters found in Sec. 4.3 up to first order in perturbation theory for arbitrary value of M F (T ). The expressions forQ(T ) andŨ (T ) are given by (51) and read where we may notice that there is no contribution to the linear polarization from the CM effect at the first order in perturbation theory. From (79) we get for the rotation angle of the polarization plane In Fig. 9 plots of the root mean square of the CMB rotation angle of the polarization plane given in (81) due to the Faraday effect are shown. In Fig. 9a the region in magenta colour between 0.36 • ≤ δψ 0 ≤ 4.3 • represents the region where experimentally do exist constraints on δψ 0 . The black points represent the values of these constraints found by different experiments at different frequencies. For example, the first black point from the top in Fig. 9a is the constraint found by BOOM3 experiment [22] at the frequency ν 0 = 145 GHz where |δψ 0 | = 4.3 • 4 . The second point from the top is the constraint found by BICEP 1 [23] at the frequency ν 0 = 129 GHz where |δψ 0 | = 2.77 • . The third and fourth points from the top are the constraint found by the QUaD collaboration [24] at the frequencies 100 GHz and 150 GHz where the constraint on δψ 0 are respectively |δψ 0 | = 1.89 • and |δψ 0 | = 0.83 • . The fifth point from the top is the constraint found by the WMAP9 collaboration [25] at the frequency ν 0 = 53 GHz where |δψ 0 | = 0.36 • , see also Ref. [26] for a general discussion on the constraints on δψ 0 . It is worth to stress that in Fig. 9 we have chosen the values of the parameters in such a way that the conditions |M C (T 0 )| < 1 and |∆M(T 0 )| < 1 are satisfied, see Fig. 1b.
We may observe from Fig. 9a that most of experimental constraints on |δψ 0 |, represented by the black points, are within the grey region between the magnetic field values 10 −8 G ≤ B e0 ≤ 8 × 10 −8 G. The only exception is the constraint found by WMPA9 where the magnetic field amplitude corresponding to |δψ 0 | = 0.36 • is by equation (81) B e0 = 7.47 × 10 −10 G. In Fig. 9b plots of the root mean square δψ 0 as a function of the CMB frequency  are shown. Similarly to the Fig. 9a, the black points represent the constraints on |δψ 0 |. The first point from the bottom correspond to the WMPA9 constraint [25], the second point from the bottom correspond to the QUaD constraint [24] and the third one correspond to the BICEP 1 [23]. Each of the points lie on the corresponding frequency line where the measurements are done. For example, the QUaD constraint on |δψ 0 | = 0.83 • is consistent with B e0 = 1.38 × 10 −8 G, while the BICEP 1 constraint on |δψ 0 | = 2.77 • is consistent with B e0 = 3.4 × 10 −8 G.

Conclusions
In this work, we have studied the generation of the CMB circular polarization and rotation angle of the CMB polarization plane due to the CM effect in a large-scale cosmic magnetic field. We worked with the Stokes parameters and derived a system of differential equations for their evolution in an expanding universe. In the equations governing the evolution of the Stokes vector, we included all standard magneto-optic effects which manifest in a magnetized plasma which are the CM and Faraday effects. Then we looked to solutions of the equations of motion of the Stokes parameters in different regimes by using several perturbative approaches such as the regular perturbation theory and the Neumann series expansion. The equations of motion that we found in (18) are a generalization to the equations of motion found in Ref. [14] in the case of an arbitrary direction of B e with respect to the photon direction of propagation and for arbitrary magnetic field profile. For an arbitrary direction of B e , the equations of motion (18) include two additional terms proportional to M C (T ), which, would be absent in the particular case when the magnetic field B e is in the same plane with the wave-vector k. These two terms proportional to M C (T ) make possible the mixing of the Q(T ) and V (T ) parameters with each other. The magnitude of the degree of circular polarization for the CM effect depends on several parameters where the most important ones are the CMB frequency ν 0 and the magnetic field amplitude B e (x, t 0 ). In addition, other parameters which play also an important role are the angles Θ and Φ. Consequently, depending on the values of these parameters, in this work, we divided our analysis of the CM effect in three major regimes. In the regime where |M C, F (T 0 )| < 1 and |∆M(T 0 )| < 1, the degree of circular polarization assumes the lowest values as shown in Figs. 2 and 3, where at best its value reaches P C (T 0 ) 10 −17 . The reason for such low values of P C (T 0 ) stands from the fact that the condition |M F (T 0 )| < 1, drastically restricts the values of the parameters ν 0 to very high frequencies and the values of B e0 to very low ones.
In the case when the Faraday effect is completely absent, which happens when the direction of B e is perpendicular to the direction of propagation of the CMB photons, we essentially have that M F (T ) = 0 and the generation of circular polarization is maximal. The absence of the Faraday effect for such specific configuration results in an enhancement of the generation of the CMB circular polarization. For such case, we have found in Sec. 5.2 that the degree of circular polarization can reach values close to the CMB degree of linear polarization in the CMB low-frequency part of the spectrum. The maximum values of the degree of circular polarization are reached in the case when we concentrate at the frequency ν 10 8 Hz, where depending on the magnetic field amplitude, the degree of circular polarization is in the range 1.19 × 10 −10 P C (T 0 ) 7.65 × 10 −7 for magnetic field values 10 −10 G ≤ B e0 ≤ 8 × 10 −8 G. These results are plotted in Figs. 4 and 5 for different values of the parameters.
In the case when the Faraday effect is present and in particular when |M F (T )| ≥ 1, the generation of the CMB circular polarizartion is strongly suppressed with respect to the case of absent Faraday effect where M F (T ) = 0. However, the generation of the CMB circular polarization in the case when |M F (T )| ≥ 1 is usually much efficient than that in the case when 0 < |M F (T )| < 1 which we studied in Sec. 5.1. Even in the case |M F (T )| ≥ 1 the degree of circular polarization depends on B e0 and ν 0 , where in some specific range of these parameters, the degree of circular polarization scales with the frequency as P C (T 0 ) ∝ ν −1 0 and with the magnetic field amplitude as P C (T 0 ) ∝ B e0 , see for example the expression (67). As shown in Fig. 6b, the degree of circular polarization can reach values in the range 3 × 10 −15 P C (T 0 ) 2.5 × 10 −13 for magnetic field values 10 −9 G ≤ B e0 ≤ 8 × 10 −8 G at the frequency ν 0 10 8 Hz, |Q i | = 10 −6 and |r| = 1. At the frequency ν 0 10 9 Hz, the values of P C (T 0 ) decrease exactly by an order of magnitude since P C (T 0 ) ∝ ν −1 0 in the frequency range considered. On the other hand, P C (T 0 ) ∝ |r|, so, higher values of |r| give higher values of P C (T 0 ) and vice-versa for smaller values of |r|.
Apart from generating circular polarization, the CM effect also generates linear polarization and this fact is evident in all expressions of the Stokes parameters that we found in Sec. 4. In connection with linear polarization, in this work, we have studied the rotation angle of the CMB polarization plane due to the CM effect in the case when the Faraday effect is absent and in combination with the Faraday effect when it is present. In the case when it is present only the CM effect, the rotation angle is δψ(T 0 ) ∝ ν −6 0 B 4 e0 and consequently, significant rotation of the polarization plane occurs in the low-frequency part of the CMB spectrum and for higher values of the magnetic field amplitude. We have found in Sec. 6 that at ν 0 10 8 Hz, the rotation angle in units of degrees is in the range 10 −3 ≤ δψ(T 0 ) ≤ 1 for magnetic field amplitude in the range 10 −9 G ≤ B e0 ≤ 8 × 10 −8 G, |r| = 1 and |Q i | = 10 −6 , see Fig. 8. For higher frequencies, |δψ(T 0 )| acquires extremely smaller values which are uninteresting for any practical purpose.
In the case when the rotation angle δψ(T 0 ) is due to a combination of the CM and Faraday effects the situation slightly changes with respect to the case of absent Faraday effect. If we are interested in taking the average value of δψ(T 0 ), the Faraday effect gives null contribution while the CM effect gives in average the same contribution as it does in the case of absent Faraday effect. If we do not take the average value of δψ(T 0 ) or we take the root mean square, the Faraday effect usually dominates over the CM effect in the case when it is present. One important aspect is that in case we take the root mean square of δψ(T 0 ), the Faraday effect generates significant rotation of the polarization plane depending on the CMB frequency and magnetic field amplitude. As shown in Fig. 9, the Faraday effect can generate substantial rotation of the polarization plane especially in the low-frequency part especially for ν 0 10 10 Hz. In the high-frequency part of the spectrum, namely for frequencies above 10 GHz, the rotation angle is still large depending on the magnetic field amplitude. An interesting fact is that most of the constraints on δψ(T 0 ) experimentally found correspond to magnetic field amplitudes in the range 10 −8 G B e0 8 × 10 −8 G.
If we have to consider currents limits on δψ(T 0 ) as a potential indicator of the existence of the large-scale magnetic field and consequently a non zero rotation angle of the polarization, these limits would allow us to make some predictions on the signal of the circular polarization due to the CM effect. Indeed, if we consider the hypothesis that the rotation angle is due to the Faraday effect only (root mean square value) and that most experimental constraints on δψ(T 0 ) would suggest a magnetic field with amplitude approximately 10 −8 G B e0 8 × 10 −8 G, we would have that the signal of circular polarization for these values of B e0 would be quite substantial. For these values of the magnetic field, in the case when the field is perpendicular to the photon direction of propagation, we would have a circular polarization signal at present in the range 3 × 10 −8 K V (T 0 ) 2 × 10 −6 K at ν 0 10 8 Hz and a signal of 3 × 10 −11 K V (T 0 ) 2 × 10 −9 K at ν 0 10 9 Hz, see Fig.  4. In the case when the magnetic field is not perpendicular, the signal of the circular polarization is reduced by many orders of magnitude and is in the range 10 −14 K V (T 0 ) 10 −13 K depending on the angles Θ and Φ, see Fig. 5.
Based on the arguments presented so far, it seems quite plausible that the CM effect is probably the most substantial effect in generating CMB circular polarization. However, the strongest signal of the circular polarization is located in the CMB frequency range 10 8 Hz ν 0 10 9 Hz. In the high-frequency range the signal of circular polarization due to the CM effect is much smaller than in the low-frequency part of the spectrum, but still, the signal is not negligible and can be comparable with the vacuum polarization circular polarization signal in a magnetic field and to that due to free photon-photon scattering. If we assume that there is not any major difficulty in arranging an experiment aiming to detect the circular polarization in the low-frequency part of the CMB spectrum, then it is quite logical to concentrate our attention in this frequency part of the spectrum where the signal is the strongest and more likely to be detected in a relatively short time.