Back-in-time dynamics of the cluster IE 0657-56 (the Bullet System)

We present a simplified dynamical model of the ``Bullet'' system of two colliding clusters. The model constrains the masses of the system by requiring that the orbits of the main and sub components satisfy the cosmological initial conditions of vanishing physical separation a Hubble time ago. This is also known as the timing argument. The model considers a system embedded in an over-dense region. We argue that a relative speed of $4500 \rm km/s$ between the two components is consistent with cosmological conditions if the system is of a total mass of $2.8\times 10^{15}h^{-1} M_\odot$ is embedded in a region of a (mild) over-density of 10 times the cosmological background density. Combining this with the lensing measurements of the projected mass, the model yields a ratio of 3:1 for the mass of the main relative to that of the subcomponent. The effect of the background weakens as the relative speed between the two components is decreased. For relative speeds lower than $\sim 3700\rm km/s$, the timing argument yields masses which are too low to be consistent with lensing.Comment: 5 pages, 3 figures, submitted to MNRA

The importance of the merging activity for the kinetic polarization of the Sunyaev-Zel'dovich signal from galaxy clusters

The polarization sensitivity of the upcoming millimetric observatories will open new possibilities for studying the properties of galaxy clusters and for using them as powerful cosmological probes. For this reason it is necessary to investigate in detail the characteristics of the polarization signals produced by their highly ionized intra-cluster medium (ICM). This work is focussed on the polarization effect induced by the ICM bulk motions, the so-called kpSZ signal, which has an amplitude proportional to the optical depth and to the square of the tangential velocity. In particular we study how this polarization signal is affected by the internal dynamics of galaxy clusters and what is its dependence on the physical modelling adopted to describe the baryonic component. This is done by producing realistic kpSZ maps starting from the outputs of two different sets of high-resolution hydrodynamical N-body simulations. The first set (17 objects) follows only non-radiative hydrodynamics, while for each of 9 objects of the second set we implement four different kinds of physical processes. Our results shows that the kpSZ signal turns out to be a very sensitive probe of the dynamical status of galaxy clusters. We find that major merger events can amplify the signal up to one order of magnitude with respect to relaxed clusters, reaching amplitude up to about 100 nuK. This result implies that the internal ICM dynamics must be taken into account when evaluating this signal because simplicistic models, based on spherical rigid bodies, may provide wrong estimates. Finally we find that the dependence on the physical modelling of the baryonic component is relevant only in the very inner regions of clusters.Comment: 13 pages, 7 figures, submitted to A&

Tracing early structure formation with massive starburst galaxies and their implications for reionization

Cosmological hydrodynamic simulations have significantly improved over the past several years, and we have already shown that the observed properties of Lyman-break galaxies (LBGs) at z=3 can be explained well by the massive galaxies in the simulations. Here we extend our study to z=6 and show that we obtain good agreement for the LBGs at the bright-end of the luminosity function (LF). Our simulations also suggest that the cosmic star formation rate density has a peak at z= 5-6, and that the current LBG surveys at z=6 are missing a significant number of faint galaxies that are dimmer than the current magnitude limit. Together, our results suggest that the universe could be reionized at z=6 by the Pop II stars in ordinary galaxies. We also estimate the LF of Lyman-alpha emitters (LAEs) at z=6 by relating the star formation rate in the simulation to the Ly-alpha luminosity. We find that the simulated LAE LFs agree with the observed data provided that the net escape fraction of Ly-alpha photon is f_{Ly-alpha} <= 0.1. We investigate two possible scenarios for this effect: (1) all sources in the simulation are uniformly dimmer by a factor of 10 through attenuation, and (2) one out of ten LAEs randomly lights up at a given moment. We show that the correlation strength of the LAE spatial distribution can possibly distinguish the two scenarios.Comment: 9 pages, 4 figures. Summary of the talk given at the "First Light & Reionization" workshop at UC Irvine, May 2005. The published article is available from http://dx.doi.org/10.1016/j.newar.2005.11.00

Hydrodynamical Simulations of the IGM at High Mach Numbers

We present a new approach to doing Eulerian computational fluid dynamics that is designed to work at high Mach numbers encountered in hydrodynamical simulations of the IGM. In conventional Eulerian CFD, the thermal energy is poorly tracked in supersonic bulk flows where local fluid variables cannot be accurately separated from the much larger bulk flow components. We described a method in which local fluid quantities can be directly tracked and the Eulerian fluid equations solved in a local frame moving with the flow. The new algorithm has been used to run large hydrodynamical simulations on a 1024^3 grid to study the kinetic SZ effect. The KSZ power spectrum is broadly peaked at l~10^4 with temperature fluctuations on micro Kelvin levels.Comment: 6 pages, to appear in the Proc. from the IGM/Galaxy Connection conferenc

Past and Future of CG J1720-67.8: Constraints from Observations and Models

We discuss the evolution of the peculiar, nearby (z = 0.045), compact galaxy group CG J1720-67.8, by interpreting a large amount of observational information on the basis of our recent results from spectrophotometric evolutionary synthesis models and new N-body/SPH simulations. The group, that is composed of two spiral galaxies with a mass ratio approximately 4:1 and an S0 galaxy in a particularly compact configuration, is undergoing an active pre-merging phase. Several tidal features are signposts of the complex dynamics of the system. We suggest that the observed structure of the tidal features can be explained only if all three galaxies are involved in a strong interaction process.Comment: 5 pages, 3 (degraded) figures. Proc. ESO Workshop "Groups of galaxies in the nearby Universe", Santiago, Chile, 5-9 Dec. 2005, ESO Astrophysics Symposia, eds. I. Saviane, V. Ivanov & J. Borissova, Springer-Verla

Galaxy formation in the Planck cosmology - II. Star-formation histories and post-processing magnitude reconstruction

We adapt the L-Galaxies semi-analytic model to follow the star-formation histories (SFH) of galaxies -- by which we mean a record of the formation time and metallicities of the stars that are present in each galaxy at a given time. We use these to construct stellar spectra in post-processing, which offers large efficiency savings and allows user-defined spectral bands and dust models to be applied to data stored in the Millennium data repository. We contrast model SFHs from the Millennium Simulation with observed ones from the VESPA algorithm as applied to the SDSS-7 catalogue. The overall agreement is good, with both simulated and SDSS galaxies showing a steeper SFH with increased stellar mass. The SFHs of blue and red galaxies, however, show poor agreement between data and simulations, which may indicate that the termination of star formation is too abrupt in the models. The mean star-formation rate (SFR) of model galaxies is well-defined and is accurately modelled by a double power law at all redshifts: SFR proportional to $1/(x^{-1.39}+x^{1.33})$, where $x=(t_a-t)/3.0\,$Gyr, $t$ is the age of the stars and $t_a$ is the loopback time to the onset of galaxy formation; above a redshift of unity, this is well approximated by a gamma function: SFR proportional to $x^{1.5}e^{-x}$, where $x=(t_a-t)/2.0\,$Gyr. Individual galaxies, however, show a wide dispersion about this mean. When split by mass, the SFR peaks earlier for high-mass galaxies than for lower-mass ones, and we interpret this downsizing as a mass-dependence in the evolution of the quenched fraction: the SFHs of star-forming galaxies show only a weak mass dependence.Comment: Accepted version of the paper, to appear in MNRAS. Compared to the original version, contains more detail on the post-processing of magnitudes, including a table of rms magnitude errors. SFHs available on Millennium database http://gavo.mpa-garching.mpg.de/MyMillennium

Resolving Subhaloes' Lives with the Hierarchical Bound-Tracing Algorithm

We develop a new code, the Hierarchical Bound-Tracing (HBT for short) code, to find and trace dark matter subhaloes in simulations based on the merger hierarchy of dark matter haloes. Application of this code to a recent benchmark test of finding subhaloes demonstrates that HBT stands as one of the best codes to trace the evolutionary history of subhaloes. The success of the code lies in its careful treatment of the complex physical processes associated with the evolution of subhaloes and in its robust unbinding algorithm with an adaptive source subhalo management. We keep a full record of the merger hierarchy of haloes and subhaloes, and allow growth of satellite subhaloes through accretion from its "satellite-of-satellites", hence allowing mergers among satellites. Local accretion of background mass is omitted, while rebinding of stripped mass is allowed. The justification of these treatments is provided by case studies of the lives of individual subhaloes and by the success in finding the complete subhalo catalogue. We compare our result to other popular subhalo finders and show that HBT is able to well resolve subhaloes in high density environment and keep strict physical track of subhaloes' merger history. This code is fully parallelized and freely available upon request to the authors.Comment: 15 pages; accepted for publication by MNRA