Revisiting cosmic microwave background radiation using blackbody radiation inversion

Blackbody radiation inversion is a mathematical process for the determination of probability distribution of temperature from measured radiated power spectrum. In this paper a simple and stable blackbody radiation inversion is achieved by using an analytical function with three determinable parameters for temperature distribution. This inversion technique is used to invert the blackbody radiation field of the cosmic microwave background, the remnant radiation of the hot big bang, to infer the temperature distribution of the generating medium. The salient features of this distribution are investigated and analysis of this distribution predicts the presence of distortion in the cosmic microwave background spectrum.

A blackbody is an ideal object which can absorb all of the incident radiation of all frequency. The total power radiated per unit frequency per unit solid angle by a unit area of a blackbody emitter can be expressed by Planck's law 1,2 where ν is frequency, T is the absolute temperature, h is Planck's constant, k is Boltzmann's constant and c is the speed of light. Usually telescopes are used to measure this power spectrum of any celestial object. But due to its finite field of view a telescope can observe a small portion of the sky at any time. These small portions consist of different blackbody radiators with different temperature T and each of them are in thermal equilibrium. When a collection of blackbodies with probability distribution α(T) and temperature T is considered, the total radiated power per unit area is given by the integration over the distribution as 3 where W(ν) is the radiated power per unit frequency per unit area and per unit solid angle and α(T) is the probability distribution of temperature of the blackbody. The dimension of α(T) is 1 K . The blackbody radiation inversion problem aims to find the probability distribution of temperature from the radiated power spectrum.
In practice, a set of discrete values of W(ν) are available experimentally. By using this set of data, α(T) can be calculated by blackbody inversion method.
Equation (3) is the first kind of Fredholm integral equation and is an ill-posed problem. Bojarski was the first to propose a solution to this problem using Laplace transform with an iterative process 4 . Since then various other methods have been proposed for solving this problem like Tikonov regularization method 5 , universal function set method 6 , Mellin transform method 7 , modified Mobius inverse formula 8 , variational expectation www.nature.com/scientificreports/ maximization method 9 , maximum entropy method 10 , regularised GMRES method 11 . There are also several other methods available in literature as solutions to this problem [12][13][14][15] . However, the required number of input data is large in the existing method. The number of data points required for successful inversion is 50 in 6 , 50 in 11 and 32 in 15 . In this paper, a simple and robust method for blackbody radiation inversion is developed which uses 3 input data. The size of programming is also small in comparison to the previous methods. The present method for blackbody radiation inversion reduces the complexity of the overall program significantly.
This method is applied to obtain the probability distribution of temperature of the universe using cosmic microwave background radiation (CMB) from COBE, FIRAS data 16 . "Method and validation" section describes the method and its validation and "Cosmic microwave background radiation" section describes the application of this method in CMB radiation.

Method and validation
Equation (3) takes all possible values of temperature into consideration; hence the limit runs from zero to infinity. Here it is assumed that the temperature of black body radiators in a collection of blackbodies vary from T 1 to T 2 and they have a finite frequency range of ν. Therefore, Eq. (3) can be written as Using change of variable T = T 1 + (T 2 -T 1 )t, Eq. (4) becomes 15 And a(t) = α(T 1 + (T 2 -T 1 ) t). The required interval of a(t) is [0,1]. Equation (6) informs that the problem of solving α(T) is equivalent to solving a(t). In the present article, an analytical function represented by Eq. (8) is proposed as a(t).

Equation (8) can be expanded as
In this method we are trying to obtain the probability distribution of temperature. The nature of the probability distribution is expected to be close to gaussian. So, Eq. (8) is chosen in such a way that for large vales of k 3 , the e −k 2 3 t part in the sine hyperbolic function is very small. When k 2 ≈ k 2 3 , Eq. (8) represents a gaussian distribution provided the value of k 1 is small.
The lower and upper limits of temperatures (T 1 , T 2 ) are taken as 1 K and 6 K respectively. The motive behind this choice is that we will use this method to analyse the CMB spectrum and it closely resembles a blackbody radiation at a temperature range similar to T 1 and T 2 . Then, t = T−1 5 or, where k 1 , k 2 and k 3 are three determinable parameters such that the interval of α(t) is [0,1] and T is absolute temperature.
The data have been simulated by using model function Eq. (11), Equation (11) is used in Eq. (4) in place of α(T) and the values of G(ν) are calculated. This process is repeated with different frequencies ν . These simulated data are put in the left-hand side of Eq. (6). Three of such equations for three different values of frequency ν are obtained. These three equations with three unknowns k 1 , k 2 and k 3 are then solved. Thus, the function α(T) is obtained with these three parameters.
, we calculate α(T) for three different frequencies of 5 × 10 11 Hz, 6 × 10 11 Hz and 7 × 10 11 Hz. Since we will be using this method in the CMB spectrum, the range of frequency is chosen such a way www.nature.com/scientificreports/ that it resembles the frequencies in the data we have 16 . It is observed in Fig. 1a that we have reconstructed the temperature distribution that resembles the model temperature distribution. The difference between b(T) and α(T) is expressed as d 1 (T) = b(T) − α(T) and it is plotted in Fig. 1b against absolute temperature. The I I value obtained from Fig. 1b is 0.0119 for T = 3 K. Here I is the value of b(T) and I is the value of d 1 (T) for a specific temperature T. The method is sensitive to the chosen frequency. To quantify this sensitivity, we choose sets of frequencies as i, where i includes three frequencies with ν 1 = i × 10 11 Hz, ν 2 = (i + 1)× 10 11 Hz and ν 3 = (i + 2)× 10 11 Hz. This set is then used in Eq. (5) to calculate α(T). The standard deviation from the model function is calculated by Eq. (12) where N is the number of data used for the calculation of standard deviation and we have taken N = 51.
It is inferred from Fig. 2 that the standard deviation is less for the sets of i = 1, 2, 5 and 6. Therefore, it is expected to use either of these set of frequencies. All the values of k 1 , k 2 and k 3 that are calculated during the validation, are listed in Table 1.

Cosmic microwave background radiation
Cosmic microwave background radiation is the afterglow as predicted by the hot big bang model. The presence of such radiation in the universe was first suggested in the late 1940s 17 . It was only in 1965 when a signal was first detected which was reported to be coming from every direction of the observed sky 18 . This was the first detection of the radiation which later became to be known as cosmic microwave background radiation. The study of CMB can unravel the mysteries of the initial stage of the universe and its evolution for the last 13.7 billion years. Right after its first detection a lot of work has been done on CMB [19][20][21][22][23][24][25] .
The first detection showed the radiation to be isotropic, i.e. similar in every direction. But subsequent studies showed that the radiation is in fact anisotropic in nature 26,27 . More recent studies focus on the different types of www.nature.com/scientificreports/ distortions in the CMB spectrum [28][29][30][31][32] . It suggests that the radiation is not of a blackbody with single temperature, rather it is a superposition of different blackbodies that are at different temperatures. When several blackbodies of different temperatures are mixed together it creates y and µ type distortions 32,33 . In this paper we have calculated the distortions present in the CMB spectrum. The blackbody radiation inversion (BRI), as discussed in "Method and validation" section, is applied to the cosmic microwave background radiation for obtaining the temperature of the universe and the probability distribution of temperature. We have used the data of COBE FIRAS to calculate intensity 16 . For the input of the BRI, the spectral irradiance I(λ) is transformed to the power spectrum W(ν) according to the relation where λ (= c ν ) is the wavelength. For each value of W(ν) corresponding to a particular frequency ν, we have an integral equation in Eq. (6). Three of such equations are taken to calculate k 1 , k 2 and k 3 of Eq. (10). Table 2 shows the values of k 1 , k 2 and k 3 that we have calculated.    The moments of different order of α(T) are calculated by using Eq. (16) where n is the order of the moment and µ is the mean value.
First order moment or mean value is calculated as The mean temperature is 2.69 K, which is close to the average value 2.725 K 23 . Second order moment or variance is calculated as So, standard deviation σ indicates the uncertainty in temperature which is 0.195 K 23 . Third order standardised moment or Skewness is calculated as An ideal normal distribution has a skewness of 0. We get a positive skewness which describes its deviation from ideal behaviour. A positive skewness suggests that the tail of the distribution right to the mean is more extended than the left-hand side tail 34,35 .
Fourth order standardised moment about mean Kurtosis represents the peakedness and tailedness of a distribution. An ideal normal distribution has a kurtosis of 3, so 3 is subtracted from β 2 to measure the deviation from ideal normal behaviour. In our calculation γ 2 (= β 2 − 3) yields 0.0563, a positive number. A distribution with positive kurtosis is called Leptokurtic. A Excess kurtosis γ 2 = β 2 − 3, where β 2 = µ 4 σ 4 = 3.056. Table 2. k 1 , k 2 and k 3 values for different probability distribution function are listed. Three wavelengths are used to obtain three integral equations. The probability distributions are denoted by a 1 (T), a 2 (T), a 3 (T) and a 4 (T) corresponding to a set of frequency. www.nature.com/scientificreports/ positive kurtosis means that the peak of the curve is slightly higher than the normal distribution while the tail and shoulder portion is slightly pushed towards the mean value 34,35 .
From the standard deviation σ and mean value µ, a Gaussian function ( 1 2σ 2 ) is constructed in Eq. (24). α(T) and s(T) are plotted against absolute temperature in Fig. 3a. And the difference between α(T) and s(T), expressed as d 2 (T) = α(T) -s(T) is plotted against absolute temperature in Fig. 3b.
A deviation from the ideal gaussian behaviour is observed. The I I value obtained from Fig. 3b is 0.0194 for T = 2.5 K. Here I is the value of s(T) and I is the value of d 2 (T) for a specific temperature T. The deviation in Fig. 3b (0.0194) is larger than the deviation in Fig. 1b (0.0119). So, this deviation in Fig. 3b is not due to the error in the inversion method we have used. A small deviation from ideal Gaussian behaviour is also predicted when non-extensive case is considered 36 . The temperature distribution of CMB is found to be primarily between 2 and 3.5 K.
To verify the accuracy of our method to obtain probability distribution of temperature, we reconstructed the intensity of the radiation by using the calculated α(T) in Eq. (2) for different frequencies ν. W(ν ) is then converted to I(λ) by using Eq. (13). Figure 4 displays the overlay of reconstructed data on the original data of COBE FIRAS.  Table 3 for the original and reconstructed data. The order of the error in the reconstructed spectrum (~ 10 −13 W/m 2 × μm × sr) is larger than the error in the original spectrum (~ 10 −14 W/m 2 × μm × sr). It is evident from Fig. 4 that the present method can faithfully reconstruct the original data.
In this paper, the original data of COBE/FIRAS are used as the input in the blackbody radiation inversion problem. These data are thus mathematically processed to obtain the distortion of the CMB spectrum. The standard deviation between the original and reconstructed data is 0.142 × 10 −10 W/m 2 × μm × sr or 5.209 × 10 −20 W/ m 2 × Hz× sr for the wavelength of 1049 μm. The deviation is the distortion present in the CMB spectrum. The spread in the probability distribution of the temperature (Fig. 3a) suggests that there are multiple blackbodies with different temperatures ( T = σ ). Due to this mixing of blackbodies, the original spectrum becomes distorted. So, the calculated deviation is interpreted as the distortion of the CMB spectrum.
The COBE data shows a spectrum similar to a perfect blackbody 21 . But the possible distortions are limited by the maximum sensitivity of the instrument. It has y distortion of |y|< 1.5 × 10 −5 and µ distortion of | µ|< 9.0 × The present set of data, collected by the COBE/FIRAS satellite is not sensitive enough to detect the distortions beyond the 10 −5 order. More precise datasets are required to study these distortions. The TRIS, used between 1996 to 2000, had the limit of µ < 6.0 × 10 −537 . A balloon borne instrument ARCADE (Absolute Radiometer for Cosmology, Astrophysics, and Diffuse Emission) used in 2006, had the upper limit of µ < 6.0 × 10 −438 . Two new projects, PIXIE 39 and PRISM 40 aim to find the distortions with 10 3 -10 4 times better sensitivity than COBE/ FIRAS. PIXIE has I = 5 × 10 −26 W/m 2 srHz and detection of |y|= 1 × 10 −8 and | µ|= 5.0 × 10 −8 is possible. PRISM is better than PIXIE with I = 6 × 10 −27 W/m 2 srHz and sensitive to y and µ distortion of ~ 10 −9 .  www.nature.com/scientificreports/ The distortions obtained in our calculation are limited by the sensitivity of available measured data. The experiments planned in the future 39,40 are expected to provide data with better precision that will help in carrying out more precise calculations and lead to a better understanding of CMB.

Discussion
A novel method of blackbody radiation inversion is studied. This technique is then applied to study cosmic microwave background radiation and some of its most important features. We have described the deviation of the temperature probability distribution from ideal gaussian distribution. The distortion in the spectrum, caused due to mixing of blackbodies are mathematically described as well. Our approach is much simpler than the existing techniques and the computational bulkiness is significantly reduced. While we can obtain the probability distribution of the temperature effectively, the present method is not completely general in nature. The frequency range needs to be selected to minimise the error in the calculation. www.nature.com/scientificreports/