text:
# XMM-Newton observation of the distant (𝑧=0.6) galaxy cluster RX J1120.1+4318

## 1 Introduction

Galaxy clusters occupy a unique position in the scenario of hierarchical structure formation, as they are still forming today. The statistical properties of this evolving population (e.g. the mass distribution function and the correlations between physical quantities, at various redshifts) therefore provide unique constraints on cosmological scenarios.

In the simplest model of structure formation, purely based on gravitation, galaxy clusters constitute a homologous family. Clusters are self-similar in shape, and predictable scaling laws relate each physical property to the cluster total mass (or temperature) and redshift $`z`$ (e.g. Bryan & Norman bryan (1998); Eke et al. eke (1998)). From observations with ROSAT and ASCA satellites, it is now well established that this simple model fails to explain all the observed structural and scaling properties of the nearby cluster population (Tozzi & Norman tozzi (2001) and reference therein). The evolution with redshift of these properties is an essential piece of information, still largely missing, to reconstruct the physics of the formation process. This information is also important for accurate $`\mathrm{\Omega }_0`$ estimates based on the evolution of the cluster mass function (e.g. Oukbir & Blanchard oukbir (1997)). To compare the results of flux-limited X-ray surveys with the prediction of the various theoretical models, it is necessary to understand the relation between observed quantities, like the cluster temperature and luminosity, and the virial mass. Precise $`\mathrm{\Omega }_0`$ estimate also requires a good understanding of the survey selection function, which further depends sensitively on cluster morphology (e.g. Adami et al. adami (2000)).

Following the pioneering work of the Einstein medium Sensitivity Survey (Gioia et al. gioia90 (1990)), several large, well controlled, X–ray samples of distant clusters have been assembled in the last years, using ROSAT observations (see Gioia gioia (2000) for a review). The exceptional sensitivity of XMM, associated with good spectroscopic and imaging capabilities, now allows the detailed analysis of these clusters, down to what was only a detection limit with ROSAT. We present here the XMM observation of RX J1120.1+4318, a cluster at $`z=0.6`$ detected in the bright Serendipitous High-redshift Archival ROSAT Cluster (SHARC) survey (Romer et al. romer (2000)). This observation was made in the framework of the XMM-Newton $`\mathrm{\Omega }`$ project, a systematic XMM Guaranteed Time follow-up of the most distant ($`z>0.45`$) SHARC clusters (Bartlett et al. bartlett (2001)).

The paper is organized as follows. In Sec. 2, we describe the data analysis performed to derive the surface brightness profile, global temperature and temperature profiles, which are presented in Sec. 3. In Sec. 4, we compare the physical properties of RX J1120.1+4318 with the predictions of the self-similar model of cluster formation. The scaled emission measure and temperature profiles, the gas mass fraction and $`L_\mathrm{X}`$$`T`$ relation are compared to those of nearby clusters. This study is made for two cosmological models: an Einstein-de Sitter Universe (EdS, $`\mathrm{\Omega }_0=1`$) and a flat low density Universe ($`\mathrm{\Lambda }\mathrm{CDM}`$ ) with $`\mathrm{\Omega }_0=0.3`$ and $`\mathrm{\Lambda }=0.7`$. In Sec. 5, we study the thermodynamical state of this cluster - possible presence of cooling gas, entropy content - to further assess the role of non-gravitational processes in cluster formation. Section 6 contains our conclusions.

Unless otherwise stated all errors on the cluster parameters are at the $`90\%`$ confidence level. A Hubble constant of $`\mathrm{H}_0=50\mathrm{km}/\mathrm{s}/\mathrm{Mpc}`$ is assumed. At the cluster redshift, 1 arcmin corresponds to $`456\mathrm{kpc}`$ and $`561\mathrm{kpc}`$ for the EdS model and the $`\mathrm{\Lambda }\mathrm{CDM}`$ model, respectively.

## 2 Data Analysis

### 2.1 Data preparation

RX J1120.1+4318 was observed for $`20`$ ksec on May 8, 2001 with the EPIC/MOS and pn camera (using the THIN optical blocking filter) in Full Frame mode (Turner et al. turner (2001); Strüder et al. struder (2001)). We generated calibrated event files using the tasks emchain and epchain of the SAS V5.1.

We discarded the data corresponding to the periods of high background induced by solar flares (e.g. see Arnaud et al. arnaud01 (2001)). We extracted the light curves in the energy band $`[1012]\mathrm{keV}`$ and $`[1214]\mathrm{keV}`$ for the MOS and pn data respectively. In these energy bands, the effective area of XMM is negligible and the emission is dominated by the particle induced background. We removed all frames corresponding to a count rate greater than $`15\mathrm{ct}/100\mathrm{s}`$ (MOS data) and $`22\mathrm{ct}/100\mathrm{s}`$ (pn data). After this selection, the remaining exposure time is $`17.6`$ ksec, $`17.9`$ ksec and $`14.2`$ ksec for the MOS1, MOS2 and pn observations.

### 2.2 Vignetting correction

The effective area of the XMM mirrors is a function of off-axis angle and this vignetting effect depends on energy. An additional vignetting effect is due to the RGA obscuration for the MOS camera. The vignetting calibration data were those available at time of release of SAS V5.1.

To correct for vignetting effects, we used the method proposed by Arnaud et al. (arnaud01 (2001)) for spectra, which can be generalized in a straightforward way to profiles or images (e.g. Majerowicz et al. majerowicz (2002)). For each event, we computed the corresponding weight coefficient, defined as the ratio of the effective area at the photon position and energy to the central effective area at that energy. When extracting spectra, image or surface brightness profile, each event is weighted by this coefficient111This vignetting correction can now be done with task EVIGWEIGHT of SAS V5.2. These ‘corrected’ products correspond to those we would obtain if the detector response were flat. The on-axis response can then be used for spectral fitting and to estimate the physical parameters of the cluster. We used the spectral response files m1$`\mathrm{\_}`$thin1v9q20t5r6$`\mathrm{\_}`$all$`\mathrm{\_}`$15.rsp, m2$`\mathrm{\_}`$thin1v9q20t5r6$`\mathrm{\_}`$all$`\mathrm{\_}`$15.rsp for EPOC/MOS1,2 and epn$`\mathrm{\_}`$ff20$`\mathrm{\_}`$sY9$`\mathrm{\_}`$thin.rmf V6.1 for EPIC/pn respectively222see http://xmm.vilspa.esa.es/calibration/. Those files include both the effective area and the redistribution matrix.

### 2.3 Background subtraction

After cleaning for flare events, the XMM background is dominated by the Cosmic X–ray background (CXB) and the non X–ray background (NXB) induced by high energy particles. The former component, which dominates at low energy, depends on the observed position in the sky and is vignetted by the optics. The latter component, which dominates at high energy, is not vignetted but varies slightly across the detector and with time.

To subtract the total XMM background, we used EPIC blank sky event files (one for each camera) obtained by combining several high galactic latitude pointings. The data are cleaned for background flares and bright point sources are excluded (Lumb lumb (2002)). The sky coordinates in the event files were modified using the aspect solution of the RX J1120.1+4318 observation, so that extraction can be done in sky coordinates, while insuring that the same detector region is considered for both the blank field and the source observation. The background level, estimated from the total count rate in the whole FOV in the high energy bands defined above, was found to be 0.85 times smaller for the source observation than for the blank sky data. This is typical of the variations observed since launch (Lumb lumb (2002)), which are explainable by variations in Cosmic Ray shielding as the magnetosphere is pumped up and down by solar activity.

The background subtraction for each source product (spectrum or profile) is done in two steps. The method is fully described in Appendix A. We first subtract the corresponding blank field product, obtained using the same spatial and energy selection, and normalized by the 0.85 factor defined above. For consistency, the blank field products were obtained using the same vignetting correction method as that used for the source. The correction factor is thus formally wrong for the NXB component. However, as we subtract blank field and source count rates estimated in the same region in detector coordinates this correction factor is the same and does not introduce bias. This first step thus allows us to subtract properly the NXB. However, we are left with a residual CXB component, which is the difference between the CXB in the source region and in the blank fields (multiplied by the normalization factor above). This residual is corrected for vignetting and is expected to be uniform all over the field of view. In a second step, we thus subtract this residual component, using data in the outer part of the FOV, outside the cluster region.

## 3 Results

The vignetting corrected image in the $`[0.35]\mathrm{keV}`$ energy band is displayed in Fig. 1. The data of all EPIC cameras are combined. The cluster has a regular spherical morphology, suggesting it is in a relaxed state.

### 3.1 Surface brightness profile

We extracted the surface brightness profile of the cluster in the $`[0.33]\mathrm{keV}`$ energy band. This band was chosen to optimize the signal-to-noise $`(S/N)`$ ratio. We cut out serendipitous sources in the field of view and binned the photons into concentric annuli with a width of 3.3” (3 pixels of the MOS camera) centered on the maximum of the X-ray emission for each camera. Since the cluster is regular, this center corresponds to the centroid of the emission. The three profiles were then summed. The vignetting correction and background subtraction was performed as described above. After subtraction of the corresponding blank field profile, only a contribution from the CXB remains (see Appendix A). Since the data are corrected for vignetting effects, the CXB surface brightness profile should be constant with radius. The profile was indeed found to be flat beyond $`4^{}`$, where we can thus consider that the cluster emission is negligible. The residual background was thus estimated from the data within $`4^{}7^{}`$ and subtracted from the profile. Starting from the central annulus, we re-binned the data in adjacent annulii so that i) at least a $`S/N`$ ratio of $`3\sigma `$ is reached after background subtraction and ii) the width of the bin increases with radius, with $`\mathrm{\Delta }(\theta )>0.1\theta `$. Such a logarithmic radial binning insures a $`S/N`$ ratio in each bin roughly constant in the outer part of the profile, when the background can still be neglected.

The resulting surface brightness profile, $`S(\theta )`$, is shown in Fig. 2. The cluster emission is significantly detected up to $`R_{\mathrm{det}}=3^{}`$ or $`1.37`$ Mpc for a critical density Universe. Beyond that radius it was not possible to create $`S/N>3`$ annulus of any width. The total count rate within the $`R<3\mathrm{}`$ region and in the considered energy band is $`0.30\pm 0.006`$ ct/s.

We fitted $`S(\theta )`$ with a $`\beta `$–model convolved with the XMM PSF, and binned as the observed profile. The PSF of each camera is modeled by a normalized King profile, with parameters depending on energy and off-axis angle (Ghizzardi ghizzardi (2001), Griffiths & Saxton griffiths02b (2002)). The overall PSF at each radius is obtained by summing the PSF of each camera, estimated at an energy of $`1\mathrm{keV}`$, weighted by the respective cluster count rate in the $`[0.33]\mathrm{keV}`$ energy band. The on-axis overall PSF is plotted in Fig. 2 (dashed line). The convolution with the PSF takes into account the (small) PSF variation across the region considered. However, in practice only the MOS PSF variations are taken into account. We used the on-axis pn PSF, due to the lack of available parametrical fit of the off-axis data.

The fit of the cluster profile gives $`\beta =0.78_{0.04}^{+0.06}`$ and a core radius of $`\theta _\mathrm{c}=0.44_{0.04^{}}^{+0.06^{}}`$. The reduced $`\chi ^2`$ is $`1`$ ($`\chi ^2=20`$ for $`21`$ d.o.f). The convolved best fit model is plotted in Fig. 2 (solid line), together with the corresponding (unbinned and non-convolved) $`\beta `$–model (dotted line). As can be seen, the PSF mostly affects the core of the profile.

We then estimated the cluster surface brightness profile corrected for the effect of the PSF, that we will use in the following (Sec. 4). In principle, we should have deconvolved the observed profile. However since the cluster profile is well fitted by a $`\beta `$–model , we used a simpler procedure. For each radial bin, we estimated the ratio of the surface brightness corresponding to the non-convolved and convolved best fit $`\beta `$–model . The observed surface brightness profile was then corrected for the effect of the PSF, by simply multiplying the observed $`S(\theta )`$ value in each bin by the model ratio obtained for that bin.

### 3.2 Mean temperature

The overall MOS1, MOS2 and pn spectra, extracted from the event file, are shown in Fig. 3. The spectra are corrected for vignetting and background subtracted (see Sect. 2.3). To optimize the $`S/N`$ ratio, the integration region for the cluster was restricted to $`2^{}`$ from the cluster center, and the residual CXB spectrum was estimated from the $`4^{}<\theta <11^{}`$ region. The spectra are binned so that the $`S/N`$ ratio is greater than 3 $`\sigma `$ in each energy bin after background subtraction.

The spectra are jointly fitted with XSPEC using a redshifted MEKAL model (Mewe et al. mewe1 (1985),mewe2 (1986); Kaastra kaastra (1992); Liedahl et al. liedahl (1995)). We let the relative normalization between the various instruments be free but assumed a common temperature and abundance. When letting the hydrogen column density $`N_\mathrm{H}`$ and redshift $`z`$ free, we obtained $`N_\mathrm{H}=2.2\pm 1.2\times 10^{20}\mathrm{cm}^2`$, in agreement with the 21 cm value ($`N_\mathrm{H}=2.26\times 10^{20}\mathrm{cm}^2`$ from Dickey & Lockman dickey (1990)) and $`z=0.61\pm 0.03`$, consistent with the optical value ($`z=0.60`$). We then fixed these parameters to the 21 cm and optical values. The best fit gives $`\mathrm{k}T=5.3\pm 0.5\mathrm{keV}`$ and an abundance of $`0.47\pm 0.19`$. The reduced $`\chi ^2`$ is $`0.94`$ ($`\chi ^2=258`$ for $`276`$ d.o.f).

The temperature values estimated separately from the MOS and pn spectra are consistent within the error bars ($`\mathrm{k}T=5.8_{0.7}^{+1.0}\mathrm{keV}`$ for MOS data and $`\mathrm{k}T=4.5_{0.5}^{+0.8}\mathrm{keV}`$ for pn data). We note, however that a lower temperature is obtained with the pn data than with the MOS data. This is likely to be partly due to the remaining calibration systematic uncertainties in the XMM spectral responses. Indeed, it is known that the EPIC-pn and MOS cameras show a relative flux difference of $`4\%`$ at low energies which increases with energy above $`4.5\mathrm{keV}`$, resulting in a MOS spectral slope flatter than the pn (Saxton saxton (2002); Griffiths et al. griffiths02a (2002)).

### 3.3 Temperature profile

We then extracted the spectra in three concentric annuli, centered on the cluster X–ray emission peak and fitted the data as described above. The corresponding temperature profile is shown in Fig. 4. It appears flat up to $`2^{}`$, within the error bars.

The emission from the central annulus is affected by the PSF and may contaminate the other annulii, blurring out gradients if they exist. Moreover the energy dependence of the PSF, if not taken into account, might in principle bias the temperature estimate, since the photon redistribution is energy dependent. However, for the MOS1 instrument for instance, the Encircled Energy Fraction (EEF) within $`15\mathrm{"}`$, the size of the first bin, is already about $`70\%`$ (Ghizzardi ghizzardi (2001)). The flux in the first bin is about 1/3 of the second bin flux. We thus expect a contamination of only $`10\%`$ in the second bin and we probably do not significantly underestimate possible gradients in the central part. Furthermore, the EEF varies only by $`4.7\%`$ between $`0.3\mathrm{keV}`$ and $`5\mathrm{keV}`$, the minimum and maximum energy for the spectra. Neglecting the PSF energy dependence results in a negligible bias ($`0.5\%`$) at high energies. We emphasize that the XMM telescope PSF has a very weak energy dependence, so that indeed the correction is tiny and we are not facing the problem experienced in spectro-imagery with ASCA, for example.

## 4 Scaling Properties

### 4.1 The self-similar model

The self-similar model is based on simple assumptions for cluster formation, derived from the top-hat spherical collapse model. The virialised part of a cluster, present at a given redshift, corresponds to a fixed density contrast as compared to the critical density of the Universe at that redshift and the internal shape of clusters of different masses and $`z`$ are similar. Self similarity applies to both the dark matter component and the hot intra-cluster medium, whose mass fraction is assumed to be constant. Consequently simple scaling laws relate each physical X–ray property, $`Q`$, to the cluster total mass (or temperature $`T`$) and redshift, in the form $`QA(z)T^\alpha `$.

This simple model thus makes a definitive prediction in terms of the evolution of cluster properties. First, the normalization (but not the slope) of the $`Q`$$`T`$ scaling relations should evolve with $`z`$. In particular, for a critical density Universe, the virial radius at fixed $`T`$ should decrease with $`z`$ as $`(z+1)^{3/2}`$, while the luminosity and central emission measure increase as $`(z+1)^{3/2}`$ and $`(z+1)^{9/2}`$ respectively. Furthermore once expressed in scaled coordinates333The considered quantity is normalized according to the scaling relation estimated at the cluster temperature and redshift and the radius is expressed in units of the virial radius, the radial profile of any physical quantity (e.g density, temperature) should be the same at all redshifts.

In this section we will compare the structural and scaling properties of RX J1120.1+4318 to the properties of nearby clusters. We will consider the scaled emission measure and temperature profiles, the gas mass fraction and the $`L_\mathrm{X}`$$`T`$ relation. Results are given for both the EdS model and the $`\mathrm{\Lambda }\mathrm{CDM}`$ model444The comparison of distant and nearby cluster properties does not depend on the assumed $`\mathrm{H}_0`$ value: scaled and physical quantities depend on $`\mathrm{H}_0`$ via a multiplying factor, the same at all redshifts.. The physical properties of RX J1120.1+4318 are summarized In Table 1.

### 4.2 The scaled emission measure profile

Recently, Arnaud et al. (arnaud02 (2002), hereafter AAN) studied the surface brightness profiles of a sample of 25 distant $`(z=0.30.83)`$ hot $`(\mathrm{k}T>3.5\mathrm{keV})`$ clusters observed with ROSAT, with published temperature from ASCA. They found that the scaled profiles of distant clusters are perfectly consistent with the average scaled profile of nearby clusters, for a flat low density Universe. For RX J1120.1+4318 , we will make the same comparison, and we only briefly summarize the method.

#### 4.2.1 The cluster profile

We used the surface brightness profile, corrected for the effect of the PSF (see Sect. 3.1). It is converted to emission measure ($`EM`$) profile:

$$EM(r)S(\theta )(1+z)^4/ϵ(T,z)$$
(1)
where $`ϵ(T,z)`$ is the emissivity in the energy band considered, taking into account the instrument response and interstellar absorption, and $`r=\theta d_\mathrm{A}`$, where $`d_\mathrm{A}`$ is the angular distance.

The $`EM`$ profile is then scaled according to the self-similar model. We use the standard scaling relations of cluster properties with redshift and temperature, with the empirical slope of the $`M_{\mathrm{gas}}`$$`T`$ relation derived by Neumann & Arnaud (neumann01 (2001)). This empirical relation is consistent with the observed steepening of the $`L_\mathrm{X}`$$`T`$ relation and reduces significantly the scatter in the scaled profile of nearby clusters. The physical radius is thus scaled to the virial radius, $`R_\mathrm{V}`$, with $`R_\mathrm{V}\mathrm{\Delta }_\mathrm{z}^{1/2}\left(1+z\right)^{3/2}T^{1/2}`$ and the emission measures by $`\mathrm{\Delta }_\mathrm{z}^{3/2}(1+z)^{9/2}T^{1.38}`$ (see AAN for details). This scaling depends on the density parameter $`\mathrm{\Omega }_0`$ and the cosmological constant $`\mathrm{\Lambda }`$, via the factor $`\mathrm{\Delta }_\mathrm{z}=(\mathrm{\Delta }_\mathrm{c}(\mathrm{\Omega }_\mathrm{z},\mathrm{\Lambda })\mathrm{\Omega }_0)/(18\pi ^2\mathrm{\Omega }_\mathrm{z})`$, where $`\mathrm{\Delta }_\mathrm{c}(\mathrm{\Omega }_\mathrm{z},\mathrm{\Lambda })`$ is the cluster density contrast at redshift $`z`$ and $`\mathrm{\Omega }_\mathrm{z}`$ the corresponding Universe density parameter. Analytical expression of $`\mathrm{\Delta }_\mathrm{c}(\mathrm{\Omega }_\mathrm{z},\mathrm{\Lambda })`$ can be found in Bryan & Norman (1998). The scaled profiles also depend on these cosmological parameters via the angular distance used to convert angular radius to physical radius.

The scaled emission measure profiles of RX J1120.1+4318 derived for the EdS model and the $`\mathrm{\Lambda }\mathrm{CDM}`$ model are displayed in the left and right panel of Fig. 5, respectively. The corresponding virial radii are given in Table 1. As in AAN, the normalization of the $`R_\mathrm{V}`$$`T`$ relation is taken from the simulation Evrard et al. (evrard (1996)).

#### 4.2.2 Comparison with local data

For each cosmological model, the scaled profile of RX J1120.1+4318 is compared in Fig. 5 with the corresponding average scaled profile of nearby clusters. This profile (full line), and the typical dispersion around it (dotted lines) is derived, as in AAN, from a sample of 15 hot nearby clusters observed with ROSAT (Neumann & Arnaud neumann99 (1999)). The comparison is performed at radii beyond $`0.1R_\mathrm{V}`$, where nearby clusters were found to obey self-similarity555A large dispersion was observed in the cluster core properties, which are dominated by non-gravitational physics..

To be meaningful, this comparison requires that the relative error in the calibration of XMM/EPIC and ROSAT/PSPC can be neglected. Our previous study of A1795 (Arnaud et al. arnaud01 (2001)), where we performed a combined fit of the EPIC/MOS and ROSAT/PSPC spectra, indicates that the fluxes derived by the two instruments match at the $`\pm 5\%`$ level. Similar relative calibration uncertainties were derived by Snowden (snowden (2002)).

For both cosmological models, the scaled profile of RX J1120.1+4318 is consistent with the local reference profile, taking into account its intrinsic scatter. The RX J1120.1+4318 data thus reinforce the validity of the self-similar model. However, some significant discrepancy in shape can be noted. At large radii the cluster profile falls off more rapidly than the reference profile. This corresponds to its higher $`\beta `$ value, $`\beta =0.78`$, as compared to the universal value of $`\beta =2/3`$, which fits well the reference curve (Neumann & Arnaud neumann99 (1999)). The scaled core radius, $`x_\mathrm{c}=r_\mathrm{c}/R_\mathrm{V}`$, on the other hand, is similar to the value ($`x_\mathrm{c}=0.12`$) of the reference profile: we obtain $`x_\mathrm{c}0.11`$ (EdS model) and $`x_\mathrm{c}0.14`$ ($`\mathrm{\Lambda }\mathrm{CDM}`$ model). This might simply be due to intrinsic scatter in the properties of distant clusters. A typical scatter of $`20\%`$ is indeed observed in the slopes of nearby cluster profiles (Neumann & Arnaud neumann99 (1999)). A scatter in cluster properties is observed as well in numerical simulations of cluster formation (Navarro, Frenk & White navarro97 (1997)). However, this discrepancy could also indicate a systematic departures from the self-similar model considered. With data on a single cluster, we obviously cannot distinguish between these two possibilities.

It is also interesting to compare in more detail the results obtained for the two cosmological models. The relative position of the scaled profile of RX J1120.1+4318 with respect to the reference profile is higher (in y axis) for the $`\mathrm{\Lambda }\mathrm{CDM}`$ model than for the EdS model. As discussed by AAN, the derived scaled profiles of distant clusters depend on the cosmological parameters, mostly via the angular distance $`d_\mathrm{A}`$. The typical dependence is $`d_\mathrm{A}^3`$ at fixed scaled radius. At $`z=0.6`$, this $`d_\mathrm{A}^3`$ factor is about $`85\%`$ higher for a $`\mathrm{\Lambda }\mathrm{CDM}`$ model than for an EdS model. The $`d_\mathrm{A}^3`$ dependence is strictly exact for a profile of logarithmic slope $`3`$ (e.g. in the cluster external region for a $`\beta `$–model with $`\beta =2/3`$). Due to the steepest slope of the RX J1120.1+4318 profile at large radii, the effect of varying the cosmological parameters is slightly smaller, of the order of $`60\%`$ (Fig. 5). The effect is smaller than the intrinsic dispersion in the nearby profiles at radius below typically $`0.2R_\mathrm{V}`$, and becomes only marginally larger than the $`\pm 25\%`$ dispersion at higher radii. As emphasised by AAN, no definitive conclusion on the cosmological parameters can then be drawn from the observation of a single distant cluster. We simply note that the profile of RX J1120.1+4318 is in better agreement with the local reference profile at large radii for a $`\mathrm{\Lambda }\mathrm{CDM}`$ model than for an EdS model. In the latter model, the profile appears to deviate more and more from the reference profile above $`0.4R_\mathrm{V}`$. The better agreement observed for the $`\mathrm{\Lambda }\mathrm{CDM}`$ model can simply be an artifact, due to the intrinsic dispersion in the cluster properties, or again due to systematic departure from the self-similar model. However, it is consistent with the statistical analysis of the scaled profiles of distant clusters, performed by AAN, which clearly favors a low density Universe.

Finally, we note the statistical quality of the scaled emission measure profile, as compared to ROSAT/ASCA data. The profile of RX J1120.1+4318 is traced nearly up to the virial radius, further out than the mean profile of nearby clusters observed with ROSAT and the errors bars on most of the data points are smaller than the typical dispersion in the reference profile. Note also the very good sampling of the profile. Using the same criteria (at least $`3\sigma `$ detection in each bin), most of the ROSAT profiles of distant clusters (AAN, Figure 4) are much more sparsely sampled.

### 4.3 The scaled temperature profile

Cluster temperature profiles can now be measured with high precision with XMM. First observations of nearby clusters indicate that the profiles are isothermal (within typically 10%) up to half the virial radius (Arnaud et al. arnaud01 (2001); Arnaud arnaud01b (2001)).

With XMM we can now also study the evolution of the temperature profiles. The cluster temperature profile, normalized to the mean temperature, is plotted in Fig. 6 as a function of scaled radius ($`\mathrm{\Lambda }\mathrm{CDM}`$ model). The temperature profile of RX J1120 is mapped up to about $`0.5R_\mathrm{V}`$. This is the first time that a temperature profile can be measured at such high redshift. Although the errors bars are large, this distant cluster appears isothermal, as observed for nearby clusters. This is consistent with the self-similarity in shape expected in the simple model of cluster formation.

### 4.4 The gas mass fraction

#### 4.4.1 VT method

We first estimated the virial radius ($`R_\mathrm{V}`$) and mass ($`M_\mathrm{V}`$) from the best fit temperature and the theoretical $`R_\mathrm{V}`$$`T`$ and $`M_\mathrm{V}`$$`T`$ scaling relations at the cluster redshift, as described in Sec. 4.1. We recall that these relations correspond to a fixed density contrast at redshift $`z`$ and are derived using the virial theorem (VT):

$$\frac{\mathrm{G}\mu m_\mathrm{p}M_\mathrm{V}}{2R_\mathrm{V}}=\beta _\mathrm{T}\mathrm{k}T$$
(2)
where the normalization factor is assumed to be $`\beta _\mathrm{T}=1.05`$, from the simulation of Evrard et al. (evrard (1996)). In that case:

$`M_\mathrm{V}`$ $`=`$ $`\mathrm{2.98\; 10}^{15}\mathrm{\Delta }_\mathrm{z}^{1/2}(1+z)^{3/2}\left({\displaystyle \frac{\mathrm{k}T}{10\mathrm{keV}}}\right)^{3/2}\mathrm{M}_{}`$ (3)
The best fit temperature corresponds to the emission weighted temperature within $`2^{}`$, which is about half of the virial radius (see Table 1). Since $`90\%`$ of the emission lies within that radius, this temperature is a good estimate of the overall emission weighted value. Within that radius, there is no significant evidence of strong gradients. Nevertheless, if the temperature profile decreases beyond $`2^{}0.5R_\mathrm{V}`$, the measured temperature would be an overestimate of the mass weighted temperature (considered in the simulations). However, if the shape of the temperature profile does not evolve with $`z`$, the corresponding bias in the VT mass estimate would be the same at all redshifts and will not affect the following comparison with nearby cluster properties.

The central hydrogen density, $`n_{\mathrm{H},0}`$666$`n_{\mathrm{H},0}=0.85n_{\mathrm{e},0}`$, where $`n_{\mathrm{e},0}`$ is the central electronic density for an ionised plasma with the observed abundance and the corresponding gas mass within the virial radius, $`M_{\mathrm{gas}}`$, are derived from the best fit $`\beta `$–model parameters (central surface brightness, $`\beta `$ and core radius), assuming a cluster extent equal to the virial radius. The errors on these quantities, due to that on the surface brightness profile, are derived as in Elbaz et al. (elbaz (1995)). The error due to the uncertainty on the temperature, which appears in the emissivity factor, is negligible. We emphasize that the estimate of the gas mass does not require severe extrapolation of the data and is therefore robust. The cluster emission is detected virtually up to the virial radius for an EdS model ($`R_{\mathrm{det}}R_\mathrm{V}`$) and up to about $`0.75R_\mathrm{V}`$ for a $`\mathrm{\Lambda }\mathrm{CDM}`$ model (Table 1). The gas mass within that radius is already $`70\%`$ of the gas mass within $`R_\mathrm{V}`$.

The central hydrogen density, gas mass, virial radius, virial mass and gas mass fraction, $`f_{\mathrm{gas}}=M_{\mathrm{gas}}/M_\mathrm{V}`$, are given in Table 1 for the EdS and $`\mathrm{\Lambda }\mathrm{CDM}`$ models. The uncertainty on this last quantity is dominated by the uncertainty on the temperature, through the estimate of the virial mass.

#### 4.4.2 BM method

We also estimated the total mass using the hydrostatic equilibrium (HE) equation and the isothermal $`\beta `$–model (BM method):

$$\frac{G\mu m_\mathrm{p}M_\mathrm{V}}{2R_\mathrm{V}}=\frac{3\beta }{2}\frac{(R_\mathrm{V}/r_\mathrm{c})^2}{1+(R_\mathrm{V}/r_\mathrm{c})^2}\mathrm{k}T$$
(4)
where $`\beta `$ and $`r_\mathrm{c}`$ are the slope and core radius parameters of the gas distribution (derived Sec. 3.1).

The corresponding virial mass and radius, corresponding to the same density contrast as in the VT approach, are given in Table 1. Comparing Eq. 2 and Eq. 4 and neglecting the terms $`(r_\mathrm{c}/R_\mathrm{V})^2`$, the ratio of the virial masses (or radius) derived from the VT and BM method is:

$`{\displaystyle \frac{R_\mathrm{V}(BM)}{R_\mathrm{V}(VT)}}\left({\displaystyle \frac{3\beta }{2\beta _\mathrm{T}}}\right)^{1/2};{\displaystyle \frac{M_\mathrm{V}(BM)}{M_\mathrm{V}(VT)}}\left({\displaystyle \frac{3\beta }{2\beta _\mathrm{T}}}\right)^{3/2}`$ (5)
The VT and BM methods give very similar results (see Table 1). As compared to the VT estimate, the virial radius and mass, estimated using the BM method, are $`5\%`$ and $`15\%`$ higher, respectively. The central hydrogen density is unchanged. The gas mass within the virial radius is about $`5\%`$ higher (slightly larger integration region) and the gas mass fraction is decreased by about $`10\%`$.

#### 4.4.3 Comparison with the local gas mass fraction

The gas mass fraction of RX J1120.1+4318 can be compared to the $`90\%`$ confidence region derived by Arnaud & Evrard (arnaud99 (1999), hereafter AE) for hot ($`\mathrm{k}T>4\mathrm{keV}`$) nearby clusters using the same total mass estimate methods. AE found777These values were derived for an EdS model, but the local value is not sensitive to the cosmological parameters. The definition of the virial region is also slightly different but close enough to have no impact on our conclusion: our definition corresponds to a density contrast of 178 at z=0 for an EdS model, while the value quoted are estimated at a density contrast of 200: $`f_{\mathrm{gas}}=20.1[\mathtt{2.5},+\mathtt{3.0}]\%`$ using the VT method and $`f_{\mathrm{gas}}=21.5[\mathtt{3.5},+\mathtt{4.5}]\%`$ with the BM method. For both the VT and BM methods, the value derived for RX J1120.1+4318 is perfectly consistent with the AE local value for a $`\mathrm{\Lambda }\mathrm{CDM}`$ model, but is significantly lower for an EdS model.

However, the temperature of RX J1120.1+4318 ($`5.3\mathrm{keV}`$) is lower than the median temperature of the AE hot cluster sample ($`8\mathrm{keV}`$). We are thus not exactly comparing clusters of same temperature and there is some indication that the gas mass fraction varies with $`T`$. The empirical slope of the $`M_{\mathrm{gas}}`$$`T`$ relation derived by Neumann & Arnaud (neumann01 (2001)) corresponds to $`f_{\mathrm{gas}}T^{0.44}`$, if the classical $`M_\mathrm{V}`$$`T`$ relation is not modified. This variation is consistent with the variation, $`f_{\mathrm{gas}}T^{0.41\pm 0.16}`$, found by Mohr et al. (mohr99 (1999)) using the virial mass estimate. A variation of $`f_{\mathrm{gas}}`$ with $`T`$ is also readily apparent in Fig.3 of AE. Correcting the AE $`f_{\mathrm{gas}}`$ values for a possible bias of $`(5.3/8)^{0.44}`$, gives a local value of $`f_{\mathrm{gas}}=16.8\%`$ (VT approach), now in agreement with the RX J1120.1+4318 estimate for both a $`\mathrm{\Lambda }\mathrm{CDM}`$ model and an EdS model, although marginally in the latter case.

The gas mass fraction has been proposed from some time (Pen pen (1997); Rines et al. rines (1999)) as a novel distance indicator. The method is based on the assumed constancy of this quantity with $`z`$, while its estimate from X–ray data depends on the assumed angular distance. As Mohr et al. (mohr00 (2000)) pointed out, the region considered to compare $`f_{\mathrm{gas}}`$ at different redshifts must be fixed in terms of scaled radius, with the evolution of the virial radius properly taken into account, to avoid biases due to the variation of the gas mass fraction with physical radius. The comparison performed above, where the gas mass fractions are estimated within the virial radius, does not suffer from this bias888Note that the dependence of $`f_{\mathrm{gas}}`$ on $`\mathrm{\Omega }_0`$ and $`\mathrm{\Lambda }`$ in the present approach (given in Neumann & Arnaud (neumann00 (2000))) is indeed different than the $`d_\mathrm{A}^{3/2}`$ dependence derived by Pen (pen (1997)) and Rines et al. (rines (1999)). However, as illustrated above, the comparison is further complicated by possible variation of $`f_{\mathrm{gas}}`$ with cluster temperature. A much better understanding of the $`f_{\mathrm{gas}}`$$`T`$ relation, both in the local and distant Universe is required, before any conclusion can be drawn on the cosmological parameters.

### 4.5 The $`L_\mathrm{X}`$$`T`$ relation

#### 4.5.1 Cluster luminosity

We computed the bolometric luminosity, $`L_\mathrm{X}`$, within the virial radius (given Table 1). The observed count rate in the \[0.3-3\] keV band, obtained by integrating the surface brightness profile up to the detection radius, is converted to bolometric luminosity using the best fit MEKAL spectral model and the instrument response. The contribution beyond the detection region is estimated from the best fit $`\beta `$–model , but is totally negligible ($`0.7\%`$ for the $`\mathrm{\Lambda }\mathrm{CDM}`$ model). The error on $`L_\mathrm{X}`$ includes both the statistical error on the count rate and on the temperature.

#### 4.5.2 Evolution of the $`L_\mathrm{X}`$$`T`$ relation

As in AAN and in Sadat et al. (sadat (1998)), we divided this luminosity by the luminosity estimated from the local $`L_\mathrm{X}`$$`T`$ relation of AE and the cluster temperature: $`L_\mathrm{X}=\mathrm{1.15\; 10}^{45}(T/6\mathrm{keV})^{2.88}`$. The resulting factor $`C`$ is given in Table 1 for the two cosmological models considered. The contribution to the error due to the uncertainties on $`T`$ and $`L_\mathrm{X}`$ are summed quadratically in the log space.

This $`C`$ factor can be compared to the evolution of the normalization of the $`L_\mathrm{X}`$$`T`$ relation, $`C_{\mathrm{mod}}(z)=\mathrm{\Delta }_\mathrm{z}^{1/2}\left(1+z\right)^{3/2}`$, expected in the self-similar model (AAN). This theoretical value, estimated at the cluster redshift, is given in Table 1 for the EdS and $`\mathrm{\Lambda }\mathrm{CDM}`$ models. The ‘error’ bars corresponds to plus or minus the intrinsic scatter in $`L_\mathrm{X}`$ estimated by AE. The luminosity of RX J1120.1+4318 is in good agreement with the expected evolution for an EdS Universe, and in the upper range of the expected evolution for a $`\mathrm{\Lambda }\mathrm{CDM}`$ Universe. For an EdS model, the luminosity is also marginally consistent with no evolution.

#### 4.5.3 Comparison with recent works

Several groups quantified the evolution of the $`L_\mathrm{X}`$$`T`$ relation using ROSAT/ASCA data (Sadat et al. sadat (1998); Reichart et al. reichart (1999); Fairley et al. fairley (2000); AAN). There is a general consensus that no significant evolution is observed for an EdS model. However, as emphasized by AAN, both the luminosity estimate and the theoretical evolution depend on the assumed cosmology and a different conclusion is reached for a $`\mathrm{\Lambda }\mathrm{CDM}`$ model. In that case, AAN found a significant evolution, consistent with the self-similar model. Our results agree with this finding. They are also consistent with the results of Reichart et al. (reichart (1999)): they found an evolution of $`C=(1+z)^{(0.911.12q_0)_{1.22}^{+0.54}}`$ in the redshift range $`z<0.5`$ or $`C=2.0_{0.9}^{+0.6}`$ extrapolated at $`z=0.6`$ for the $`\mathrm{\Lambda }\mathrm{CDM}`$ model considered here. This is similar to the value, $`C=2.4_{0.5}^{+0.8}`$, derived for RX J1120.1+4318 .

On the other hand, Borgani et al. (2001b ), using recent CHANDRA data up to $`z=1.26`$, did not find any evidence of significant evolution, for the EdS model but also for the $`\mathrm{\Lambda }\mathrm{CDM}`$ model. We note first the large uncertainties on the CHANDRA data, especially above $`z=1`$. Moreover, there is a large dispersion, with data points above and below the local $`L_\mathrm{X}`$$`T`$ curve for the $`\mathrm{\Lambda }\mathrm{CDM}`$ model. To illustrate this point, let us consider the two clusters RX J0848+4456 ($`z=0.57,\mathrm{k}T=3.6\pm 0.5\mathrm{keV}`$) and MS1137.5+6625 ($`z=0.78,\mathrm{k}T=5.7_{0.6}^{+0.8}\mathrm{keV}`$). They are both relaxed clusters, with precise $`\mathrm{k}T`$ measurements, and are in the same redshift and temperature range than RX J1120.1+4318 . The $`C`$ value of MS1137.5+6625, $`C=2.0_{0.8}^{+0.7}`$, is perfectly consistent with the expected evolution: $`C_{\mathrm{mod}}(z)=1.88`$ at the cluster redshift and is similar to the value obtained for RX J1120.1+4318 . On the other hand, the luminosity of RX J0848+4456 is particularly low for its measured temperature: $`C=0.6_{0.2}^{+0.4}`$, corresponding to a negative evolution and of course inconsistent with the expected positive evolution of $`C_{\mathrm{mod}}(z)=1.6`$. The measured evolution is sensitive to systematic uncertainties on the temperature ($`CT^3`$). RX J0848+4456 was observed with ACIS-I (Holden et al. holden (2001)), for which there are still large calibration uncertainties, especially at low energies. This is a potential worry for temperature measurements of high $`z`$ clusters with this instrument. The temperature estimate of MS1137.5+6625, observed with ACIS-S (Borgani et al. 2001b ), is a priori more secure. MS1137.5+6625 was also in the AAN sample and the Chandra/ACIS-S data are consistent with the ASCA data. First cross calibration studies of ACIS-S and XMM/EPIC also show a good agreement between these instruments (Snowden snowden (2002)). Further cross-calibration studies and accurate measurements on a larger cluster sample are definitively required to assess if the presently observed dispersion is real.

### 4.6 Consistency of the results

In our comparisons with the nearby cluster properties, we found that the scaled $`EM`$ profile of RX J1120.1+4318 , as well as its gas mass fraction and luminosity, are consistent with the predictions of the self-similar model of cluster formation. Since all three observables are not independent - they are directly related to the gas density - the consistency of our results is a priori not surprising. It is, however not entirely trivial, as we discuss now.

A good match of the scaled $`EM`$ profile of RX J1120.1+4318 with the local profile, would imply that i) the shape of the gas profile is similar to that of nearby clusters ii) the $`EM`$ and thus gas density scales with $`z`$ as expected. Therefore its luminosity (integrated emission measure profile) and gas mass (integrated density profile) would naturally be found to follow the standard evolution. Since the standard evolution of the gas mass is the same than the evolution assumed to compute the virial mass, we would also find a gas mass fraction consistent with the local value. However, the emission measure, luminosity and possibly gas mass fraction also depend on the temperature. We emphasize that independent studies of the evolution of these quantities are thus expected to yield consistent results, but only if one considers a consistent scaling with $`T`$ for the three quantities. The present analysis is consistent in that sense. As mentioned above, the theoretical model used to scale the $`EM`$ profiles ($`M_{\mathrm{gas}}T^{1.94}`$) is consistent with the slope of the local $`L_\mathrm{X}`$$`T`$ relation (used to normalize the cluster luminosity) and corresponds to $`f_{\mathrm{gas}}T^{0.44}`$ (used to estimate the local gas mass fraction at the cluster temperature).

In spite of the general agreement with the self-similar model, some differences appear when comparing the results of the EdS and $`\mathrm{\Lambda }\mathrm{CDM}`$ models. The analysis of the gas mass fraction and scaled profiles would rather favor a $`\mathrm{\Lambda }\mathrm{CDM}`$ model. However for such a model, the luminosity is in the upper range of the expected evolution and a better agreement is observed for an EdS model. This is actually a direct consequence of the specific shape of RX J1120.1+4318 . The concentration of the gas distribution of RX J1120.1+4318 is slightly larger than the reference profile, yielding a higher luminosity as compared to the gas mass. These differences can also be understood by looking at the scaled profiles. As outlined above, the scaled profile of RX J1120.1+4318 is in a good agreement with the nearby profile at high radii for a $`\mathrm{\Lambda }\mathrm{CDM}`$ model. For an EdS model, the scaled profile lies more and more below the reference curve at $`r0.4R_\mathrm{V}`$, indicating a lower gas density than expected. Since the external regions contribute most to the gas mass ($`75\%`$ of the mass lies beyond $`0.3R_\mathrm{V}`$), this explains why the gas mass fraction for an EdS model is marginally too low, even as compared to the corrected AE local value, while the agreement is better for a $`\mathrm{\Lambda }\mathrm{CDM}`$ model. On the other hand, for an EdS model, the cluster scaled profile matches better the local profile in the central region (see Fig.5), which contributes most to the X–ray luminosity.

## 5 Core properties and non-gravitational effects

### 5.1 Radiative cooling

The XMM data suggest that RX J1120.1+4318 does not host a cooling flow. There is indeed no indication of a temperature drop in the central bin and a good fit is obtained with a $`\beta `$–model , down to the center, with no indication of central excess.

We computed the central cooling time, $`t_{\mathrm{cool}}`$ corresponding to the measured central density: $`t_{\mathrm{cool}}=\mathrm{2.9\; 10}^{10}\mathrm{yrs}\sqrt{T}/n_\mathrm{H}`$, where $`T`$ is in $`\mathrm{keV}`$ and $`n_\mathrm{H}`$ in $`\mathrm{cm}^3`$ (Sarazin sarazin (1986)). We obtained $`t_{\mathrm{cool}}\mathrm{1.\; 10}^{10}\mathrm{yrs}`$. The age of the Universe at the cluster redshift is $`t=\mathrm{6.5\; 10}^9\mathrm{yrs}`$ for the EdS model and $`t=\mathrm{1.1\; 10}^{10}\mathrm{yrs}`$ for the $`\mathrm{\Lambda }\mathrm{CDM}`$ model. The cooling time, even at the cluster center where it is close to its minimum value, is thus similar to or larger than the age of the Universe. We thus do not expect a strong cooling flow in the center, in spite of the apparent relaxed state of the cluster.

### 5.2 Entropy and non-gravitational heating

As mentioned in the introduction, the simplest self-similar model of cluster formation fails to explain all the observed properties of the nearby cluster population. A definitive evidence of breaking of self-similarity is the entropy excess (the “entropy floor”) detected in the core of cool systems with a baseline entropy of about $`70140\mathrm{keV}\mathrm{cm}^2`$ (Ponman et al. ponman (1999); Lloyd-Davies et al. lloyd (2000)). The origin of this break of similarity is not yet understood.

With the present data, we can for the first time study the entropy profile of a $`z=0.6`$ cluster, ie look at cluster entropy evolution. Note that RX J1120.1+4318 appears to be in a relaxed state and the entropy distribution is not likely to be affected by shocks induced during recent mergers. The cluster entropy profile, as a function of scaled radius, is plotted in Fig. 7 for the EdS and $`\mathrm{\Lambda }\mathrm{CDM}`$ cosmologies. The entropy $`S=T/n_\mathrm{e}^{2/3}`$ is estimated from the best fit $`\beta `$–model and the cluster mean temperature. The uncertainty on the profile is dominated by the uncertainty on the temperature distribution. Typical errors, corresponding to the temperature error in each bin of the temperature profile are indicated in the figure.

One observes a nearly constant entropy level of about $`140\mathrm{keV}\mathrm{cm}^2`$ within the central $`r<0.1R_\mathrm{V}`$ region and a rising entropy profile beyond that region. This shape is not surprising, the entropy profile shape reflects the shape of the inverse of the density profile, which has a core radius of about $`0.1R_\mathrm{V}`$. The core entropy is, however, in remarkable agreement with the entropy “floor” measured in nearby clusters. The same agreement was noted by Arabadjis et al. (arabadjis (2002)) in the case of EMSS 1358+6245 at $`z=0.33`$.

If confirmed on more clusters, this coincidence has potentially substantial implications for the physics of cluster formation. We first note that it is consistent with the expectation of the early pre-heating scenario, where the gas is pre-heated at a given entropy level before collapse (e.g. Borgani et al. 2001a ; Tozzi & Norman tozzi (2001)). At the beginning of the cluster formation, the entropy floor prevents shock heating and the collapse is adiabatic. The initial entropy is thus preserved in the core. When the mass of the system increases, the infall of incoming shells becomes supersonic and a shock regime begins (gravitational heating). The external profile follows the classical ’self-similar’ rising profile, as we observe. The core entropy of RX J1120.1+4318 suggests that the entropy floor was already established at high $`z`$: at least $`z0.6`$ and probably much before, since in a hierarchical scenario the core of the cluster must have collapsed earlier. Note also that the metal abundance found for this distant cluster is similar to the abundance observed in nearby clusters (e.g. de Grandi & Molendi degrandi (2001) and reference therein). If early galactic winds are responsible for the gas pre-heating (e.g. Kaiser kaiser91 (1991); Evrard & Henry evrard91 (1991)), we do expect that this is accompanied by an early enrichment. Finally, our data suggest that pre-heating in hot clusters has mostly an effect on the core properties up to high $`z`$. This might explain why AAN verified self-similarity of the $`EM`$ profile of hot clusters above $`0.1R_\mathrm{V}`$ up to z=0.8.

Radiative cooling has been also proposed to explain the entropy “floor”: it can remove low entropy gas from near the cluster center, triggering the inflow of higher entropy material (Pearce et al. pearce (2000)). RX J1120.1+4318 , contrary to most nearby relaxed clusters (and EMSS 1358+6245), is a cluster, for which radiative cooling has probably not yet affected the core properties (see above). This suggests that radiative cooling cannot be the dominant process in the establishment of the entropy floor.

## 6 Conclusion

XMM-Newton data allow us to measure, with unprecedented accuracy at a redshift of $`z=0.6`$, the gas and temperature distribution of the distant cluster RX J1120.1+4318 . The cluster has a regular spherical morphology, suggesting it is in a relaxed state. The cluster temperature is $`\mathrm{k}T=5.3\pm 0.5\mathrm{keV}`$ and there is no significant evidence of temperature radial gradient up to half the virial radius. The surface brightness profile, measured nearly up to the virial radius, is well fitted by a $`\beta `$–model with $`\beta =0.78_{0.04}^{+0.06}`$.

The RX J1120.1+4318 data reinforce the validity of the self-similar cluster formation model. For both an EdS and a $`\mathrm{\Lambda }\mathrm{CDM}`$ model, the scaled emission measure profile beyond $`0.1R_\mathrm{V}`$ is consistent with the nearby cluster reference profile, taking into account its intrinsic scatter. Consistently, the gas mass fraction is in agreement with the local value (although marginally for an EdS Universe), and the luminosity of the cluster, taking into account its temperature, is consistent with the expected evolution of the $`L_\mathrm{X}`$$`T`$ relation.

There is no evidence of a cooling flow at the cluster center, in spite of its apparent relaxed state. This is consistent with the estimated cooling time, larger than the age of the Universe at the cluster redshift, indicating that radiative cooling has not yet affected the cluster properties. The entropy profile shows a flat core with a central entropy of $`140\mathrm{keV}\mathrm{cm}^2`$, remarkably similar to the entropy floor observed in nearby clusters. This favors early pre-heating models for the establishment of this entropy floor.

A statistical study of large sample, as the SHARC cluster sample, covering a wide range of redshift and luminosity (thus temperature) is required to better constrain the physics of cluster formation, from the evolution of cluster properties. However, the present data already demonstrate that XMM provides the statistical quality required by such study.

###### Acknowledgements.

The present work is based on observations obtained with XMM-Newton, an ESA science mission with instruments and contributions directly funded by ESA Member States and the USA (NASA). We acknowledge financial support from the ”Programme National de Cosmologie” funded by INSU, CEA and CNES. AKR acknowledges partial financial support from the NASA XMM-Newton GO program. DJB acknowledges the support for NASA contract NAS8-39073

## Appendix A Background subtraction method

In this section we detail the background subtraction method. This method applies after screening for flares. We suppose we have a ‘template’ event file, for the background estimate, obtained by collecting several high latitude observations. In the following, the data are supposed to be corrected for vignetting effect, as described in Arnaud et al. (arnaud01 (2001)) and in section 2.2: each event detected with energy $`E`$ at location $`x,y`$ is weighted by a coefficient $`w(x,y,E)`$ which is the ratio of the effective area at position $`(x,y)`$ to the central effective area, for the energy $`E`$. The Cosmic X-ray Background (CXB) varies across the sky, but can be considered as uniform at the scale of $`30^{}`$, the size of the field of view. The non X–ray background (NXB) is not uniform in the FOV, but is not vignetted.

For the template file, the corrected count rate, measured at a given location $`(x,y)`$ and energy $`E`$, $`T(x,y,E)`$, is the sum of the CXB and NXB contributions:

$$T(x,y,E)=T_{\mathrm{CXB}}(x,y,E)+T_{\mathrm{NXB}}(x,y,E)w(x,y,E)$$
(6)
Since the data are corrected for vignetting effects, the CXB component is, apart from statistical fluctuations, uniform over the FOV and corresponds to the average CXB for the blank field observations:

$$T_{\mathrm{CXB}}(x,y,E)T_{\mathrm{CXB}}(E)$$
(7)
The NXB component is the NXB for the observations, $`T_{\mathrm{NXB}}(x,y,E)`$, multiplied by the weight factor.

Similarly for the observation data set, the count rate is:

$`O_(x,y,E)`$ $`=`$ $`S(x,y,E)+O_{\mathrm{CXB}}(x,y,E)`$ (8)
$`+O_{\mathrm{NXB}}(x,y,E)w(x,y,E)`$
with

$$O_{\mathrm{CXB}}(x,y,E)O_{\mathrm{CXB}}(E)$$
(9)
where $`S(x,y,E)`$ is the source contribution and $`O_{\mathrm{CXB}}(E)`$ corresponds to the CXB at the pointing position of the considered observation. Since the CXB varies across the sky, a priori $`O_{\mathrm{CXB}}(E)T_{\mathrm{CXB}}(E)`$. The quiescent NXB can be reasonably considered to have the same spectral and spatial characteristics for all observations but there is evidence of long term variations (D.Lumb lumb (2002)). We can thus write:

$$O_{\mathrm{NXB}}(x,y,E)QT_{\mathrm{NXB}}(x,y,E)$$
(10)
The normalization factor, $`Q`$, is estimated by considering the total count rate in the whole FOV in the high energy band, where any cosmic X–ray emission is negligible.

The background subtraction for each source product (spectrum or profile) is done in two steps. In a first step, we subtract the corresponding product of the blank field observation, obtained using the same energy and spatial selection, normalized by the factor $`Q`$. From Eq. 6, Eq. 7, Eq. 8, Eq. 9 and Eq.10:

$`O(x,y,E)QT(x,y,E)`$ $``$ $`S(x,y,E)+O_{\mathrm{CXB}}(E)`$ (11)
$`QT_{\mathrm{CXB}}(E)`$
and the NXB contribution is thus removed. The remaining CXB component, which can be considered uniform, can thus be estimated, using data outside the source region, where:

$`O(x,y,E)QT(x,y,E)`$ $``$ $`O_{\mathrm{CXB}}(E)QT_{\mathrm{CXB}}(E)`$ (12)
It can then subtracted in a second step.

To illustrate the method, let us consider that we want to estimate the source spectrum in a particular region, $`Reg`$:

$$S_{\mathrm{Reg}}(E)=\underset{\mathrm{Reg}}{}S(x,y,E)$$
(13)
In a first step, we extract the region spectrum from the observation data set. The spectrum extracted in the same region from the template file is then subtracted, after normalization by the factor $`A_{\mathrm{norm}}`$ above. The resulting spectrum is:

$`I_{\mathrm{Reg}}(E)`$ $`=`$ $`{\displaystyle \underset{\mathrm{Reg}}{}}O(x,y,E)Q{\displaystyle \underset{\mathrm{Reg}}{}}T(x,y,E)`$ (15)
$``$ $`S_{\mathrm{Reg}}(E)`$
$`+A_{Reg}(O_{\mathrm{CXB}}(E)QT_{\mathrm{CXB}}(E))`$
where $`A_{Reg}`$ is the surface of the extraction region. We are thus left with a residual X-ray background component, which is the estimate of the difference between the CXB spectrum in the observation and in the template file (multiplied by the normalization factor).

In a second step, we do the same for a region $`Reg^{}`$ outside the source. The resulting spectrum is:

$$I_{\mathrm{Reg}^{}}(E)A_{\mathrm{Reg}^{}}(O_{\mathrm{CXB}}(E)QT_{\mathrm{CXB}}(E))$$
(16)
The source spectrum can thus be estimated by subtracting this spectrum to the spectrum obtained in Step 1, after normalization to the size of the extraction region. From Eq. 15 and Eq. 16 one get:

$$I_{\mathrm{Reg}}(E)\frac{A_{\mathrm{Reg}}}{A_{\mathrm{Reg}^{}}}I_{\mathrm{Reg}^{}}(E)S_{\mathrm{Reg}}(E)$$
(17)
The same method can be applied to extract the source surface brightness profile, in a given energy band. We extract the surface brightness profile from the observation data set. The surface brightness profile extracted in the same energy band from the template file is then subtracted, after normalization by the factor $`Q`$ above. The resulting surface brightness profile is an estimate of:

$$I(r)S(r)+(S_{O_{\mathrm{CXB}}}QS_{T_{\mathrm{CXB}}})$$
(18)
where $`S(r)`$ is the source surface brightness profile and $`S_{O_{\mathrm{CXB}}}`$ and $`S_{T_{\mathrm{CXB}}}`$ are the CXB surface brightness in the energy band considered for the observation and template file respectively. The residual CXB contribution can be estimated from the measured profile, $`I(r)`$, in the region outside the source and subtracted.


format: text

text:
# An alternative approach to 𝑏→𝑠𝛾in the unconstrained MSSM

## 1 Introduction

The precision measurements of the inclusive radiative decay $`BX_s\gamma `$ provides an important benchmark for the Standard Model (SM) and New Physics (NP) models at the weak scale, such as low-energy supersymmetric (SUSY) models. In the SM, flavour changing neutral currents (FCNC) are forbidden at tree level. The first SM contribution to the $`bs\gamma `$transition appears at one loop level from the CKM flavour changing structure, showing the characteristic Cabibbo suppression. NP contributions to $`bs\gamma `$typically also arise at one loop, and in general can be much larger than the SM contributions if no mechanisms for suppressing the new sources of flavour violation exist (see $`^\mathrm{?}`$ for a complete set of references).

Experimentally, the inclusive $`BX_s\gamma `$ Branching Ratio (BR) has been measured by ALEPH, BELLE and CLEO, resulting in the current experimental weighted average $`BR(BX_s\gamma )_{exp}=(3.23\pm 0.41)\times 10^4`$, with new results expected shortly from BABAR and BELLE which could further reduce the experimental errors. Squeezing the theoretical uncertainties down to the 10% level has been (and still is) a crucial task. The SM theoretical prediction has been the subject of intensive theoretical investigation in the past several years, leading to the completion of the NLO QCD calculations. The original SM NLO calculation $`^\mathrm{?}`$ gives, for $`\sqrt{z}=m_c/m_b=0.29`$, the following result: $`BR(BX_s\gamma )_{SM}=(3.28\pm 0.33)\times 10^4`$. The main source of theoretical uncertainty is due to NNLO QCD ambiguities. In $`^\mathrm{?}`$ it was shown that using $`\sqrt{z}=0.22`$ (i.e. the running charm mass instead of the pole mass) is more justifiable and causes an enhancement of about 10% of the $`bs\gamma `$BR, leading to the current preferred value: $`BR(BX_s\gamma )_{SM}=(3.73\pm 0.30)\times 10^4`$. Although these theoretical uncertainties can be addressed only with a complete NNLO calculation, the SM value for the BR is in agreement with the experimental measurement within the $`12\sigma `$ level.

The general agreement between the SM theoretical prediction and the experimental results has provided useful guidelines for constraining the parameter space of models with NP present at the electroweak scale, such as the 2HDM and the minimal supersymmetric standard model (MSSM). In SUSY models superpartners and charged Higgs loops contribute to $`bs\gamma `$, with contributions that typically rival or exceed the SM one in size. For calculational ease only simplified MSSM scenarios (like cMSSM or MSSM with minimal flavour violation (MFV)) have usually been assumed. Netherveless, as the origin and dynamical mechanism of SUSY breaking are unknown, there is no reason a priori to expect that the soft parameters will be flavour-blind (or violate flavour in the same way as the SM). Of course, the kaon system has provided strong FCNC constraints for the mixing of the first and second generations, which severely limit the possibility of flavour violation in that sector. However the constraints for third generation mixings are significantly weaker, with $`bs\gamma `$usually providing the most stringent constraints.

A discussion of the $`bs\gamma `$process in the general unconstrained MSSM is in principle possible, but it is necessary to deal with two unavoidable problems: (i) a large number of free, essentially unconstrained parameters; (ii) the need to achieve a quite accurate cancellation between the sizeable different contributions (SM, Higgs, chargino, neutralino and gluino). In the following we’ll provide a particularly interesting and simple analysis of $`bs\gamma `$in the uMSSM.

## 2 Alternative solution to $`bs\gamma `$branching ratio

The low-energy effective Hamiltonian, at the bottom mass scale $`\mu _b`$, is defined as $`_{eff}=(4G_F/\sqrt{2})V_{tb}V_{ts}^{}_iC_i(\mu _b)Q_i(\mu _b)`$. The operators relevant to the $`bs\gamma `$process are:
$`Q_2`$ $`=`$ $`\overline{s}_L\gamma _\mu c_L\overline{c}_L\gamma ^\mu b_L,`$    $`Q_7`$ $`=`$ $`{\displaystyle \frac{e}{16\pi ^2}}m_b\overline{s}_L\sigma ^{\mu \nu }b_RF_{\mu \nu },`$    $`Q_8`$ $`=`$ $`{\displaystyle \frac{g_s}{16\pi ^2}}m_b\overline{s}_L\sigma ^{\mu \nu }G_{\mu \nu }^aT_ab_R,`$        $`Q_2^{}`$ $`=`$ $`\overline{s}_R\gamma _\mu c_R\overline{c}_R\gamma ^\mu b_R,`$    $`Q_7^{}`$ $`=`$ $`{\displaystyle \frac{e}{16\pi ^2}}m_b\overline{s}_R\sigma ^{\mu \nu }b_LF_{\mu \nu },`$    $`Q_8^{}`$ $`=`$ $`{\displaystyle \frac{e}{16\pi ^2}}m_b\overline{s}_R\sigma ^{\mu \nu }G_{\mu \nu }^aT_ab_L.`$  (1)   Effects of NP generally appear as modifications of the Wilson Coefficients (WC) $`C_{7,8}^{(_{^{}})}`$ associated to the operators $`Q_{7,8}`$ and their chirality conjugated ones. In the majority of the previous studies of the $`bs\gamma `$process, the main focus was to calculate the SM or NP contributions to $`C_{7,8}`$. The contributions coming from $`C_{7,8}^{}`$ have usually been neglected on the assumption that they are suppressed with respect to $`C_{7,8}`$ by the ratio $`m_s/m_b`$. While this is always valid in the SM, in the 2HDM or within specific MSSM scenarios (with MFV), this mass suppression can be absent in the uMSSM where the gluino contributions to $`C_{7,8}`$ and $`C_{7,8}^{}`$ are naturally of the same order$`^\mathrm{?}`$.

Therefore, in the following we present an alternative approach to $`bs\gamma `$in the uMSSM. We assume a particular scenario in which the total contribution to $`C_{7,8}`$ is negligible and the main contribution to the $`bs\gamma `$BR is given by $`C_{7,8}^{}`$. This “$`C_7^{}`$-dominated” scenario is realized when the chargino, neutralino, and gluino contributions to $`C_{7,8}`$ sum up in such a way as to cancel the W and Higgs contributions almost completely. In our opinion this situation does not require substantially more fine tuning than what is required in the usual MFV scenario, where conversely the NP contributions to $`C_{7,8}`$ essentially cancel between themselves.

In the following we will focus on the gluino contribution to $`C_{7,8}`$ and $`C_{7,8}^{}`$. There is only one gluino diagram that contributes to $`C_7`$ and $`C_7^{}`$, with the external photon line attached to the down-squark line, while two diagrams can contribute to the $`C_8`$ and $`C_8^{}`$ coefficients, as the gluon external line can be attached to the squark or the gluino lines. The one-loop gluino contributions to the $`C_{7,8}^{}`$ coefficients are given, at first and second order in the MI, respectively by:

$`C_7^{{}_{}{}^{}\stackrel{~}{g}}(1)`$ $`=`$ $`{\displaystyle \frac{8g_s^2}{3g^2}}{\displaystyle \frac{Q_d}{V_{tb}V_{ts}^{}}}{\displaystyle \frac{m_W^2}{\stackrel{~}{m}_D^2}}\left\{\delta _{23}^{RR}F_2^{(1)}(x_D^g){\displaystyle \frac{\stackrel{~}{m}_{\stackrel{~}{g}}}{m_b}}\delta _{23}^{RL}F_4^{(1)}(x_D^g)\right\},`$ (2)
$`C_8^{{}_{}{}^{}\stackrel{~}{g}}(1)`$ $`=`$ $`{\displaystyle \frac{g_s^2}{3g^2}}{\displaystyle \frac{Q_d}{V_{tb}V_{ts}^{}}}{\displaystyle \frac{m_W^2}{\stackrel{~}{m}_D^2}}\left\{\delta _{23}^{RR}F_{21}^{(1)}(x_D^g){\displaystyle \frac{\stackrel{~}{m}_{\stackrel{~}{g}}}{m_b}}\delta _{23}^{RL}F_{43}^{(1)}(x_D^g)\right\}`$ (3)
and

$`C_7^{{}_{}{}^{}\stackrel{~}{g}}(2)`$ $`=`$ $`{\displaystyle \frac{4g_s^2}{3g^2}}{\displaystyle \frac{Q_d}{V_{tb}V_{ts}^{}}}{\displaystyle \frac{m_W^2}{\stackrel{~}{m}_D^2}}{\displaystyle \frac{m_b(A_b\mu \mathrm{tan}\beta )}{\stackrel{~}{m}_D^2}}\left\{\delta _{23}^{RL}F_2^{(2)}(x_D^g){\displaystyle \frac{\stackrel{~}{m}_{\stackrel{~}{g}}}{m_b}}\delta _{23}^{RR}F_4^{(2)}(x_D^g)\right\},`$ (4)
$`C_8^{{}_{}{}^{}\stackrel{~}{g}}(2)`$ $`=`$ $`{\displaystyle \frac{g_s^2}{6g^2}}{\displaystyle \frac{Q_d}{V_{tb}V_{ts}^{}}}{\displaystyle \frac{m_W^2}{\stackrel{~}{m}_D^2}}{\displaystyle \frac{m_b(A_b\mu \mathrm{tan}\beta )}{\stackrel{~}{m}_D^2}}\left\{\delta _{23}^{RL}F_{21}^{(2)}(x_D^g){\displaystyle \frac{\stackrel{~}{m}_{\stackrel{~}{g}}}{m_b}}\delta _{23}^{RR}F_{43}^{(2)}(x_D^g)\right\}.`$ (5)
The gluino contributio to $`C_{7,8}`$ can be obtained exchanging L $``$ R in eqs. (25). In deriving eqs.(25) to the second order in the MI parameters, we have kept only the dominant term proportional to $`\mathrm{tan}\beta `$(the $`A_b`$ term is retained in the above expression for defining our convention for the $`\mu `$ term; see later) and neglected all of the other off-diagonal MIsaaaIn $`^\mathrm{?}`$ a complete derivation of the general results and conventions used in eqs. (2-5) is presented.. Clearly the dominant terms in eqs. (2-5) are those proportional to the gluino chirality flip, so that the gluino contribution to $`C_7`$ ($`C_7^{}`$) depends, at first order, only on the MI term $`\delta _{23}^{LR}`$($`\delta _{23}^{RL}`$). However, for large $`\mathrm{tan}\beta `$and $`\mu \stackrel{~}{m}_A`$, the second order MI terms in eqs. (4,5) can become comparable in size with the first order mass insertions. Thus, two different MI parameters are relevant in the L/R sectors: ($`\delta _{23}^{LR}`$, $`\delta _{23}^{LL}`$) and ($`\delta _{23}^{RL}`$, $`\delta _{23}^{RR}`$), contrary to common wisdom. To which extent the LL and RR MIs are relevant depends of course on the values chosen for $`\mu `$ and $`\mathrm{tan}\beta `$, but in a large part of the allowed SUSY parameter space they cannot in general be neglected. Moreover, the fact that the gluino WCs depend on two different MI parameters will have important consequences in the study of the $`bs\gamma `$CP asymmetry.

### 2.1 Single and multiple MI-dominance analysis

From eqs.(25), one can read (in MI language) the off-diagonal entries that are relevant to the gluino contribution to the $`C_{7,8}`$ and $`C_{7,8}^{}`$ WCs. Note that limits on $`\delta _{23}^{LR}`$$`O(10^2)`$ have previously been obtained in $`^\mathrm{?}`$. No stringent bound has been derived there for $`\delta _{23}^{LL}`$, as this term does not come, at lowest order, with the $`\stackrel{~}{m}_{\stackrel{~}{g}}/m_b`$ enhancement (see eqs.(2)). No limits were showed on $`\delta _{23}^{RL}`$and $`\delta _{23}^{RR}`$because the MI formula are symmetric in the $`LR`$ exchange and in the scenarios generally adopted in the literature the “opposite chirality” MIs are suppressed by a factor $`m_s/m_b`$ and so negligible. An analysis of the $`\delta _{23}^{RL}`$dependence has been performed in $`^\mathrm{?}`$, in which the $`W`$ contribution to $`C_{7,8}`$ was not set to zero (sometimes also Higgs and MFV chargino contributions to $`C_{7,8}`$ were included). Consequently their bounds on the down-squark off-diagonal MIs contributing to $`C_{7,8}^{}`$ are more stringent than the bounds we derive in our scenario, for which the total contribution to $`C_{7,8}`$ is assumed to be negligible. It is clearly only in the scenario we study that an absolute constraint on these MIs can be derived. Moreover no analysis on $`\delta _{23}^{LL}`$and $`\delta _{23}^{RR}`$was performed in $`^\mathrm{?}`$ as these contributions are not relevant in the small $`\mathrm{tan}\beta `$region, as can be seen from eqs.(4,5).

In fig. 1 we show the dependence of the $`bs\gamma `$BR on the MI terms $`\delta _{23}^{RL}`$and $`\delta _{23}^{RR}`$for different values of $`x_D^g=\stackrel{~}{m}_{\stackrel{~}{g}}^2/\stackrel{~}{m}_D^2`$ and for $`\mathrm{tan}\beta =20`$ and $`\mu =350`$ GeV. All the other off-diagonal entries in the down-squark mass matrix are assumed to vanish. “Individual” limits $`\delta _{23}^{RL}`$$`<10^2`$ and $`\delta _{23}^{RR}`$$`<1.5\times 10^1`$ can be obtained respectively from the left and right side plot of fig. 1. Of course, the required cancellation of the total $`C_{7,8}`$ contribution may in general need non-vanishing off-diagonal entries of the up and down squark mass matrices. However, the specific values of these entries do not significantly affect the absolute limits on the $`\delta _{23}^{RL}`$and $`\delta _{23}^{RR}`$MIs shown in fig. 1. As expected from eqs.(4,5), the bounds obtained for $`\delta _{23}^{RR}`$are strongly dependent on the product $`\mu \mathrm{tan}\beta `$. More stringent bounds on $`\delta _{23}^{RR}`$can be obtained for larger $`\mathrm{tan}\beta `$. For $`\mathrm{tan}\beta `$$`>35`$ the bounds on $`\delta _{23}^{RR}`$can become as stringent as the $`\delta _{23}^{RL}`$bounds.

A general analysis of the gluino contribution to $`C_{7,8}^{}`$ depends simultaneously on both the $`\delta _{23}^{RL}`$and $`\delta _{23}^{RR}`$MIs. For a complete specification of our scenario the only other free parameters that need be fixed are the ratio between the gluino mass and the common down-squark mass, $`\stackrel{~}{m}_{\stackrel{~}{g}}/\stackrel{~}{m}_D`$, the product $`\mu `$ $`\mathrm{tan}\beta `$, and the relative phase between $`\delta _{23}^{RL}`$and $`\delta _{23}^{RR}`$. The influence of all the other down-sector squark matrix off-diagonal entries and MSSM parameters in the $`C_{7,8}^{}`$ sector can safely be neglected. In fig. 2 (left) we show the $`1\sigma `$ experimentally allowed region in the ($`\delta _{23}^{RL}`$, $`\delta _{23}^{RR}`$) parameter space for a specific choice of $`\stackrel{~}{m}_{\stackrel{~}{g}}/\stackrel{~}{m}_D=350/500`$, $`\mu =350`$ GeV, and for three different values of $`\mathrm{tan}\beta `$=$`3,20`$ and $`35`$. For $`\delta _{23}^{RL}`$or $`\delta _{23}^{RR}`$vanishing, one obtains the regions depicted in fig. 1. Larger regions in the ($`\delta _{23}^{RL}`$, $`\delta _{23}^{RR}`$) parameter space are obtained when both the MIs take non-vanishing values. It is clear that no absolute limit can be derived for the two MIs simultaneously. The values ($`\delta _{23}^{RR}`$, $`\delta _{23}^{RR}`$) $`(1,0.1)`$ are, for example, possible for $`\mathrm{tan}\beta `$$`=35`$. In fact, as can be seen in fig. 2 (left), there is always a “flat direction” where large values of $`\delta _{23}^{RL}`$and $`\delta _{23}^{RR}`$can be tuned in such a way that the gluino contribution to $`C_{7,8}^{}`$ is consistent with the experimental bound. This flat direction clearly depends on the chosen values for $`\stackrel{~}{m}_{\stackrel{~}{g}}/\stackrel{~}{m}_{\stackrel{~}{q}}`$ and $`\mu `$ $`\mathrm{tan}\beta `$. The presence of this particular direction is explained by the fact that we are allowing complex off-diagonal entries. Hence the relative phase between $`\delta _{23}^{RL}`$and $`\delta _{23}^{RR}`$can be fixed in such a way that the needed amount of cancellation can be obtained between the first and second order MI contribution. In the notation used in eqs.(25) the line of maximal cancellation is obtained for $`\phi =\mathrm{arg}[`$$`\delta _{23}^{RL}`$$`\delta _{23}^{RR}`$$`]=\pm \pi `$.

### 2.2 CP asymmetry and LR asymmetry

In addition to the $`bs\gamma `$BR, the experimental collaborations will provide in the coming years more precise measurements of the $`bs\gamma `$CP asymmetry. The present experimental value gives, at $`90\%`$ CL, the range $`0.27<𝒜_{CP}(bs\gamma )<0.10`$ which is still too imprecise to provide useful tests for NP, although the measurement is expected to be upgraded soon.

The only flavour-violating and CP-violating source in the SM (and MFV scenarios) is given by the CKM matrix, which results in a very small prediction for the CP asymmetry. In the SM an asymmetry of approximatively $`0.5\%`$. If other sources of CP violation are present, a much larger CP asymmetry could be produced. In our $`C_7^{}`$-dominated scenario, one can derive the following approximate relation for the CP asymmetry, in terms of the $`\delta _{23}^{RL}`$and $`\delta _{23}^{RR}`$MIs:

$$𝒜_{CP}(bs\gamma )=\frac{4}{9}\alpha _s(\mu _b)\frac{\mathrm{Im}\left[C_7^{}C_8^{}\right]}{|C_7^{}|^2}k(x_D^g)\left(\frac{m_b\mu \mathrm{tan}\beta }{\stackrel{~}{m}_D^2}\right)|\delta _{23}^{RL}\delta _{23}^{RR}|\mathrm{sin}\phi ,$$
(6)
where $`\phi `$ is the relative phase between $`\delta _{23}^{RL}`$and $`\delta _{23}^{RR}`$as previously defined. One can immediately note that if only one MI is considered, the CP asymmetry is automatically zero. In addition a non-vanishing phase in the off-diagonal down-squark mass matrix is necessary. No sensitive bounds on this phase can be extracted from EDMs in a general flavour-violating scenario.

In fig. 2 (right), we show the results obtained for the BR and CP asymmetry in which $`\delta _{23}^{RL}`$, $`\delta _{23}^{RR}`$and the relative phase $`\phi `$ are varied arbitrarily for fixed value of $`\stackrel{~}{m}_{\stackrel{~}{g}}/\stackrel{~}{m}_{\stackrel{~}{q}}=350/500`$ and $`\mathrm{tan}\beta `$$`=35`$. The full vertical lines represent the $`1\sigma `$ region experimentally allowed by the $`bs\gamma `$BR measurements. It is possible, using $`C_{7,8}^{}`$ alone, to saturate the $`bs\gamma `$measured BR and at the same time have a CP asymmetry even larger than $`\pm 10\%`$, the sign of the asymmetry being determined by the sign of $`\mathrm{sin}\phi `$. As shown in fig. 2 (right), in the relevant BR range the CP asymmetry range is constant. No strong dependence from $`\mathrm{tan}\beta `$, in the large $`\mathrm{tan}\beta `$region, is present. The points with large asymmetry ($`>5\%`$) lie in the “flat direction” observed in fig. 2 and they have almost $`\phi \pm \pi `$ (obviously for $`\phi =\pm \pi `$ the CP asymmetry vanishes). The explanation of this fact is the following. The numerator is proportional to $`\mathrm{sin}\phi `$ and so goes to 0 as $`\phi `$ approaches $`\pm \pi `$. However, at the same time it is enhanced for large MI values. This happens when the flat direction condition is (almost) satisfied. Here, in fact, a cancellation between the two (large) MI terms takes place, providing the enhancement of the CP asymmetry as the denominator remains practically constant, fixed by the allowed experimental measurement on the BR. Note also that for parameter values outside the flat direction condition a CP asymmetry of a few per cent can still be observed, about ten times bigger than the SM prediction. In our scenario even smaller values of the CP asymmetry can be obtained, e.g. if one of the two off-diagonal entries is negligible, or the two MIs are “aligned”.

A possible method for disentangling the relative contributions to the $`bs\gamma `$BR from the $`Q_7`$ and $`Q_7^{}`$ operators utilizes an analysis of the photon polarization. For simplicity, let us define the following “theoretical” LR asymmetry at LO:

$$𝒜_{LR}(bs\gamma )=\frac{BR(bs\gamma _L)BR(bs\gamma _R)}{BR(bs\gamma _L)+BR(bs\gamma _R)}=\frac{|C_7(\mu _b)|^2|C_7^{}(\mu _b)|^2}{|C_7(\mu _b)|^2+|C_7^{}(\mu _b)|^2},$$
(7)
which could in principle disinguish between $`C_7`$ and $`C_7^{}`$ dominated scenarios. Here L,R is the polarization of the external photon. This quantity is related to the quark chiralities of the $`Q_7,Q_7^{}`$ operators. Such a measurement is not yet available, as only the average quantity $`BR(bs\gamma _L)+BR(bs\gamma _R)`$ is reported experimentally. In the SM case, and in general in all the MFV and mSUGRA scenarios, only the $`C_7`$ coefficient gives a non-negligible contribution to the $`bs\gamma `$BR, in such a way that $`𝒜_{LR}(bs\gamma )=1`$. In our scenario, where the total contribution to $`C_7`$ is negligible, one obtains $`𝒜_{LR}(bs\gamma )=1`$. In any other uMSSM scenario, any LR asymmetry between $`1`$ and $`1`$ is allowed. Consequently, a measurement of $`A_{LR}(bs\gamma )`$ different from one will be a clear indication of physics beyond the SM with a non-minimal flavour structure.

## 3 Conclusions

We have discussed an alternative explanation of the $`bs\gamma `$BR in the unconstrained MSSM. We analyzed in particular the gluino contribution to the WC $`C_7^{}`$ associated with the chirality operator $`Q_7^{}`$. We show that this coefficient arises mainly from two off-diagonal entries: $`\delta _{23}^{RL}`$and $`\delta _{23}^{RR}`$. For scenarios where the $`C_{7,8}`$ contributions to $`bs\gamma `$are small (i.e. for regions in the MSSM parameter space where Ws, Higgs, chargino and gluino contributions to $`C_{7,8}`$ tend to cancel each other), $`C_{7,8}^{}`$ provides the dominant effect. We derived absolute bounds separately on each of these coefficients. We then described the allowed region of ($`\delta _{23}^{RL}`$, $`\delta _{23}^{RR}`$) parameter space, as a function of $`\mathrm{tan}\beta `$. We observed that (for a fixed ratio $`\stackrel{~}{m}_{\stackrel{~}{g}}/\stackrel{~}{m}_{\stackrel{~}{q}}`$ and for each chosen value of $`\mu \mathrm{tan}\beta `$), there exists a “flat direction” where large (even $`O(1)`$) off-diagonal entries are allowed. For the majority of parameter space the CP asymmetry is less than $`5\%`$, even if asymmetries as large as $`20\%`$ can be obtained along these “flat directions”. Finally, we suggested that the measure of the LR asymmetry could help to disentangle the $`C_7`$ from the $`C_7^{}`$ contribution to the $`bs\gamma `$BR. Any $`𝒜_{LR}(bs\gamma )1`$ would be an irrefutable proof, not only of physics beyond the SM, but also it would indicate the existence of a non-minimal flavour violation structure of the down-squark mass matrix. It would be very interesting if such a quantity could be measured. One implication of our analysis is that previous results on MSSM parameters, including constraints on the “sign of $`\mu `$bbbThe relative sign between the parameters $`\mu `$ and $`A_t`$. See $`^\mathrm{?}`$ for a discussion about the $`g2`$ “sign of $`\mu `$”. which are more model-dependent than has generally been assumed.

## References


text:
# Constraining Fourth Generation with 𝐵→𝑋_𝑠𝛾

## 1 Introduction

Today, despite the success of the Standard Model (SM), from the theoretical point of view, it is incomplete. Number of generations of fermions can be mentioned as one of the open problems of SM, for which we do not have a clear argument to restrict the SM to three known generations. Mass of the extra generations, if ever exists, can be extracted from the measurements of neutrino experiments, which set a lower bound for extra generations ($`m_{\nu _4}>45GeV`$) .

The idea of generalizing SM is not a new one. Probable effects of extra generations was studied in many works . Generalizations of the SM can be used to introduce a new family, which was performed previously . Using similar techniques, one can search fourth generation effects in B meson decays. The existing electroweak data on the $`Z`$–boson parameters, the $`W`$ boson and the top quark masses excluded the existence of the new generations with all fermions heavier than the $`Z`$ boson mass , nevertheless, the same data allows few extra generations, if neutral leptons have masses close to $`50GeV`$.

$`BX_s\gamma `$ is one of the most promising areas in search of the fourth generation, via its indirect loop effects, which was performed previously . This decay is one of the most appropriate candidates to be searched in the extensions of SM, since we have solid experimental and theoretical background for the process under consideration.

In this work we study the contribution of the fourth generation in the rare $`BX_s\gamma `$ decay, to obtain constrains on the parameter space of the fourth generation. Our basic assumption is to fill the gap between theoretical and experimental results of $`BX_s\gamma `$, with the fourth generation. Of course, due to the mentioned assumption, decay width will change at the order of difference between theoretical and experimental results, however, predicted CP asymmetry is interesting when SM contribution is neglected. As it is well known, new physics effects can manifest themselves through the Wilson coefficients and their values can be different from the ones in the SM , as well as through the new operators . Note that the inclusive $`BX_s\gamma `$ decay have already been studied with the inclusion of the fourth generation  to constrain $`V_{t^{}b}^{}V_{t^{}s}`$. The restrictions of the parameter space of nonstandard models based on LO analysis are not as sensitive as in the case of NLO analysis. Therefore we preferred to work at NLO, for the decay under consideration.

On the experimental side, values related with the $`BX_s\gamma `$ are well known. First measurement of the $`BX_s\gamma `$ was performed by CLEO collaboration, leading to CLEO branching ratio

$$BX_s\gamma =(2.32\pm 0.57\pm 0.35)\times 10^4.$$
(1)
In 1999, CLEO has presented an improved result

$$BX_s\gamma =(3.15\pm 0.35\pm 0.32\pm 0.26)\times 10^4.$$
(2)
The errors are statistical, systematic, and model dependent respectively. The rate measured by ALEPH  is consistent with the CLEO measurement. There exists also results of BELLE with a larger central value :

$$BX_s\gamma =(3.37\pm 0.53\pm 0.42\pm 0.54)\times 10^4.$$
(3)
Observing CP asymmetry in the decay $`BX_s\gamma `$ is interesting, presented by CLEO collaboration recently

$$A_{CP}(BX_s\gamma )=(0.079\pm 0.108\pm 0.022)\times (1.0\pm 0.030),$$
(4)
for which prediction of the SM is $`0.6\%`$ .

On the theoretical side, situation within and beyond the SM is well settled. A collective theoretical effort has led to the practical determination of $`BX_s\gamma `$ at the NLO, which was completed recently, as a joint effort of many different groups (,, , , ,). For a recent review, to complete the computation of NLO QCD corrections, we refer to ref.  and references therein. It is necessary to have precise calculations also in the extensions of the SM, which was performed for certain models . With the appearance of more accurate data we will be able to provide stringent constraints on the free parameters of the models beyond SM. We can state that, the aim of the present paper is to obtain such constraints when the fourth generation is considered.

The paper is organized as follows. In section 2, we present the necessary theoretical expressions for the $`BX_s\gamma `$ decay in the SM with four generations, where we investigated the effect of introducing fourth generation at different scales upon branching ratio and CP asymmetry. Section 3 is devoted to the numerical analysis and our conclusion.

## 2 Theoretical results

We use the framework of an effective low-energy theory, obtained by integrating out heavy degrees of freedoms, which in our case W-boson and top quark and an additional $`t^{}`$. Mass of the $`t^{}`$ is at the order of $`m_W`$. In this approximation the effective Hamiltonian relevant for $`bs\gamma `$ decay reads

$`_{eff}`$ $`=`$ $`{\displaystyle \frac{4G_F}{\sqrt{2}}}V_{ts}^{}V_{tb}{\displaystyle \underset{i=1}{\overset{8}{}}}C_i(\mu )O_i(\mu ),`$ (5)
where $`G_F`$ is the Fermi coupling constant $`V`$ is the Cabibbo-Kobayashi-Maskawa (CKM) quark mixing matrix, the the full set of the operators $`O_i(\mu )`$ and the corresponding expressions for the Wilson coefficients $`C_i(\mu )`$ in the SM can be found in ().

In the model under consideration, the fourth generation is introduced in a similar way the three generations are introduced in the SM, no new operators appear and clearly the full operator set is exactly the same as in SM . The fourth generation changes values of the Wilson coefficients $`C_7(\mu _W),C_8(\mu _W)`$, via virtual exchange of the fourth generation up quark $`t^{}`$ at matching scale. With the definition $`\lambda _i=V_{is}^{}V_{ib},i=\{u,c,t,t^{}\}`$, the above mentioned Wilson coefficients, can be written in the following form

$`C_7^{eff}(\mu _W)`$ $`=`$ $`C_7^{SM}(\mu _W)+{\displaystyle \frac{\lambda _t^{}}{\lambda _t}}C_7^{New}(\mu _W),`$
$`C_8^{eff}(\mu _W)`$ $`=`$ $`C_8^{SM}(\mu _W)+{\displaystyle \frac{\lambda _t^{}}{\lambda _t}}C_8^{New}(\mu _W),`$ (6)
where the last terms in these expressions describe the contributions of the $`t^{}`$ quark to the Wilson coefficients and $`V_{t^{}s}^{}`$ and $`V_{t^{}b}`$ are the two elements of the $`4\times 4`$ Cabibbo–Kobayashi–Maskawa (CKM) matrix. The explicit forms of the $`C_i^{New}`$ can easily be obtained from the corresponding Wilson coefficient expressions in SM by simply substituting $`m_tm_t^{}`$ (see ). Neglecting the $`s`$ quark mass we can define the Wilson coefficients at the matching scale, where the LO functions are :

$`C_7^{SM}`$ $`=`$ $`{\displaystyle \frac{x}{24}}{\displaystyle \frac{8x^3+3x^2+12x7+(18x^212x)\mathrm{ln}x}{(x1)^4}}`$
$`C_8^{SM}`$ $`=`$ $`{\displaystyle \frac{x}{8}}{\displaystyle \frac{x^3+6x^23x26x\mathrm{ln}x}{(x1)^4}},`$ (7)
where ($`x=m_t^2/m_W^2`$).

In the calculations we used the NLO theoretical expressions, and different experimental values to constraint the $`\lambda _t^{}`$ paramater. Since extended models are very sensitive to NLO corrections, we used the NLO expression for the branching ratio of the radiative decay $`BX_s\gamma `$, which has been presented in ref. :

$`BR(BX_s\gamma )`$ $`=`$ $`BR(BX_ce\overline{\nu }_e)\left|{\displaystyle \frac{V_{ts}^{}V_{tb}}{V_{cb}}}\right|^2{\displaystyle \frac{6\alpha _e}{\pi f(z)\kappa (z)}}{\displaystyle \frac{\overline{m}_b^2(\mu _b)}{m_b^2}}`$ (8)
$`\times \left(|D|^2+A\right)\left(1{\displaystyle \frac{\delta _{SL}^{NP}}{m_b^2}}+{\displaystyle \frac{\delta _\gamma ^{NP}}{m_b^2}}+{\displaystyle \frac{\delta _c^{NP}}{m_c^2}}\right).`$
Explicit forms of virtual, bremsstrahlung and non-perturbative parts of Eq. (8) can be found in  and references therein. In the numerical analysis we obtained $`BX_s\gamma `$ branching ratio in the Standard Model $`BR(BX_s\gamma )=(3.48\pm 0.33)\times 10^4`$, which remains in agreement with the previous literature. But we considered only the central value in our analysis, with the expectation of absorbing errors into the different experimental values.

To obtain quantitative results we need the value of the fourth generation CKM matrix element $`\lambda _t^{}`$. For this aim following , we will use the experimental results of the decays $`BR(BX_s\gamma )`$ and determine the fourth generation CKM factor $`\lambda _t^{}`$ When we consider the possible effects of the fourth generation, we demanded the theoretical value to be equal to the experimental values presented in the previous section. Which can be expressed as

$`BR(BX_s\gamma )_{4th}`$ $`=`$ $`\{2.66,3.15,3.37\}.`$ (9)
Theoretical results of the branching ratio for $`m_t^{}=75,\mathrm{},500GeV`$ values are obtained as as function of $`\lambda _t^{}`$. Notice that in the expressions related with $`BR(BX_s\gamma )_{4th}`$, theoretical and experimental results are multiplied by a factor of $`10^4`$ . For instance when we chose $`m_t^{}=75`$ GeV, and use the approach of Eq. (6):

$`BR(BX_s\gamma )_{4th}`$ $`=`$ $`0.654502+\mathrm{\hspace{0.17em}6.69962}\lambda _t^{}+\mathrm{\hspace{0.17em}20.3501}\lambda _t^{^{}}^2+`$
$`0.396254|0.305738\mathrm{\hspace{0.17em}1.87828}\lambda _t^{}|^2+`$
$`23.9926|(0.340878\mathrm{\hspace{0.17em}0.0154077}i)(\mathrm{\hspace{0.17em}1.64283}+\mathrm{\hspace{0.17em}0.05443}i)\lambda _t^{})|^2.`$
When $`\lambda _t^{}`$ is neglected branching ratio reduces to the re-scaled central value 3.48 of SM prediction. During the calculations we obtained similar expressions for different $`m_t^{}`$ values. It suffices to present the case of a very heavy quark, for $`m_t^{}=500GeV`$ :

$`BR(BX_s\gamma )_{4th}`$ $`=`$ $`0.654502+\mathrm{\hspace{0.17em}20.9868}\lambda _t^{}+\mathrm{\hspace{0.17em}198.863}\lambda _{t^{}}^{}{}_{}{}^{2}+`$
$`0.396254|0.305738\mathrm{\hspace{0.17em}5.95621}\lambda _t^{}|^2+`$
$`23.9926|(0.340878\mathrm{\hspace{0.17em}0.0154077}i)(\mathrm{\hspace{0.17em}5.1664}+\mathrm{\hspace{0.17em}0.11899}i)\lambda _t^{})|^2.`$
In the numerical analysis, as a first step, $`\lambda _t^{}`$ is assumed real and constraints are obtained as a function of mass of the extra generation top-quark $`m_t^{}`$, and the values are presented in tab. (1) and can be obtained from fig.(1). Those values can also be extracted from the figures $`(3,4,and5)`$ (a) where the solution is the intersection point on the $`BR_\gamma =1`$ line . Notice that in the figures we normalized branching ratio to 1, using the experimental values 2.66, 3.15 and 3.37 respectively, hence $`\lambda _t^{}`$ values can be obtained from the intersection points and this is true for all $`figures`$ except the ones related with $`A_{CP}`$ , presented in the following subsections.

We also performed a very similar analysis for introducing the fourth generation effects at the $`\mu _b`$ scale to see the difference between the previous results. Following  it can be written as follows:

$`C_7^{eff}(\mu _b)`$ $`=`$ $`C_7^{SM}(\mu _b)+{\displaystyle \frac{\lambda _t^{}}{\lambda _t}}C_7^{New}(\mu _b),`$
$`C_8^{eff}(\mu _b)`$ $`=`$ $`C_8^{SM}(\mu _b)+{\displaystyle \frac{\lambda _t^{}}{\lambda _t}}C_8^{New}(\mu _b),`$ (12)
Using Eq.(9), and demanding theretical results to be equal to the experimental results again, we obtained following expression for $`m_t^{}=75`$:

$`BR(BX_s\gamma )_{4th}`$ $`=`$ $`\mathrm{\hspace{0.17em}0.691542}+\mathrm{\hspace{0.17em}23.9926}\times `$ (13)
$`|(0.3408780.0154077i)(8.13033+0.423786i)\lambda _t^{}|^2`$
As another example for $`m_t^{}=500`$ we obtained

$`BR(BX_s\gamma )_{4th}`$ $`=`$ $`\mathrm{\hspace{0.17em}0.691542}+\mathrm{\hspace{0.17em}23.9926}\times `$ (14)
$`|(0.3408780.0154077i)(12.4566+0.484564i)\lambda _t^{}|^2.`$
It is interesting to notice that, if we assume $`\lambda _t^{}`$ can have imaginary parts, experimental values can also be satisfied. This case is presented with a graphical solution in figure (2) for $`m_t=75`$ and the decomposition $`\lambda _t^{}=\lambda _t^{real}+i\lambda _t^{imaginary}`$. When $`\lambda _t^{}`$ is assumed real constraints can be extracted from figures (3,4, and 5) (b) on the normalized line. Real and imaginary parts or this approach is presented in tables (2) and (3) respectively.

In order to check the consistency of the results of present work one can demand $`\lambda _t^{}`$ values to satisfy the unitarity condition. If we impose the unitarity condition of the CKM matrix we then have

$`\lambda _u+\lambda _c+\lambda _t+\lambda _t^{}=0.`$ (15)
With the values of the CKM matrix elements in the SM , the sum of the first three terms in Eq. (15) is about $`7.6\times 10^2`$, where the error in sum of first three terms is about $`\pm 0.6\times 10^2`$. By substituting the values of $`\lambda _t^{}`$ from tables 1 and 2, we observe that the sum of the four terms on the left–hand side of Eq. (15) may get very close to zero or diverge from the prediction of SM . When $`\lambda _t^{}`$ is very close to the sum of the first three terms, but with opposite sign, this is a very desirable result. Using table 2 for $`m_t^{}=100`$ GeV and the experimental branching ratio $`3.37\times 10^{}4`$, our prediction reads $`\lambda _t^{}^{(V)}=7.56\times 10^2`$. On the other hand the same prediction contains an imaginary part $`0.19i\times 10^2`$, which may be absorbed within the error range. In other words, results presented in table (2) satisfy unitarity constrain to a good extend. Nevertheless, it is a matter of taste to accept or reject $`\lambda _t^{}`$ values, according to unitarity condition. Because, it is possible that , existence of extra generations can affect present constraints on $`V_{CKM}`$ to a certain extend, and hence, constraints may get relaxed , which is beyond the scope of this work. From this respect it is hard to claim that all results presented here can satisfy unitarity. Nevertheless, in order to give the full picture, we did not exclude regions that violates unitarity.

### 2.1 Differences in the definitions of $`\lambda _t^{}`$

In order to explain the difference, on the results of the two different approaches given in Eq. (6) and Eq. (12) or tables (1) and (2), we can perform the analysis in LO, to extract the value of the fourth generation CKM matrix element $`\lambda _t^{}`$. Following , one can use the experimental results of the decays $`BR(BX_s\gamma )`$ and $`BR(BX_ce\overline{v}_e)`$, as in . In order to reduce the uncertainties arising from $`b`$ quark mass, consider the following ratio

$$R=\frac{Br(BX_s\gamma )}{Br(BX_ce\overline{v}_e)}.$$
(16)
In leading logarithmic approximation, for low energy scale approximation ratio can be written as

$$R=\alpha _m|C_7^{\mathrm{eff}}(\mu _b)|^2$$
(17)
where $`\alpha _m=\frac{|V_{ts}^{}V_{tb}|^2}{|V_{cb}|^2}\frac{6\alpha }{\pi f(\widehat{m}_c)\kappa (\widehat{m}_c)}`$ , the phase factor $`f(\widehat{m}_c)`$ and $`𝒪(\alpha _s)`$, QCD correction factor $`\kappa (\widehat{m}_c)`$ of $`bcl\overline{\nu }`$ are given in ref.. Using the LO definition of $`C_7^{eff}(\mu _b)`$ one can write ,

$`C_7^{eff}(\mu _b)=\eta ^{16/23}C_7^{eff}(\mu _W)+{\displaystyle \frac{8}{3}}(\eta ^{14/23}\eta ^{16/23})C_8^{eff}(\mu _W)+C_2^{eff}(\mu _W){\displaystyle \underset{i=1}{\overset{8}{}}}h_i\eta ^{a_i}`$ (18)
for the present purpose, which can be written as

$`C_7^{eff}(\mu _b)=\eta _1C_7^{eff}(\mu _W)+\eta _2C_8^{eff}(\mu _W)+\eta _3C_2^{eff}(\mu _W)`$ (19)
When the effect of 4-generation it is defined as Eq. (12)

$`C_{7,8}^{eff}(\mu _b)`$ $`=`$ $`C_{7,8}^{SM}(\mu _b)+{\displaystyle \frac{\lambda _t^{}}{\lambda _t}}C_{7,8}^{New}(\mu _b),`$
solution of Eq. (17) for $`\lambda _t^{}`$ can be written as follows

$$\lambda _t^{}^\pm =\left[\pm \sqrt{\frac{R}{\alpha _m}}C_7^{\mathrm{SM}}(\mu _b)\right]\frac{\lambda _t}{C_7^{\mathrm{New}}(\mu _b)}.$$
(21)
whereas in the case of the following approach ( Eq. (12))

$`C_{7,8}^{eff}(\mu _W)`$ $`=`$ $`C_{7,8}^{SM}(\mu _W)+{\displaystyle \frac{\lambda _t^{}}{\lambda _t}}C_{7,8}^{New}(\mu _W),`$
Eq. (21) is modified into the following form

$$\lambda _t^{}^\pm =\left[\pm \sqrt{\frac{R}{\alpha _m}}C_7^{\mathrm{SM}}(\mu _b)\right]\frac{\lambda _t}{\left[\eta _1C_7^{\mathrm{New}}(\mu _b)+\eta _2C_8^{\mathrm{New}}(\mu _b)\right]}.$$
(23)
This analysis can also be performed for NLO expressions. By comparing Eq. (21) and Eq. (23) the difference in tables (1) and and (2) can be inferred. It should be stressed that, for Eq.(17), possibility of a complex solution for $`\lambda _t^{}`$ should not be excluded.

### 2.2 Direct CP violation in $`BX_s\gamma `$

Observation of CP violation in $`BX_s\gamma `$ is attractive, because it could lead to an evidence related with the new physics. Theoretical predictions for $`BX_s\gamma `$ can be written as

$$A_{CP}(BX_s\gamma )=\frac{\mathrm{\Gamma }(\overline{B}X_s\gamma )\mathrm{\Gamma }(BX_{\overline{s}}\gamma )}{\mathrm{\Gamma }(\overline{B}X_s\gamma )+\mathrm{\Gamma }(BX_{\overline{s}}\gamma )}.$$

Numerically, prediction of the SM is

$$A_{CP}(BX_s\gamma )0.6\%,$$
(24)
when the best-fit values for the CKM parameters  are used. From the experimental side, we have the CLEO measurement of the CP asymmetry in the $`bs\gamma `$ decays ,

$$A_{CP}(BX_s\gamma )=(0.079\pm 0.108\pm 0.022)\times (1.0\pm 0.030),$$
(25)
We used the CP asymmetry formulae to look for 4 generation effects ,

$`A_{CP}(BX_s\gamma ){\displaystyle \frac{10^2}{|C_7|^2}}(1.17\times \mathrm{Im}\left[C_2C_7^{}\right]9.51\times \mathrm{Im}\left[C_8C_7^{}\right]`$ (26)
$`+0.12\times \mathrm{Im}\left[C_2C_8^{}\right]9.40\times \mathrm{Im}\left[ϵ_sC_2(C_7^{}0.013C_8^{})\right]);`$
$`ϵ_s={\displaystyle \frac{V_{us}^{}V_{ub}}{V_{ts}^{}V_{tb}}}\lambda ^2(\rho i\eta ).`$
As it is stated in the same reference, the large coefficient of the second term in (26) is very attractive.

We observed that, enhanced chromomagnetic dipole contribution, $`C_8`$, induces a large direct CP violation in the decay $`BX_s\gamma `$. This is due to complex phases of $`\lambda _t^{}`$, which in result affects $`C_7,C_8`$. Such an enhancement of the chromomagnetic contribution may lead to a natural explanation of the phenomenology of semileptonic $`B`$ decays and charm production in $`B`$ decays .

Notice that in $`A_{CP}`$ figures, when the real values of $`\lambda _t^{}`$ is around $`6\times 10^2`$, even for very small imaginary parts, peak values of $`A_{CP}`$ can be observed. Evolution of $`A_{CP}(BX_s\gamma )`$ is presented in figures $`\{7,8,9,10\}`$. CP asymmetry is not sensitive to very heavy $`m_t^{}`$ quark masses.

## 3 Conclusion

To summarize, the $`BX_s\gamma `$ decay has a clean experimental and theoretical base, very sensitive to the various extensions of the Standard Model, can be used to constrain the fourth generation model. In the present work, this decay is studied in the SM with the four generation model. The solutions of the fourth generation CKM factor $`\lambda _t^{}`$ have been obtained. It is observed that different choices of the factor $`\lambda _t^{}`$, could be very informative, especially due to new CP violation effects, in searching new physics.

CP asymmetry in the $`BX_s\gamma `$ decay can be enhanced up to $`5\%`$, which is ten times larger compared to the SM prediction. Hence it could be mentioned among the probes of new physics, especially in the case of fourth generation.