Numerical convergence of physical variables in hydrodynamical simulations of cooling clusters

被引:20
作者
Valdarnini, R [1 ]
机构
[1] SISSA, I-34014 Trieste, Italy
关键词
cooling flows; cosmology : theory; galaxies : clusters : general; intergalactic medium; methods : numerical; X-rays : galaxies : clusters;
D O I
10.1086/338548
中图分类号
P1 [天文学];
学科分类号
0704 ;
摘要
Results from smoothed particle hydrodynamical (SPH) simulations of galaxy clusters are used to investigate the dependence of the final cluster X-ray properties on the numerical resolution and the assumed models for the physical gas processes. Two different spatially flat cosmological models have been considered : a low-density cold dark matter universe with a vacuum energy density Omega(Lambda) = 0.7 (LambdaCDM) and a cold+hot dark matter model (CHDM). For each of these models two different clusters have been extracted from a cosmological N-body simulation. A series of hydrodynamical simulations has then been performed for each of them using a TREESPH code. These simulations first include radiative cooling and then also conversion of cold gas particles into stars; because of supernova explosions these particles can release energy in the form of thermal energy to the surrounding intracluster gas. For a specific treatment for the thermal state of the gas, simulation runs have been performed with different numerical resolutions. This is in order to disentangle in the final results for the cluster profiles, the effects of the resolution from those due to the assumed model for the gas thermal evolution. The numerical resolution of the simulation is controlled by the number of gas particles N-g and the chosen value for the gas gravitational softening parameter epsilon(g). The latter is proportional to the minimum SPH smoothing length and therefore sets a maximum spatial resolution for the simulations. For the cooling runs, final X-ray luminosities have been found to be diverging according to L-X proportional to 1/epsilon(g)(similar to5). The gas density profiles are also diverging at the cluster center. This is in agreement with previous findings. When cold gas particles are allowed to convert into stars, the divergences are removed. The final gas profiles show a well-defined core radius, and the temperature profiles are nearly flat. For the most massive test cluster in the LambdaCDM model, these simulations show a prominent cooling flow in the cluster core. This cluster was analyzed in detail, running simulations with different star formation methods and increasing numerical resolution. A comparison between different runs shows that the results of simulations, based on star formation methods in which gas conversion into stars is controlled by an efficiency parameter c(*), are sensitive to the numerical resolution of the simulation. In this respect star formation methods based instead on a local density threshold, as in Navarro & White, are shown to give more stable results. Final X-ray luminosities are found to be numerically stable, with uncertainties of a factor similar to2. These simulations are also in good agreement with observational data when the final results are compared with the observed star formation rate and the luminosity-temperature relation from cooling flow clusters. Therefore I find that hydrodynamical simulations of cooling clusters can be used to give reliabe predictions of the cluster X-ray properties. For a given numerical resolution, the conversion of cool gas particles into stars as in Navarro & White should be preferred.
引用
收藏
页码:741 / 761
页数:21
相关论文
共 103 条
[1]   The relationship between cooling flows and metallicity measurements for X-ray-luminous clusters [J].
Allen, SW ;
Fabian, AC .
MONTHLY NOTICES OF THE ROYAL ASTRONOMICAL SOCIETY, 1998, 297 (03) :L63-L68
[2]   The impact of cooling flows on the TX-LBol relation for the most luminous clusters [J].
Allen, SW ;
Fabian, AC .
MONTHLY NOTICES OF THE ROYAL ASTRONOMICAL SOCIETY, 1998, 297 (03) :L57-L62
[3]   ABUNDANCES OF THE ELEMENTS - METEORITIC AND SOLAR [J].
ANDERS, E ;
GREVESSE, N .
GEOCHIMICA ET COSMOCHIMICA ACTA, 1989, 53 (01) :197-214
[4]   Hierarchical numerical cosmology with hydrodynamics: Resolving X-ray clusters [J].
Anninos, P ;
Norman, ML .
ASTROPHYSICAL JOURNAL, 1996, 459 (01) :12-26
[5]   The most massive distant clusters:: Determining Ω and σ8 [J].
Bahcall, NA ;
Fan, XH .
ASTROPHYSICAL JOURNAL, 1998, 504 (01) :1-6
[6]  
BAHCALL NA, 1997, APJ, V485, P53
[7]   The dependence of cluster galaxy star formation rates on the global environment [J].
Balogh, ML ;
Schade, D ;
Morris, SL ;
Yee, HKC ;
Carlberg, RG ;
Ellingson, E .
ASTROPHYSICAL JOURNAL, 1998, 504 (02) :L75-L78
[8]   Statistical properties of X-ray clusters: Analytic and numerical comparisons [J].
Bryan, GL ;
Norman, ML .
ASTROPHYSICAL JOURNAL, 1998, 495 (01) :80-99
[9]   X-RAY-CLUSTERS FROM A HIGH-RESOLUTION HYDRODYNAMIC PPM SIMULATION OF THE COLD DARK-MATTER UNIVERSE [J].
BRYAN, GL ;
CEN, RY ;
NORMAN, ML ;
OSTRIKER, JP ;
STONE, JM .
ASTROPHYSICAL JOURNAL, 1994, 428 (02) :405-418
[10]   Galaxy formation and evolution - II. Energy balance, star formation and feedback [J].
Buonomo, F ;
Carraro, G ;
Chiosi, C ;
Lia, C .
MONTHLY NOTICES OF THE ROYAL ASTRONOMICAL SOCIETY, 2000, 312 (02) :371-379