Subscribe to the gold package and get unlimited access to Shamra Academy
Register a new userStudy of galaxy cluster properties from high-resolution SPH simulations
55
0
0.0
(
0
)
Ask ChatGPT about the research
No Arabic abstract
We present some of the results of an ongoing collaboration to sudy the dynamical properties of galaxy clusters by means of high resolution adiabatic SPH cosmological simulations. Results from our numerical clusters have been tested against analytical models often used in X-ray observations: $beta$ model (isothermal and polytropic) and those based on universal dark matter profiles. We find a universal temperature profile, in agreement with AMR gasdynamical simulations of galaxy clusters. Temperature decreases by a factor 2-3 from the center to virial radius. Therefore, isothermal models (e.g. $beta$ model) give a very poor fit to simulated data. Moreover, gas entropy profiles deviate from a power law near the center, which is also in very good agreement with independent AMR simulations. Thus, if future X-ray observations confirm that gas in clusters has an extended isothermal core, then non-adiabatic physics would be required in order to explain it.
rate research
Read More
We present results from a set of high (512^3 effective resolution), and ultra-high (1024^3) SPH adiabatic cosmological simulations of cluster formation aimed at studying the internal structure of the intracluster medium (ICM). We derive a self-consistent analytical model of the structure of the intracluster medium (ICM). We discuss the radial structure and scaling relations expected from purely gravitational collapse, and show that the choice of a particular halo model can have important consequences on the interpretation of observational data. The validity of the approximations of hydrostatic equilibrium and a polytropic equation of state are checked against results of our simulations. The properties of the ICM are fully specified when a universal profile is assumed for either the dark or the baryonic component. We also show the first results from an unprecedented large-scale simulation of 500 Mpc/h and 2 times 512^3 gas and dark matter particles. This experiment will make possible a detailed study of the large-scale distribution of clusters as a function of their X-ray properties.
We perform N-Body/SPH simulations of disk galaxy formation inside equilibrium spherical and triaxial cuspy dark matter halos. We systematically study the disk properties and morphology as we increase the numbers of dark matter and gas particles from 10^4 to 10^6 and change the force resolution. The force resolution influences the morphological evolution of the disk quite dramatically. Unless the baryon fraction is significantly lower than the universal value, with high force resolution a gaseous bar always forms within a billion years after allowing cooling to begin. The bar interacts with the disk, transferring angular momentum and increasing its scale length. In none of the simulations does the final mass distribution of the baryons obey a single exponential profile. Indeed within a few hundred parsecs to a kiloparsec from the center the density rises much more steeply than in the rest of the disk, and this is true irrespective of the presence of the bar.
We investigate, by means of numerical simulations, the phenomenology of star formation triggered by low-velocity collisions between low-mass molecular clumps. The simulations are performed using an SPH code which satisfies the Jeans condition by invoking On-the-Fly Particle Splitting. Clumps are modelled as stable truncated (non-singular) isothermal, i.e. Bonnor-Ebert, spheres. Collisions are characterised by M_0 (clump mass), b (offset parameter, i.e. ratio of impact parameter to clump radius), and M (Mach Number, i.e. ratio of collision velocity to effective post-shock sound speed). The gas subscribes to a barotropic equation of state, which is intended to capture (i) the scaling of pre-collision internal velocity dispersion with clump mass, (ii) post-shock radiative cooling, and (iii) adiabatic heating in optically thick protostellar fragments. The efficiency of star formation is found to vary between 10% and 30% in the different collisions studied and it appears to increase with decreasing M_0, and/or decreasing b, and/or increasing M. For b<0.5 collisions produce shock compressed layers which fragment into filaments. Protostellar objects then condense out of the filaments and accrete from them. The resulting accretion rates are high, 1 to 5 x 10^{-5} M_sun yr^{-1}, for the first 1 to 3 x 10^4 yrs. The densities in the filaments, n >~ 5 x 10^5 cm^{-3}, are sufficient that they could be mapped in NH_3 or CS line radiation, in nearby star formation regions.
We present the McMaster Unbiased Galaxy Simulations (MUGS), the first 9 galaxies of an unbiased selection ranging in total mass from 5$times10^{11}$ M$_odot$ to 2$times10^{12}$ M$_odot$ simulated using n-body smoothed particle hydrodynamics (SPH) at high resolution. The simulations include a treatment of low temperature metal cooling, UV background radiation, star formation, and physically motivated stellar feedback. Mock images of the simulations show that the simulations lie within the observed range of relations such as that between color and magnitude and that between brightness and circular velocity (Tully-Fisher). The greatest discrepancy between the simulated galaxies and observed galaxies is the high concentration of material at the center of the galaxies as represented by the centrally peaked rotation curves and the high bulge-to-total ratios of the simulations determined both kinematically and photometrically. This central concentration represents the excess of low angular momentum material that long has plagued morphological studies of simulated galaxies and suggests that higher resolutions and a more accurate description of feedback will be required to simulate more realistic galaxies. Even with the excess central mass concentrations, the simulations suggest the important role merger history and halo spin play in the formation of disks.
Abridged: We simulate a massive galaxy cluster in a LCDM Universe using three different approaches to solving the equations of non-radiative hydrodynamics: `classic Smoothed Particle Hydrodynamics (SPH); a novel SPH with a higher order dissipation switch (SPHS); and adaptive mesh refinement (AMR). We find that SPHS and AMR are in excellent agreement, with both forming a well-defined entropy core that rapidly converges with increasing mass and force resolution. By contrast, SPH exhibits rather different behaviour. At low redshift, entropy decreases systematically with decreasing cluster-centric radius, converging on ever lower central values with increasing resolution. At higher redshift, SPH is in better agreement with SPHS and AMR but shows much poorer numerical convergence. We trace these discrepancies to artificial surface tension in SPH at phase boundaries. At early times, the passage of massive substructures close to the cluster centre stirs and shocks gas to build an entropy core. At later times, artificial surface tension causes low entropy gas to sink artificially to the centre of the cluster. We use SPHS to study the contribution of numerical versus physical dissipation on the entropy core, and argue that numerical dissipation is required to ensure single-valued fluid quantities in converging flows. However, provided this dissipation occurs only at the resolution limit, and provided that it does not propagate errors to larger scales, its effect is benign. There is no requirement to build `sub-grid models of unresolved turbulence for galaxy cluster simulations. We conclude that entropy cores in non-radiative simulations of galaxy clusters are physical, resulting from entropy generation in shocked gas during cluster assembly, putting to rest the long-standing puzzle of cluster entropy cores in AMR simulations versus their apparent absence in classic SPH simulations.
Log in to be able to interact and post comments
comments
Fetching comments
Sign in to be able to follow your search criteria
