## Examining the local Universe isotropy with galaxy cluster velocity dispersion scaling relations

A. Pandya 1 K. Migkas 2 3 1 A. Stanford 4 F. Pacaud 1 G. Schellenberger 5 6 5

, , , T. H. Reiprich , , , , L. Lovisari , , M. E. Ramos-Ceja 7 , N. T. Nguyen-Dang 8 , and S. Park 9 10 ,

- 1 Argelander-Institut für Astronomie (AIfA), Universität Bonn, Auf dem Hügel 71, 53121 Bonn, Germany e-mail: apandya@astro.uni-bonn.de
- 2 Leiden Observatory, Leiden University, PO Box 9513, 2300 RA Leiden, the Netherlands
- 3 SRON Netherlands Institute for Space Research, Niels Bohrweg 4, NL-2333 CA Leiden, the Netherlands
- 4 University of California, Davis, CA 95616, USA
- 5 Center for Astrophysics | Harvard &amp; Smithsonian, 60 Garden St., Cambridge, MA 02138, USA
- 6 INAF, Istituto di Astrofisica Spaziale e Fisica Cosmica di Milano, via A. Corti 12, 20133 Milano, Italy
- 7 Max Planck Institute for Extraterrestrial Physics, Giessenbachstrasse 1, 85748 Garching, Germany
- 8 Institute for Astronomy and Astrophysics of Tübingen (IAAT), Germany.
- 9 Korea Astronomy and Space Science Institute, Daejeon 34055, Republic of Korea
- 10 University of Science and Technology, Daejeon 34113, Republic of Korea

Received ...; accepted ...

## ABSTRACT

Context. In standard cosmology, the late Universe is assumed to be statistically homogeneous and isotropic. This assumption suggests that the expansion rate of the Universe, as measured by the Hubble parameter, should be the same in all directions. However, a recent study based on galaxy clusters by Migkas et al. (2021) found an apparent spatial variation of approximately ∼ 9% in the Hubble constant, H 0 , across the sky. In the study, the authors utilised galaxy cluster scaling relations between various cosmologydependent cluster properties and a cosmology-independent property, i.e., the temperature of the intracluster gas ( T ). A positiondependent systematic bias of T measurements can, in principle, result in an overestimation of apparent H 0 variations. Therefore, it is crucial to confirm or exclude this possibility.

Aims. In this work, we search for directional T measurement biases by examining the relationship between the member galaxy velocity dispersion and gas temperature ( σ v -T ) of galaxy clusters. Both measurements are independent of any cosmological assumptions and do not su ff er from the same potential systematic biases. Additionally, we search for apparent H 0 angular variations independently of T by analysing the relations between the X-ray luminosity and Sunyaev-Zeldovich signal with the velocity dispersion, L X -σ v and Y SZ -σ v .

Methods. To study the angular variation of scaling relation parameters, we determine the latter for di ff erent sky patches across the extra-galactic sky. We constrain the possible directional T bias using the σ v -T relation, as well as the apparent H 0 variations using the L X -σ v and Y SZ -σ v relations. We utilise Monte Carlo simulations of isotropic cluster samples to quantify the statistical significance of any observed anisotropies. We find and rigorously take into account a correlation of L X and Y SZ residuals.

Results. No significant directional T measurement biases are found from the σ v -T anisotropy study. The probability that the previously observed H 0 anisotropy is caused by a directional T bias is only 0 002%. On the other hand, from the joint analysis of the . L X -σ v and Y SZ -σ v relations, the maximum variation of H 0 is found in the direction of (295 ◦ ± 71 ◦ , -30 ◦ ± 71 ) with a statistical ◦ significance of 3 64 . σ , fully consistent with Migkas et al. (2021).

Conclusions. Our findings, based on the analysis of new scaling relations utilising a completely independent cluster property, σ v , strongly corroborate the previously detected spatial anisotropy of galaxy cluster scaling relations. The underlying cause, e.g., a spatial H 0 variation or large-scale bulk flows of matter, remains to be identified.

Key words. Galaxies: clusters: general - X-rays: galaxies: clusters - Cosmology: observations - large-scale structure of Universe

## 1. Introduction

Galaxy clusters are the largest gravitationally virialised systems in the Universe. They are crucial tools for astrophysical and cosmological studies as they provide valuable information about the large-scale structure and evolution of the Universe (Pratt et al. 2019). Galaxy clusters can be observed in multiple wavelengths, providing insights into di ff erent cluster properties. Galaxy cluster scaling relations describe the correlation between various important cluster properties using simple power laws (Giodini et al. 2013; Lovisari &amp; Maughan 2022). These relations were predicted first by Kaiser (1986), and the model predicts that the ob- jects of di ff erent sizes are merely scaled versions of each other. Due to this, it is also referred to as the self-similar model.

The evaluation of certain cluster observable properties involves cosmological assumptions. For instance, the X-ray luminosity ( L X) and the total integrated Compton parameter ( Y SZ) rely on cosmological distances inferred from the estimated cluster redshift. The relation between these two quantities is a function of the assumed cosmological model. However, the measurement of some properties, such as intracluster gas temperature, T , and galaxy velocity dispersion, σ v , depend very weakly on such cosmological assumptions. Valuable insights about di ff er-

ent aspects of cosmology can be obtained by utilising the scaling relations between these two types of cluster properties (for a recent review, see Migkas (2024)).

The cosmological principle is a fundamental assumption of the standard cosmological model, and it states that the Universe is homogeneous and isotropic on su ffi ciently large scales. The Friedmann equations - the basic set of equations that underpins most cosmological models - are derived under this assumption. The Λ CDM model assumes that matter should converge to an isotropic behaviour for cosmic volumes greater than 150 Mpc. This implies that the cosmic expansion rate, cosmological parameters, and the distances to extra-galactic objects depend solely on redshift, z , regardless of direction. Therefore, given the fundamental importance of this assumption, it is necessary to assess the validity of it.

The cosmic microwave background (CMB) observations support the cosmological principle within the CMB rest frame, exhibiting remarkable isotropy at small angular scales (Bennett et al. 1994; Bennett et al. 2013; Planck Collaboration et al. 2020b). However, the CMB provides us with the rest frame of cosmic radiation in the early Universe, whereas examining the (an)isotropic behaviour of the matter rest frame in the local Universe (e.g., galaxies, clusters, and supernovae at z ≲ 0 2 . -0 3) . through CMB data is particularly challenging and not extensively explored, with very few exceptions (e.g., Yeung &amp; Chu 2022). Therefore, assessing if cosmic matter behaves isotropically in the late Universe is crucial. Thus, additional cosmological probes are required to test the cosmological principle in the late Universe. Some of these probes include the use of Type Ia supernovae (SNIa) (Appleby et al. 2015), infrared quasars (Secrest et al. 2022), the distribution of optical (Javanmardi &amp; Kroupa 2017; Sarkar et al. 2019) and infrared galaxies (Yoon et al. 2014; Rameez et al. 2018) and the distribution of distant radio sources (Singal 2011; Rubart &amp; Schwarz 2013; Tiwari &amp; Nusser 2016; Colin et al. 2019) and gamma-ray bursts ( ˇ ípa &amp; R Shafieloo 2017; Andrade et al. 2019).

Galaxy clusters have become a useful tool to test the cosmological principle in recent years. A new method based on cluster scaling relations was introduced by Migkas &amp; Reiprich (2018). Migkas et al. (2020) (hereafter M20) and Migkas et al. (2021) (hereafter M21) then applied various scaling relations for galaxy clusters to examine the isotropy of the local Universe quantitatively. The key principle behind these studies involved pairing a cosmology-dependent quantity (e.g., L X and Y SY) with the cosmology-independent T . This pairing allowed them to draw conclusions on the directionality of cosmological parameters. They found a (9 0 . ± 1 7)% variation in the Hub-. ble constant ( H 0 ) across the sky which could be alternatively attributed to a ∼ 900 km s -1 cluster bulk flow extending up to ∼ 500 Mpc ( z ∼ 0 12). These results strongly disagree with the . isotropy assumption underlying the standard model of cosmology ( Λ CDM) (for a recent review, see Migkas (2024)).

Galaxy cluster T is a measure of the average kinetic energy of the electrons in the hot intra-cluster medium (ICM) and is measured using X-ray spectroscopy. However, there is a possibility that such T measurements could be a ff ected by systematic biases and, eventually, bias cosmological conclusions. Indeed, it has long been known that X-ray T measurements may be biased in various ways (see, e.g., Chapter 4 in the review by Reiprich et al. 2013). However, a directionality or spatial variation of such a bias is not expected. Possible reasons for directional-dependent systematic biases include under or overestimation of absorption corrections in T in certain directions. However, this e ff ect would likely be very small, as we use a range of 0 7 . -7 keV for the

T measurement. Another potential factor could be the presence of a hot, di ff use cloud that is not accounted for in the X-ray background. Regardless, if present, such a direction-dependent bias may lead to variations in best-fit parameters of scaling relations that depend on sky position. If the variations observed by M20 and M21 were indeed due to such systematic biases, a ∼ 10 -13% overestimation of T towards the most anisotropic direction identified by M21 could alleviate the tension between their findings and the Λ CDM model 1 . Therefore, it is essential to carefully investigate systematic T biases in the data used by M21 throughout the whole sky.

In this work, we look for potential T measurement biases across the sky by pairing the cluster T with the σ v since its measurement also does not involve any cosmological assumptions. We convert the variations obtained in the σ v -T relation's bestfit parameters into T over- or underestimations across the sky. σ v has completely di ff erent systematics as compared to T but is una ff ected by Galactic absorption. Therefore, if the latter is the reason for a possible T bias, this will, in principle, show up in this scaling relation's anisotropy. A significant variation in the best-fit normalisation across the sky that results in a noticeable overestimation of T in a particular region would suggest that systematic biases influence the results reported in M21. Conversely, if the best-fit parameters are uniform across the sky, it implies that a systematic T bias is a highly unlikely explanation for the observed anisotropies in M21.

We also explore cosmic isotropy by using σ v in combination with cluster properties, whose measurements rely on cosmological parameters. Two scaling relations, the X-ray luminosityvelocity dispersion ( L X -σ v ) and the total integrated Comptonization parameter-velocity dispersion ( Y SZ -σ v ), are studied for this purpose. The variations in the best-fit parameters of these relations are converted to variations in H 0 . This test provides new insights into a potential H 0 angular variation and acts as a crosscheck of the M21 results. For this work, a flat Λ CDM model is assumed with H 0 = 70 km s -1 Mpc -1 , Ω m = 0 3 and . ΩΛ = 0 7. .

The paper is structured as follows: Section 2 describes the various data samples and cluster properties used in this work. Section 3 explains the fitting procedure for scaling relations and the techniques used to study their variations across di ff erent parts of the sky. In Sect. 4, we present the general behaviour of the scaling relations used in this work. Section 5 provides detailed information on the variations in the best-fit parameters of the σ v -T relation, T bias across the sky, and the isotropic Monte Carlo simulation results. Section 6 explores the variations of L X -σ v and Y SZ -σ v scaling relations, H 0 variations across the sky from their joint analysis, and their comparison with isotropic Monte Carlo simulations. In Sect. 7, we discuss possible systematics and compare our results with M21. Finally, in Sect. 8, the conclusions of this work are given.

## 2. Data samples

This work utilises four key cluster parameters: T , σ v , L X, and Y SZ. These parameters are obtained for clusters present in three di ff erent data catalogues: the Meta Catalogue of X-ray detected Clusters of galaxies (MCXC, Pi ff aretti et al. 2011), the extremely expanded HIghest X-ray FLUx Galaxy Cluster Sample (eeHIFLUGCS, Pacaud et al., in preparation), and the Euclid velocity dispersion metacatalog (Euclid Collaboration: Melin et al. 2024, in preparation).

1 However, this scenario was tested in M20 and M21 and strongly disfavoured

MCXC is a comprehensive catalogue of X-ray-detected galaxy clusters, which is based on the ROSAT All-Sky Survey (RASS). The RASS mapped the entire sky in the X-ray band with the ROSAT PSPC detector in the energy range of 0 1 . -2 4 keV (Voges et al. 1999). The catalogue contains 1743 clus-. ters with virtually no duplicate entries. For each cluster, the catalogue provides the position, redshift, standardised 0 1 . -2 4 keV . band luminosity L X = L 500 , total mass M 500 2 , and radius R 500 3 .

eeHIFLUGCS is a homogeneously selected sample of galaxy clusters that has been created by imposing several selection criteria on the MCXC. The selection criteria include a flux limit on the unabsorbed X-ray flux of f 0 1 . -2 4 keV . ≥ 5 × 10 -12 erg s -1 cm -2 . Clusters in the Galactic plane ( | b | ≤ 20 ), ◦ the Magellanic clouds, and the Virgo cluster area are excluded. Another selection criterion used here is the availability of a high quality Chandra (Weisskopf et al. 2000) or XMM-Newton (Jansen et al. 2001) observation.

The Euclid velocity dispersion metacatalog is a collection of homogeneous velocity dispersion measurements from various sources compiled by the Euclid Collaboration. The catalogue includes 614 clusters from previous catalogues such as the Planck cluster sample (Aguado-Barahona et al. 2022), Abell (Girardi et al. 1998), ACO (Mazure et al. 1996), XDCP (Nastasi et al. 2014), SDSS Abell (Popesso et al. 2007), and other samples. /

## 2.1. X-ray luminosity L X and Temperature T

We utilise the L X from the eeHIFLUGCS catalogue, which contains the luminosity measurements from the RASS (Voges et al. 1999) covering the entire R 500 of a cluster. This is di ff erent from most of the XMM-Newton and Chandra clusters available in this catalogue. The luminosities of the MCXC have been corrected for the Galactic absorption based on the Willingale et al. (2013) N Htot values as described in M21.

The cluster T data is sourced from the eeHIFLUGCS catalogue, as provided in Migkas et al. (2020). The T values are measured within the range of 0 2 . -0 5 . R 500 for each cluster and are derived from observations made by Chandra and XMM-Newton . In total, T values are available for 313 eeHIFLUGCS clusters.

## 2.2. Total integrated Compton parameter Y SZ

The Sunyaev-Zeldovich (SZ) e ff ect provides information about the total gas pressure of the ICM and is characterised by the Compton parameter denoted by the symbol y . The integrated Compton parameter Y is obtained by integrating y over the cluster's solid angle. This quantity is multiplied by the square of the angular diameter distance to obtain the total integrated Compton parameter Y SZ, which has units of kpc , and scales with other 2 cluster properties (detailed description available in M21).

We determine Y SZ values for all MCXC clusters using the final data release of Planck Collaboration et al. (2020a) as explained in M21. We apply several signal-to-noise (S N) cuts to / the Y SZ measurements to see how our results depend on the Y SZ selection. For the default analysis, we use an S N cut of / ≥ 2, which results in 1093 clusters being included. Results in Sect. 6.5 show that no systematic trends are observed when we vary the S N cut. We are already aware of the existence of clusters in / these areas. Therefore, there is no necessity to increase the S N / cut for detection purposes.

2 M 500 obtained from L X -M 500 scaling relation.

3 The su ffi x 500 refers to the radius within which the average density is 500 times the critical density of the Universe.

## 2.3. Velocity dispersion σ v

The velocity dispersion of a galaxy cluster is a measure of the spread of velocities among the member galaxies of the cluster. It is obtained by measuring the line-of-sight velocities of the member galaxies in optical bands. We utilise the Euclid velocity dispersion metacatalog for the σ v measurements. This catalogue has several important properties: Firstly, all the clusters have at least ten member galaxies with spectroscopic redshift measurements, and interlopers have been adequately rejected while defining cluster members. Secondly, the velocity dispersion calculations are based on the methodology of Beers et al. (1990) with a S N / ≥ 4 and a minimum aperture of 0 5 Mpc. In addi-. tion to the velocity dispersion errors from the parent catalogues, dispersion uncertainties are also calculated in a standard way (Sereno &amp; Ettori (2015)).

## 2.4. Matching different catalogues

We cross-matched the σ v catalogue to eeHIFLUGCS and the Y SZ data (MCXC clusters) to study various scaling relations. Our primary matching criteria were that the cluster positions should be within 5 ′ of each other, and the di ff erence in redshifts ∆ z of the matching clusters should not exceed 0 01. We discovered . 165 matching clusters with both σ v and T measurements and 200 clusters with both L X and σ v measurements for the eeHIFLUGCS sample. We identified 284 matching clusters between the σ v catalogue and the Y SZ data (S N / ≥ 2).

## 3. Analysing the scaling relations

To study the three distinct scaling relations, a general form of the scaling relation is adopted. The scaling relation between two variables, Y and X , is expressed as

$$\stackrel { \sim } { \text{$are\,mea$} } } _ { \text{$\colon and\,are} } \quad \frac { Y } { C _ { Y } } \times E ( z ) ^ { \gamma _ { \gamma _ { X } } } = A _ { Y X } \times \left ( \frac { X } { C _ { X } } \right ) ^ { B _ { Y X } }.$$

Here CY and CX are the pivot points for the Y and X quantity respectively, E z ( ) is the redshift evolution factor GLYPH&lt;16&gt; E z ( ) = p Ω m(1 + z ) 3 + ΩΛ GLYPH&lt;17&gt; with γ YX being its power law index. AYX and BYX are the normalisation and slope of the scaling relation, respectively. The pivot points CY and CX are chosen close to the median of the Y and X , respectively. This is done to minimise the correlation between the best-fit parameters. These terms and the γ YX for all three scaling relations are mentioned in the Table 1.

## 3.1. Bayesian linear regression

Alinear regression is performed using Bayesian statistics to find the best-fit parameters for scaling relations in the logarithmic space. We constrain the best-fit parameters by maximising the posterior probability distribution of the parameters. The Markov Chain Monte Carlo (MCMC) method is used for sampling these distributions. We take its logarithm to convert the power law relation to a linear relation. The general form of scaling relations in terms of a linear model is given as

$$y = m \times x + c,$$

$$\wp e \text{ density } \quad \text{where } y \, = \, \log _ { 1 0 } \left ( \frac { Y } { C _ { Y } } \times E ( z ) ^ { \gamma _ { Y X } } \right ), \, x \, = \, \log _ { 1 0 } \left ( \frac { X } { C _ { X } } \right ), \, m \, = \, B _ { Y X }, \, \text{and}$$

Fig. 1. Sky distribution of the galaxy clusters used in the scaling relations σ v -T ( left ), L X -σ v ( middle ), and Y SZ -σ v ( right ). The clusters' redshifts are colour-coded.

![image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00726v1-figure-e5619023-3721-4f3c-bdcd-aa9af1992d8b.jpg)

To speed up the calculation process, the logarithm of the posterior probability distribution is used. The log-likelihood function is given as as the independent variable. The log-likelihood function for X Y | is given as

$$\log \left ( \mathcal { L } \right ) = - \frac { 1 } { 2 } \sum _ { i = 1 } ^ { N } \left ( \frac { [ y _ { i } - m x _ { i } - c ] ^ { 2 } } { \sigma _ { i } ^ { 2 } } + \log \left ( 2 \pi \sigma _ { i } ^ { 2 } \right ) \right ), \quad \quad ( 3 )$$

where N is the number of data points, yi and xi are the i th data points of Y and X respectively, and σ i is the total uncertainty given by:

$$\sigma _ { i } = \sqrt { \sigma _ { y, i } ^ { 2 } + m ^ { 2 } \sigma _ { x, i } ^ { 2 } + \sigma _ { i n t r } ^ { 2 } }. & & ( 4 ) & & \sigma \sigma _ { \cdot \cdot \cdot }$$

Here σ , y i and σ , x i are the uncertainties 4 for the yi and xi , respectively, and σ intr is the intrinsic scatter of the scaling relation measured in Y direction. Thus, the three parameters to be constrained are m c , , and σ intr . Flat priors are used for m and c with upper and lower bounds of + 10 and -10, respectively. Since σ intr is a positive definite quantity, a flat prior with a lower bound of 0 and an upper bound of + 10 is selected. Even though the prior choices are uninformative, the resulting posterior distribution always converges to a normal distribution.

The chain is initialised with a random set of parameters. It is run for at least 20,000 iterations in 4 chains to ensure it has converged 5 . The best-fit parameters are obtained by taking the median of the parameter space distribution. Lower and upper bounds are then determined by using the 16 th and 84 th percentile values of the distribution, respectively. The implementation of Bayesian linear regression is coded in Python 6 and the numba (Lam et al. 2015) package is used to improve computational time. The validity of the code is verified by comparing the results with the LINMIX (Kelly 2007) and PyMC (Abril-Pla et al. 2023) packages, and we find results that are consistent.

This method for Bayesian linear regression is Y X since the | y -axis distance of the data points from the best-fit line is minimised, and x is treated as the independent variable. For certain scaling relations, X Y best-fit implementation was used (more | details in Sect. 4). In this method, the x -axis distance of the data points from the best-fit line is minimised, and y is treated

4 Uncertainties in the linear space are converted to the logarithmic space using σ log 10 a = log 10 ( e ) × a + - -a 2 a . Here, a + and a -are the upper and lower limits of the parameter a , respectively.

5 The convergence is tested by discarding burn-in samples and checking the acceptance fraction.

6 The code is available at https://github.com/AdiPandya/Fast\_ Bayesian\_Regression

$$\log \left ( \mathcal { L } \right ) _ { X | Y } = - \frac { 1 } { 2 } \sum _ { i = 1 } ^ { N } \left ( \frac { \left [ x _ { i } - \left ( \frac { y _ { i } - c } { m } \right ) \right ] ^ { 2 } } { \sigma _ { i, X | Y } ^ { 2 } } + \log \left ( 2 \pi \sigma _ { i, X | Y } ^ { 2 } \right ) \right ), \quad \ \ ( 5 )$$

where σ, i XY is | q σ , 2 x i + GLYPH&lt;16&gt; σ , y i m GLYPH&lt;17&gt; 2 + σ intrX Y | 2 . The σ intrX Y | is the intrinsic scatter of the scaling relation measured in X direction. Note that the slope m and the normalisation (10 ) refer to the c best-fit parameters of the Y X form of the relation ( | L X -σ v ), and not to the X Y form ( | σ v -L X).

## 3.2. Removing outliers

We use an iterative 3 σ clipping method to detect and remove outliers in the data. This method removes outliers based on their residual distance from the best-fit line. When using the Y X fit-| ting method, we consider residuals in the Y direction, and for the X Y method, we consider residuals in the X direction. The | method assumes that residuals follow a normal distribution, and points lying outside the 3 σ of this distribution are removed as outliers. Using the 3 σ cut, we expect approximately one out of every 370 clusters to fall outside this range purely due to chance. Since the number of clusters in all relations is less than 370, we can safely use the 3 σ cut. This approach does not eliminate any clusters that adhere to normal scaling relation behaviour and helps eliminate problematic measurements.

It is important to note that this method may remove additional outliers when repeated because the new sample will have di ff erent best-fit parameters, resulting in di ff erent residual distributions. The method is repeated several times until no new outliers are found. Typically, this method removes 2-5 outliers from each scaling relation sample. Table 1 lists the final sample sizes for the scaling relations. Fig. 1 displays the sky distribution of the clusters used in the scaling relations σ v -T , L X -σ v , and Y SZ -σ v , along with their redshifts.

## 3.3. Scanning the sky

We study the consistency of the scaling relations in di ff erent sky directions by considering sections of the sky and calculating the best-fit parameters of the clusters inside those regions. We use the two-dimensional scanning method adopted from M21. In this method, a cone of radius θ is constructed for a given galactic longitude l and latitude b , and the scaling relation best-fit parameters of the clusters inside this cone are constrained. To cover the entire sky, the central longitudes and latitudes of the cones are

varied with a step size of 5 ◦ for their whole range. This results in 72 × 37 = 2664 di ff erent cones covering the entire sky 7 . Ideally, the cone sizes should be as small as possible to minimise overlap. However, the cone size should also be large enough to have su ffi cient clusters 8 inside the cone to e ffi ciently constrain the best-fit parameters. We tested various cone sizes and found a cone size of radius θ = 75 ◦ to be optimal for this analysis (additional details in Appendix A). An example of such a cone, along with the statistical weights of clusters inside this cone for the σ v -T sample, is shown in the bottom panel of Fig. 2.

Fig. 2. Example of a region used for 2-dimensional scanning centred at ( l , b ) = (60 ◦ , 15 ) with cone size, ◦ θ , of radius 75 ◦ for σ v -T relation. Di ff erent colours represent the statistical weights assigned to the clusters inside these regions. P represents the number of clusters inside this region.

![image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00726v1-figure-37d162c8-49b1-40ce-8c0f-48be7f8fa224.jpg)

We also apply statistical weights to clusters inside the scanning regions. This is done by dividing the statistical uncertainties of the observables ( σ y and σ x ) by the cosine of the angular separation of the cluster from the cone centre 9 . This results in a lower statistical weight for clusters that are further from the cone centre. 2-D maps are created to visualise the results of 2-D scanning. These maps are a grid plot with a box size of 5 ◦ on the Aito ff projection of the sky using the Python package desiutil 10 . The colour of each box in these maps represents the parameter value of interest for the cone centred at a given l and b .

## 3.4. Constraining H 0 angular variations

In the analysis of L X -σ v and Y SZ -σ v scaling relations, we assume that the intracluster physics remain unchanged for di ff erent directions, and variations of cosmological parameters cause the apparent anisotropy of scaling relations. According to this assumption, H 0 can vary across the sky while the true normalisation of the scaling relations remains constant. Alternatively, the A variations can also be interpreted in other ways, such as the presence of large-scale bulk flows, as shown in M21. However, due to the large scatter in the scaling relations, we restrict our analysis to constraining H 0 angular variations.

7 The final statistical directional uncertainties are much larger than the step size; thus, the step size of 5 ◦ is su ffi ciently small.

8 We ensure that there are at least 30 clusters in each cone.

9 Use of this method underestimates the intrinsic scatter and the statistical significance since the statistical uncertainties of the observables are increased. However, the isotropic MC simulations capture this information, and its significance is unbiased.

10 Package available at https://github.com/desihub/desiutil

We use best-fit A A / all of every cone to constrain H 0 angular variations using the relation

$$\underset { \text{thumb} } { \text{thumb} } \text{to} \\ \text{constrain} \quad H _ { 0 } = H _ { 7 0 } \times \left ( \frac { A } { A _ { \text{all} } } \right ) ^ { \frac { 1 } { 2 } }.$$

Here H 70 = 70 km s -1 Mpc -1 and A all is the best-fit A of the entire sample. A and H 0 are degenerate, and absolute constraints can only be put on the quantity A × H 2 0 . One of these parameters is assumed to be fixed to constrain the other (more details in Appendix F).

This method cannot be used to put absolute constraints on H 0 because a value of H 0 was assumed to calibrate the relation initially. The calculated value of H 0 expresses relative di ff erences between the regions; thus, the variations are always around the assumed value of H 0 .

## 3.5. Statistical significance of the variations

We use the best-fit parameters and their uncertainties in di ff erent parts of the sky to quote the statistical significance of the variations between the two regions. For two independent sub-samples i and j , the statistical significance of their deviation in terms of number of sigma is given by

$$\text{centred at} \quad \text{No. of } \sigma = \frac { p _ { i } - p _ { j } } { \sqrt { \sigma _ { p _ { i } } ^ { 2 } + \sigma _ { p _ { j } } ^ { 2 } } },$$

where p i and p j are the best-fit parameters of the two subsamples, and σ p i and σ p j are their uncertainties. For a cone centred at ( , l b ), all the clusters inside it are considered one subsample, and the clusters outside the cone are considered the other sub-sample. We find the number of sigma or significance values for each cone to create a significance (sigma) map. The colour of each box in this map represents the statistical significance of deviation compared to the rest of the sky at that location.

## 3.6. Isotropic Monte Carlo simulations

The statistical methods used above provide a way to study the anisotropies in the scaling relations. However, unknown statistical biases could still be present in the analysis methods, which could lead to underor overestimation of the significance of anisotropies. We perform isotropic Monte Carlo simulations to test the method's validity and provide robust significance. To create an isotropic sample for a scaling relation, we perform the following steps:

- 1. We start by fixing the cluster coordinates, redshifts, and σ v (including the uncertainties) to their actual values. This is done to encompass potential e ff ects that may give rise to anisotropies, including the spatial arrangement of real clusters.
- 2. The cluster T , L X, and Y SZ are simulated based on the best-fit parameters and scatter of the real data. First, their predicted values corresponding to the fixed σ v are calculated using the respective best-fit parameters.
- 3. Once all the points lie on the best-fit line, a random o ff set is added to these predicted values to obtain a simulated sample. This is done by drawing random values from a log-normal distribution centred at this value with a standard deviation equal to the total scatter (intrinsic + statistical) of the respective scaling relation and given cluster.

Fig. 3. From left to right: Best-fit plots for the scaling relation σ v -T (Y X), | L X -σ v (X Y), and | Y SZ -σ v (X Y). Shaded regions represent 1 | σ uncertainties of the best-fit parameters.

![image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00726v1-figure-785dd8f5-7b8b-47de-a67f-9eb4690b9a06.jpg)

![image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00726v1-figure-663d1660-f311-4ac5-bf87-0b5530924909.jpg)

![image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00726v1-figure-f153b18a-c38a-4359-8fb4-7fbf9e8f4af8.jpg)

Table 1. Best-fit parameters for the three scaling relations for the full sample. The table contains number of matching clusters ( N ), pivot points CY and CX , redshift evolution power γ of the relation, best-fit normalisation ( A ), best-fit slope ( B ) and the intrinsic scatter ( σ intr ). For the L X -σ v and Y SZ -σ v relations, results are shown using the X Y method. The intrinsic scatter mentioned here is measured in the direction in which the scatter | is minimised.

| Relation   |   N | C Y             | C X          | γ   | Method   | A                           | B                           | σ intr (dex)                |
|------------|-----|-----------------|--------------|-----|----------|-----------------------------|-----------------------------|-----------------------------|
| σ v - T    | 160 | 750 km s - 1    | 5 keV        | -   | YX |     | 1 . 016 + 0 . 014 - 0 . 014 | 0 . 394 + 0 . 025 - 0 . 026 | 0 . 068 + 0 . 005 - 0 . 005 |
| L X - σ v  | 195 | 10 44 erg s - 1 | 750 km s - 1 | - 1 | XY |     | 1 . 515 + 0 . 134 - 0 . 122 | 5 . 988 + 0 . 452 - 0 . 394 | 0 . 080 + 0 . 005 - 0 . 005 |
| Y SZ - σ v | 284 | 30 kpc 2        | 800 km s - 1 | + 1 | XY |     | 1 . 490 + 0 . 112 - 0 . 106 | 6 . 419 + 0 . 423 - 0 . 384 | 0 . 073 + 0 . 004 - 0 . 004 |

- 4. Step 3 is repeated 1000 times, each time adding a random o ff set to the predicted values to create 1000 simulated samples.

## 4.2. L X -σ v relation

We compare the maximum variations in scaling relations across the sky for these samples with those obtained from the observations to quantify the statistical significance of the observed di ff erences against random chance.

## 4. General behaviour of the three scaling relations

The three scaling relations used are σ v -T , L X -σ v , and Y SZ -σ v . This section presents best-fit parameters for these scaling relations for the entire sample. These relations are plotted in Fig. 3, and Table 1 shows an overview of the best-fit results.

## 4.1. σ v -T relation

To study the σ v -T relation, 160 galaxy clusters with measured velocity dispersion and temperatures are used. We obtain a bestfit slope of B = 0 394 . + 0 025 . -0 026 . , which is lower than the predicted self-similar value of 0 5 (Lovisari et al. 2021). Recent studies . (Xue &amp; Wu 2000; Ortiz-Gil et al. 2004; Nastasi et al. 2014; Wilson et al. 2016) obtain a slightly higher slope than the predicted self-similar scaling relation. One important thing to note is that most previous studies use the Orthogonal Distance Regression (ODR) method to fit the data. In contrast, Y X is used for this | analysis, which usually returns a flatter slope than ODR. Using the ODR method for our sample, a slope of B = 0 434 . + 0 030 . -0 029 . is obtained, similar to the results from Y X but slightly larger. |

195 clusters are used to study L X -σ v relation. The X Y method | is used to get best-fit parameters since the Y X method results in | residual trends when plotted against the cluster redshifts (refer to Appendix B for more details). However, as shown in Sect. 6.5, the choice of fitting method does not significantly influence the final results.

Using the X Y method we find | B = 5 988 . + 0 452 . -0 394 . while the Y X | leads to a shallower slope B = 2 722 . + 0 214 . -0 213 . . The slope from the YX method is similar to the self-similar slope | B self = 2 7 . -2 8 . when considering 0 1 . -2 4 keV luminosities for clusters (Lovis-. ari et al. 2021). Previous studies (Mahdavi &amp; Geller 2001; OrtizGil et al. 2004; Zhang et al. 2011; Nastasi et al. 2014; Sohn et al. 2019) using bolometric luminosities for clusters have obtained slopes consistent with the predicted slope of B self = 4. Most of these studies use the ODR method to fit the data. We obtain similar results using X Y and ODR methods. |

## 4.3. Y SZ -σ v relation

For the relation Y SZ -σ v , 284 clusters are used. These clusters are selected based on S N / ≥ 2 for the Y SZ parameter. The theoretical slope using the self-similar model is B self = 5 (derived from Lovisari et al. (2021)). This is di ff erent from what is observed as, using X Y, we find a steeper slope | B = 6 419 . + 0 423 . -0 384 . , and Y X leads | to a shallow slope B = 2 915 . + 0 194 . -0 195 . . These trends are similar to what was observed in the L X -σ v relation. Previous results from Rines et al. (2016) are quite similar to those obtained for X Y and | YXwith a slope of 5 68 | . ± 0 64 and 2 36 . . ± 0 30, respectively. .

Fig. 4. Left panel: Maps of best-fit normalisation ( top ) and the slope ( bottom ) compared to the full sample for the σ v -T relation. Note that the colour scale is di ff erent for both maps. This highlights the small variations in the A A / all map. Right panel: Significance maps for the normalisation ( top ) and slope ( bottom ). Both maps have the same colour scale ( -3 σ to + 3 σ ). A negative sigma value refers to a value lower than the rest of the sky.

![image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00726v1-figure-4d461a19-0d4d-457a-a159-d98b574f19b9.jpg)

## 5. Measuring the systematic temperature bias

tistical insignificance of the σ v -T angular variations is further confirmed in Sect. 5.3.

To perform systematic checks on the T measurements using σ v -T , we study the 2-D variations, and the best-fit results are converted to T over- or underestimation. The statistical significance of the analysis is checked using the isotropic Monte Carlo simulations.

## 5.1. 2-dimensional variations

The 2-D variations of best-fit parameters are studied by shifting a cone of radius 75 ◦ across the sky. The left panel of Fig. 4 shows maps of the best-fit A and B in di ff erent regions compared to the full sample ( A all and B all respectively). Maps for A A / all show values within a few per cent of 1, suggesting no substantial variations. To comment on the statistical significance of the variations, significance maps for A and B are created (right panel of Fig. 4). The significance maps for A show that the most deviating region is ( l , b ) = (195 ◦ , -20 ) with a significance of 2 27 ◦ . σ . This region deviates from the rest of the sky by 7 20 . ± 3 49%. .

The sigma maps for B show that the most deviating region is ( l , b ) = (65 ◦ , -25 ) with a significance of 2 42 ◦ . σ . The variations obtained in the best-fit slopes have a slightly higher significance, but the most deviating region is still well below 3 σ . These simple significance estimates already lead us to conclude that no significant variations are detected in the σ v -T relation. The sta-

## 5.2. Temperature bias

We perform the following steps to convert the variations in bestfit normalisation into T bias. To quantify the change in T for a given σ v = A × T B across the sky, we use

$$\ y \text{ shift} \text{-} \frac { A _ { 2 } } { A _ { 1 } } = \left ( \frac { T _ { 1 } } { T _ { 2 } } \right ) ^ { B }.$$

We take Region One as the region of interest and Region Two as the rest of the sky. The value of B is the best-fit slope for the rest of the sky and is assumed to be the same for both regions. The quantity ∆ T = T 1 -T 2 is the temperature bias between a region and the rest of the sky. We quantify this in terms of percentage bias using the relation

$$\stackrel { \dots } { \text{if $2.27\sigma$.} } \\ 3. 4 9 \%. \\ \text{region is}$$

A 2-D map of % T bias is created to identify regions with over- or underestimation of T (Fig. 5). Since ∆ T is calculated from A , the region of most bias remains the same. This region shows a (16 3 . ± 7 1)% underestimation in . T compared to the rest of the sky. There are notable T biases across the sky in the

Fig. 5. Map of percentage T bias between a region and rest of the sky. Negative T variations refer to an underestimation of T compared to the rest of the sky. To alleviate the observed anisotropies in M21, the behaviour required would have to be opposite to what is observed.

![image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00726v1-figure-9133ea6c-164a-48f7-8834-400d0641bb24.jpg)

2-D map, indicating anisotropy. However, these are also accompanied by high uncertainty and already the simple significance analyses above have shown they are not significant. Note that here, the region with the largest deviation from isotropy in M21 has an insignificant but nonetheless negative T bias. In contrast, a positive bias would be needed to alleviate the anisotropy in M21.

## 5.3. Isotropic Monte Carlo simulations

Isotropic Monte Carlo simulations are performed for the σ v -T scaling relation to check the validity of the methodology used to obtain the results from the 2-D analysis and to obtain improved significance estimates. Isotropic samples are created using the method mentioned in Sect. 3.6. For these simulated samples, significance maps are created using the same procedure as the real data, and the maximum significance value is noted. We also take note of its corresponding percentage T bias to understand the expected T bias solely due to the statistical scatter of the relation. It is important to note that an isotropic sample with a large scatter around the best-fit line will exhibit a T bias. We can assess the statistical significance of the obtained bias by comparing the expected T bias from the scatter and the real data. The distribution of these quantities obtained from the simulations is shown in Fig. 6. Since the maximum variation in the observed data corresponds to a negative T bias, we plot only the negative T bias values from the simulations for their comparison.

The figure shows that the probability of obtaining a sigma value of 2 27 . σ or higher is 72% ( p = 0 72), and . ∼ 18% bias in T is expected purely due to the scatter of the relations. The 2-D distribution of these simulated results demonstrates that real data lies well within the 68% confidence interval, indicating a high probability of such results occurring due to statistical fluctuations.

These results further confirm that the variations in the bestfit normalisation (and, therefore, T biases) are not statistically significant.

## 6. Probing cosmic anisotropies

In this section, we present the constraints of the apparent H 0 anisotropies using the L X -σ v and Y SZ -σ v relations. We cross-

Fig. 6. Distribution of maximum sigma values and its % T bias obtained from the 1000 isotropic Monte Carlo simulated samples of the σ v -T relation. Shaded contours represent 68%, 95%, and 99% confidence levels. The value obtained in the real data is shown by a black dashed line on the histograms and by a point on the contour plot. The p -value in the histogram of maximum sigma values represents the probability of getting higher significance than the real data.

![image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00726v1-figure-f4fef524-67c5-406b-9249-d701fd76261d.jpg)

check the results of M21 with an independent cluster property that is di ff erent from that used in the original study.

## 6.1. L X -σ v relation

Using 2-D scanning, we find A variations across the sky and compare them with A of the full sample ( A all ). The 2-D maps of A A / all along with significance maps are shown in Fig. 7. The sigma maps for A show that the most deviating region is ( l , b ) = (225 ◦ ± 41 ◦ , -25 ◦ ± 41 ) with a significance of ◦ -3 66 . σ . Some variations in the best-fit slopes are also observed in certain regions but at a lower significance. However, slope variations are not a focal point of this test as they do not strongly a ff ect the inferred H 0 variations (no strong correlation exists between slope and H 0 constraints).

Fig. 7. Map of best-fit A compared to the full sample ( left ) and its significance map ( right ) for the L X -σ v relation.

![image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00726v1-figure-c42fc5d9-5793-41c4-bb74-69e9c7e4e948.jpg)

The H 0 angular variation map obtained from the L X -σ v relation is shown in Fig. 8. The region exhibiting the most statistically significant deviation remains consistent with the A maps. In terms of H 0 , this region di ff ers from the rest of the sky by ∆ H 0 = 28 3 . ± 6 6%. .

Fig. 8. H 0 angular variation map created from the L X -σ v relation.

![image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00726v1-figure-4fb88a90-c834-4043-b8a2-6044b7092e2b.jpg)

In order to better understand the variation in data points (i.e., individual clusters) between the region of maximum anisotropy ( l = 225 ◦ , b = -25 ) and the rest of the sky, we created a distri-◦ bution of H 0 values corresponding to the data points in these regions and compared them with each other (see Fig. 9). We assigned a value of H 0 to each data point based on its vertical distance from the best-fit line. Points on the line were assigned H 0 = 70 km s -1 Mpc -1 . The H 0 value for a point with a y value of yi and the best-fit y = y fit is calculated using the formula H 0 , i = 70 × GLYPH&lt;16&gt; yi y fit GLYPH&lt;17&gt; 1 2 . On average, data points in the region ( l = 225 ◦ , b = -25 ) tend to have lower ◦ H 0 and thus lie below the best-fit line. The rest of the sky contains a larger number of data points that lie above the best-fit line, as observed from the distribution tail.

Fig. 9. Distribution of H 0 values corresponding to each data point in the region of maximum anisotropy and rest of the sky for the relations L X -σ v ( top ) and Y SZ -σ v ( bottom ). The median of both distributions is shown by dashed lines along with H 0 = 70 km s -1 Mpc -1 . Note that both regions have di ff erent numbers of clusters, and thus, to compare both, the density of the distribution is plotted.

![image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00726v1-figure-679f1ff7-eb94-4752-b721-6f900fd6319f.jpg)

## 6.2. Y SZ -σ v relation

Fig. 10 shows the maps of best-fit A of di ff erent regions compared to the full sample and their sigma maps. The direction showing the largest A deviation is ( l , b ) ∼ (295 ◦ ± 58 ◦ , -35 ◦ ± 58 ) ◦ with a statistical significance of -4 13 . σ . This is in a similar direction to the most deviating region in the L X -σ v analysis.

Fig. 10. Map of best-fit A compared to the full sample ( left ) and its significance map ( right ) for the Y SZ -σ v relation. Both maps have the same colour scale as Fig. 7

![image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00726v1-figure-0dae2f00-f6da-4ff7-8d21-14d5b4599f07.jpg)

One interesting feature of this map is that, on average, the Northern Galactic Hemisphere shows higher A values than the Southern one. The best-fit A of all the clusters in the two Galactic Hemispheres di ff ers by 3 35 . σ . We were unable to identify a systematic bias causing this phenomenon; therefore, it is plausible that this discrepancy may be a random occurrence. We test the possibility of assuming an incorrect redshift evolution factor and its e ff ects on detecting anisotropy and this feature in Sect. 7.4

The H 0 angular variation maps that are derived from the Y SZ -σ v relation are shown in the Fig. 11. The region exhibiting the most statistically significant deviation di ff ers from the rest of the sky by ∆ H 0 = 28 5 . ± 6 4%. The . H 0 distributions in the region of maximum anisotropy ( l = 295 ◦ , b = -35 ) and the rest of ◦ the sky (Fig. 9) reveal a similar distribution to the ones obtained from L X -σ v analysis.

Fig. 11. H 0 angular variation map created from the Y SZ -σ v relation. The map has the same colour scale as Fig. 8

![image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00726v1-figure-f36e3b72-73d2-4f3f-a2de-e9b9e0ba014c.jpg)

## 6.3. Joint analysis of L X -σ v and Y SZ -σ v relations

The apparent anisotropy information of the L X -σ v and Y SZ -σ v scaling relations can be combined into one map that shows H 0 angular variations. By multiplying the H 0 posterior probability distributions of the two relations, the combined H 0 posterior probability distribution can be obtained. By using the combined

posterior probability distribution to constrain the best-fit H 0 for each cone, we generate a combined map showing the percentage change in H 0 and its significance.

## 6.3.1. Limitations of joint analysis

There are a couple of things to note here that could overestimate the significance of the detected anisotropies. First, we assume there is no bias in σ v measurements. Since it is used in both the scaling relations, a bias could overestimate the significance. Second, this method of combining the likelihoods assumes there is no correlation between the scatter of L X and Y SZ at fixed σ v .

Previous studies like Nagarajan et al. (2019) have shown a positive correlation between the scatter. We have 180 clusters with both L X and Y SZ measurements, and we use them to find the correlation between their scatter from their respective scaling relation. Fig. 12 shows a plot of residuals of the two relations along with their best-fit line. The correlation coe ffi cient between the residuals is 0 957, and the best-fit slope of 1 01 . . ± 0 04. The strong . correlation could be the result of a large scatter uncertainty in / σ v , causing the cluster to shift in the same direction for both the L X -σ v and Y SZ -σ v relations.

Fig. 12. Correlation between the residuals of the relations L X -σ v and Y SZ -σ v . The best-fit line is shown in orange, along with the correlation coe ffi cient.

![image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00726v1-figure-fd738140-497c-4d6b-8d87-558a9d02e892.jpg)

Despite a correlation, ∼ 37% clusters present in Y SZ -σ v are unique to this relation and are completely independent of the L X -σ v relation.

## 6.3.2. Apparent H 0 anisotropy from joint analysis

The joint H 0 map and its significance maps are shown in Fig. 13. The region with the most deviation is found to be in the direction of ( l , b ) ∼ (295 ◦ ± 71 ◦ , -30 ◦ ± 71 ) with ◦ ∆ H 0 = 27 6 . ± 4 4% and . a statistical significance of -5 27 . σ compared to the rest of the sky. Notice that the significance of the variation has increased compared to the individual analysis, which is a result of smaller H 0 uncertainties. The strong correlation between the residuals suggests that the two relations are clearly not independent and that simply combining the likelihoods will undoubtedly result in an overestimation of the significance.

Fig. 13. H 0 angular variation maps created from the joint analysis of relations L X -σ v and Y SZ -σ v ( top ) and its significance map ( bottom ). Colour scales for the H 0 angular variation map are the same as before, but the colour scales for the sigma map have been increased to account for increased significance.

![image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00726v1-figure-01225eb4-1912-4750-b6d3-93e11f8b28bb.jpg)

## 6.4. Isotropic Monte Carlo simulations

We perform isotropic Monte Carlo simulations by considering the residual correlation to calculate realistic significance. The method for creating an isotropic MC simulated sample is described in the Sect. 3.6.

First, we generate simulated L X and Y SZ values based on the best-fit line, taking into account the total scatter of their respective relations. However, simulated Y SZ values are calculated only for clusters that are unique in the Y SZ -σ v dataset. For the 180 clusters that have both L X and Y SZ values, we obtain simulated Y SZ by taking into account the correlated scatter of L X and Y SZ with σ v . To do so, we start by determining the residuals of Y SZ that correspond to residuals of the simulated L X using the correlation between the two quantities (Sect. 6.3.1). These residuals are then used to predict Y SZ by applying the best-fit parameters of the Y SZ -σ v relation. Finally, a random o ff set is added to these values, which is drawn from a log-normal distribution centred at

$$\int \infty _ { \substack {. 4 \% \text{ and } \quad Y _ { S Z, \text{pred} \\ \dots \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \colon \\ \int \infty _ { \substack {.. \dots \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \dots \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \colon \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \ \end{array} \lambda }$$

This is done 1000 times to create 1000 simulated L X and Y SZ samples. The distribution of maximum sigma values and ∆ H 0 variations obtained from the simulations are compared with the real data in Fig. 14 11 . We examine the isotropic sample's ∆ H 0 % corresponding to maximum sigma values to know the expected

11 Absolute values for both these quantities are considered.

Fig. 14. Distribution of maximum sigma values and its ∆ H 0 variations obtained from the 1000 isotropic Monte Carlo simulated samples of the L X -σ v ( left ) and Y SZ -σ v ( right ) relations. Shaded contours represent 68%, 95%, and 99% confidence levels. The value obtained in the real data is shown by a black dashed line on the histograms and by a point on the contour plot. The p -value in the histogram of maximum sigma values represents the probability of getting higher significance than the real data.

![image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00726v1-figure-f47d019e-53b0-48b8-967b-0c427956fac7.jpg)

∆ H 0 caused by the spread in the scaling relation. A large scatter in the scaling relation leads to significant ∆ H 0 when considering its sub-samples. Therefore, by comparing ∆ H 0 obtained from isotropic samples and the real data, we can assess the significance of these variations given the scatter in the relation.

have higher anisotropy in both relations. Thus, the probability is given as p = 0 001. .

For the L X -σ v , 20 out of the 1000 simulated samples have higher maximum deviation than the real data. Thus, using this method, there is 2% probability ( p = 0 02) . of observing a ≥ 3 66 . σ anisotropy in an isotropic Universe. The observed variations in H 0 are higher than the expected ∆ H 0 of 22 0 . + 4 3 . -4 2 . %due to scatter in the relation. The 2-D distribution of these quantities shows that the real data lies within the 95% confidence level. This suggests that significance estimates from Eq. 7 lead to an overestimation of significance due to the large scatter in the scaling relation.

For the Y SZ -σ v relation, only 11 out of the 1000 simulated samples had higher maximum deviation than the real data. This gives a probability of 1 1% ( . p = 0 011) of observing a . ≥ 4 13 . σ anisotropy in an isotropic Universe. The expected ∆ H 0 due to scatter in the relation is 19 2 . + 3 7 . -3 5 . %, and the real data lies between the 95% and 99% confidence levels in the 2-D distribution of these quantities. Again, significance estimates from Eq. 7 lead to overestimating the significance. However, due to lower scatter, the ∆ H 0 for the Y SZ -σ v relation is more significant than that obtained from the L X -σ v relation.

If the two relations were independent, the joint probability would be the product of the two probabilities. However, as the σ v data in both the relations are the same and their residuals are correlated, the two relations are not entirely independent. The joint probability of observing L X -σ v and Y SZ -σ v anisotropies simultaneously is calculated by noting down number of simulated samples that simultaneously satisfy the criteria of having ≥ 3 66 . σ and ≥ 4 13 . σ anisotropies in their respective relations. Using this method, only 1 out of 1000 samples simultaneously

In real data, the anisotropies in the two relations are separated by 60 2 . From the simulated samples, roughly 27 6% . ◦ . samples show lower separation in detected anisotropies than the real data. The probability of observing higher anisotropies with lower separation is obtained by multiplying the probabilities of the two events. Thus, the probability is p = 2 76 . × 10 -4 and it corresponds to a Gaussian σ of 3 64 . σ .

The results of the simulations suggest that even in isotropic samples, some variation in ∆ H 0 is expected due to the scatter of the relations being used, which biases the measured ∆ H 0 . It is essential to correct this bias before determining the actual ∆ H 0 . As a result, the observed ∆ H 0 in M21 is also somewhat overestimated. However, this overestimation would not a ff ect the statistical significance of the ∆ H 0 they detected (at a level of 5 4 . σ ) because this e ff ect is considered in the Monte Carlo simulations they performed.

A summary of results obtained from the relations L X -σ v and Y SZ -σ v along with their joint analysis is presented in Table 2.

## 6.5. Comparison with different data samples and best-fit methods

We compared our results with di ff erent data sets to check for apparent cosmic anisotropy. For L X, we used the MCXC catalogue data and compared it with the results obtained from the eeHIFLUGCS catalogue. For the Y SZ data, we compared the results with di ff erent S N cuts such as S N / / ≥ 3 and S N / ≥ 4 5. . The latter is chosen since it is the S N cut chosen in the o / ffi -cial PLANCKSZ2 cluster catalogue (Planck Collaboration et al. 2016). To check if the results are consistent between the fitting methods, the results from both the X Y and Y X methods are | |

Table 2. Maximum anisotropy direction and its absolute amplitude for the L X -σ v and Y SZ -σ v relations along with absolute H 0 variation compared to the rest of the sky. Results from isotropic MC simulations are shown in terms of the p -value, their corresponding Gaussian σ values, and the expected H 0 variation due to scatter in the relation.

| Scaling relation   | Max. anisotropy direction ( l , b )   | H 0 variation (%)   | Maximum significance ( σ )   | Isotropic MC p -value ( σ )     | Expected H 0 variation (%)   |
|--------------------|---------------------------------------|---------------------|------------------------------|---------------------------------|------------------------------|
| L X - σ v          | (225 ◦ ± 41 ◦ , - 25 ◦ ± 41 ◦ )       | 28 . 3 ± 6 . 6%     | 3 . 66 σ                     | p = 0 . 020 (2 . 33 σ )         | 22 . 0 + 4 . 3 - 4 . 2 %     |
| Y SZ - σ v         | (295 ◦ ± 58 ◦ , - 35 ◦ ± 58 ◦ )       | 28 . 5 ± 6 . 4%     | 4 . 13 σ                     | p = 0 . 011 (2 . 54 σ )         | 19 . 2 + 3 . 7 - 3 . 5 %     |
| Joint results      | (295 ◦ ± 71 ◦ , - 30 ◦ ± 71 ◦ )       | 27 . 6 ± 4 . 4%     | 5 . 27 σ                     | p = 2 . 76 × 10 - 4 (3 . 64 σ ) | 19 . 7 + 2 . 9 - 2 . 5 %     |

Table 3. Maximum anisotropy direction for the L X -σ v and Y SZ -σ v relations along with H 0 angular variation compared to the rest of the sky and the absolute significance of the variations. Results are shown for the X Y and Y X methods and di | | ff erent datasets for the two relations.

| Scaling relation   | Dataset used   | Best-fit method   | Max. anisotropy direction ( l , b )   | H 0 variation (%)                | Maximum significance ( σ )   |
|--------------------|----------------|-------------------|---------------------------------------|----------------------------------|------------------------------|
| L X - σ v          | eeHIFLUGCS     | XY | YX |         | (225 ◦ , - 25 ◦ ) (305 ◦ , - 20 ◦ )   | 28 . 3 ± 6 . 6% 22 . 1 ± 5 . 2%  | 3 . 66 σ 3 . 79 σ            |
| L X - σ v          | MCXC           | XY | YX |         | (280 ◦ , - 25 ◦ ) (300 ◦ , - 30 ◦ )   | 27 . 0 ± 6 . 2% 21 . 8 ± 4 . 6%  | 3 . 71 σ 4 . 13 σ            |
| Y SZ - v           | S / N ≥ 2      | XY | YX |         | (295 ◦ , - 35 ◦ ) (310 ◦ , - 35 ◦ )   | 28 . 5 ± 6 . 4% 21 . 8 ± 4 . 6%  | 4 . 13 σ 4 . 51 σ            |
| σ                  | S / N ≥ 3      | XY | YX |         | (185 ◦ , - 40 ◦ ) (315 ◦ , - 30 ◦ )   | 23 . 2 ± 5 . 1% 21 . 1 ± 4 . 9%  | 3 . 86 σ 4 . 19 σ            |
| Y SZ - v           | S / N ≥ 4 . 5  | XY | YX |         | (185 ◦ , - 35 ◦ ) (185 ◦ , - 40 ◦ )   | 22 . 5 ± 6 . 1% 17 . 9 ± 4 . 41% | 3 . 24 σ 3 . 66 σ            |

presented. The general behaviour of these datasets for di ff erent best-fit analysis methods is shown in Table C.1.

Table 3 shows the direction of maximum anisotropy, its amplitude, and ∆ H 0 %for these datasets using both fitting methods. The results from di ff erent fitting methods and datasets return consistent results with each other. These results indicate that the general direction of the most deviating regions is similar for different datasets and analysis methods. The direction of anisotropy is found in the general direction of l &gt; 180 ◦ and b &lt; 0 ◦ for all datasets. The Y X fitting method consistently yields less angular | H 0 variation than the X Y fitting method, possibly due to less | scatter in the Y direction. The significances are calculated using the Eq. 7 and are mentioned only for a consistency check. As shown in Sect. 6.4, these significances are overestimated due to a large scatter in the relation and isotropic MC simulations are required to quantify the true statistical significance of the observed anisotropy.

![image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00726v1-figure-97910636-d277-4a4e-9415-2bd390ee2f2b.jpg)

## 7. Discussion

7.1. Comparison with M21

## 7.1.1. Probability of T bias

In the work done by M21, the most deviating region found for the relations L X -T and Y SZ -T are ( l = 270 , ◦ b = -9 ) and ◦ ( l = 268 , ◦ b = -16 ) respectively. If the anisotropies in these ◦ relations are due to T biases, then for a fixed L X or Y SZ, T in these regions should be biased high compared to the rest of the sky (10 6 . ± 4 6% and 13 2 . . ± 4 3% respectively). These are op-. posite to the results obtained in this work since we obtain T bias of ∼ -9 3% and . ∼ -13 1% in the same regions. For ease of . notation, the regions of maximum anisotropy for both relations are referred to as the RLT and RYT regions, respectively. Their T

Fig. 15. Distribution of T bias in the RLT ( top ) and RYT ( bottom ) regions for 1 million simulated samples. The vertical black dashed line indicates the T bias needed to account for the results of M21 ( ∆ TLT and ∆ TYT ), and the shaded region represents its 1 σ uncertainties. The p -value displayed next to it shows the probability of obtaining samples with a higher T bias than ∆ TLT and ∆ TYT , respectively.

![image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00726v1-figure-82c7a61b-1cd4-4287-b5a0-2ad441ddcb50.jpg)

biases are referred to as ∆ TLT and ∆ TYT , respectively. We perform additional simulations to quantify how strongly the results

from the analysis of σ v -T confirm or reject the possibility of a systematic bias in the cluster T measurements.

We draw random values of T from a log-normal distribution centred at the real T measurement, and its standard deviation is given as the total scatter of the scaling relation. These simulations are similar to the isotropic simulations with the main di ff erence being that the T are drawn around their real values instead of those predicted from the scaling relation. The T bias for all simulated samples from the most deviating regions RLT and RYT are noted and are compared with ∆ TLT and ∆ TYT . Fig. 15 shows the distribution of T bias in these regions obtained from one million simulated samples. Negative T bias refers to an underestimation in T compared to the rest of the sky.

The probability that the observed L X -T and Y SZ -T anisotropies in M21 are caused by an overestimation of T is 2 2 . × 10 -5 (4 24 . σ ) and &lt; 10 -6 ( &gt; 4 89 . σ ) respectively. The presence of a T decrease in our work - albeit non-significant - strongly disfavours a T bias as the cause of the M21 observed anisotropies.

## 7.1.2. Joint result of cosmic anisotropy with M21

The M21 study found a (9 0 . ± 1 7)% angular variation in . H 0 in the direction of ( l , b ) = GLYPH&lt;16&gt; 280 ◦ + 35 ◦ -35 ◦ , -15 ◦ + 20 ◦ -20 ◦ GLYPH&lt;17&gt; . There is a strong agreement in the direction of maximum anisotropy between these findings and the results from L X -σ v and Y SZ -σ v analysis. One of the main distinctions between the two studies is the (27 6 . ± 4 4)% angular variation in . H 0 from the joint analysis in this work, which is significantly higher than the variation found by M21 due to a large scatter.

Fig. 16. Map of H 0 angular variations ( left ) and its significance ( right ) for the joint analysis of relations L X -σ v &amp; Y SZ -T ( top ) and Y SZ -σ v &amp; L X -T ( bottom ).

![image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00726v1-figure-8005e385-baab-4a00-83be-cfbdecc04356.jpg)

We perform a joint analysis with the results of M21 by pairing the datasets that are completely independent, i.e., pairing L X -σ v and Y SZ -σ v with M21's Y SZ -T and L X -T respectively. By multiplying the posteriors of H 0 for both relations, we create maps showing the H 0 angular variation and its significance (see Fig. 16). The maximum H 0 angular variation decreased to ∼ 11% in both joint results, which is similar to the results of M21. The maximum significance in the joint analysis of L X -σ v and Y SZ -T remains the same when compared to the results of L X -σ v analysis. However, there is a decrease in significance in the joint analysis of other relations when compared to the results of Y SZ -σ v analysis. A common trend in both maps is that the results of M21 dominate them due to having a low scatter in relations compared to our results.

We also check for any correlation between the residuals of the two relations that are being used in the joint analysis. We do not find any correlation between the L X -σ v and Y SZ -T relations. We can, therefore, treat the two relations as independent while analysing their joint results. However, we find a correlation of ρ = -0 66 between the residuals of the . Y SZ -σ v and L X -T relations. Thus, the joint result of these relations could underestimate the significance.

## 7.2. Isotropic samples based on random cluster positions

Aside from using isotropic MC samples, we employ a commonly used method for generating isotropic samples, which involves shu ffl ing the cluster positions randomly and calculating the probability of obtaining higher anisotropy in these samples compared to the observed data. We create 1000 such isotropic samples for both the relations, find their maximum significance values, and compare their distributions to the observed data (see Fig. 17). Our results indicate that the probability of observing higher significance in an isotropic sample (based on random cluster positions) compared to the real data for the L X -σ v and Y SZ -σ v relations are ∼ 5% and ∼ 1% respectively.

Fig. 17. Distribution of maximum significance obtained from the 1000 isotropic samples created from the random shu ffl ing of cluster positions for the L X -σ v ( top ), and Y SZ -σ v ( bottom ) relations. The dashed line represents the maximum significance obtained in the observed data, and the p -value next to it shows the probability of obtaining higher significance than the observed data.

![image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00726v1-figure-5b6f4161-ecf4-43ca-87e7-80fe63886183.jpg)

When we compare these to the isotropic MC simulation results in Sect. 6.4, we find that the probability remains unchanged for the Y SZ -σ v relation. However, for the L X -σ v relation, we note an increase in the p -value compared to isotropic MC simulation results. This decrease in significance can be attributed to our use of large cone sizes for the analysis. When we shu ffl e the cluster positions, it is likely that, at times, the clusters a ff ected by the anisotropy might fall within the same cone or close to one an-

other. As a result, the mock sample exhibits higher significance variations.

## 7.3. Effects of best-fit slope variations

For all three relations, we treat the slope as a free parameter, and we find some variations in best-fit slopes across the sky (see Fig. 4 for σ v -T and Appendix C for L X -σ v and Y SZ -T ). We do not expect a significant correlation between the best-fit parameters because the pivot points are chosen close to the data median. However, some cones may have a strong correlation that could bias the observed results. To examine this, we calculate the correlation between A and B in each region and create a 2-D map based on this (Fig. 18). For σ v -T , the correlation between the best-fit parameters in all the regions is less than ± 0 4. Similar . results are also found for the Y SZ -σ v relation. In the analysis of L X -σ v , we observe certain regions with a moderate correlation of + 0 5, but these regions are distant from the anisotropy region, . so we don't have any reason to suspect a bias due to parameter correlation.

Fig. 18. 2-D Maps of the correlation coe ffi cient between the best-fit parameters for the σ v -T ( top left ), L X -σ v ( top right ), and Y SZ -σ v ( bottom ) relations. All maps have the same colour scale ( -0 5 to . + 0 5). .

![image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00726v1-figure-e3678cbe-fd01-45dd-8d02-59022fca27c0.jpg)

We test our findings with the slope fixed to that of the full data to see the impact of the slight correlation between parameters (refer to Appendix C for more details). We find no significant di ff erence in the results when the slope is fixed as compared to the free slope analysis.

## 7.4. Anisotropy at different z scales and with different z evolutions

We examine the distribution of redshift across di ff erent regions of the sky. If any region shows significantly higher or lower average redshift ( z ) values, it could indicate that we are comparing di ff erent scales, which is not ideal. We create the redshift distribution for clusters present in the region of maximum anisotropy and compare it with the rest of the sky for the L X -σ v and Y SZ -σ v relations (see Fig. E.1). We find that, on average, the region with maximum anisotropy has lower z , whereas the rest of the sky has a large population of clusters at higher z . To ensure that our results are not biased by the population of clusters at di ff erent z scales, we apply various z cuts and compare these results with our initial findings.

We apply several lower cuts on z , such as z &gt; 0 03, . z &gt; 0 05, . and z &gt; 0 07 and an uppercut of . z &lt; 0 15 to our data individually, . and repeat our analysis. Redshift cuts higher than z &gt; 0 07 are . not considered since we are limited by the number of clusters in di ff erent cones. The L X -σ v relation has fewer clusters, and several cones for the cut of z &gt; 0 07 had fewer than 20 clusters, . and thus, we do not include this cut for the relation.

Fig. 19. Positions and amplitude of maximum anisotropy detected for several z cuts and z evolution power for L X -σ v ( top ), and Y SZ -σ v ( bottom ) relations. The shaded regions represent the directional uncertainties of the full sample.

![image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00726v1-figure-4febaf97-c6fb-4395-ba9e-abaebcc5e692.jpg)

We plot the positions and amplitude of regions with maximum anisotropy obtained from these cuts for both relations, including the results obtained with the full dataset in Fig. 19. These cuts return consistent results with the full sample, and we see that our results are not strongly biased due to cluster populations at di ff erent z scales. Along with these results, we also plot the results with di ff erent z evolution power ( γ ).

We observe distinct trends in the 2-D map of median z used within each cone, with the northern galactic sky showing clusters at higher z compared to the southern half (see Fig. E.2). A similar galactic north-south divide is also present in the 2-D maps of both relations (see Fig. 7 and Fig. 10), due to the northern galactic sky having higher A compared to the southern half. Regions with higher median z are strongly a ff ected by incorrect redshift evolution assumptions. Therefore, it is crucial to test if the assumed γ for the two relations results in a bias and the observed galactic north-south divide in our significance maps. As mentioned in the Table 1, γ for the L X -σ v and Y SZ -σ v relations are -1 and + 1 respectively. We change γ to -2 and + 2, respectively, to see the e ff ect of choosing a strong redshift evolution on the anisotropy detected.

The results for L X -σ v indicate almost no di ff erence in significance and in the direction of anisotropy. In the Y SZ -σ v relations, we observe minor di ff erences from γ = + 1 due to clusters with higher z in the data sample. Even when we consider the most extreme value of γ for the two relations ( -10 and + 10, respectively), we are unable to eliminate the north-south divide in the significance maps. Therefore, we can conclude that our results are not biased due to incorrect z evolution.

## 8. Summary

In this work, we investigate the isotropy of the local universe using galaxy cluster scaling relations between cosmologydependent cluster properties and a cosmology-independent property. We utilise a cosmology-independent variable that has not been used before: the velocity dispersion of a galaxy cluster. Previous studies by M21 used cluster T as their cosmologyindependent quantity and found an apparent spatial variation of approximately 9% in the Hubble constant, H 0 , across the sky.

We examined the σ v -T relation across the sky to check if a position-dependent systematic bias of T measurements causes an overestimation of apparent H 0 variations in the work of M21. Weobtain no significant anisotropies across the sky in the σ v -T relation. The probability of obtaining the observed variations for σ v -T relation in an isotropic Universe is 0 72. The region with . the most variations in σ v -T relation is found in a similar direction as M21 but with a T underestimation of (16 3 . ± 7 1)%, . hinting at the possibility that the significance of the M21 results might even be underestimated. A positive T bias needed to explain the results of M21 can be rejected at a high probability ( p ∼ 10 -5 ).

We utilised the L X -σ v and the Y SZ -σ v relations to probe the (an)isotropy of the universe. From the joint analysis, we find the most anisotropic region in the direction of (295 ◦ ± 71 ◦ , -30 ◦ ± 71 ) at 3 64 ◦ . σ . The statistical significance is obtained using the isotropic Monte Carlo simulations. The direction of maximum anisotropy is similar to the results of M21 ( l , b ) = GLYPH&lt;16&gt; 280 ◦ + 35 ◦ -35 ◦ , -15 ◦ + 20 ◦ -20 ◦ GLYPH&lt;17&gt; . These results are in disagreement with the standard cosmological model. The results obtained from both these analyses further strengthen the results of M21.

Acknowledgements. L.L. acknowledges the financial contribution from the INAF grant 1.05.12.04.01.

## References

Abril-Pla, O., Andreani, V., Carroll, C., et al. 2023, PeerJ Computer Science, 9, e1516

- Aguado-Barahona, A., Rubiño-Martín, J. A., Ferragamo, A., et al. 2022, Astronomy and Astrophysics, 659, A126
- Andrade, U., Bengaly, C. A. P., Alcaniz, J. S., &amp; Capozziello, S. 2019, MNRAS, 490, 4481
- Appleby, S., Shafieloo, A., &amp; Johnson, A. 2015, ApJ, 801, 76

Beers, T. C., Flynn, K., &amp; Gebhardt, K. 1990, Astronomical Journal, 100, 32

- Bennett, C. L., Kogut, A., Hinshaw, G., et al. 1994, The Astrophysical Journal, 436, 423
- Bennett, C. L., Larson, D., Weiland, J. L., et al. 2013, The Astrophysical Journal, 208, 20

Colin, J., Mohayaee, R., Rameez, M., &amp; Sarkar, S. 2019, A&amp;A, 631, L13 Giodini, S., Lovisari, L., Pointecouteau, E., et al. 2013, Space Sci. Rev., 177, 247 Girardi, M., Giuricin, G., Mardirossian, F., Mezzetti, M., &amp; Boschin, W. 1998, The Astrophysical Journal, 505, 74

- Jansen, F., Lumb, D., Altieri, B., et al. 2001, Astronomy and Astrophysics, 365, L1

Javanmardi, B. &amp; Kroupa, P. 2017, A&amp;A, 597, A120

Kaiser, N. 1986, J Monthly Notices of the Royal Astronomical Society, 222, 323 Kelly, B. C. 2007, Astrophysical Journal, 665, 1489

- Lam, S. K., Pitrou, A., &amp; Seibert, S. 2015, in Proceedings of the Second Workshop on the LLVM Compiler Infrastructure in HPC, 1-6

Lovisari, L., Ettori, S., Gaspari, M., &amp; Giles, P. A. 2021, Universe, 7, 139

Lovisari, L. &amp; Maughan, B. J. 2022, in Handbook of X-ray and Gamma-ray Astrophysics, 65

Mahdavi, A. &amp; Geller, M. J. 2001, Astrophysical Journal, Letters, 554, L129 Mazure, A., Katgert, P., den Hartog, R., et al. 1996, Astronomy and Astrophysics, 310, 31

- Migkas, K. 2024, arXiv e-prints, arXiv:2406.01752
- Migkas, K., Pacaud, F., Schellenberger, G., et al. 2021, A&amp;A, 649, A151
- Migkas, K. &amp; Reiprich, T. H. 2018, A&amp;A, 611, A50
- Migkas et al. 2020, A&amp;A, 636, A15
- Nagarajan, A., Pacaud, F., Sommer, M., et al. 2019, Monthly Notices of the RAS, 488, 1728
- Nastasi, A., Böhringer, H., Fassbender, R., et al. 2014, Astronomy and Astrophysics, 564, A17
- Ortiz-Gil, A., Guzzo, L., Schuecker, P., Böhringer, H., &amp; Collins, C. A. 2004, Monthly Notices of the RAS, 348, 325
- Pi ff aretti, R., Arnaud, M., Pratt, G. W., Pointecouteau, E., &amp; Melin, J. B. 2011, Astronomy and Astrophysics, 534, A109

Planck Collaboration, Ade, P. A. R., Aghanim, N., et al. 2016, A&amp;A, 594, A27 Planck Collaboration, Aghanim, N., Akrami, Y., et al. 2020a, A&amp;A, 641, A1 Planck Collaboration, Aghanim, N., Akrami, Y., et al. 2020b, J Astronomy and Astrophysics, 641, A6

- Popesso, P., Biviano, A., Boehringer, H., &amp; Romaniello, M. 2007, VizieR Online Data Catalog, J A / + A 461 397 / /

Pratt, G. W., Arnaud, M., Biviano, A., et al. 2019, Space Sci. Rev., 215, 25 Rameez, M., Mohayaee, R., Sarkar, S., &amp; Colin, J. 2018, MNRAS, 477, 1772 Reiprich, T. H., Basu, K., Ettori, S., et al. 2013, Space Science Reviews, 177, 195

- Rines, K. J., Geller, M. J., Diaferio, A., &amp; Hwang, H. S. 2016, Astrophysical Journal, 819, 63
- Rubart, M. &amp; Schwarz, D. J. 2013, A&amp;A, 555, A117
- Sarkar, S., Pandey, B., &amp; Khatri, R. 2019, MNRAS, 483, 2453
- Secrest, N. J., von Hausegger, S., Rameez, M., Mohayaee, R., &amp; Sarkar, S. 2022, ApJ, 937, L31
- Sereno, M. &amp; Ettori, S. 2015, VizieR Online Data Catalog, J MNRAS 450 3675 / / / Singal, A. K. 2011, ApJ, 742, L23

Sohn, J., Geller, M. J., &amp; Zahid, H. J. 2019, Astrophysical Journal, 880, 142

Tiwari, P. &amp; Nusser, A. 2016, J. Cosmology Astropart. Phys., 2016, 062

- Voges, W., Aschenbach, B., Boller, T., et al. 1999, Astronomy and Astrophysics, 349, 389
- ˇ ípa, J. &amp; Shafieloo, A. 2017, ApJ, 851, 15 R

Weisskopf, M. C., Tananbaum, H. D., Van Speybroeck, L. P., &amp; O'Dell, S. L. 2000, in Society of Photo-Optical Instrumentation Engineers (SPIE) Conference Series, Vol. 4012, X-Ray Optics, Instruments, and Missions III, ed. J. E. Truemper &amp; B. Aschenbach, 2-16

- Willingale, R., Starling, R. L. C., Beardmore, A. P., Tanvir, N. R., &amp; O'Brien, P. T. 2013, MNRAS, 431, 394
- Wilson, S., Hilton, M., Rooney, P. J., et al. 2016, Monthly Notices of the RAS, 463, 413
- Xue, Y.-J. &amp; Wu, X.-P. 2000, Astrophysical Journal, 538, 65
- Yeung, S. &amp; Chu, M. C. 2022, Phys. Rev. D, 105, 083508
- Yoon, M., Huterer, D., Gibelyou, C., Kovács, A., &amp; Szapudi, I. 2014, MNRAS, 445, L60
- Zhang, Y. Y., Andernach, H., Caretta, C. A., et al. 2011, Astronomy and Astrophysics, 526, A105

## Appendix A: Choice for the cone and interval sizes

The cone size is chosen such that the number of clusters inside the cone is as large as possible without covering too large an area. Four cone sizes are considered to decide the best size for the analysis. The radii of these cones are 90 ◦ , 75 , 60 , and 45 . ◦ ◦ ◦ Fig. A.1 shows the number of clusters inside a cone for the relation σ v -T for di ff erent cone sizes.

Fig. A.1. Maps of the number of clusters for di ff erent cone sizes in the relation σ v -T . The cone sizes are 90 ◦ ( top left ), 75 ◦ ( top right ), 60 ◦ ( bottom left ), and 45 ◦ ( bottom right ). Note that the colour scale is di ff erent for the four plots.

![image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00726v1-figure-726bab65-4969-4d19-8a4a-12bbcac9d0b4.jpg)

The cone size of 45 ◦ covers a small area, and the number of clusters inside the cone is also small. There are less than 10 clusters in the region of the Galactic belt. The cone size of 90 ◦ covers a large area, and the number of clusters inside the cone is also large. This creates a significant overlap between the cones. The cone size of 60 ◦ has a high number of clusters for most regions, but certain regions have less than 20 clusters. The cone size of 75 ◦ strikes a good balance between the area covered and the number of clusters inside the cone. Thus, the cone size of 75 ◦ is chosen for the analysis.

Fig. A.2. Sigma map of the A from the σ v -T relation for di ff erent step sizes. The step sizes are 10 ◦ ( top left ), 5 ◦ ( top right ), and 1 ◦ ( bottom ).

![image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00726v1-figure-e5e60f27-a348-4845-b0d7-5fbd98cdc0cd.jpg)

The centre for these cones is taken at a step size of 5 ◦ . This is mainly done to reduce the computation time. Fig. A.2 shows the significance of the best-fit A for the relation σ v -T for di ff erent step sizes. The results show that the choice of step size does not a ff ect the results significantly.

In the significance map with a step size of 1 ◦ , some artefacts are seen as curved lines. These artefacts can be seen in the significance map with a step size of 5 , but they are not as prominent. ◦ These artefacts are due to strong outliers present in a particular region. These clusters are not considered outliers when considering the entire sample. All regions with this strong outlier will produce similar best-fit results, creating the artefacts.

## Appendix B: Redshift trends in the residuals

Using the Y X best-fit minimisation method, strong | L X (and Y SZ) residual trends are observed as a function of cluster redshift. Fig. B.1 shows the y -axis residuals for the relations L X -σ v and Y SZ -σ v as a function of z using the Y X and X Y minimisa-| | tion methods.

It is generally known that when there are significant uncertainties on the x -axis while minimising on the y -axis, the best-fit results can show strong biases. Therefore, we adopt the X Y min-| imisation for this scaling relation as these trends are not observed in this method.

Fig. B.1. Y-axis residual trends with respect to cluster redshift for the L X -σ v ( left ) and Y SZ -σ v ( right ) relations. The top panel displays the residuals using the Y X minimisation, and the bottom panel displays the | residuals using the X Y minimisation method. |

![image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00726v1-figure-e5d1d1b1-1769-48b3-9322-47084a0e409d.jpg)

## Appendix C: Comparison between free slope and fixed slope

We observe variations in the best-fit slope in certain regions for all three relations (see Fig. 4 and Fig. C.1). Although there are no strong correlations between the best-fit parameters, additional tests are conducted by keeping the slope fixed at the value obtained from the full sample.

Fig. C.1. Significance maps of best-fit B for L X -σ v ( left ), and Y SZ -σ v ( right ) relations.

![image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00726v1-figure-10596916-278b-4c36-81a8-1324b2ed84b3.jpg)

Figure C.2 compares the results from free and fixed slope analysis for all three relations. We calculate ∆ σ A as the di ff erence of σ A obtained from the fixed and free slope analysis. For the σ v -T and L X -σ v relations, these di ff erences are less than ± 0 5 . σ . However, for the Y SZ -σ v , there are somewhat higher

Table C.1. Comparison of best-fit parameters for the L X -σ v and Y SZ -σ v scaling relations with di ff erent datasets. Results are shown for both XYand YX methods. The intrinsic scatter mentioned here is measured in the direction in which the scatter is minimised. | |

| Dataset               | C Y                   | C X                   | Method                | A                                                     | B                                                     | σ int (dex)                                           | N                     |
|-----------------------|-----------------------|-----------------------|-----------------------|-------------------------------------------------------|-------------------------------------------------------|-------------------------------------------------------|-----------------------|
| L X - σ v ( γ = - 1)  | L X - σ v ( γ = - 1)  | L X - σ v ( γ = - 1)  | L X - σ v ( γ = - 1)  | L X - σ v ( γ = - 1)                                  | L X - σ v ( γ = - 1)                                  | L X - σ v ( γ = - 1)                                  | L X - σ v ( γ = - 1)  |
| eeHIFLUGCS            | 10 44 erg s - 1       | 750 km s - 1          | XY | YX |             | 1 . 515 + 0 . 134 - 0 . 122 1 . 538 + 0 . 103 - 0 097 | 5 . 988 + 0 . 452 - 0 . 394 2 . 722 + 0 . 214 - 0 213 | 0 . 080 + 0 . 005 - 0 . 005 0 . 381 + 0 . 022 - 0 021 | 195 200               |
| MCXC                  | 10 44 erg s - 1       | 750 km s - 1          | XY | YX |             | 1 . 383 + 0 . 112 - 0 . 105 1 . 434 + 0 . 081 - 0 077 | 6 . 536 + 0 . 410 - 0 . 374 2 . 910 + 0 . 186 - 0 189 | 0 . 082 + 0 . 004 - 0 . 004 0 . 384 + 0 . 018 - 0 017 | 316 319               |
| Y SZ - σ v ( γ = + 1) | Y SZ - σ v ( γ = + 1) | Y SZ - σ v ( γ = + 1) | Y SZ - σ v ( γ = + 1) | Y SZ - σ v ( γ = + 1)                                 | Y SZ - σ v ( γ = + 1)                                 | Y SZ - σ v ( γ = + 1)                                 | Y SZ - σ v ( γ = + 1) |
| S / N ≥ 2             | 30 kpc 2              | 800 km s - 1          | XY | YX |             | 1 . 490 + 0 . 112 - 0 . 106 1 . 341 + 0 . 068 - 0 069 | 6 . 419 + 0 . 423 - 0 . 384 2 . 915 + 0 . 194 - 0 195 | 0 . 073 + 0 . 004 - 0 . 004 0 . 342 + 0 . 018 - 0 017 | 284 285               |
| S / N ≥ 3             | 30 kpc 2              | 800 km s - 1          | XY | YX |             | 1 . 413 + 0 . 118 - 0 . 106 1 . 391 + 0 . 075 - 0 070 | 6 . 643 + 0 . 499 - 0 . 427 2 . 892 + 0 . 208 - 0 211 | 0 . 071 + 0 . 004 - 0 . 004 0 . 334 + 0 . 019 - 0 017 | 256 256               |
| S / N ≥ 4 . 5         | 30 kpc 2              | 800 km s - 1          | XY | YX |             | 1 . 343 + 0 . 126 - 0 . 119 1 . 520 + 0 . 089 - 0 083 | 6 . 727 + 0 . 593 - 0 . 519 2 . 702 + 0 . 233 - 0 229 | 0 . 071 + 0 . 005 - 0 . 004 0 . 325 + 0 . 020 - 0 019 | 201 201               |

Fig. C.2. Di ff erence between the significance of A obtained from fixed slope and free slope analysis for the σ v -T ( top left ), L X -σ v ( top right ), and Y SZ -σ v ( bottom ) relations. Top panel maps have the same colour scale ( -0 5 to . + 0 5), whereas the bottom panel map has a colour scale . of -0 9 to . + 0 9 to visualise higher variations. .

![image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00726v1-figure-c209d0f6-d2e0-4087-a377-f171249487ec.jpg)

di ff erences ( ∼ 0 8 . σ ) in certain regions. Nonetheless, the direction and magnitude of the maximum significance obtained for this relation remains unchanged. Therefore, we can infer that the variations in the best-fit slope are not a major concern as they do not impact the conclusions drawn from the analysis as expected, given the rather weak correlation observed between slope and normalisation.

## Appendix D: General behaviour of scaling relations for different datasets and fitting methods

The Table C.1 presents the best-fit parameters for the relations L X -σ v and Y SZ -σ v obtained using di ff erent datasets and fitting methods. The best-fit parameters are consistent across the di ff erent datasets. Notably, the X Y fitting method consistently |

yields a higher slope than the Y X minimisation method for the | various datasets.

## Appendix E: z distributions

The Fig. E.1 shows the distribution of cluster redshift in the region of maximum anisotropy and the rest of the sky for the two relations L X -σ v and Y SZ -σ v . We observe some di ff erences between the distributions of the two regions. On average, clusters in the most anisotropic regions have lower z .

Fig. E.1. Distribution of redshifts for the clusters inside the cone centred at maximum anisotropy and rest of the sky for L X -σ v ( left ), and Y SZ -σ v ( right ) relations

![image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00726v1-figure-448b932b-0f5f-4ab0-b108-65a2e9477f24.jpg)

Fig. E.2. 2-D map of median redshifts of the clusters inside a cone for L X -σ v ( left ), and Y SZ -σ v ( right ) relations

![image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00726v1-figure-0b34d94e-d9f2-4ea5-b768-1ae636c0612b.jpg)

We also create a 2-D map of the median z within each cone to identify other z trends in our data. The most significant trend

observed is the galactic north-south divide for the L X -σ v relation. On average, clusters are found at higher z in the northern galactic sky compared to the southern half.

## Appendix F: Relation between A and H 0

The best-fit normalisation is related to the Hubble parameter, and variations in the best-fit normalisation can be converted to variations in the Hubble parameter. For the relation L X -σ v ( L X = A × σ B v ),

$$f _ { X } K _ { \text{cor} } \times 4 \pi ( D _ { \text{L} } ( z ) ) ^ { 2 } = A \times \sigma _ { \text{v} } ^ { B } \,.$$

Here, f X is the observed X-ray flux, K corr is the K correction used to convert L X from the rest frame of the observer to the cluster's rest frame, and D L( ) is the luminosity distance, which is given z by

$$D _ { \mathbb { L } } ( z ) = \frac { c } { H _ { 0 } } ( 1 + z ) \int _ { 0 } ^ { z } \frac { d z ^ { \prime } } { E ( z ^ { \prime } ) } \,.$$

Thus, the equation becomes

$$f _ { X } K _ { \text{cor} } \times 4 \pi \left ( c ( 1 + z ) \int _ { 0 } ^ { z } \frac { d z ^ { \prime } } { E ( z ^ { \prime } ) } \right ) ^ { 2 } = H _ { 0 } ^ { 2 } A \times \sigma _ { \text{v} } ^ { B } \,. \quad \quad \text{(F.)}$$

This is also true for the Y SZ -σ v relation since the quantity Y SZ depends on the integrated Compton parameter Y and square of the angular diameter distance ( D A( )). Thus, the equaz tion Y SZ = A × σ B v transforms to

$$Y \left ( \frac { c } { ( 1 + z ) } \int _ { 0 } ^ { z } \frac { d z ^ { \prime } } { E ( z ^ { \prime } ) } \right ) ^ { 2 } = H _ { 0 } ^ { 2 } A \times \sigma _ { \text{v} } ^ { B }$$

This shows that absolute constraints can only be put on the quantity H A 2 0 . To obtain the value of H 0 in di ff erent regions of the sky from A , equation 6 is used.
A. Pandya K. Migkas T. H. Reiprich A. Stanford F. Pacaud G. Schellenberger L. Lovisari M. E. Ramos-Ceja N. T. Nguyen-Dang S. Park
Examining the local Universe isotropy with galaxy cluster velocity dispersion scaling relations
In standard cosmology, the late Universe is assumed to be statistically
homogeneous and isotropic. However, a recent study based on galaxy clusters by
Migkas et al. (2021, arXiv:2103.13904) found an apparent spatial variation of
approximately $9\%$ in the Hubble constant, $H_0$, across the sky. The authors
utilised galaxy cluster scaling relations between various cosmology-dependent
cluster properties and a cosmology-independent property, i.e., the temperature
of the intracluster gas $(T)$. A position-dependent systematic bias of $T$
measurements can, in principle, result in an overestimation of apparent $H_0$
variations. In this study, we search for directional $T$ measurement biases by
examining the scaling relation between the member galaxy velocity dispersion
and the gas temperature $(\sigma_\mathrm{v}-T)$. Additionally, we search for
apparent $H_0$ angular variations independently of $T$ by analysing the
relations between the X-ray luminosity and Sunyaev-Zeldovich signal with the
velocity dispersion, $L_\mathrm{X}-\sigma_\mathrm{v}$ and
$Y_\mathrm{SZ}-\sigma_\mathrm{v}$. We utilise Monte Carlo simulations of
isotropic cluster samples to quantify the statistical significance of any
observed anisotropies. We find no significant directional $T$ measurement
biases, and the probability that a directional $T$ bias causes the previously
observed $H_0$ anisotropy is only $0.002\%$. On the other hand, from the joint
analysis of the $L_\mathrm{X}-\sigma_\mathrm{v}$ and
$Y_\mathrm{SZ}-\sigma_\mathrm{v}$ relations, the maximum variation of $H_0$ is
found in the direction of $(295^\circ\pm71^\circ, -30^\circ\pm71^\circ)$ with a
statistical significance of $3.64\sigma$, fully consistent with
arXiv:2103.13904. Our findings strongly corroborate the previously detected
spatial anisotropy of galaxy cluster scaling relations using a new independent
cluster property, $\sigma_\mathrm{v}$.