N-body realizations of cuspy dark matter haloes

We describe an algorithm for generating equilibrium initial conditions for numerical experiments with dark matter haloes. Our haloes are modelled using a general form for the mass density p{r), making it possible to represent most of the popular density profiles in the literature. The finite mass 7-models and the cuspy density profiles found in recent high-resolution cosmological TV-body simulations having a density power-law fall-off at large distances proportional to are included as special cases. The algorithm calculates the phase-space distribution function of each model assuming spherical symmetry and either an isotropic velocity dispersion tensor or an anisotropic velocity dispersion tensor of the type proposed by Osipkov and Merritt. The particle velocities are assigned according to the exact velocity distribution, making this method ideal for experiments requiring a high degree of stability. Numerical tests confirm that the resulting models are highly stable. This approach is motivated by the instabilities that arise when a local Maxwellian velocity distribution is adopted. For example, after approximating the velocity distribution by a Gaussian we show that a Hernquist halo with an initial r(^-1) density cusp immediately develops a constant density core. Moreover, after a single crossing time the orbital anisotropy has evolved over the entire system. Previous studies that use this approximation to construct halo or galaxy models could be compromised by this behaviour. Using the derived distribution functions we show the exact 1-d velocity distributions and we compare them with the Gaussian velocity distributions with the same second moment for different distances from the halo centre. We show that instabilities arise because a Gaussian velocity distribution is a very poor approximation to the true velocity distribution of particles. We also perform a series of numerical simulations evolving several dark matter halo models in isolation, with the intention of checking the stability of the initialization procedure in both configuration and velocity space. A subset of the models are evolved under the assumption that the velocity distribution at any given point is a Gaussian and the time evolution of the density profiles and velocity structure is monitored. Finally, a number of applications are discussed, including issues of relaxation in dark matter haloes as well as mergers of haloes in scattering experiments

Cold Dark Matter Substructure and Galactic Disks

We perform a set of high-resolution, dissipationless N-body simulations to investigate the influence of cold dark matter (CDM) substructure on the dynamical evolution of thin galactic disks. Our method combines cosmological simulations of galaxy-sized CDM halos to derive the properties of substructure populations and controlled numerical experiments of consecutive subhalo impacts onto initially-thin, fully-formed disk galaxies. We demonstrate that close encounters between massive subhalos and galactic disks since z~1 should be common occurrences in LCDM models. In contrast, extremely few satellites in present-day CDM halos are likely to have a significant impact on the disk structure. One typical host halo merger history is used to seed controlled N-body experiments of subhalo-disk encounters. As a result of these accretion events, the disk thickens considerably at all radii with the disk scale height increasing in excess of a factor of 2 in the solar neighborhood. We show that interactions with the subhalo population produce a wealth of distinctive morphological signatures in the disk stars including: conspicuous flares; bars; low-lived, ring-like features in the outskirts; and low-density, filamentary structures above the disk plane. We compare a resulting dynamically-cold, ring-like feature in our simulations to the Monoceros ring stellar structure in the MW. The comparison shows quantitative agreement in both spatial distribution and kinematics, suggesting that such observed complex stellar components may arise naturally as disk stars are excited by encounters with subhalos. These findings highlight the significant role of CDM substructure in setting the structure of disk galaxies and driving galaxy evolution.Comment: 10 pages, 4 figures. To appear in the proceedings of the IAU Symposium No. 254 "The Galaxy Disk in Cosmological Context", Copenhagen 9-13 June 2008, Denmark, (Eds.) J. Andersen, J. Bland-Hawthorn & B. Nordstrom, Cambridge University Pres

Biases in mass estimates of dSph galaxies

Using a high resolution N-body simulation of a two-component dwarf galaxy orbiting in the potential of the Milky Way, we study two effects that lead to significant biases in mass estimates of dwarf spheroidal galaxies. Both are due to the strong tidal interaction of initially disky dwarfs with their host. The tidal stripping of dwarf stars leads to the formation of strong tidal tails that are typically aligned with the line of sight of an observer positioned close to the host center. The stars from the tails contaminate the kinematic samples leading to a velocity dispersion profile increasing with the projected radius and resulting in an overestimate of mass. The tidal stirring of the dwarf leads to the morphological transformation of the initial stellar disk into a bar and then a spheroid. The distribution of stars in the dwarf remains non-spherical for a long time leading to an overestimate of its mass if it is observed along the major axis and an underestimate if it seen in the perpendicular direction.Comment: 5 pages, 3 figures, contribution to the proceedings of "Hunting for the Dark: The Hidden Side of Galaxy Formation", Malta, 19-23 Oct. 2009, eds. V.P. Debattista & C.C. Popescu, AIP Conference Series, in pres

The stellar structure and kinematics of dwarf spheroidal galaxies formed by tidal stirring

Using high-resolution N-body simulations we study the stellar properties of dwarf spheroidal galaxies resulting from the tidally induced morphological transformation of disky dwarfs on a cosmologically motivated eccentric orbit around the Milky Way. Dwarf galaxy models initially consist of an exponential stellar disk embedded in an extended spherical dark matter halo. Depending on the initial orientation of the disk with respect to the orbital plane, different final configurations are obtained. The least evolved dwarf is triaxial and retains a significant amount of rotation. The more evolved dwarfs are prolate spheroids with little rotation. We show that the final density distribution of stars can be approximated by a simple modification of the Plummer law. The kinematics of the dwarfs is significantly different depending on the line of sight which has important implications for mapping the observed stellar velocity dispersions of dwarfs to subhalo circular velocities. When the dwarfs are observed along the long axis, the measured velocity dispersion is higher and decreases faster with radius. In the case where rotation is significant, when viewed perpendicular to the long axis, the effect of minor axis rotation is detected, as expected for triaxial systems. We model the velocity dispersion profiles and rotation curves of the dwarfs by solving the Jeans equations for spherical and axisymmetric systems and adjusting different sets of free parameters. We find that the mass is typically overestimated when the dwarf is seen along the long axis and underestimated when the observation is along the short or intermediate axis. The effect of non-sphericity cannot however bias the inferred mass by more than 60 percent in either direction, even for the most strongly stripped dwarf which is close to disruption.Comment: 17 pages, 15 figures, revised version accepted for publication in Ap

Black Hole Starvation and Bulge Evolution in a Milky Way-like Galaxy

We present a new zoom-in hydrodynamical simulation, "Erisbh", which follows the cosmological evolution and feedback effects of a supermassive black hole at the center of a Milky Way-type galaxy. ErisBH shares the same initial conditions, resolution, recipes of gas cooling, star formation and feedback, as the close Milky Way-analog "Eris", but it also includes prescriptions for the formation, growth and feedback of supermassive black holes. We find that the galaxy's central black hole grows mainly through mergers with other black holes coming from infalling satellite galaxies. The growth by gas accretion is minimal because very little gas reaches the sub-kiloparsec scales. The final black hole is, at z=0, about 2.6 million solar masses and it sits closely to the position of SgrA* on the MBH-MBulge and MBH-sigma planes, in a location consistent with what observed for pseudobulges. Given the limited growth due to gas accretion, we argue that the mass of the central black hole should be above 10^5 solar masses already at z~8. The effect of AGN feedback on the host galaxy is limited to the very central few hundreds of parsecs. Despite being weak, AGN feedback seems to be responsible for the limited growth of the central bulge with respect to the original Eris, which results in a significantly flatter rotation curve in the inner few kiloparsecs. Moreover, the disk of ErisBH is more prone to instabilities, as its bulge is smaller and its disk larger then Eris. As a result, the disk of ErisBH undergoes a stronger dynamical evolution relative to Eris and around z=0.3 a weak bar grows into a strong bar of a few disk scale lengths in size. The bar triggers a burst of star formation in the inner few hundred parsecs, provides a modest amount of new fuel to the central black hole, and causes the bulge of ErisBH to have, by z=0, a box/peanut morphology.(Abridged)Comment: 16 pages, 16 figures. Submitted to MNRA