Studying galaxy clusters through their Sunyaev-Zeldovich (SZ) imprint on the Cosmic Microwave Background has many important advantages. The total SZ signal is an accurate and precise tracer of the total pressure in the intra-cluster medium and of clu
ster mass, the key observable for using clusters as cosmological probes. Band 5 observations with SKA-MID towards cluster surveys from the next generation of X-ray telescopes such as e-ROSITA and from Euclid will provide the robust mass estimates required to exploit these samples. This will be especially important for high redshift systems, arising from the SZs unique independence to redshift. In addition, galaxy clusters are very interesting astrophysical systems in their own right, and the SKAs excellent surface brightness sensitivity down to small angular scales will allow us to explore the detailed gas physics of the intra-cluster medium.
We present results of cosmological simulations of disk galaxies carried out with the GADGET-3 TreePM+SPH code, where star formation and stellar feedback are described using our MUlti Phase Particle Integrator (MUPPI) model. This description is based
on simple multi-phase model of the interstellar medium at unresolved scales, where mass and energy flows among the components are explicitly followed by solving a system of ordinary differential equations. Thermal energy from SNe is injected into the local hot phase, so as to avoid that it is promptly radiated away. A kinetic feedback prescription generates the massive outflows needed to avoid the over-production of stars. We use two sets of zoomed-in initial conditions of isolated cosmological halos with masses (2-3) * 10^{12} Msun, both available at several resolution levels. In all cases we obtain spiral galaxies with small bulge-over-total stellar mass ratios (B/T approx 0.2), extended stellar and gas disks, flat rotation curves and realistic values of stellar masses. Gas profiles are relatively flat, molecular gas is found to dominate at the centre of galaxies, with star formation rates following the observed Schmidt-Kennicutt relation. Stars kinematically belonging to the bulge form early, while disk stars show a clear inside-out formation pattern and mostly form after redshift z=2. However, the baryon conversion efficiencies in our simulations differ from the relation given by Moster et al. (2010) at a 3 sigma level, thus indicating that our stellar disks are still too massive for the Dark Matter halo in which they reside. Results are found to be remarkably stable against resolution. This further demonstrates the feasibility of carrying out simulations producing a realistic population of galaxies within representative cosmological volumes, at a relatively modest resolution.
Forthcoming experiments will enable us to determine tomographic shear spectra at a high precision level. Most predictions about them have until now been biased on algorithms yielding the expected linear and non-linear spectrum of density fluctuations
. Even when simulations have been used, so-called Halofit (Smith et al 2003) predictions on fairly large scales have been needed. We wish to go beyond this limitation. We perform N-body and hydrodynamical simulations within a sufficiently large cosmological volume to allow a direct connection between simulations and linear spectra. While covering large length-scales, the simulation resolution is good enough to allow us to explore the high-l harmonics of the cosmic shear (up to l ~ 50000), well into the domain where baryon physics becomes important. We then compare shear spectra in the absence and in presence of various kinds of baryon physics, such as radiative cooling, star formation, and supernova feedback in the form of galactic winds. We distinguish several typical properties of matter fluctuation spectra in the different simulations and test their impact on shear spectra. We compare our outputs with those obtainable using approximate expressions for non--linear spectra, and identify substantial discrepancies even between our results and those of purely N-body results. Our simulations and the treatment of their outputs however enable us, for the first time, to obtain shear results taht are fully independent of any approximate expression, also in the high-l range, where we need to incorporate a non-linear power spectrum of density perturbations, and the effects of baryon physics. This will allow us to fully exploit the cosmological information contained in future high--sensitivity cosmic shear surveys, exploring the physics of cosmic shears via weak lensing measurements.
We present the results of a series of numerical simulations aimed to study the evolution of a disc galaxy within the global tidal field of a group environment. Both the disc galaxy and the group are modelled as multi-component, collision-less, N-body
systems, composed by both dark matter and stars. In our simulations, the evolution of disc galaxies is followed as their orbits sink towards the group centre, under the effect of dynamical friction. We explore a broad parameter space, covering several aspects of the galaxy-group interaction that are potentially relevant to galaxy evolution. Namely, prograde and retrograde orbits, orbital eccentricities, disc inclination, role of a central bulge in discs, internal disc kinematics, and galaxy-to-group mass ratios. We find that significant disc transformations occur only after the mean density of the group, measured within the orbit of the galaxy, exceeds ~0.3-1 times the central mean density of the galaxy. The morphological evolution of discs is found to be strongly dependent on the initial inclination of the disc with respect to its orbital plane. That is, discs on face-on and retrograde orbits are shown to retain longer their disc structures and kinematics, in comparison to prograde discs. This suggests that after interacting with the global tidal field alone, a significant fraction of disc galaxies should be found in the central regions of groups. Prominent central bulges are not produced, and pre-existing bulges are not enhanced in discs after the interaction with the group. Assuming that most S0 are formed in group environments, this implies that prominent bulges should be formed mostly by young stars, created only after a galaxy has been accreted by a group. Finally, contrary to some current implementations of tidal stripping in semi-analytical models of galaxy evolution, we find that more massive galaxies suffer more tidal stripping.
Dark matter halo merger trees are now routinely extracted from cosmological simulations of structure formation. These trees are frequently used as inputs to semi-analytic models of galaxy formation to provide the backbone within which galaxy formatio
n takes place. By necessity, these merger trees are constructed from a finite set of discrete snapshots of the N-body simulation and so have a limited temporal resolution. To date, there has been little consideration of how this temporal resolution affects the properties of galaxies formed within these trees. In particular, the question of how many snapshots are needed to achieve convergence in galaxy properties has not be answered. Therefore, we study the convergence in the stellar and total baryonic masses of galaxies, distribution of merger times, stellar mass functions and star formation rates in the Galacticus model of galaxy formation as a function of the number of snapshot times used to represent dark matter halo merger trees. When utilizing snapshots between z=20 and z=0, we find that at least 128 snapshots are required to achieve convergence to within 5% for galaxy masses. This convergence is obtained for mean quantities averaged over large samples of galaxies - significant variance for individual galaxies remains even when using very large numbers of snapshots. We find only weak dependence of the rate of convergence on the distribution of snapshots in time - snapshots spaced uniformly in the expansion factor, uniformly in the logarithm of expansion factor or uniformly in the logarithm of critical overdensity for collapse work equally well in almost all cases. We provide input parameters to Galacticus which allow this type of convergence study to be tuned to other simulations and to be carried out for other galaxy properties.
We present results based on an implementation of the Godunov Smoothed Particle Hydrodynamics (GSPH), originally developed by Inutsuka (2002), in the GADGET-3 hydrodynamic code. We first review the derivation of the GSPH discretization of the equation
s of moment and energy conservation, starting from the convolution of these equations with the interpolating kernel. The two most important aspects of the numerical implementation of these equations are (a) the appearance of fluid velocity and pressure obtained from the solution of the Riemann problem between each pair of particles, and (b the absence of an artificial viscosity term. We carry out three different controlled hydrodynamical three-dimensional tests, namely the Sod shock tube, the development of Kelvin-Helmholtz instabilities in a shear flow test, and the blob test describing the evolution of a cold cloud moving against a hot wind. The results of our tests confirm and extend in a number of aspects those recently obtained by Cha (2010): (i) GSPH provides a much improved description of contact discontinuities, with respect to SPH, thus avoiding the appearance of spurious pressure forces; (ii) GSPH is able to follow the development of gas-dynamical instabilities, such as the Kevin--Helmholtz and the Rayleigh-Taylor ones; (iii) as a result, GSPH describes the development of curl structures in the shear-flow test and the dissolution of the cold cloud in the blob test. We also discuss in detail the effect on the performances of GSPH of changing different aspects of its implementation. The results of our tests demonstrate that GSPH is in fact a highly promising hydrodynamic scheme, also to be coupled to an N-body solver, for astrophysical and cosmological applications. [abridged]
We list here the contents of the Proceedings of the Wide Field X-ray Telescope conference held in Bologna, Italy on 25-26 Nov 2009. The conference highlighted the scientific potential and discovery space provided by an X-ray mission concept character
ized by a wide field-of-view (1 sq.deg.), large effective area (1 sq.mt.) and approximately constant PSF (~5 arcsec HEW) across the whole FOV. The index is in html form with clickable links to the individual contributions.
We review recent progress in the description of the formation and evolution of galaxy clusters in a cosmological context by using numerical simulations. We focus our presentation on the comparison between simulated and observed X-ray properties, whil
e we will also discuss numerical predictions on properties of the galaxy population in clusters. Many of the salient observed properties of clusters, such as X-ray scaling relations, radial profiles of entropy and density of the intracluster gas, and radial distribution of galaxies are reproduced quite well. In particular, the outer regions of cluster at radii beyond about 10 per cent of the virial radius are quite regular and exhibit scaling with mass remarkably close to that expected in the simplest case in which only the action of gravity determines the evolution of the intra-cluster gas. However, simulations generally fail at reproducing the observed cool-core structure of clusters: simulated clusters generally exhibit a significant excess of gas cooling in their central regions, which causes an overestimate of the star formation and incorrect temperature and entropy profiles. The total baryon fraction in clusters is below the mean universal value, by an amount which depends on the cluster-centric distance and the physics included in the simulations, with interesting tensions between observed stellar and gas fractions in clusters and predictions of simulations. Besides their important implications for the cosmological application of clusters, these puzzles also point towards the important role played by additional physical processes, beyond those already included in the simulations. We review the role played by these processes, along with the difficulty for their implementation, and discuss the outlook for the future progress in numerical modeling of clusters.
