New Constraints on Spatial Variations of the Fine Structure Constant from Clusters of Galaxies

We have constrained the spatial variation of the fine structure constant using multi-frequency measurements of the thermal Sunyaev-Zeldovich effect of 618 X-ray selected clusters. Although our results are not competitive with the ones from quasar absorption lines, we improved by a factor 10 and ~2.5 previous results from Cosmic Microwave Background power spectrum and from galaxy clusters, respectively.

comparable to the direction of other so-called anomalies and/or dipoles [6,7]. However, complementary analyses carried out on the Cosmic Microwave Background (CMB) power spectrum [8,9], and studying the effect of such modulation on the orbital motion of the major bodies of the solar system [10], do not report any detection. Since their sensitivity is, however, one thousand times worse than the one achieved with QSO data, there is a need to carry out tests with other observational datasets such as multi-frequency measurements of the thermal Sunyaev Zeldovich effect (TSZ, [11]) in cluster of galaxies.
The TSZ effect induces secondary anisotropies on the CMB power spectrum. Such anisotropies are usually expressed in terms of the Comptonization parameter (Y c ) as ∆T T SZ = T 0 G(ν)Y c , where G(ν) is the TSZ spectral dependence, Y c is proportional to the integral along the line of sight of the pressure profile of the intra-cluster medium, and T 0 is the current value of the CMB black-body temperature. Moreover,ν = hν(z)/k B T CM B (z) is the reduced frequency, k B is the Boltzmann constant, h is the Planck constant and T CM B (z) is the CMB black-body temperature at the cluster location.
In the concordance ΛCDM model, the CMB black-body temperature is T CM B,std (z) = T 0 (1 + z). Any departure from such behavior should imply an extension of the concordance model. In particular classes of models, where a scalar field is coupled to the Maxwell F 2 term in the matter Lagrangian, the photon number conservation is violated. Thus, the violation of the standard T CM B (z) law is coupled to the variation of the fine-structure as and their common origin implies that the coefficient ǫ is expected to be of order unity [12]; for example, in a somewhat simplistic adiabatic limit ǫ = 1/4. Thus, if one is able to determine T CM B (z) at the cluster location using multi-frequency measurements of the TSZ effect, then Equation (1) can be used as a phenomenological relation to observationally test the spatial variation of α.
We apply the techniques developed in [13] to measure the TSZ emission at the location of 618 X-ray selected clusters. Then, we use such measurements to test the spatial variation of the fine structure constant. The article is organized as follows: in Section 2 we will describe the data; in Section 3 we will explain the methodology; in Section 4 we will illustrate and discuss our results; and finally, in Section 5, we will give our main conclusions.

Data
Galaxy clusters were selected from the ROSAT-ESO Flux Limited X-ray catalog (REFLEX, [14]), the extended Brightest Cluster Sample (eBCS, [15,16]) and the Clusters in the Zone of Avoidance (CIZA, [17]). The catalog lists 618 galaxy clusters for which all quantities of interest for our analysis are directly measured or derived: positions, spectroscopic redshifts, the central Comptonization parameter (y c,0 ), the scale radius at which the mean overdensity of the cluster is 500 times the critical density (r 500 ), and the corresponding angular size θ 500 = r 500 /d A (z). Here, d A (z) is the angular diameter distance of each cluster. At the location of each cluster we measure the TSZ anisotropies using the Planck 2013 maps (Maps were originally released in a Healpix format with resolution N side = 2048 [18]. Data are available at http://www.cosmos.esa.int/web/planck).
In March 2013, the Planck Collaboration made publicly available its nine Nominal maps covering the frequency channels from 30 to 857 GHz. Since the TSZ effect has a peculiar spectral dependence it can be distinguished from other components and reliably detected once the foreground emissions and the cosmological CMB signal are reduced. For that purpose, we apply a cleaning procedure, fully described in [13,19,20], to the High Frequency Instrument data (100-857 GHz), since they have better resolution and lower instrumental noise than the Low Frequency Instrument data . Before applying the cleaning procedure, we remove the intrinsic CMB monopole and dipole from the Nominal maps (We use the remove_dipole.pro facility of Healpix. It simultaneously performs a least square fit of the monopole and dipole on all the valid pixels and removes them from the map. More details can be found at http://healpix.jpl.nasa.gov/html/subroutinesnode86.htm). The cleaning procedure returns, for each Planck frequency channel, a foreground cleaned patch centered at the position of each cluster in our catalog. Then, in each patch and for each frequency, we measure the average TSZ temperature (δT /T 0 (ν)) over discs with angular extent equal to θ 500 . Hereafter, the averaged quantities are always evaluated on discs of radius θ 500 and indicated with a bar. To assign an error bar to each of the measurements, we carry out 1000 random simulations. In each one, we evaluate the mean temperature fluctuations in patches randomly placed out of the cluster positions and clean using the same procedure adopted for the real cluster population. Finally, we compute the correlation matrix between different frequencies (C ij ) where the average was computed over all simulations, The square roots of the diagonal elements of the correlation matrix are the error bars associated to the TSZ temperature anisotropies of each cluster.

Methodology
To estimate the CMB temperature at the cluster location, as well as to constrain the spatial variations of the fine structure constant, we use a Monte Carlo Markov Chain (MCMC) technique to explore the corresponding relevant parameter space. Our pipelines employ the Metropolis-Hastings sampling algorithm [21,22] with different (randomly set) starting points, and the Gelman-Rubin criteria [23] to test the convergence and the mixing of each chain. Specifically, we compute the ratio of the variances in the target distribution, and require it to be R < 1.03 (mathematical definitions are given in Section 3.2 of [24]). We run four chains that stop when they satisfy the criteria and contain at least 30,000 steps. Finally the chains are merged, and the expectation value of the 1D marginalized likelihood distribution is computed using all the steps.
First we need to estimate the CMB temperature at the cluster location. Having measured for each galaxy cluster the averaged TSZ emission as a function of frequency, we predict the theoretical counterpart at the same aperture ∆T (p, ν i )/T 0 = G(ν i , T CM B (z))Ȳ c , where p = [T CM B (z),Ȳ c ] are the free parameters of the model, andȲ c is the averaged Comptonization parameter. Finally, we compute the likelihood −2 log L = χ 2 (p) as where , and N = 4 is the number of data points (the four frequencies).
Once our simulations reached the convergence criteria, we estimated the CMB temperature for each individual cluster, and its 68% of confidence level error bar, from the 1D marginalized likelihood distribution. See [25] for a full description of the methodology. These measurements are used to estimate the α-anisotropies as ∆α Finally, we probe 3 different models allowing for a monopole (m) amplitude plus dipole (d) variation. Specifically, the first model assumes the functional form defined in [1] to fit quasar data (Model 1) where Θ is the angle on the sky between the line of sight of each cluster and the best fit dipole direction. Then, we fit the model used in [26] to the cluster data. Such a model (Model 2) allows a variation with the ΛCDM look-back time, Finally, to explore the possibility of an intrinsic dipole, we also try to fit the following relation to the CMB temperature measurements (Model 3) For each model we analyze four different configurations of the parameter space: (A) m = 0 for the Model 1 and 2, and unity for Model 3, the direction of dipole is fixed to the one from QSO [1], while the amplitude vary in the interval [−1, 1]; (B) the direction is still fixed, but we also vary the monopole in the intervals [−1, 1] for Model 1 and 2, and [0, 2] for Model 3; (C) and (D) we repeat the configurations of (A) and (B) but we also vary the dipole direction over the whole sky.

Results and Discussion
Our results show that the monopole and dipole amplitudes are always compatible with their expected valued in standard cosmology within 1. and (B), respectively. The direction of the "possible" dipole is always compatible with the one predicted from QSO in [1], but when they are varied the constraints on the monopole and dipole amplitudes are degraded since the number of free parameters is increased. Nevertheless, our best fit directions are always compatible with the one from QSO data [1,2] and other CMB anomalies [6,7], while they are several degrees away from the directions of the intrinsic CMB and the Dark Flow dipoles, from the kinetic CMB asymmetry, and from the CMB cold spot location [27][28][29][30][31][32][33][34][35][36][37][38][39]. All results are summarized in the Table 1 and, for easier comparison, the anomalies directions are plotted in Figure 1.
Since our analysis shows that the possible dipole direction coincides with the one obtained from QSO [1] and dark energy models [6,7], it could be interpreted as the fact that such dipoles, if confirmed with more extended and independent analysis, could represent signature of the failure of the ΛCDM model in favor of varying constant models or extended theories of gravity. Finally let us remark that, although we are not competitive with the QSO constraint from [1], our best fits for Model 2 provide a precision improvement of a factor ∼2.5, while being compatible with the results in [26]. For Model 1, we improve on the previous constraints in [8,9] by a factor ∼10. The improvement is mainly due to the fact that previous analyses, based on the CMB power spectrum, need to fit not only all "standard" parameters of the cosmological model, but also the additional ones related to the variation of fundamental constants; as a consequence, the resulting constraints were less tight.

Figure 1.
Directions in galactic coordinates for the spatial variation of the fine structure constant from [1,2] in blue and cyan, respectively, for the Dark Energy dipoles from [6] in green, Dark Flow direction from [27,[30][31][32][33] in red, for the Cosmic Microwave Background (CMB) asymmetry from [7] in yellow, and for our results from the analysis (C) of models 1,2 and 3 in orange, magenta and brown, respectively. Finally, we also indicate in blue and green dashed circles the direction of the Cold Spot anomaly and the intrinsic CMB dipole, respectively [28,29,34].

Conclusions
We have considered a model in which the variation of the fine structure constant is coupled to the variation of the CMB temperature. Using Equation (1), we have constrained the spatial variation of the fine structure constant using the Planck 2013 Nominal maps and a X-ray selected cluster catalog. First, we have cleaned the Planck maps and extracted the CMB temperature at the location of 618 X-ray selected clusters using the TSZ multi-frequency measurements. Then, we have carried out a statistical analysis to test three models allowing for monopole and dipole amplitudes, and describing both spatial variations of ∆α/α and of T CM B itself. All results of our analysis are summarized in Table 1. Table 1. Results of the analysis for Models 1, 2 and 3 from [25]. Although at the present time galaxy clusters are not competitive with high-resolution spectroscopic measurements in absorption systems along the line-of-sight of bright quasars, our analysis was improved by an order of magnitude previous constraints by the Planck Collaboration [8] and by a factor ∼2.5 the results from Galli (2013) [26].