Halo Concentrations and the Fundamental Plane of Galaxy Clusters

TL;DR

This paper reveals a fundamental plane relating halo structure and ICM temperature in galaxy clusters, showing that cluster formation history influences their properties and provides a new way to calibrate cluster masses.

Contribution

It demonstrates that galaxy clusters form a tight fundamental plane linking halo parameters and ICM temperature, explained by structure formation models, challenging the assumption of virial equilibrium.

Findings

Clusters evolve along the fundamental plane.

The plane explains the mass-temperature relation without virial equilibrium.

Cluster scatter on the plane relates to formation age diversity.

Abstract

According to the standard cold dark matter (CDM) cosmology, the structure of dark halos including those of galaxy clusters reflects their mass accretion history. Older clusters tend to be more concentrated than younger clusters. Their structure, represented by the characteristic radius $r_s$ and mass $M_s$ of the Navarro--Frenk--White (NFW) density profile, is related to their formation time. In~this study, we showed that $r_s$, $M_s$, and the X-ray temperature of the intracluster medium (ICM), $T_X$, form a thin plane in the space of $(\log r_s, \log M_s, \log T_X)$. This tight correlation indicates that the ICM temperature is also determined by the formation time of individual clusters. Numerical simulations showed that clusters move along the fundamental plane as they evolve. The plane and the cluster evolution within the plane could be explained by a similarity solution of structure formation of the universe. The angle of the plane shows that clusters have not achieved "virial equilibrium" in the sense that mass/size growth and pressure at the boundaries cannot be ignored. The distribution of clusters on the plane was related to the intrinsic scatter in the halo concentration--mass relation, which originated from the variety of cluster ages. The well-known mass--temperature relation of clusters ($M_\Delta\propto T_X^{3/2}$) can be explained by the fundamental plane and the mass dependence of the halo concentration without the assumption of virial equilibrium. The fundamental plane could also be used for calibration of cluster masses.

Equations

$${\rho }_{\mathrm{DM}}(r)=\frac{{\delta }_c{\rho }_c}{(r/r_s)(1+r/r_s{)}^2},$$

$$M(r)=4\pi {\delta }_c{\rho }_c{r}_s^3\left[\ln \left(1+\frac{r}{r_s}\right)-\frac{r}{r+r_s}\right].$$

$$M_{\Delta }=\frac{4\pi }{3}{\rho }_{\Delta }{r}_{\Delta }^3.$$

$$c_{\Delta }=r_{\Delta }/r_s$$

$$P(\lambda )D(\lambda {)}^{-\gamma }M(\lambda {)}^{10/3-3\gamma }=\mathrm{const},$$

$$\rho (r,t)={\rho }_{\mathrm{H}}D(\lambda ),p(r,t)={\rho }_{\mathrm{H}}(r_{\mathrm{ta}}/t{)}^{2}P(\lambda ),m(r,t)=(4\pi /3){\rho }_{\mathrm{H}}{r}_{\mathrm{ta}}^{3}M(\lambda ),$$

$$r_{\mathrm{ita}}\propto m_{\mathrm{ita}}^{(n+5)/6},t_{\mathrm{ita}}\propto m_{\mathrm{ita}}^{(n+3)/4},$$

$$r_s^2{M}_s^{-7/3}{T}_X\propto m_{\mathrm{ita}}^{(n-3)/6}.$$

$$r_s^2{M}_s^{-(n+11)/6}{T}_X=\mathrm{const},$$

$$M_{\Delta }\propto T_X^{3/2}{\Delta }^{-1/2}E(z{)}^{-1},$$

$$T_X=T_{X0}{\left(\frac{r_s}{r_{s0}}\right)}^{-2}{\left(\frac{M_s}{M_{s0}}\right)}^{(n+11)/6},$$

$$c_{200}(M_{200},z)=6.71{\left(\frac{M_{200}}{2\times 10^{12}{h}^{-1}{M}_{\odot }}\right)}^{-0.091}(1+z{)}^{-0.44}$$

$$r_{\Delta }={\left(\frac{3M_{\Delta }}{4\pi \Delta {\rho }_c(z)}\right)}^{1/3}.$$

$$M_s=M_{\Delta }\frac{\ln 2-1/2}{\ln (1+c_{\Delta })-c_{\Delta }/(1+c_{\Delta })}.$$

Full text

\pubvolume

6 \issuenum1 \articlenumber1

\historyReceived: 25 October 2018; Accepted:21 December 2018; Published: date \TitleHalo Concentrations and the Fundamental Plane of Galaxy Clusters \AuthorYutaka Fujita *1,**\orcidA, Megan Donahue 2, Stefano Ettori 3,4, Keiichi Umetsu 5, Elena Rasia 6, Massimo Meneghetti 3,4, Elinor Medezinski 7, Nobuhiro Okabe 8 and Marc Postman 9

\AuthorNamesYutaka Fujita

\corresCorrespondence: fujita@astro-osaka.jp; Tel.: +81-6-6850-5484

1 Introduction

Clusters of galaxies are the most massive objects in the Universe. Since the fraction of baryons in clusters is not much different from the cosmic mean value, dark matter accounts for most of the mass of clusters ($\sim$84%) 2011ApJS..192…18K ; 2018arXiv180706209P . Thus, the structure of the clusters is mainly determined by the halos of dark matter, or the dark halos. Cold dark matter (CDM) cosmology predicts that more massive halos form later. Thus, clusters form after galaxies do. However, the definition of the formation is not obvious, because halos are continuously growing through mergers and accretion from their environments. A current trend may be associating the formation time with the internal structure of dark halos.

The density distribution of dark halos is well-represented by the Navarro–Frenk–White (NFW) density profile nav96a ; nav97a :

[TABLE]

where $r$ is the clustercentric radius, $r_{s}$ is the characteristic radius, and $\rho_{c}$ is the critical density of the universe. The normalization of the profile is given by $\delta_{c}$. The characteristic mass $M_{s}$ is defined as the mass inside $r_{s}$ and the characteristic density is written as $\rho_{s}\equiv 3\>M_{s}/(4\pi r_{s}^{3})$. The mass profile of the NFW profile is written as

[TABLE]

Another commonly defined characteristic radius of clusters is that based on the critical density $\rho_{c}$; it is represented by $r_{\Delta}$, which is the radius of a sphere of mean interior density $\rho_{\Delta}\equiv\Delta\rho_{c}$, where $\Delta$ is the constant. The mass within $r_{\Delta}$ is written as

[TABLE]

The radius when $\Delta=200$ or $r_{200}$ is often called the “virial radius”. Since it is generally difficult to observationally study cluster properties out to $r\sim r_{200}$, $\Delta=500$ is also used as a representative value. The ratio

[TABLE]

is called the halo concentration parameter and $c_{\Delta}>1$ for $\Delta=200$ and 500 for clusters.

Navarro et al. nav97a pointed out that the characteristic parameters of the NFW profile (e.g., $\rho_{s}$ and $c_{\Delta}$) reflect the density of the background universe when the halo was formed. This means that, since older halos formed when the density of the universe was higher, they tend to have larger characteristic densities $\rho_{s}$ and become more concentrated with larger $c_{\Delta}$. This issue has been addressed in many studies, especially by $N$-body simulations nav97a ; bul01a ; 2001ApJ…554..114E ; wec02a ; zha03a ; 2006ApJ…646..815S ; 2007MNRAS.381.1450N ; 2008MNRAS.387..536G ; zha09a ; 2012MNRAS.423.3018P ; lud13a ; cor15c ; cor15b ; mor15b ; 2018arXiv181008212R . These studies have indicated that the halo structure is determined by their mass-growth history. The inner region ($r\ ^{\displaystyle<}_{\displaystyle\sim}\ r_{s}$)111More precisely, the boundary radius between the inner and outer regions is a few times $r_{s}$ (e.g. zha03a ). of current halos develops in the early “fast-rate growth” phase when the halos grow rapidly through matter accumulation. Their outer region ($r\ ^{\displaystyle>}_{\displaystyle\sim}\ r_{s}$) is formed in the subsequent “slow-rate growth” phase in which halos grow slowly through moderate matter accumulation. During this phase, the inner region is almost preserved. Thus, halos form “inside-out”. The formation time of a halo can be defined as the transition time from the fast-rate growth phase to the slow-rate growth phase. This shift of the growth phase is largely associated with the decrease in the average density of the universe in the $\Lambda$CDM cosmology. There are a few specific definitions of the formation time that well represent the transition time. One is the time at which the mass of the main progenitor was equal to the characteristic mass $M_{s}$ of the halo at its observed redshift $z_{\rm obs}$ lud13a ; cor15c . The formation redshift ($z_{f}$) corresponding to the formation time should be larger than $z_{\rm obs}$, or $z_{f}\geq z_{\rm obs}$. For a given $z_{\rm obs}$, clusters with a larger $z_{f}$ have a larger $\rho_{s}$ and $c_{\Delta}$.

Moreover, numerical simulations have shown that clusters are dynamically evolving systems and such evidence is often found in their outskirts. In fact, the ambient material is continuously falling toward clusters, which creates “surfaces” around clusters. For example, the outskirt profiles of dark matter halos can become extremely steep over a narrow range of radius (“splashback radius”). This features in the density profiles are caused by splashback of collisionless dark matter on its first apocentric passage after accretion 2014JCAP…11..019A ; 2014ApJ…789….1D . Accretion of collisional gas toward clusters also creates discontinuities in the form of shock fronts in their outskirts 2000ApJ…542..608M ; 2003ApJ…593..599R . These discontinuities mean that clusters are neither isolated nor in an equilibrium state.

In this review, we first introduce the fundamental plane we discovered, and its implications for structure formation of the universe. In particular, we show that clusters in general have not achieved virial equilibrium in contrast with conventional views. Then, using the fundamental plane, we discuss that a scaling relation (mass-temperature relation) can be explained without assuming virial equilibrium. We also show that the fundamental plane can be used for mass calibration. We assume a spatially-flat $\Lambda$CDM cosmology with $\Omega_{\mathrm{m}}=0.27$, $\Omega_{\Lambda}=0.73$, and the Hubble constant of $H_{0}=70$ km s*-1* Mpc*-1*.

2 Fundamental Plane

The hot intracluster medium (ICM) is distributed in the potential well of dark halos. Since the X-ray emission from the ICM is proportional to the square of the density, it mainly comes from the central region of the cluster where the density is high. Thus, the observed X-ray temperature $T_{X}$ represents that of the central region and should reflect the gravitational potential there. Since the potential is determined by $r_{s}$ and $M_{s}$, we can expect some relation among $T_{X}$, $r_{s}$, and $M_{s}$.

Based on this motivation, Fujita et al. fuj18a analyzed 20 massive clusters in the Cluster Lensing And Supernova survey with Hubble (CLASH) observational sample pos12a ; men14a . For these clusters, $r_{s}$ and $M_{s}$ had been obtained from the joint analysis ume16a of strong lensing observations with 16-band Hubble Space Telescope observations 2015ApJ…801…44Z and weak-lensing observations mainly with Suprime-Cam on the Subaru Telescope ume14a . The X-ray temperature had been obtained with Chandra pos12a ; don14a . Temperatures are estimated for a cylindrical volume defined by the projected radii $r=50$–500 kpc to avoid the influence of cool cores. Figure 1a shows the data distribution in the $(\log r_{s},\log M_{s},\log T_{X})$ space. As can be seen, the data are distributed on a plane. The plane is described by $r_{s}^{a}M_{s}^{b}T_{X}^{c}=\mathrm{const.}$, with $a=0.76^{+0.03}_{-0.05}$, $b=-0.56^{+0.02}_{-0.02}$, and $c=0.32^{+0.10}_{-0.09}$. Figure 1b shows the cross-section of the plane; the dispersion of the data around the plane is very small and is only $0.045^{+0.008}_{-0.007}$ dex (all uncertainties are quoted at the $1\>\sigma$ confidence level unless otherwise mentioned). In Figure 1c, error bars for individual clusters are shown. In the vertical direction ($T_{X}$), we show the temperature errors of individual clusters. In the horizontal direction, the errors of $r_{s}$ and $M_{s}$ are strongly correlated, and we display them as a single bar. This means that we draw a bar connecting ($r_{s}^{u}$, $M_{s}^{u}$) and ($r_{s}^{l}$, $M_{s}^{l}$) for each cluster, where the superscripts $u$ and $l$ are the upper and the lower limits, respectively. We note that, when we calculate the plane parameters, we properly account for the correlation for each cluster using the joint posterior probability distribution of the NFW parameters (mass and concentration) fuj18a . Thus, the actual error is not represented by a single bar in a precise sense. Figure 2 shows the direction of the plane normal $P_{3}$ in the $(\log r_{s},\log M_{s},\log T_{X})$ space fuj18a . The contours are drawn considering the errors and show that the direction is inconsistent with the prediction of a simplified dimensional analysis or $T_{X}\propto M_{s}/r_{s}$. We note that it is meaningless to discuss cluster distribution in the space of $(\log r_{200},\log M_{200},\log T_{X})$, where $r_{200}$, $M_{200}$, and $T_{X}$ are their current values. Although the clusters form a plane in that space, it is just the obvious relation of $M_{200}=4\pi\rho_{200}r_{200}^{3}/3\propto r_{200}^{3}$ regardless of $T_{X}$ (Section 5.3 in fuj18a ).

The “fundamental plane” 222Other fundamental planes of clusters with different combinations of three parameters have also been studied 1993MNRAS.263L..21S ; 1998A&A…331..493A ; 1999ApJ…519L..51F ; 2002ApJ…581….5V ; 2004ApJ…600..640L ; 2006ApJ…640..673O ; 2009MNRAS.400.1317A .has been reproduced by numerical simulations. Figure 3 shows the results of MUSIC N-body/hydrodynamical simulations (see details in men14a ; fuj18a ). These simulations do not include radiative cooling or non-gravitational feedback by active galactic nuclei (AGNs) and supernovae (SNe). The mass resolutions for the dark-matter particles that for the gas particles are $m_{\rm DM}=9.01\times 10^{8}\>h^{-1}\>M_{\odot}$ and $m_{\rm SPH}=1.9\times 10^{8}\>h^{-1}\>M_{\odot}$, respectively, where the Hubble constant is written as $H_{0}=100\>h$ km s*-1* Mpc*-1* and $h=0.7$. The gravitational softening is set to be 6 $h^{-1}$ kpc for the both gas and dark-matter particles in high-resolution regions. We chose all 402 clusters at $z=0.25$ with $M_{200}>2\times 10^{14}\>h^{-1}\>M_{\odot}$ regardless of dynamical state. In this analysis, we included the core because these simulations are non-radiative and thus do not present cool-core features. The absolute position of the plane is very close to that of the CLASH observational data (Figure 1b). Figure 2 shows that the plane angle for the MUSIC sample is consistent with the CLASH data at the 90% confidence level.

Figure 4 presents the results of other numerical simulations including phenomena such as heating by AGNs and SNe in addition to radiative cooling (see details in ras15a ; fuj18a ). Sample FB0 (blue dots) consists of the clusters at $z=0$, while sample FB1 (green dots) refers to the runs at $z=1$. These runs simulate 29 Lagrangian regions around massive clusters with $M_{\rm 200}\sim 1$–$30\times 10^{14}\>h^{-1}\>M_{\odot}$ at $z=0$. The mass resolution for the dark-matter particles and the initial gas particles are $m_{\rm DM}=8.3~{}\times~{}10^{8}\>h^{-1}\>M_{\odot}$ and $m_{\rm SPH}=1.5\times 10^{8}\>h^{-1}\>M_{\odot}$, respectively. In high-resolution regions the gravitational softening is set to be $3.75\>h^{-1}$ kpc (Biffi et al., 2017). For the gas particles, this is always in comoving units, while for the DM particles it changes to physical units from $z=2$ to $z=0$. For these samples, the temperature is estimated for $r=50$–500 kpc, and thus the influence of cool cores is removed. Both groups of dots are located on almost the same fundamental plane, and the plane angles for the two samples are almost the same (Figure 2). This means that clusters evolve along the unique plane in the direction of $P_{1}$ in Figure 4a. The plane angles for FB0 and FB1 are not much different from those for the CLASH data and the MUSIC adiabatic simulations (Figure 2). In Figure 2, NF0 is the result of a simulation that is the same as FB0 but not including radiative cooling and feedback. Since their angles are almost the same, this indicates that radiative cooling and feedback do not much affect the fundamental plane. This is because we are discussing cluster properties on a scale of $r\sim r_{s}\ ^{\displaystyle>}{}_{\displaystyle\sim}\ 300$ kpc, and the influences of cool cores, where radiative cooling and feedback are especially important, are ignorable.

3 Origin of the Fundamental Plane and Cluster Growth

Fujita et al. fuj18a explained the origin of the fundamental plane using an analytic similarity solution developed by Bertschinger ber85a (see also shi16b ). This solution treats spherical collapse of an overdense region in the Einstein–de Sitter universe ($\Omega_{0}=1$) and subsequent matter accretion onto the collapsed object. In the solution, matter profiles are represented by non-dimensional radius, $\lambda$, density $D(\lambda)$, pressure $P(\lambda)$, and mass $M(\lambda)$. The solution has a constant called the “entropy constant”.

[TABLE]

where $\gamma=5/3$ is the adiabatic index. The non-dimensional parameters are related to dimensional density $\rho$, pressure $p$, and mass $m$:

[TABLE]

where $r_{\rm ta}(t)$ is the maximum radius that a mass shell reaches (the turnaround radius), $\rho_{\rm H}\propto t^{-2}$ is the density of the background universe, and $t$ is the cosmological time. The non-dimensional radius is given by $\lambda=r/r_{\rm ta}$. We note that the similarity solution describes the matter profile in the region where the matter is later accreted (say, $r\ ^{\displaystyle>}_{\displaystyle\sim}\ r_{s}$). The solution was originally developed for objects totally composed of baryons, thus $\rho$, $p$, and $m$ are for the gas. However, the non-dimensional profiles ($D$, $p$, and $M$) are not much changed, even if objects are mostly composed of dark matter ber85a . Thus, the profiles $\rho$, $p$, and $m$ can be regarded as the values for the gas as long as we do not discuss the normalizations. Although the solution is constructed for the Einstein–de Sitter universe, it well-reproduces the structure of objects except for the outermost region even for a $\Lambda$CDM universe, because the inner region was established when the background density of the universe was large ber85a . From Equations (5) and (6), we obtain $p\rho^{-5/3}m^{-5/3}\propto A_{\rm ita}^{-3}$, where $A_{\rm ita}=r_{\rm ita}/t_{\rm ita}^{8/9}$ and is time-independent. Here, $r_{\rm ita}$ and $t_{\rm ita}$ are the turnaround radius (the maximum radius before the collapse) and the turnaround time (the time when the radius reaches turnaround radius) of the overdense region, respectively. The evolution of the overdense region is described by the conventional spherical collapse model. Thus, it should follow the spectrum of initial density perturbations of the universe and the mass of the overdense region $m_{\rm ita}$ has scaling relations of

[TABLE]

where $n$, the spectral index, is of the initial density perturbations kai86a and $n\sim-2$ is expected around the mass scales of clusters eis98a ; die15a . Here, we emphasize that Equation (7) is applied to the overdense region, and not to the entire cluster, because we separately treat the initial collapse of the overdense region and the subsequent matter accretion. Thus, we obtain $p\rho^{-5/3}m^{-5/3}\propto m_{\rm ita}^{(n-3)/6}$. Assuming that $p\propto\rho T_{X}$, $\rho\propto M_{s}/r_{s}^{3}$ and $m\sim M_{s}$ at $r\sim r_{s}$, it is written as

[TABLE]

The radius $r_{\rm ita}$333Note that, although the radius $r_{\rm ita}$ is the turnaround radius of the overdense region, it is proportional to the radius of the region after the collapse because the solution is similar. and the mass $m_{\rm ita}$ of the overdense region can be connected to the characteristic radius $r_{s}$ and mass $M_{s}$ of the NFW profile. This is because the evolution of both overdense region in the similarity solution and the inner region of the NFW profile ($r\ ^{\displaystyle<}_{\displaystyle\sim}\ r_{s}$) is related to the background universe, and they evolve in a similar way. In fact, the evolution of the former is described by the conventional spherical collapse of an overdense region ber85a , and thus the typical density is proportional to that of the background universe at the collapse. The same applies to the latter because the characteristic density $\rho_{s}$ is always $\sim$900 times as large as that of the background universe at the formation redshift $z_{f}$ cor15c . Thus, we can assume that $r_{s}\propto r_{\rm ita}$ and $M_{s}\propto m_{\rm ita}$, and that the collapse time of the overdense region ($\sim$2 $t_{\rm ita}$; see, e.g., pee80 , section 19) corresponds to the formation redshift $z_{f}$. In summary, the similarity solution has scales such as $r_{\rm ita}$ and $m_{\rm ita}$. Since $r_{\rm ita}$ and $m_{\rm ita}$ are proportional to $r_{s}$ and $M_{s}$, respectively, the solution has scales of $r_{s}$ and $M_{s}$. These scales represent the border between the initially-collapsed overdense region and the later-accreted region, although the actual transition is gradual. This also means that the initial collapse and the later accretion correspond to the fast-rate growth phase and the slow-rate growth phase, respectively. Finally, from Equation (8), we obtain

[TABLE]

or $T_{X}\propto M_{s}^{(n+11)/6}/r_{s}^{2}$. Equation (9) forms a plane in the $(\log r_{s},\log M_{s},\log T_{X})$ space. The direction of the normal when $n=-2$ is shown in Figure 2 as “SSol”, and it is consistent with the CLASH observations and the results of numerical simulations. Note that this relation (Equation (9)) is independent of redshift $z$ at least $z\ ^{\displaystyle<}_{\displaystyle\sim}\ 1$ (Figure 4), because $r_{s}$ and $M_{s}$ are physical values that have already reflected the high density of the background universe in the past.

The similarity solution indicates that clusters are not in virial equilibrium, because clusters are growing through matter accretion from their outer environments ber85a ; shi16b . That is one reason clusters follow Equation (9) instead of $T_{X}\propto M_{s}/r_{s}$, which could be realized if clusters are in virial equilibrium at their formation. The condition of the virial “equilibrium” is represented by $2\>K+W=0$, where $K$ is the kinetic and/or thermal energy and $W$ is the gravitational energy. However, according to the virial “theorem”, which is a higher-level concept of the virial “equilibrium”, additional terms are required when clusters are growing ber85a ; shi16b . One is the term representing the increase of mass and size of clusters and another is the boundary term originating from the flux of inertia through the boundary and the pressure at the boundary. The boundary corresponds to the splashback radius for dark matter and the shock front for gas (see Section 1). Note that the similarity solution shows that clusters are almost in hydrostatic equilibrium in the sense that gas motion is negligible within clusters even if they are not in virial equilibrium ber85a . The relation between matter accretion and the cluster structure has also been numerically studied (e.g., Ref. 2018arXiv181008212R ).

Figure 5 shows the projection of the fundamental plane shown in Figure 3 on the $\log r_{s}$–$\log M_{s}$ plane. The solid arrow is parallel to the line of $r_{s}\propto M_{s}^{1/2}$ along which the distribution of the MUSIC clusters (red points) is elongated. This direction is also close to that of cluster evolution ($P_{1}$) in Figure 4 projected on the $\log r_{s}$–$\log M_{s}$ plane fuj18a . Since we assumed that $r_{s}\propto r_{\rm ita}$ and $M_{s}\propto m_{\rm ita}$, the line corresponds to the first relation of Equation (7) when $n=-2$. Considering the derivation of the relations in Equation (7) (see kai86a ), this indicates that the evolution of clusters on the fundamental plane reflects the spectrum of the initial density perturbations of the universe and follows $M_{s}\propto m_{\rm ita}\propto(1+z_{f})^{-6/(n+3)}$ kai86a . Figure 5a also shows that the characteristic density $\rho_{s}$ decreases as a cluster moves in the direction of the solid arrow. While the formation redshift $z_{f}$ is formally related to the collapse time of the overdense region, in reality, it is often related to the time of major cluster mergers. That is, the formation redshift $z_{f}$ is reset when the cluster experiences a major merger, and $z_{f}$ estimated from $\rho_{s}$ for a given cluster at $z_{\rm obs}$ often corresponds to the time when the cluster underwent its last major merger. In fact, numerical simulations have shown that an individual cluster intermittently moves in the direction of the solid arrow in Figure 5 every time it undergoes mergers fuj18a . While the cluster temporarily deviates the general motion in the middle of a major merger, the deviation is confined in the fundamental plane and thus mergers do not much affect the thinness of the plane fuj18a . In other words, the effect of major cluster mergers introduces some random history that could be different for clusters of the same mass, but since the mergers move cluster properties within the limits of the plane, the scatter of the plane does not increase very much.

We would like to point out that in Figure 5 simulated clusters are not tightly distributed along the line of $r_{s}\propto M_{s}^{1/2}$ (solid arrow), and there is a scatter about the line. This reflects the fact that the density perturbations of the universe are described by a Gaussian random field (see, e.g., bar01a ). Thus, while the variance of the perturbation field $\sigma(M)$ is a decreasing function of mass scale $M$, the amplitudes of the perturbations that collapse into objects with a given mass $M$ are not always $\sigma(M)$. Owing to this, for example, $\rho_{s}$ and $M_{s}$ are not perfectly in one-to-one correspondence, and $\rho_{s}$ has some range for a given $M_{s}$, which produces the band-like distribution of clusters in Figure 5 and on the fundamental plane (Figures 1a, 3a and 4a). In other words, clusters form a two-parameter family. Thus, a correlation between two physical quantities is generally represented by a band rather than a line unless some special combination of quantities is chosen. In that sense, it is natural that the relation between $c_{\Delta}$ and $M_{\Delta}$ has a large dispersion bul01a ; lud13a ; 2008MNRAS.390L..64D ; men14a ; cor15c , which will be discussed in Section 4. On the fundamental plane, different clusters move along nearly parallel but different tracks each of which approximately follows the relation of $r_{s}\propto M_{s}^{1/2}$ fuj18a . While the temperature of each cluster $T_{X}$ is affected by its formation time, it also depends on the track the cluster chooses.

4 Mass–Temperature Relation and the Concentration Parameter

The fundamental plane can be used to relate the cluster structure to the temperature. As an application, we discuss the mass–temperature relation in this section. It is well-known that the mass of clusters and the X-ray temperature has a relation of $M_{\Delta}\propto T_{\Delta}^{3/2}$. This relation is obtained by both observations and numerical simulations bry98a ; ett02b ; 2016A&A…592A…4L ; tru18a . Conventionally, this relation has been explained based on the following three assumptions: (i) the typical density of a cluster is $\rho_{\Delta}=\Delta\rho_{c}$ (not $\rho_{s}$); (ii) clusters are well-relaxed or virialized, and they are almost isothermal within $r_{\Delta}$; and (iii) the X-ray temperature is determined on a scale of $r_{\Delta}$ (not $r_{s}$). Here, we consider cluster temperature outside cool cores.

The density $\rho_{\Delta}$ is represented by $\rho_{\Delta}\propto\Delta E(z)^{2}$, where $E(z)$ is the Hubble parameter at redshift $z$ normalized by the current value $H_{0}$. Equation (3) is associated with Assumption (i). From Assumptions (ii) and (iii), we obtain $T_{X}\propto M_{\Delta}/r_{\Delta}\propto\rho_{\Delta}r_{\Delta}^{2}\propto\Delta E(z)^{2}r_{\Delta}^{2}$. Eliminating $r_{\Delta}$ by using the relation $r_{\Delta}\propto M_{\Delta}/T_{X}$, the mass–temperature relation is obtained:

[TABLE]

which well reproduces the results of observations and simulations bry98a ; 2011ASL…..4..204B ; 2015SSRv..188…93P . However, the assumptions are clearly inconsistent with the inside-out scenario of cluster formation and the fundamental plane. For example, the inside-out scenario indicates that clusters are not well relaxed and keep the memory of their formation in their structure. The angle of the fundamental plane shows that clusters are not virialized, as discussed in Section 3. The NFW profile (Equation (1)) is not an isothermal profile ($\rho_{\rm DM}\propto r^{-2}$). These are inconsistent with Assumption (ii). Moreover, the tight correlation of the fundamental plane shows that $T_{X}$ is determined by $r_{s}$ and $M_{s}$, which contradicts Assumption (iii).

In fuj18b , Fujita et al. showed that the relation in Equation (10) can be derived using the fundamental plane and the halo concentration–mass ($c_{\Delta}$–$M_{\Delta}$) relation. The fundamental plane relation in Equation (9) is rewritten as

[TABLE]

where $(r_{s0},M_{s0},T_{X0})$ corresponds to a representative point on the fundamental plane, and we adopt $r_{s0}=414$ kpc, $M_{s0}=1.4\times 10^{14}\>M_{\odot}$, and $T_{X0}=3.7$ keV based on the results of the MUSIC simulations (men14a, ; fuj18a, ). Based on the inside-out scenario, there are analytical forms of the concentration parameter represented as a function of $M_{\Delta}$ and $z$. One example is

[TABLE]

for $\Delta=200$ that was obtained by Duffy et al. 2008MNRAS.390L..64D (see also 2013ApJ…766…32B ; dut14b ; men14a ; die15a ; 2018arXiv180907326D ). From Equation (3), we obtain

[TABLE]

Equations (4), (12) and (13) indicate that $r_{s}$ is a function of $M_{\Delta}$ for a given $z$. Moreover, Equation (2) suggests that $M_{s}$ is also a function of $M_{\Delta}$:

[TABLE]

Thus, using Equation (11), $T_{X}$ can be represented as a function of $M_{\Delta}$ for a given $z$. Figure 6a shows the results for $n=-2$ using a general formula of $c_{200}(M_{200},z)$ developed by Correa et al. cor15c instead of Equation (12). The slope is $\alpha=1.33$ for $z=0$ and $\alpha=1.28$ for $z=1$ ($M_{\Delta}\propto T_{X}^{\alpha}$). The slope is close to but slightly smaller than $\alpha=1.5$. However, the derivation of the fundamental plane in Section 3 might be too simplified, and there might be some minor uncertainties on $n$ fuj18b . In fact, if we take $n=-2.5$, the slope becomes $\alpha=1.53$ for $z=0$ and $\alpha=1.45$ for $z=1$. Note that, even if we assume $n=-2.5$, the direction of the fundamental plane (Equation (9)) is consistent with observations and simulations fuj18b 444Here, we see $n$ as a parameter of the direction of the fundamental plane, and we do not intend to claim that the spectral index of the initial density perturbations is exactly $-$2.5.. Thus, the relation of $M_{\Delta}\propto T_{X}^{3/2}$ can be reproduced without the virial assumption or $T_{X}\propto M_{\Delta}/r_{\Delta}$. Note that Figure 6 indicates that the red lines ($z=1$) are slightly below the black lines ($z=0$). This may cause some bias about the slope index $\alpha$ if clusters with various redshifts are plotted at the same time. For example, if higher-redshift clusters ($z\sim 1$) tend to have smaller masses and lower temperatures than lower-redshift clusters ($z\sim 0$), the slope is slightly steepened (larger $\alpha$). We note that Voit 2000ApJ…543..113V (see also 1998ApJ…500L.111V ) already addressed this issue. He considered accretion history of clusters and the effects of cluster surfaces as we do. While we focused on the inner structure of clusters, he studied the evolution of global properties of clusters. He concluded that the approximate agreement between the $M_{\Delta}$–$T_{X}$ relation derived via the traditional collapse model (Equation (10)) and those of simulations and observations is largely coincidental. Although our approach is different, our results support the conclusion.

The relation of $c_{\Delta}$–$M_{\Delta}$ or the function $c_{\Delta}=c_{\Delta}(M_{\Delta},z)$ can be converted into the relation between $r_{s}$ and $M_{s}$ using Equations (3), (4), and (14), and the result is shown by the solid black line in Figure 5b. The black dotted and dashed lines correspond to the dispersion of $c_{\Delta}$–$M_{\Delta}$ relation indicated by numerical simulations. The three black lines in Figure 5b are almost parallel to the lines of $\rho_{s}=$ const or the three black dashed lines in Figure 5a. This means that the dispersion of $c_{\Delta}$–$M_{\Delta}$ relation is almost the same as that of $\rho_{s}$ or the dispersion of cluster formation time $z_{f}$. Figure 5b also indicates that the minor axis of the cluster distribution (red points) corresponds to the dispersion of the $c_{\Delta}$–$M_{\Delta}$ relation. The dispersion of the $c_{\Delta}$–$M_{\Delta}$ relation is also associated with that of the $M_{\Delta}$–$T_{X}$ relation, which is indicated by the black dotted and dashed lines in Figure 6. In Figure 7, we present the evolution of simulated clusters along the $M_{\Delta}$–$T_{X}$ relation. As expected, the clusters move along the bands enclosed by the dotted and dashed lines. The clusters frequently move in the horizontal direction, which corresponds to temporal temperature increase during cluster mergers. However, even during the mergers, the clusters are located within the bands, which means that the $M_{\Delta}$–$T_{X}$ relation is not much affected by mergers.

5 Cluster Mass Calibration

The thinness and solidity of the fundamental plane inspires applications in cosmology. Here, we show that the plane can be used to calibrate cluster mass fuj18b . Precise estimation of cluster mass is important. For example, when cosmological parameters are derived from cluster number counts, scaling relations among observables are used and they are affected by the calibration of cluster mass pla14b .

Figure 8a shows the cross sections of the fundamental plane. The red open circles are the clusters of the CLASH sample pos12a , for which $r_{s}$ and $M_{s}$ are derived through gravitational lensing. The black dots are those of an X-ray sample ett10a , for which $r_{s}$ and $M_{s}$ are derived through X-ray observations assuming that the ICM is in hydrostatic equilibrium. We discuss the fundamental plane formed by the CLASH sample (CFP) and the one formed by the X-ray sample (XFP) separately. Fixing the direction of the plane normals at the one shown by SSol in Figure 2, the distance between the two fundamental planes is estimated to be $d_{\rm FP}=0.031^{+0.027}_{-0.039}$ dex in the space of $(\log r_{s},\log M_{s},\log T_{X})$. Thus, the position of the fundamental planes are consistent with each other. However, the XFP seems to be located slightly above the CFP in Figure 8a. The shift $d_{\rm FP}$ may be caused by a possible systematic difference of observed $r_{s}$ or $M_{s}$ between CFP and XFP because they are obtained through different methods (gravitational lensing and X-ray observations). The plane shift in the direction of $r_{s}$ or $M_{s}$ can be estimated from $d_{\rm FP}$. Then, assuming the NFW profile (Equation (1) or (2)), the shift in the direction of $M_{\Delta}$ can be derived fuj18b .

Figure 8b shows the systematic difference of $M_{\Delta}$, which is defined by $f_{M\Delta}\equiv M_{\Delta X}/M_{\Delta C}$, where $M_{\Delta X}$ is the mass of a cluster on the XFP, and $M_{\Delta C}$ is the mass of the same cluster on the CFP. While the ratio $f_{M\Delta}$ depends on the concentration parameter $c_{\Delta}$, the dependence is weak. Figure 8b shows that $f_{M\Delta}\sim 0.85_{-0.2}^{+0.2}$, which means that the cluster mass estimated through X-ray observations assuming hydrostatic equilibrium (hydrostatic mass) is $\sim$15% smaller than that estimated through gravitational lensing. Since the error is rather large, the current dataset may not be accurate enough for the calibration purpose. However, the error could be reduced by using larger and more accurate datasets in the future. Assuming that gravitational lensing mass is solid, the value of $f_{M\Delta}\sim 0.85$ is consistent with the results of numerical simulations showing that hydrostatic mass tends to be smaller than the true mass nag07a ; 2008A&A…491…71P ; 2010A&A…510A..76L ; 2012NJPh…14e5018R .

6 Sparsity

Finally, we would like to make comments on the halo “sparsity”, which has been proposed recently (balmes14, ; corasaniti18, ) as a valid alternative to the full description of the dark matter profile. It measures the ratio of halo mass at two different radii (e.g., $M_{500}/M_{1000}$) and, in the case that the halo follows a NFW profile, it is directly related to the halo concentration. The advantage in using the halo sparsity is that it has an ensemble average value at a given redshift with a scatter much smaller than that associated to the distribution in mass concentration and does not require any modeling of the mass density profile, which might significantly deviate from a NFW one in particular in systems still in the process of complete relaxation, but only the integrated mass measurements within two overdensities. The use of the halo sparsity has also been proposed as new cosmological probe for galaxy clusters (corasaniti18, ) because it carries cosmological information encoded in the halo mass profile and, at given redshift, the average sparsity can be predicted from prior knowledge of the halo mass function.

Both the fundamental plane and the halo sparsity reflect the halo concentration of clusters. While the fundamental plane gives us the direct information of cluster formation time, it is generally difficult to measure $r_{s}$ and $M_{s}$ observationally, compared with the sparsity. In future study, we will discuss the relation between the fundamental plane and the halo sparsity.

7 Conclusions

It has been known that the concentration of dark halos reflects their formation history. In particular, the halo structure represented by the characteristic radius $r_{s}$, and mass $M_{s}$ is related to the formation time of the halo. In this study, we showed that $r_{s}$, $M_{s}$, and the X-ray temperature $T_{X}$ of observed clusters form a plane (fundamental plane) in the space of $(\log r_{s},\log M_{s},\log T_{X})$ with a very small orthogonal scatter. The tight correlation shows that $T_{X}$ is also affected by the formation time of individual clusters. Numerical simulations supported the results and showed that clusters evolve along the plane. The plane and its angle in the space of $(\log r_{s},\log M_{s},\log T_{X})$ can be explained by a similarity solution, which indicates that clusters are still growing and have not reached a state of virial equilibrium. In other words, when cluster formation and the internal structure is considered, matter accretion after the initial collapse cannot be ignored. The motion of clusters on the plane was determined by the spectrum of the initial density perturbations of the universe. The spread of clusters on the fundamental plane is related to the scatter of the halo concentration–mass relation.

We also discussed applications of the fundamental plane. For example, we showed that the mass–temperature relation of clusters ($M_{\Delta}\propto T_{X}^{3/2}$) can be explained by the fundamental plane and the halo concentration–mass relation without assuming virial equilibrium. We also showed that the solidity and thinness of the fundamental plane can be used to calibrate cluster mass. Since the fundamental plane associates the structure of dark halos with the gas temperature, other applications may be possible. For example, the gas temperature $T_{X}$ of a dark halo could be estimated from $r_{s}$ and $M_{s}$ obtained through $N$-body simulations without calculating gas dynamics.

\authorcontributions

Y.F. coordinated the research, wrote the manuscript, obtained the fundamental plane, and made a theoretical interpretation. S.E. provided the X-ray data and contributed in the interpretation of the findings. K.U. analyzed the data of gravitational lensing and contributed in the interpretation of the findings. E.R. and M.M. analyzed the numerical simulations, kindly provided by the computational groups in Trieste and Madrid, and fit the data with the NFW profile. M.D., E.M., N.O., and M.P. contributed in the interpretation of the optical and X-ray results and provided extensive feedback on the study.

\funding

This work was supported by MEXT KAKENHI No. 18K03647 (YF). S.E. acknowledges financial contribution from the contracts NARO15 ASI-INAF I/037/12/0, ASI 2015-046-R.0 and ASI-INAF n.2017-14-H.0. K.U. acknowledges support from the Ministry of Science and Technology of Taiwan (grant MoST 106-2628- M-001-003-MY3) and from Academia Sinica (grant AS-IA-107-M01). E.R. acknowledge support from the ExaNeSt and EuroExa projects, funded by the European Union’s Horizon 2020 research and innovation programme under grant agreements No. 671553 and No. 754337, respectively.

Acknowledgements.

We thank referees for their useful comments, which were very helpful to clarify explanations. \conflictsofinterestThe authors declare no conflict of interest.