No Arabic abstract
Stellar tidal streams are sensitive tracers of the properties of the gravitational potential in which they orbit and detailed observations of their density structure can be used to place stringent constraints on fluctuations in the potential caused by, e.g., the expected populations of dark matter subhalos in the standard cold dark matter paradigm (CDM). Simulations of the evolution of stellar streams in live $N$-body halos without low-mass dark-matter subhalos, however, indicate that streams exhibit significant perturbations on small scales even in the absence of substructure. Here we demonstrate, using high-resolution $N$-body simulations combined with sophisticated semi-analytic and simple analytic models, that the mass resolutions of $10^4$--$10^5,rm{M}_{odot}$ commonly used to perform such simulations cause spurious stream density variations with a similar magnitude on large scales as those expected from a CDM-like subhalo population and an order of magnitude larger on small, yet observable, scales. We estimate that mass resolutions of $approx100,rm{M}_{odot}$ ($approx1,rm{M}_{odot}$) are necessary for spurious, numerical density variations to be well below the CDM subhalo expectation on large (small) scales. That streams are sensitive to a simulations particle mass down to such small masses indicates that streams are sensitive to dark matter clustering down to these low masses if a significant fraction of the dark matter is clustered or concentrated in this way, for example, in MACHO models with masses of $10$--$100,rm{M}_{odot}$.
We describe a major upgrade of a Monte Carlo code which has previously been used for many studies of dense star clusters. We outline the steps needed in order to calibrate the results of the new Monte Carlo code against $N$-body simulations for large $N$ systems, up to $N=200000$. The new version of the Monte Carlo code (called MOCCA), in addition to the features of the old version, incorporates the direct Fewbody integrator (Fregeau et al. 2004) for three- and four-body interactions, and a new treatment of the escape process based on Fukushige & Heggie (2000). Now stars which fulfil the escape criterion are not removed immediately, but can stay in the system for a certain time which depends on the excess of the energy of a star above the escape energy. They are called potential escapers. With the addition of the Fewbody integrator the code can follow all interaction channels which are important for the rate of creation of various types of objects observed in star clusters, and ensures that the energy generation by binaries is treated in a manner similar to the $N$-body model. There are at most three new parameters which have to be adjusted against $N$-body simulations for large $N$: two (or one, depending on the chosen approach) connected with the escape process, and one responsible for the determination of the interaction probabilities. The values adopted for the free parameters have at most a weak dependence on $N$. They allow MOCCA to reproduce $N$-body results with reasonable precision, not only for the rate of cluster evolution and the cluster mass distribution, but also for the detailed distributions of mass and binding energy of binaries. Additionally, the code can follow the rate of formation of blue stragglers and black hole - black hole binaries.
We use direct $N$-body calculations to study the evolution of the unusually extended outer halo globular cluster Palomar 4 (Pal~4) over its entire lifetime in order to reproduce its observed mass, half-light radius, velocity dispersion and mass function slope at different radii. We find that models evolving on circular orbits, and starting from a non-mass segregated, canonical initial mass function (IMF) can reproduce neither Pal 4s overall mass function slope nor the observed amount of mass segregation. Including either primordial mass segregation or initially flattened IMFs does not reproduce the observed amount of mass segregation and mass function flattening simultaneously. Unresolved binaries cannot reconcile this discrepancy either. We find that only models with both a flattened IMF and primordial segregation are able to fit the observations. The initial (i.e. after gas expulsion) mass and half-mass radius of Pal~4 in this case are about 57000 M${odot}$ and 10 pc, respectively. This configuration is more extended than most globular clusters we observe, showing that the conditions under which Pal~4 formed must have been significantly different from that of the majority of globular clusters. We discuss possible scenarios for such an unusual configuration of Pal~4 in its early years.
We present a study of the effects of collisional dynamics on the formation and detectability of cold tidal streams. A semi-analytical model for the evolution of the stellar mass function was implemented and coupled to a fast stellar stream simulation code, as well as the synthetic cluster evolution code EMACSS for the mass evolution as a function of a globular cluster orbit. We find that the increase in the average mass of the escaping stars for clusters close to dissolution has a major effect on the observable stream surface density. As an example, we show that Palomar 5 would have undetectable streams (in an SDSS-like survey) if it was currently three times more massive, despite the fact that a more massive cluster loses stars at a higher rate. This bias due to the preferential escape of low-mass stars is an alternative explanation for the absence of tails near massive clusters, than a dark matter halo associated with the cluster. We explore the orbits of a large sample of Milky Way globular clusters and derive their initial masses and remaining mass fraction. Using properties of known tidal tails we explore regions of parameter space that favour the detectability of a stream. A list of high probability candidates is discussed
We investigate clustering properties of quasars using a new version of our semi-analytic model of galaxy and quasar formation with state-of-the-art cosmological N-body simulations. In this study, we assume that a major merger of galaxies triggers cold gas accretion on to a supermassive black hole and quasar activity. Our model can reproduce the downsizing trend of the evolution of quasars. We find that the median mass of quasar host dark matter haloes increases with cosmic time by an order of magnitude from z=4 (a few 1e+11 Msun) to z=1 (a few 1e+12 Msun), and depends only weakly on the quasar luminosity. Deriving the quasar bias through the quasar--galaxy cross-correlation function in the model, we find that the quasar bias does not depend on the quasar luminosity, similar to observed trends. This result reflects the fact that quasars with a fixed luminosity have various Eddington ratios and thus have various host halo masses that primarily determine the quasar bias. We also show that the quasar bias increases with redshift, which is in qualitative agreement with observations. Our bias value is lower than the observed values at high redshifts, implying that we need some mechanisms that make quasars inactive in low-mass haloes and/or that make them more active in high-mass haloes.
Hierarchical clustering represents the favoured paradigm for galaxy formation throughout the Universe; due to its proximity, the Magellanic system offers one of the few opportunities for astrophysicists to decompose the full six-dimensional phase-space history of a satellite in the midst of being cannibalised by its host galaxy. The availability of improved observational data for the Magellanic Stream and parallel advances in computational power has led us to revisit the canonical tidal model describing the disruption of the Small Magellanic Cloud and the consequent formation of the Stream. We suggest improvements to the tidal model in light of these recent advances.