The effect of radiative cooling on scaling laws of X-ray groups and clusters

We have performed cosmological simulations in a ΛCDM cosmology with and without radiative cooling in order to study the effect of cooling on the cluster scaling laws. Our simulations consist of 4.1 million particles each of gas and dark matter within a box size of 100 h-1 Mpc, and the run with cooling is the largest of its kind to have been evolved to z = 0. Our cluster catalogs both consist of over 400 objects and are complete in mass down to ~1013 h-1 M☉. We contrast the emission-weighted temperature-mass (Tew-M) and bolometric luminosity-temperature (Lbol-Tew) relations for the simulations at z = 0. We find that radiative cooling increases the temperature of intracluster gas and decreases its total luminosity, in agreement with the results of Pearce et al. Furthermore, the temperature dependence of these effects flattens the slope of the Tew-M relation and steepens the slope of the Lbol-Tew relation. Inclusion of radiative cooling in the simulations is sufficient to reproduce the observed X-ray scaling relations without requiring excessive nongravitational energy injection

Simulation of primordial object formation

We have included the chemical rate network responsible for the formation of molecular Hydrogen in the N-body hydrodynamic code, Hydra, in order to study the formation of the first cosmological at redshifts between 10 and 50. We have tested our implementation of the chemical and cooling processes by comparing N-body top hat simulations with theoretical predictions from a semi-analytic model and found them to be in good agreement. We find that post-virialization properties are insensitive to the initial abundance of molecular hydrogen. Our main objective was to determine the minimum mass ($M_{SG}(z)$) of perturbations that could become self gravitating (a prerequisite for star formation), and the redshift at which this occurred. We have developed a robust indicator for detecting the presence of a self-gravitating cloud in our simulations and find that we can do so with a baryonic particle mass-resolution of 40 solar masses. We have performed cosmological simulations of primordial objects and find that the object's mass and redshift at which they become self gravitating agree well with the $M_{SG}(z)$ results from the top hat simulations. Once a critical molecular hydrogen fractional abundance of about 0.0005 has formed in an object, the cooling time drops below the dynamical time at the centre of the cloud and the gas free falls in the dark matter potential wells, becoming self gravitating a dynamical time later.Comment: 45 pages, 17 figures, submitted to Ap

Constraining global properties of the Draco dwarf spheroidal galaxy

By fitting a flexible stellar anisotropy model to the observed surface brightness and line-of-sight velocity dispersion profiles of Draco we derive a sequence of cosmologically plausible two-component (stars + dark matter) models for this galaxy. The models are consistent with all the available observations and can have either cuspy Navarro-Frenk-White or flat-cored dark matter density profiles. The dark matter halos either formed relatively recently (at z~2...7) and are massive (up to ~5x10^9 M_Sun), or formed before the end of the reionization of the universe (z~7...11) and are less massive (down to ~7x10^7 M_Sun). Our results thus support either of the two popular solutions of the "missing satellites" problem of Lambda cold dark matter cosmology - that dwarf spheroidals are either very massive, or very old. We carry out high-resolution simulations of the tidal evolution of our two-component Draco models in the potential of the Milky Way. The results of our simulations suggest that the observable properties of Draco have not been appreciably affected by the Galactic tides after 10 Gyr of evolution. We rule out Draco being a "tidal dwarf" - a tidally disrupted dwarf galaxy. Almost radial Draco orbits (with the pericentric distance <15 kpc) are also ruled out by our analysis. The case of a harmonic dark matter core can be consistent with observations only for a very limited choice of Draco orbits (with the apocentric-to-pericentric distances ratio of <2.5).Comment: 18 pages, 14 figures; accepted by Ap

Peculiar Velocities of Galaxy Clusters

We investigate the peculiar velocities predicted for galaxy clusters by theories in the cold dark matter family. A widely used hypothesis identifies rich clusters with high peaks of a suitably smoothed version of the linear density fluctuation field. Their peculiar velocities are then obtained by extrapolating the similarly smoothed linear peculiar velocities at the positions of these peaks. We test these ideas using large high resolution N-body simulations carried out within the Virgo supercomputing consortium. We find that at early times the barycentre of the material which ends up in a rich cluster is generally very close to a high peak of the initial density field. Furthermore the mean peculiar velocity of this material agrees well with the linear value at the peak. The late-time growth of peculiar velocities is, however, systematically underestimated by linear theory. At the time clusters are identified we find their rms peculiar velocity to be about 40% larger than predicted. Nonlinear effects are particularly important in superclusters. These systematics must be borne in mind when using cluster peculiar velocities to estimate the parameter combination $\sigma_8\Omega^{0.6}$.Comment: 8 pages, 4 figures; submitted to MNRA

The Bispectrum as a Signature of Gravitational Instability in Redshift-Space

The bispectrum provides a characteristic signature of gravitational instability that can be used to probe the Gaussianity of the initial conditions and the bias of the galaxy distribution. We study how this signature is affected by redshift distortions using perturbation theory and high-resolution numerical simulations. We obtain perturbative results for the multipole expansion of the redshift-space bispectrum which provide a natural way to break the degeneracy between bias and $\Omega$ present in measurements of the redshift-space power spectrum. We propose a phenomenological model that incorporates the perturbative results and also describes the bispectrum in the transition to the non-linear regime. We stress the importance of non-linear effects and show that inaccurate treatment of these can lead to significant discrepancies in the determination of bias from galaxy redshift surveys. At small scales we find that the bispectrum monopole exhibits a strong configuration dependence that reflects the velocity dispersion of clusters. Therefore, the hierarchical model for the three-point function does not hold in redshift-space.Comment: 19 pages, 4 figures. Revised version accepted for publication in Ap

Star Formation, Supernovae Feedback and the Angular Momentum Problem in Numerical CDM Cosmogony: Half Way There?

We present a smoothed particle hydrodynamic (SPH) simulation that reproduces a galaxy that is a moderate facsimile of those observed. The primary failing point of previous simulations of disk formation, namely excessive transport of angular momentum from gas to dark matter, is ameliorated by the inclusion of a supernova feedback algorithm that allows energy to persist in the model ISM for a period corresponding to the lifetime of stellar associations. The inclusion of feedback leads to a disk at a redshift $z=0.52$, with a specific angular momentum content within 10% of the value required to fit observations. An exponential fit to the disk baryon surface density gives a scale length within 17% of the theoretical value. Runs without feedback, with or without star formation, exhibit the drastic angular momentum transport observed elsewhere.Comment: 4 pages, 3 figures, accepted for publication in ApJ Letter

Ray splitting in paraxial optical cavities

We present a numerical investigation of the ray dynamics in a paraxial optical cavity when a ray splitting mechanism is present. The cavity is a conventional two-mirror stable resonator and the ray splitting is achieved by inserting an optical beam splitter perpendicular to the cavity axis. We show that depending on the position of the beam splitter the optical resonator can become unstable and the ray dynamics displays a positive Lyapunov exponent.Comment: 13 pages, 7 figures, 1 tabl