The dearth of dwarf galaxies in the local universe is hard to reconcile with the large number of low mass haloes expected within the concordance $\Lambda$CDM paradigm. In this paper we perform a systematic evaluation of the uncertainties affecting the measurement of DM halo abundance using galaxy kinematics. Using a large sample of dwarf galaxies with spatially resolved kinematic data we derive a correction to obtain the observed abundance of galaxies as a function of their halo maximum circular velocity from the line-of-sight velocity function in the Local Volume. This estimate provides a direct means of comparing the predictions of theoretical models and simulations (including nonstandard cosmologies and novel galaxy formation physics) to the observational constraints. The new "galactic $V_{max}$" function is steeper than the line-of-sight velocity function but still shallower than the theoretical CDM expectation, showing that some unaccounted physical process is necessary to reduce the abundance of galaxies and/or drastically modify their density profiles compared to CDM haloes. Using this new galactic $V_{max}$ function, we investigate the viability of baryonic solutions such as feedback-powered outflows and photoevaporation of gas from an ionising radiation background. At the 3-$\sigma$ confidence level neither energetic feedback nor photoevaporation are effective enough to reconcile the disagreement. In the case of maximum baryonic effects, the theoretical estimate still deviates significantly from the observations for $V_{max} < 60$ km/s. CDM predicts at least 1.8 times more galaxies with $V_{max} = 50$ km/s and 2.5 times more than observed at $30$ km/s. Recent hydrodynamic simulations seem to resolve the discrepancy but disagree with the properties of observed galaxies with resolved kinematics. (abridged)

The observed number of dwarf galaxies as a function of rotation velocity is significantly smaller than predicted by the standard model of cosmology. This discrepancy cannot be simply solved by assuming strong baryonic feedback processes, since they would violate the observed relation between maximum circular velocity ($v_{\rm max}$) and baryon mass of galaxies. A speculative but tantalising possibility is that the mismatch between observation and theory points towards the existence of non-cold or non-collisionless dark matter (DM). In this paper, we investigate the effects of warm, mixed (i.e warm plus cold), and self-interacting DM scenarios on the abundance of dwarf galaxies and the relation between observed HI line-width and maximum circular velocity. Both effects have the potential to alleviate the apparent mismatch between the observed and theoretical abundance of galaxies as a function of $v_{\rm max}$. For the case of warm and mixed DM, we show that the discrepancy disappears, even for luke-warm models that evade stringent bounds from the Lyman-$\alpha$ forest. Self-interacting DM scenarios can also provide a solution as long as they lead to extended ($\gtrsim 1.5$ kpc) dark matter cores in the density profiles of dwarf galaxies. Only models with velocity-dependent cross sections can yield such cores without violating other observational constraints at larger scales.

Merging has been proposed to explain multiple populations in globular clusters (GCs) where there is a spread in iron abundance (hereafter, iron-complex GCs). By means of N-body simulations, we investigate if merging is consistent with the observations of sub-populations and rotation in iron-complex GCs. The key parameters are the initial mass and density ratios of the progenitors. When densities are similar, the more massive progenitor dominates the central part of the merger remnant and the less massive progenitor forms an extended rotating population. The low-mass progenitor can become the majority population in the central regions of the merger remnant only if its initial density is higher by roughly the mass ratio. To match the radial distribution of multiple populations in two iron-complex GCs ({\omega} Cen and NGC 1851), the less massive progenitor needs to be four times as dense as the larger one. Our merger remnants show solid-body rotation in the inner parts, becoming differential in the outer parts. Rotation velocity V and ellipticity {\epsilon} are in agreement with models for oblate rotators with isotropic dispersion. We discuss several kinematic signatures of a merger with a denser lower mass progenitor that can be tested with future observations.

Metal poor globular clusters (MPGCs) are a unique probe of the early universe, in particular the reionization era. Systems of globular clusters in galaxy clusters are particularly interesting as it is in the progenitors of galaxy clusters that the earliest reionizing sources first formed. Although the exact physical origin of globular clusters is still debated, it is generally admitted that globular clusters form in early, rare dark matter peaks (Moore et al. 2006; Boley et al. 2009). We provide a fully numerical analysis of the Virgo cluster globular cluster system by identifying the present day globular cluster system with exactly such early, rare dark matter peaks. A popular hypothesis is that that the observed truncation of blue metal poor globular cluster formation is due to reionization (Spitler et al. 2012; Boley et al. 2009; Brodie & Strader 2006); adopting this view, constraining the formation epoch of MPGCs provides a complementary constraint on the epoch of reionization. By analyzing both the line of sight velocity dispersion and the surface density distribution of the present day distribution we are able to constrain the redshift and mass of the dark matter peaks. We find and quantify a dependence on the chosen line of sight of these quantities, whose strength varies with redshift, and coupled with star formation efficiency arguments find a best fitting formation mass and redshift of $\simeq 5 \times 10^8 \rm{M}_\odot$ and $z\simeq 9$. We predict $\simeq 300$ intracluster MPGCs in the Virgo cluster. Our results confirm the techniques pioneered by Moore et al. (2006) when applied to the the Virgo cluster and extend and refine the analytic results of Spitler et al. (2012) numerically.

Tidal streams are a powerful probe of the Milky Way (MW) potential shape. In this paper, we introduce a simple test particle method to fit stream data, using a Markov Chain Monte Carlo technique to marginalise over uncertainties in the progenitor's orbit and the Milky Way halo shape parameters. Applying it to mock data of thin streams in the MW halo, we show that, even for very cold streams, stream-orbit offsets - not modelled in our simple method - introduce systematic biases in the recovered shape parameters. For the streams that we consider, and our particular choice of potential parameterisation, these errors are of order ~20% on the halo flattening parameters. However, larger systematic errors can arise for more general streams and potentials; such offsets need to be correctly modelled in order to obtain an unbiased recovery of the underlying potential. Assessing which of the known Milky Way streams are most constraining, we find NGC 5466 and Pal 5 are the most promising candidates. These form an interesting pair as their orbital planes are both approximately perpendicular to each other and to the disc, giving optimal constraints on the MW halo shape. We show that - while with current data their constraints on potential parameters are poor - good radial velocity data along the Pal 5 stream will provide constraints on qz - the flattening perpendicular to the disc. Furthermore, as discussed in a companion paper, NGC 5466 can provide rather strong constraints on the MW halo shape parameters, if the tentative evidence for a departure from the smooth orbit towards its western edge is confirmed.

Every Galactic environment is characterized by a stellar density and a velocity dispersion. With this information from literature, we simulated flyby encounters for several Galactic regions, numerically calculating stellar trajectories as well as orbits for particles in disks; our aim was to understand the effect of typical stellar flybys on planetary (debris) disks in the Milky Way Galaxy. For the Solar neighborhood, we examined nearby stars with known distance, proper motions, and radial velocities. We found occurrence of a disturbing impact to the Solar planetary disk within the next 8 Myr to be highly unlikely; perturbations to the Oort cloud seem unlikely as well. Current knowledge of the full phase space of stars in the Solar neighborhood, however, is rather poor, and thus we cannot rule out the existence of a star that is more likely to approach than those for which we have complete kinematic information. We studied the effect of stellar encounters on planetary orbits within the habitable zones of stars in more crowded stellar environments, such as stellar clusters. We found that in open clusters habitable zones are not readily disrupted; this is true if they evaporate in less than 108 years. For older clusters the results may not be the same. We specifically studied the case of Messier 67, one of the oldest open clusters known, and show the effect of this environment on debris disks. We also considered the conditions in globular clusters, the Galactic nucleus, and the Galactic bulge-bar. We calculated the probability of whether Oort clouds exist in these Galactic environments.

We apply a new method to determine the local disc matter and dark halo matter density to kinematic and position data for \sim2000 K dwarf stars taken from the literature. Our method assumes only that the disc is locally in dynamical equilibrium, and that the 'tilt' term in the Jeans equations is small up to \sim1 kpc above the plane. We present a new calculation of the photometric distances to the K dwarf stars, and use a Monte Carlo Markov Chain to marginalise over uncertainties in both the baryonic mass distribution, and the velocity and distance errors for each individual star. We perform a series of tests to demonstrate that our results are insensitive to plausible systematic errors in our distance calibration, and we show that our method recovers the correct answer from a dynamically evolved N-body simulation of the Milky Way. We find a local dark matter density of {\rho}dm = 0.025+0.014-0.013 M\odotpc^{-3} (0.95+0.53-0.49 GeV cm^{-3}) at 90% confidence assuming no correction for the non-flatness of the local rotation curve, and {\rho}dm = 0.022+0.015-0.013 M\odotpc^-3 (0.85+0.57-0.50 GeV cm^{-3}) if the correction is included. Our 90% lower bound on {\rho}dm is larger than the canonical value typically assumed in the literature, and is at mild tension with extrapolations from the rotation curve that assume a spherical halo. Our result can be explained by a larger normalisation for the local Milky Way rotation curve, an oblate dark matter halo, a local disc of dark matter, or some combination of these.

Stellar streams provide unique probes of galactic potentials, with the longer streams normally providing the cleaner measurements. In this paper, we show an example of a short tidal stream that is particularly sensitive to the shape of the Milky Way's dark matter halo: the globular cluster tidal stream NGC 5466. This stream has an interesting deviation from a smooth orbit at its western edge. We show that such a deviation favours an underlying oblate or triaxial halo (irrespective of plausible variations in the Milky Way disc properties and the specific halo parametrisation chosen); spherical or prolate halo shapes can be excluded at a high confidence level. Therefore, more extensive data sets along the NGC 5466 tidal stream promise strong constraints on the Milky Way halo shape.

We study the systematic problems in determining the local dark matter density $\rho_{dm}(R_\odot)$ from kinematics of stars in the Solar Neighbourhood, using a simulated Milky Way-like galaxy. We introduce a new unbiased method for recovering $\rho_{dm}(R_\odot)$ based on the moments of the Jeans equations, combined with a Monte Carlo Markov Chain technique and apply it to real data. ov

We revisit systematics in determining the local dark matter density (rho_dm) from the vertical motion of stars in the Solar Neighbourhood. Using a simulation of a Milky Way-like galaxy, we determine the data-quality required to detect the dark matter density at its expected local value. We introduce a new method for recovering rho_dm that uses moments of the Jeans equations, combined with a Monte Carlo Markov Chain technique to marginalise over the unknown parameters. Given sufficiently good data, we show that our method can recover the correct local dark matter density even in the face of disc inhomogeneities, non-isothermal tracers and a non-separable distribution function. We illustrate the power of our technique by applying it to Hipparcos data [Holmberg & Flynn 2000,2004]. We first make the assumption that the A and F star tracer populations are isothermal. This recovers rho_dm=0.003^{+0.009}_{-0.007}Msun/pc^3 (with 90 per cent confidence), consistent with previous determinations. However, the vertical dispersion profile of these tracers is poorly known. If we assume instead a non-isothermal profile similar to the blue disc stars from SDSS DR-7 [Abazajian et al. 2009] measured by Bond et al. (2009), we obtain a fit with a very similar chi^2 value, but with rho_dm=0.033^{+0.008}_{-0.009}Msun/pc^3 (with 90 per cent confidence). This highlights that it is vital to measure the vertical dispersion profile of the tracers to recover an unbiased estimate of the local dark matter density.

We have computed trajectories, distances and times of closest approaches to the Sun by stars in the Solar neighbourhood with known position, radial velocity and proper motions. For this purpose we have used a full potential model of the Galaxy that reproduces the local z-force, the Oort constants, the local escape velocity, and the rotation curve of the Galaxy. From our sample we constructed initial conditions, within observational uncertainties, with a Monte Carlo scheme for the twelve most suspicious candidates because of their small tangential motion. We find that the star Gliese 710 will have the closest approach to the Sun, with a distance of approximately 0.34 pc at 1.36 Myr in the future. We show that the effect of a flyby with the characteristics of Gliese 710 on a 100 AU test particle disk representing the Solar system is negligible. However, since there is a lack of 6D data for a large percentage of stars in the Solar neighbourhood, closer approaches may exist. We calculate parameters of passing stars that would cause noticeable effects on the Solar disk. Regarding the birth cloud of the Sun, we performed experiments to reproduce roughly the observed orbital parameters such as eccentricities and inclinations of the Kuiper Belt. It is known now that in Galactic environments, such as stellar formation regions, the stellar densities of new born stars, are high enough to produce close encounters within 200 AU. Moreover, in these Galactic environments, the velocity dispersion is relatively low, typically approximately 1-3 km s-1.We find that with a velocity dispersion of approximately 1 km s-1 and an approach distance of about 150 AU, typical of these regions, we obtain approximately the eccentricities and inclinations seen in the current Solar system. Simple analytical calculations of stellar encounters effects on the Oort cloud are presented.

We examine systematic problems in determining the local matter density from the vertical motion of stars, i.e. the 'Oort limit'. Using collisionless simulations and a Monte Carlo Markov Chain technique, we determine the data quality required to detect local dark matter at its expected density. We find that systematic errors are more important than observational errors and apply our technique to Hipparcos data to reassign realistic error bars to the local dark matter density.

Making robust predictions for the phase space distribution of dark matter at the solar neighbourhood is vital for dark matter direct detection experiments. To date, almost all such predictions have been based on simulations that model the dark matter alone. Here, we use three cosmological hydrodynamics simulations of bright, disc dominated galaxies to include the effects of baryonic matter self-consistently for the first time. We find that the addition of baryonic physics drastically alters the dark matter profile in the vicinity of the Solar neighbourhood. A stellar/gas disc, already in place at high redshift, causes merging satellites to be dragged preferentially towards the disc plane where they are torn apart by tides. This results in an accreted dark matter disc that contributes ~0.25 - 1.5 times the non-rotating halo density at the solar position. The dark disc, unlike dark matter streams, is an equilibrium structure that must exist in disc galaxies that form in a hierarchical cosmology. Its low rotation lag with respect to the Earth significantly boosts WIMP capture in the Earth and Sun, boosts the annual modulation signal, and leads to distinct variations in the flux as a function of recoil energy that allow the WIMP mass to be determined.

We calculate orbits for the Milky Way dwarf galaxies with proper motions, and compare these to subhalo orbits in a high resolution cosmological simulation. We use the simulation data to assess how well orbits may be recovered in the face of measurement errors, a time varying triaxial gravitational potential, and satellite-satellite interactions. For present measurement uncertainties, we recover the apocentre r_a and pericentre r_p to ~40%. With improved data from the Gaia satellite we should be able to recover r_a and r_p to ~14%, respectively. However, recovering the 3D positions and orbital phase of satellites over several orbits is more challenging. This owes primarily to the non-sphericity of the potential and satellite interactions during group infall. Dynamical friction, satellite mass loss and the mass evolution of the main halo play a more minor role in the uncertainties. We apply our technique to nine Milky Way dwarfs with observed proper motions. We show that their mean apocentre is lower than the mean of the most massive subhalos in our cosmological simulation, but consistent with the most massive subhalos that form before z=10. This lends further support to the idea that the Milky Way's dwarfs formed before reionisation.

We calculate orbits for the Milky Way dwarf galaxies with proper motions, and compare these to subhalo orbits in a high resolution cosmological simulation. We use this same simulation to assess how well are able to recover orbits in the face of measurement errors, a time varying triaxial gravitational potential, and satellite-satellite interactions. We find that, for present measurement uncertainties, we are able to recover the apocentre r_a and pericentre r_p to ~ 40%. However, even with better data the non-sphericity of the potential and satellite interactions during group infall make the orbital recovery more challenging. Dynamical friction, satellite mass loss and the mass evolution of the main halo play a more minor role. We apply our technique to nine Milky Way dwarfs with observed proper motions. We show that their mean apocentre is consistent with the most massive subhalos that form before z=10, lending support to the idea that the Milky Way dwarfs formed before reionisation.

We trace the evolution of central galaxies in three ~10^13 M_sun galaxy groups simulated at high resolution in cosmological hydrodynamical simulations. The evolution in the group potential leads, at z=0, to central galaxies that are massive, gas-poor early-type systems supported by stellar velocity dispersion resembling either elliptical or S0 galaxies. Their z~2-2.5 main progenitors are massive M* ~ 3-10 x 10^10 M_sun, star forming (20-60 M_sun/yr) galaxies which host substantial reservoirs of cold gas (~5 x 10^9 M_sun) in extended gas disks. Our simulations thus show that star forming galaxies observed at z~2 are likely the main progenitors of central galaxies in galaxy groups at z=0. Their central stellar densities stay approximately constant from z~1.5 down to z=0. Instead, the galaxies grow inside-out, by acquiring a stellar envelope outside the innermost ~2 kpc. Consequently the density within the effective radius decreases by up to two orders of magnitude. Both major and minor mergers contribute to most of the mass accreted outside the effective radius and thus drive the evolution of the half-mass radii. In one of the three simulated groups the short central cooling time leads to a dramatic rejuvenation of the central group galaxy at z<1, affecting its morphology, kinematics and colors. This episode is eventually terminated by a group-group merger. Our simulations demonstrate that, in galaxy groups, the interplay between halo mass assembly, galaxy merging and gas accretion has a substantial influence on the star formation histories and z=0 morphologies of central galaxies.[Abridged]

We place constraints on the formation redshifts for blue globular clusters (BGCs), independent of the details of hydrodynamics and population III star formation. The observed radial distribution of BGCs in the Milky Way Galaxy suggests that they formed in biased dark matter halos at high redshift. As a result, simulations of a ~1 Mpc box up to z~10 must resolve BGC formation in LCDM. We find that most halo stars could be produced from destroyed BGCs and other low-mass clusters that formed at high redshift. We present a proof-of-concept simulation that captures the formation of globular-like star clusters.

In LambdaCDM, massive satellites are dragged into the disk-plane by dynamical friction where they dissolve into a stellar thick disk and a more massive dark matter disk. The distinctive kinematics of the dark disk matches the stars that also entered in the satellites. The lower velocities of the dark disk with respect to the Earth enhances detection rates at low recoil energy. For WIMP masses > 50 GeV/c^2, the detection rate increases by up to a factor of 3 in the 5 - 20 keV recoil energy range. Comparing this with rates at higher energy will improve constraints on the WIMP mass, particularly for masses > 100 GeV/c^2. The annual modulation signal of the dark disk is significantly boosted and its modulation phase is shifted by ~3 weeks relative to the dark halo. The variation of the observed phase with recoil energy can also be used to determine the WIMP mass once the dark disk properties are fixed by future astronomical surveys. The constraints on the WIMP interaction cross section from current experiments improve by factors of 1.4 to 3.5 when a typical contribution from the dark disk is included.

As disc galaxies form in a hierarchical cosmology, massive merging satellites are preferentially dragged towards the disc plane. The material accreted from these satellites forms a dark matter disc that contributes 0.25 - 1.5 times the non-rotating halo density at the solar position. Here, we show the importance of the dark disc for indirect dark matter detection in neutrino telescopes. Previous predictions of the neutrino flux from WIMP annihilation in the Earth and the Sun have assumed that Galactic dark matter is spherically distributed with a Gaussian velocity distribution, the standard halo model. Although the dark disc has a local density comparable to the dark halo, its higher phase space density at low velocities greatly enhances capture rates in the Sun and Earth. For typical dark disc properties, the resulting muon flux from the Earth is increased by three orders of magnitude over the SHM, while for the Sun the increase is an order of magnitude. This significantly increases the sensitivity of neutrino telescopes to fix or constrain parameters in WIMP models. The flux from the Earth is extremely sensitive to the detailed properties of the dark disc, while the flux from the Sun is more robust. The enhancement of the muon flux from the dark disc puts the search for WIMP annihilation in the Earth on the same level as the Sun for WIMP masses < 100 GeV.

Cold dark matter (CDM) hierarchical structure formation models predict the existence of large-scale accretion shocks between the virial and turnaround radii of clusters of galaxies. Kocsis et al. (2005) suggest that the Sunyaev-Zel'dovich (SZ) signal associated with such shocks might be observable with the next generation radio interferometer, ALMA. We study the three--dimensional distribution of accretion shocks around individual clusters of galaxies drawn from adaptive mesh refinement (AMR) and smoothed particle hydrodynamics (SPH) simulations of LCDM (dark energy dominated CDM) models. In relaxed clusters, we find two distinct sets of shocks. One set ("virial shocks"), with Mach numbers of 2.5-4, is located at radii 0.9-1.3 Rvir, where Rvir is the spherical infall estimate of the virial radius, covering about 40-50% of the total surface area around clusters at these radii. Another set of stronger shocks ("external shocks") is located farther out, at about 3 Rvir, with large Mach numbers (~100), covering about 40-60% of the surface area. We simulate SZ surface brightness maps of relaxed massive galaxy clusters drawn from high resolution AMR runs, and conclude that ALMA should be capable of detecting the virial shocks in massive clusters of galaxies. More simulations are needed to improve estimates of astrophysical noise and to determine optimal observational strategies.

The fact that dark matter (DM), thus far, has revealed itself only on scales of galaxies and larger, again thrusts onto astrophysics the opportunity and the responsibility to confront the age old mystery "What is the nature of matter?" By deriving basic data on the nature of DM - e.g., mass of its particle(s), present mean temperature, distribution in galaxies and other structures in the universe, and capacity for dissipational collapse - we will be uncovering the properties of the dominant species of matter in the universe and significantly extending the standard models of particle physics. Determining the mass of the DM particle to an order of magnitude would help to sort out the particle family to which it (or they) belongs. Beyond mass, there are issues of stability. The DM particle may be unstable with a measurable half-life, or it may become unstable after absorbing a certain amount of energy from collisions. In both cases it would contribute to the present hot dark matter component. Some key parameters of DM can most accurately be measured in the very nearby universe because DM dominates the mass in the outer Milky Way (MW), in other galaxies in the Local Group, and in the Local Group in its entirety. The presence and amount of DM can be quantified by study of dynamical processes observable in fine detail within these entities. Precise measurements of 3-D velocities for stars, coherent star streams, and stars in satellite stellar systems out to the edge of the Galaxy can reveal "what is the shape, orientation, density law, and lumpiness of the dark matter halo" as well as "what is the total mass of the Galaxy?"

Predicting the local flux of dark matter particles is vital for dark matter direct detection experiments. To date, such predictions have been based on simulations that model the dark matter alone. Here we include the influence of the baryonic matter for the first time. We use two different approaches. Firstly, we use dark matter only simulations to estimate the expected merger history for a Milky Way mass galaxy, and then add a thin stellar disc to measure its effect. Secondly, we use three cosmological hydrodynamic simulations of Milky Way mass galaxies. In both cases, we find that a stellar/gas disc at high redshift (z~1) causes merging satellites to be preferentially dragged towards the disc plane. This results in an accreted dark matter disc that contributes ~0.25 - 1 times the non-rotating halo density at the solar position. An associated thick stellar disc forms with the dark disc and shares a similar velocity distribution. If these accreted stars can be separated from those that formed in situ, future astronomical surveys will be able to infer the properties of the dark disc from these stars. The dark disc, unlike dark matter streams, is an equilibrium structure that must exist in disc galaxies that form in a hierarchical cosmology. Its low rotation lag with respect to the Earth significantly boosts WIMP capture in the Earth and Sun, increases the likelihood of direct detection at low recoil energy, boosts the annual modulation signal, and leads to distinct variations in the flux as a function of recoil energy that allow the WIMP mass to be determined (see contribution from T. Bruch this volume).

Observations of turbulent velocity dispersions in the HI component of galactic disks show a characteristic floor in galaxies with low star formation rates and within individual galaxies the dispersion profiles decline with radius. We carry out several high resolution adaptive mesh simulations of gaseous disks embedded within dark matter haloes to explore the roles of cooling, star-formation, feedback, shearing motions and baryon fraction in driving turbulent motions. In all simulations the disk slowly cools until gravitational and thermal instabilities give rise to a multi-phase medium in which a large population of dense self-gravitating cold clouds are embedded within a warm gaseous phase that forms through shock heating. The diffuse gas is highly turbulent and is an outcome of large scale driving of global non-axisymmetric modes as well as cloud-cloud tidal interactions and merging. At low star-formation rates these processes alone can explain the observed HI velocity dispersion profiles and the characteristic value of ~10 km/s observed within a wide range of disk galaxies. Supernovae feedback creates a significant hot gaseous phase and is an important driver of turbulence in galaxies with a star-formation rate per unit area >10^-3 M_sun/yr/kpc^2.

In hierarchical structure formation models of disk galaxies, a dark matter disk forms as massive satellites are preferentially dragged into the disk-plane where they dissolve. Here, we quantify the importance of this dark disk for direct and indirect dark matter detection. The low velocity of the dark disk with respect to the Earth enhances detection rates in direct detection experiments at low recoil energy. For WIMP masses M_{WIMP} >~ 50 GeV, the detection rate increases by up to a factor of 3 in the 5 - 20 keV recoil energy range. Comparing this with rates at higher energy is sensitive to M_{WIMP}, providing stronger mass constraints particularly for M_{WIMP}>~100 GeV. The annual modulation signal is significantly boosted by the dark disk and the modulation phase is shifted by ~3 weeks relative to the dark halo. The variation of the observed phase with recoil energy determines M_{WIMP}, once the dark disk properties are fixed by future astronomical surveys. The low velocity of the particles in the dark disk with respect to the solar system significantly enhances the capture rate of WIMPs in the Sun, leading to an increased flux of neutrinos from the Sun which could be detected in current and future neutrino telescopes. The dark disk contribution to the muon flux from neutrino back conversion at the Earth is increased by a factor of ~5 compared to the SHM, for rho_d/rho_h=0.5.

We consider the possibility that the Magellanic Clouds were the largest members of a group of dwarf galaxies that entered the Milky Way (MW) halo at late times. Seven of the eleven brightest satellites of the MW may have been part of this system. The proximity of some dwarfs to the plane of the orbit of the Large Magellanic Cloud (LMC) has been used to argue that they formed from tidal debris from the LMC and Small Magellanic Cloud (SMC). Instead, they may owe to the tidal breakup of the Magellanic Group. This can explain the association of many of the dwarf galaxies in the Local Group with the LMC system. It provides a mechanism for lighting up dwarf galaxies and reproduces the bright end of the cumulative circular velocity distribution of the satellites in the MW without invoking a stripping scenario for the subhalos to match the satellite distribution expected according to CDM theory. Finally, our model predicts that other isolated dwarfs will be found to have companions. Evidence for this prediction is provided by nearby, recently discovered dwarf associations.

The Magellanic Clouds were the largest members of a group of dwarf galaxies that entered the Milky Way (MW) halo at late times. This group, dominated by the LMC, contained ~4% of the mass of the Milky Way prior to its accretion and tidal disruption, but ~70% of the known dwarfs orbiting the MW. Our theory addresses many outstanding problems in galaxy formation associated with dwarf galaxies. First, it can explain the planar orbital configuration populated by some dSphs in the MW. Second, it provides a mechanism for lighting up a subset of dwarf galaxies to reproduce the cumulative circular velocity distribution of the satellites in the MW. Finally, our model predicts that most dwarfs will be found in association with other dwarfs. The recent discovery of Leo V (Belokurov et al. 2008), a dwarf spheroidal companion of Leo IV, and the nearby dwarf associations supports our hypothesis.

In a LCDM cosmology, the Milky Way accretes satellites into the stellar disc. We use cosmological simulations to assess the frequency of near disc plane and higher inclination accretion events, and collisionless simulations of satellite mergers to quantify the final state of the accreted material and the effect on the thin disc. On average, a Milky Way-sized galaxy has 3 subhalos with vmax>80km/s; 7 with vmax>60km/s; and 15 with vmax>40km/s merge at redshift z>1. Assuming isotropic accretion, a third of these merge at an impact angle <20 degrees and are dragged into the disc plane by dynamical friction. Their accreted stars and dark matter settle into a thick disc. The stellar thick disc qualitatively reproduces the observed thick disc at the solar neighbourhood, but is less massive by a factor ~2-10. The dark matter disc contributes 0.25-1 times the halo density at the solar position. Although not likely to be dynamically interesting, the dark disc has important implications for the direct detection of dark matter because of its low velocity with respect to the Earth. Higher inclination encounters (>20 degrees) are twice as likely as low inclination ones. These lead to structures that closely resemble the recently discovered inner/outer stellar halos. They also do more damage to the Milky Way stellar disc creating a more pronounced flare, and warp; both long-lived and consistent with current observations. The most massive mergers (vmax>80km/s) heat the thin disc enough to produce a thick disc. These heated thin disc stars are essential for obtaining a thick disc as massive as that seen in the Milky Way; they likely comprise some ~50-90% of the thick disc stars. The Milky Way thin disc must reform from fresh gas after z=1 [abridged].

We construct mock galaxy catalogues to analyse clustering properties of a Lambda cold dark matter (LCDM) universe within a cosmological dark matter simulation of sufficient resolution to resolve structure down to the scale of dwarfs. We show that there is a strong age-clustering correlation for objects likely to host luminous galaxies, which includes the satellite halo (subhalo) population. Older mock galaxies are significantly more clustered in our catalog, which consists of satellite haloes as well as the central peaks of discrete haloes, selected solely by peak circular velocity. This age dependence is caused mainly by the age-clustering relation for discrete haloes, recently found by Gao et al., acting mostly on field members, combined with the tendency for older mock galaxies to lie within groups and clusters, where galaxy clustering is enhanced. Our results suggest that the clustering age dependence is manifested in real galaxies. At small scales (less than ~5 Mpc/h), the very simple assumption that galaxy colour depends solely on halo age is inconsistent with the strength of the observed clustering colour trends, where red galaxies become increasingly more clustered than blue galaxies toward smaller scales, suggesting that luminosity weighted galaxy ages do not closely trace the assembly epoch of their dark matter hosts. The age dependence is present but is much weaker for satellite haloes lying within groups and clusters than for the global population.

We use large cosmological N-body simulations to study the subhalo population in galaxy group sized halos. In particular, we look for fossil group candidates with typical masses ~10-25% of Virgo cluster but with an order of magnitude less substructure. We examine recent claims that the earliest systems to form are deficient enough in substructure to explain the luminosity function found in fossil groups. Although our simulations show a correlation between the halo formation time and the number of subhalos, the maximum suppression of subhalos is a factor of 2-2.5, whereas a factor of 6 is required to match fossil groups and galaxies. While the number of subhalos depends weakly on the formation time, the slope of the halo substructure velocity function does not. The satellite population within Cold Dark Matter (CDM) halos is self-similar at scales between galaxies and galaxy clusters regardless of mass, whereas current observations show a break in self-similarity at a mass scale corresponding to group of galaxies.

Using analytic calculations and N-body simulations we show that in constant density (harmonic) cores, sinking satellites undergo an initial phase of very rapid (super-Chandrasekhar) dynamical friction, after which they experience no dynamical friction at all. For density profiles with a central power law profile of log-slope, $-\alpha$, the infalling satellite heats the background and causes $\alpha$ to decrease. For $\alpha < 0.5$ initially, the satellite generates a small central constant density core and stalls as in the $\alpha = 0$ case. We discuss some astrophysical applications of our results to decaying satellite orbits, galactic bars and mergers of supermassive black hole binaries. In a companion paper we show that a central constant density core can provide a natural solution to the timing problem for Fornax's globular clusters.

Asteroids can be eclipsed by other bodies in the Solar System, but no direct observation of an asteroid eclipse has been reported to date. We describe a statistical method to predict an eclipse for an asteroid based on the analysis of the orbital elements covariance matrix. By propagating a set of Virtual Asteroids to an epoch correspondent to a close approach with a Solar System planet or natural satellite, it is possible to estimate the probability of a partial or total eclipse. The direct observation of an eclipse can provide data useful to improve the asteroid orbit, especially for dim asteroids typically observed only for a few days. We propose two different methods: the first, based on the inclusion of the apparent magnitude residuals into the orbit's least squares minimization process, capable of improving the asteroid's nominal orbit and the related covariance matrix; the second, based on weighting different Virtual Asteroids in relation to their apparent magnitude during the eclipse, useful for recovery purposes. As an application, we have numerically investigated the possibility of a Near Earth Asteroid eclipsed by the Moon or the Earth in the 1990-2050 period. A total of 74 distinct eclipses have been found, involving 59 asteroids. In particular, the asteroid (99942) Apophis has a probability of about 74% to enter the Moon's penumbra cone and a probability of about 6% to enter the umbra cone on April 14, 2029, less than six hours after a very close approach to Earth.

We use cosmological LCDM numerical simulations to model the evolution of the substructure population in sixteen dark matter haloes with resolutions of up to seven million particles within the virial radius. The combined substructure circular velocity distribution function (VDF) for hosts of 10^11 to 10^14 Msun at redshifts from zero to two or higher has a self-similar shape, is independent of host halo mass and redshift, and follows the relation: dn/dv=(1/8)(v_cmax/v_cmax,host)^-4. Halo to halo variance in the VDF is a factor of roughly two to four. At high redshifts, we find preliminary evidence for fewer large substructure haloes (subhaloes). Specific angular momenta are significantly lower for subhaloes nearer the host halo centre where tidal stripping is more effective. The radial distribution of subhaloes is marginally consistent with the mass profile for r >~ 0.3r_vir, where the possibility of artificial numerical disruption of subhaloes can be most reliably excluded by our convergence study, although a subhalo distribution that is shallower than the mass profile is favoured. Subhalo masses but not circular velocities decrease toward the host centre. Subhalo velocity dispersions hint at a positive velocity bias at small radii. There is a weak bias toward more circular orbits at lower redshift, especially at small radii. We additionally model a cluster in several power law cosmologies of P ~ k^n, and demonstrate that a steeper spectral index, n, results in significantly less substructure.

We use numerical simulations in a Lambda CDM cosmology to model density profiles in a set of 16 dark matter haloes with resolutions of up to 7 million particles within the virial radius. These simulations allow us to follow robustly the formation and evolution of the central cusp over a large mass range of 10^11 to 10^14 M_sun, down to approximately 0.5% of the virial radius, and from redshift 5 to the present. The cusp of the density profile is set at redshifts of 2 or greater and remains remarkably stable to the present time, when considered in non-comoving coordinates. We fit our haloes to a 2 parameter profile where the steepness of the asymptotic cusp is given by gamma, and its radial extent is described by the concentration, c_gamma. In our simulations, we find gamma = 1.4 - 0.08Log(M/M_*) for haloes of 0.01M_* to 1000M_*, with a large scatter of gamma ~ +/-0.3$; and c_gamma = 8*M/M_*^{-0.15}, with a large M/M_* dependent scatter roughly equal to +/- c_gamma. Our redshift zero haloes have inner slope parameters ranging approximately from r^{-1} to r^{-1.5}, with a median of roughly r^{-1.3}. This 2 parameter profile fit works well for all our halo types, whether or not they show evidence of a steep asymptotic cusp. We also model a cluster in power law cosmologies of P ~ k^n (n=0,-1,-2,-2.7). We find larger concentration radii and shallower cusps for steeper n. The minimum resolved radius is well described by the mean interparticle separation. The trend of steeper and more concentrated cusps for smaller $M/M_*$ haloes clearly shows that dwarf sized Lambda CDM haloes have, on average, significantly steeper density profiles within the inner few percent of the virial radius than inferred from recent observations. Code to reproduce this profile can be downloaded from http://www.icc.dur.ac.uk/~reed/profile.html

The Cold Dark Matter (CDM) theory predicts a wealth of substructure within dark halos. These predictions match observations of galaxy clusters like the nearby Virgo cluster. However, CDM has a "small scale crisis" since galaxies dominate the halo with little substructure while the model predicts that galaxies should be scaled versions of galaxy clusters with abundant substructure. Compared to CDM predictions, the Milky Way and Andromeda are "missing" objects with velocity dispersions $\sigma \geq 10$ km s$^{-1}$. The energy scale of these missing satellites is low enough that stellar winds and supernovae might remove gas and suppress the formation of their luminous stellar components. Here, we show that the small scale crisis persists in fossil groups that have masses of up to 40% of the nearby Virgo cluster of galaxies. Fossil groups are missing satellites with luminosities that occur at the predicted frequency in the Virgo cluster. Moreover, the "missing galaxies" in fossil groups are nearly as luminous as the Milky Way with a velocity dispersion $\sigma \leq$ 150 km s$^{-1}$.

We simulate the formation of a realistic disk galaxy within the hierarchical scenario of structure formation and study its internal properties to the present epoch. We compare results from a LambdaCDM simulation with a LambdaWDM (2keV) simulation that forms significantly less small scale structure. We show how high mass and force resolution in both the gas and dark matter components play an important role in solving the angular momentum catastrophe claimed from previous simulations of galaxy formation within the hierarchical framework. The stellar material in the disk component has a final specific angular momentum equal to 40% and 90% of that of the dark halo in the LambdaCDM and LambdaWDM models respectively. The LambdaWDM galaxy has a drastically reduced satellite population and a negligible stellar spheroidal component. Encounters with satellites play only a minor role in disturbing the disk. Satellites possess a variety of star formation histories linked to mergers and pericentric passages along their orbit around the primary galaxy. In both cosmologies, the galactic halo retains most of the baryons accreted and builds up a hot gas phase with a substantial X-ray emission. Therefore, while we have been successful in creating a realistic stellar disk in a massive galaxy within the LambdaCDM scenario, energy injection emerges as necessary ingredient to reduce the baryon fraction in galactic halos, independent of the cosmology adopted. (abridged)

We use a high resolution $\Lambda$CDM numerical simulation to calculate the mass function of dark matter haloes down to the scale of dwarf galaxies, back to a redshift of fifteen, in a 50 $h^{-1}$Mpc volume containing 80 million particles. Our low redshift results allow us to probe low $\sigma$ density fluctuations significantly beyond the range of previous cosmological simulations. The Sheth and Tormen mass function provides an excellent match to all of our data except for redshifts of ten and higher, where it overpredicts halo numbers increasingly with redshift, reaching roughly 50 percent for the $10^{10}-10^{11} \msun$ haloes sampled at redshift 15. Our results confirm previous findings that the simulated halo mass function can be described solely by the variance of the mass distribution, and thus has no explicit redshift dependence. We provide an empirical fit to our data that corrects for the overprediction of extremely rare objects by the Sheth and Tormen mass function. This overprediction has implications for studies that use the number densities of similarly rare objects as cosmological probes. For example, the number density of high redshift (z $\simeq$ 6) QSOs, which are thought to be hosted by haloes at 5$\sigma$ peaks in the fluctuation field, are likely to be overpredicted by at least a factor of 50%. We test the sensitivity of our results to force accuracy, starting redshift, and halo finding algorithm.

We present a set of gas-dynamical simulations of galaxy groups and clusters aimed at exploring the effect of non-gravitational heating. We use GASOLINE, a parallel Tree+SPH code, to simulate the formation of four cosmic halos with temperature 0.5<T<8 keV. Non-gravitational heating is implemented in two different ways: (1) by imposing a minimum entropy floor at a given redshift, 1<z<5; (2) by gradually heating gas, proportionally to the SN rate expected from semi-analytical modeling of galaxy formation. Our main results are the following. (a) An extra heating energy of about 1 keV per gas particle is required to reproduce the observed Lx-T relation, independent of whether it is provided so as to create an entropy floor of 50-100 keV cm^2, or is modulated in redshift; our SN feedback recipe provides only 1/3 keV/part. (b) The M-T relation is almost unaffected by non-gravitational heating and follows the M T^{3/2} scaling, with a normalization ~40% higher than observed, independent of the heating scheme. The inclusion of cooling in a run of a small group has the effects of increasing T_ew by ~30%, possibly reconciling simulated and observed M-T relations, and of decreasing Lx by ~40%. In spite of the inclusion of SN feedback energy, almost 40% of the gas becomes cold, in excess of current observational estimates. (abridged)

We present results from high-resolution Tree+SPH simulations of galaxy clusters and groups, aimed at studying the effect of non-gravitational heating on the entropy of the ICM. We simulate three systems, having emission-weighted temperature T=0.6,1 and 3 keV, with spatial resolution better than 1% of the virial radius. We consider the effect of different prescriptions for non-gravitational ICM heating, such as SN energy feedback, as predicted by semi-analytical models of galaxy formation, and two different minimum entropy floors, S_fl=50 and 100 keV cm^2, imposed at z=3. Simulations with only gravitational heating nicely reproduce predictions from self-similar ICM models, while extra heating is shown to break the self-similarity, by a degree which depends on total injected energy and on cluster mass. We use observational results on the excess entropy in central regions of galaxy systems, to constrain the amount of extra-heating required. We find that setting the entropy floor S_fl=50 keV cm^2, which corresponds to an extra heating energy of about 1 keV per particle, is able to reproduce the observed excess of ICM entropy.

We use a cosmological simulation of the Local Group to make quantitative and speculative predictions for direct detection experiments. Cold dark matter (CDM) halos form via a complex series of mergers, accretion events and violent relaxation which precludes the formation of significant caustic features predicted by axially symmetric collapse. The halo density profiles are combined with observational constraints on the galactic mass distribution to constrain the local density of cold dark matter to lie in the range 0.18 <~ rho_CDM(R_solar)/GeV cm^-3 <~ 0.30. In velocity space, coherent streams of dark matter from tidally disrupted halos fill the halo and provide a tracer of the merging hierarchy. The particle velocities within triaxial CDM halos cannot be approximated by a simple Maxwellian distribution and is radially biased at the solar position. The detailed phase space structure within the solar system will depend on the early merger history of the progenitor halos and the importance of major mergers over accretion dominated growth. We follow the formation of a ``Draco'' sized dSph halo of mass 10^8M_solar with several million particles and high force accuracy. Its internal structure and substructure resembles that of galactic or cluster mass halos: the density profile has a singular central cusp and it contains thousands of sub-halos orbiting within its virial radius demonstrating a self-similar nature to collisionless dark matter sub-clustering. The singular cores of substructure halos always survive complete tidal disruption although mass loss is continuous and rapid. Extrapolating wildly to earth mass halos with velocity dispersion of 1 m s^-1 (roughly equal to the free streaming scale for neutralinos) we find that most of the dark matter may remain attached to bound subhalos. (Abridged)

We present results from high-resolution N-Body/SPH simulations of rotationally supported dwarf irregular galaxies moving on bound orbits in the massive dark matter halo of the Milky Way.The dwarf models span a range in disk surface density and the masses and sizes of their dark halos are consistent with the predictions of cold dark matter cosmogonies. We show that the strong tidal field of the Milky Way determines severe mass loss in their halos and disks and induces bar and bending instabilities that transform low surface brightness dwarfs (LSBs) into dwarf spheroidals (dSphs) and high surface brightness dwarfs (HSBs) into dwarf ellipticals (dEs) in less than 10 Gyr. The final central velocity dispersions of the remnants are in the range 8-30 km/s and their final $v/\sigma$ falls to values $< 0.5$, matching well the kinematics of early-type dwarfs. The transformation requires the orbital time of the dwarf to be $\simlt 3-4$ Gyr, which implies a halo as massive and extended as predicted by hierarchical models of galaxy formation to explain the origin of even the farthest dSph satellites of the Milky Way, Leo I and Leo II. Only dwarfs with central dark matter densities as high as those of Draco and Ursa Minor can survive for 10 Gyr in the proximity of the Milky Way: this is naturally achieved within hierarchical models, where the densest objects should have small orbital times due to their early formation epochs. Part of the gas is stripped and part is funneled to the center due to the bar, generating one strong burst of star formation in HSBs and smaller, multiple bursts in LSBs. Extended low-surface brightness stellar and gaseous streams originate from LSBs and, when projected along the line of sight, can lead to overestimate the mass-to-light ratio of the bound remnant by a factor $\simlt 2$,

N-Body/SPH simulations are used to study the evolution of dwarf irregular galaxies (dIrrs) entering the dark matter halo of the Milky Way or M31 on plunging orbits. We propose a new dynamical mechanism driving the evolution of gas rich, rotationally supported dIrrs, mostly found at the outskirts of the Local Group (LG), into gas free, pressure supported dwarf spheroidals (dSphs) or dwarf ellipticals (dEs), observed to cluster around the two giant spirals. The initial model galaxies are exponential disks embedded in massive dark matter halos and reproduce nearby dIrrs. Repeated tidal shocks at the pericenter of their orbit partially strip their halo and disk and trigger dynamical instabilities that dramatically reshape their stellar component. After only 2-3 orbits low surface brightness (LSB) dIrrs are transformed into dSphs, while high surface brightness (HSB) dIrrs evolve into dEs. This evolutionary mechanism naturally leads to the morphology-density relation observed for LG dwarfs. Dwarfs surrounded by very dense dark matter halos, like the archetypical dIrr GR8, are turned into Draco or Ursa Minor, the faintest and most dark matter dominated among LG dSphs. If disks include a gaseous component, this is both tidally stripped and consumed in periodic bursts of star formation. The resulting star formation histories are in good qualitative agreement with those derived using HST color-magnitude diagrams for local dSphs.

Many models of globular cluster formation assume the presence of cold dense clouds in early universe. Here we re-examine the Fall & Rees (1985) model for formation of proto-globular cluster clouds (PGCCs) via thermal instabilities in a protogalactic halo. We first argue, based on the previous study by others, that under the protogalactic environments only nonlinear density inhomogeneities can condense into PGCCs. We then carry out numerical simulations of the collapse of overdense clouds in one-dimensional spherical geometry, including self-gravity and radiative cooling down to T=10^4 K. Since imprinting of Jeans mass at 10^4 K is essential to this model, here we focus on the cases where external UV background radiation prevents the formation of H2 molecules and so prevent the cloud from cooling below 10^4 K. The quantitative results from these simulations can be summarized as follows: 1) Perturbations smaller than M_min ~ (10^{5.6} M_sun) (n_h/0.05 cm3)^{-2} cool isobarically, while perturbations larger than M_max ~ (10^8 M_sun) (n_h/0.05 cm3)^{-2} cool isochorically. On the other hand, intermediate size perturbations (M_min< M_pgcc < M_max) are compressed supersonically. 2) For supersonically collapsing clouds, the density compression factor after they cool to T_c=10^4 K range 10^{2.5}-10^6. 3) For supersonically collapsing clouds the Jeans mass can be reduced to as small as 10^{5.5} M_sun (n_h/0.05 cm3)^{-1/2} at the maximum compression. 4) The density profile of simulated PGCCs can be approximated by a constant core with a halo of rho ~ r^{-2} rather than a singular isothermal sphere.

Can N-body simulations reliably determine the structural properties of dark matter halos? Focussing on a Virgo-sized galaxy cluster, we increase the resolution of current ``high resolution simulations'' by almost an order of magnitude to examine the convergence of the important physical quantities. We have 4 million particles within the cluster and force resolution 0.5 kpc/h (0.05% of the virial radius). The central density profile has a logarithmic slope of -1.5, as found in lower resolution studies of the same halo, indicating that the profile has converged to the ``physical'' limit down to scales of a few kpc. Also the abundance of substructure is consistent with that derived from lower resolution runs; on the scales explored, the mass and circular velocity functions are close to power laws of exponents ~ -1.9 and -4. Overmerging appears to be globally unimportant for suhalos with circular velocities > 100 km/s. We can trace most of the cluster progenitors from z=3 to the present; the central object (the dark matter analog of a cD galaxy)is assembled between z=3 and 1 from the merging of a dozen halos with v_circ \sim 300 km/s. The mean circular velocity of the subhalos decreases by ~ 20% over 5 billion years, due to tidal mass loss. The velocity dispersions of halos and dark matter globally agree within 10%, but the halos are spatially anti-biased, and, in the very central region of the cluster, they show positive velocity bias; however, this effect appears to depend on numerical resolution.

We combine semi-analytical methods with a ultra-high resolution simulation of a galaxy cluster (of mass 2.3 10^14h-1Msolar, and 4 10^6 particles within its virial radius) formed in a standard CDM universe to study the spatial distribution and orbital properties of the present-day descendents of Lyman Break Galaxies (LBGs). At the present time only five (out of 12) of halos containing LBGs survive as separate entities inside the cluster virial radius. Their circular velocities are in the range 200 - 550 km/sec. Seven halos merged together to form the central object at the very center of the cluster. Using semi-analytical modeling of galaxy evolution we show that descendents of halos containing LBGs now host giant elliptical galaxies. Galaxy orbits are radial, with a pericenter to apocenter ratio of about 1:5. The orbital eccentricities of LBGs descendents are statistically indistinguishable from those of the average galaxy population inside the cluster, suggesting that the orbits of these galaxies are not significantly affected by dynamical friction decay after the formation of the cluster's main body. In this cluster, possibly due to its early formation time, the descendents of LBGs are contained within the central 60% of the cluster virial radius and have an orbital velocity dispersion lower than the global galaxy population, originating a mild luminosity segregation for the brightest cluster members. Mass estimates based only on LBGs descendents (especially including the central cD) reflect this bias in space and velocity and underestimate the total mass of this well virialized cluster by up to a factor of two compared to estimates using at least 20 cluster members.

We use numerical simulations to examine the substructure within galactic and cluster mass halos that form within a hierarchical universe. Clusters are easily reproduced with a steep mass spectrum of thousands of substructure clumps that closely matches observations. However, the survival of dark matter substructure also occurs on galactic scales, leading to the remarkable result that galaxy halos appear as scaled versions of galaxy clusters. The model predicts that the virialised extent of the Milky Way's halo should contain about 500 satellites with circular velocities larger than Draco and Ursa-Minor i.e. bound masses > 10^8Mo and tidally limited sizes > kpc. The substructure clumps are on orbits that take a large fraction of them through the stellar disk leading to significant resonant and impulsive heating. Their abundance and singular density profiles has important implications for the existence of old thin disks, cold stellar streams, gravitational lensing and indirect/direct detection experiments.

We show that a universe dominated by cold dark matter fails to reproduce the rotation curves of dark matter dominated galaxies, one of the key problems that it was designed to resolve. We perform numerical simulations of the formation of dark matter halos, each containing \gsim 10^6 particles and resolved to 0.003 times the virial radius, allowing an accurate comparison with rotation curve data. A good fit to both galactic and cluster sized halos can be achieved using the density profile rho(r) \propto [(r/r_s)^1.5(1+(r/r_s)^1.5)]^-1, where r_s is a scale radius. This profile has a steeper asymptotic slope, rho(r) \propto r^-1.5, and a sharper turnover than found by lower resolution studies. The central structure of relaxed halos that form within a hierarchical universe has a remarkably small scatter (unrelaxed halos would not host disks). We compare the results with a sample of dark matter dominated, low surface brightness (LSB) galaxies with circular velocities in the range 100-300 km/s. The rotation curves of disks within cold dark matter halos rise too steeply to match these data which require a constant mass density in the central regions. The same conclusion is reached if we compare the scale free shape of observed rotation curves with the simulation data. It is important to confirm these results using stellar rather than HI rotation curves for LSB galaxies. We test the effects of introducing a cut-off in the power spectrum that may occur in a universe dominated by warm dark matter. In this case halos form by a monolithic collapse but the final density profile hardly changes, demonstrating that the merger history does not play a role in determining the halo structure.

We follow the evolution of disk galaxies within a cluster that forms hierarchically in a standard cold dark matter N-body simulation. At a redshift z=0.5 we select several dark matter halos that have quiet merger histories and are about to enter the newly forming cluster environment. The halos are replaced with equilibrium high resolution model spirals that are constructed to represent luminous examples of low surface brightness (LSB) and high surface brightness (HSB) galaxies. Whilst the models have the same total luminosity, \sim L_*, they have very different internal mass profiles, core radii and disk scale lengths, however they all lie at the same place on the Tully-Fisher relation. Due to their ``soft'' central potentials, LSB galaxies evolve dramatically under the influence of rapid encounters with dark matter halos, substructure and strong tidal shocks from the global cluster potential - galaxy harassment. As much as 90% of the LSB disk stars are tidally stripped and congregate in large diffuse tails that trace the orbital path of the galaxy and form the diffuse intra-cluster light. The bound stellar remnants closely resemble the dwarf spheroidals (dE's) that populate nearby clusters, with large scale lengths and low central surface brightness.

We examine the distributions of eccentricities of orbits within mass distributions like those we see for galaxies and clusters. A comprehensive understanding of these orbital properties is essential to calculate the rates of physical processes relevant to the formation and evolution of galaxies and clusters. We derive the orbital eccentricity distributions for a number of spherical potentials. These distributions depend strongly on the velocity anisotropy, but only slightly on the shape of the potential. We also present high resolution N-body simulations of the orbital decay of satellite systems on eccentric orbits in an isothermal halo. The dynamical friction timescales are found to decrease with increasing orbital eccentricity, but the dependence is weaker than previously suggested. The orbital eccentricity stays remarkably constant throughout the decay; although the eccentricity decreases near pericenter, it increases again near apocenter, such that there is no net circularization. Finally, we discuss a number of applications. We find that the distribution of orbital eccentricities of the globular clusters in the Milky Way is consistent with a close to isotropic velocity distribution. The limited data on the galactic system of satellites, however, appears to be different and warrants further investigation as a clue to the formation and evolution of our Milky Way and its halo substructure.

We follow the evolution of disk galaxies within a cluster that forms hierarchically in a cold dark matter N-body simulation. At a redshift z=0.5 we select several dark matter halos that have quiet merger histories and are about to enter the newly forming cluster environment. The halos are replaced with equilibrium high resolution model spirals that are constructed to represent examples of low surface brightness (LSB) and high surface brightness (HSB) galaxies. Varying the disk and halo structural parameters reveals that the response of a spiral galaxy to tidal encounters depends primarily on the potential depth of the mass distribution and the disk scale length. LSB galaxies, characterised by slowly rising rotation curves and large scale lengths, evolve dramatically under the influence of rapid encounters with substructure and strong tidal shocks from the global cluster potential --- galaxy harassment. We find that up to 90% of their stars are tidally stripped and congregate in large diffuse tails that trace the orbital path of the galaxy and form the diffuse intra-cluster light. The bound stellar remnants closely resemble the dwarf spheroidals (dE's) that populate nearby clusters. HSB galaxies are stable to the chaos of cluster formation and tidal encounters. These disks lie well within the tidally limited dark matter halos and their potentials are more concentrated. Although very few stars are stripped, the scale height of the disks increases substantially and no spiral features remain, therefore we speculate that these galaxies would be identified as S0 galaxies in present day clusters.

We analyse parallel N-body simulations of three Cold Dark Matter (CDM) universes to study the abundance and clustering of galaxy clusters. The simulations cover a volume comparable to the forthcoming SDSS. We are able to make robust measurements of cluster properties to a redshift larger than unity. We extract halos using two independent, public domain group finders (FOF & HOP) and find consistent results. The correlation function of clusters is in very good agreement with a simple analytic prescription based upon a Lagrangian biasing scheme developed by Mo & White (1996) and the Press-Schechter (PS) formalism for the mass function. The R_0--D_c relation for the open CDM model is in good agreement with the results from the APM Cluster Survey. The SCDM universe shows a robust deviation in the shape and evolution of the mass function when compared with that predicted by the PS formalism. Critical models with a low sigma_8 normalization or small shape parameter Gamma show an excess of massive clusters compared with the PS prediction. When cluster normalized, the SCDM universe at z =1 contains 10 times more clusters with temperatures greater than 7keV, compared with the PS prediction. The agreement between the analytic and N-body mass functions of SCDM can be improved if the value of the delta_c (the extrapolated linear theory threshold for collapse) is revised to be $ delta_c(z) = 1.685[(0.7/sigma_8)(1+z)]^{-0.125}. Our best estimate for the amplitude of fluctuations inferred from the local cluster abundance for SCDM is sigma_{8} = 0.5 \pm 0.04. However, the discrepancy between the temperature function predicted in a critical density universe and that observed at z=0.33 (Henry et al. 1998) remains. (abridged)

Higher order statistics are investigated in ($\Omega$)CDM universes by analyzing $500\mpc$ high resolution tree N-body simulations with both $\Omega = 1$, and $\Omega < 1$. The amplitudes of the N-point correlation functions are calculated from moments of counts-in-cells determined by a pair of new algorithms especially developed for large simulations. This approach enables massive oversampling with $\simeq 10^{9-14}$ cells for accurate determination of factorial moments from up to 47 million particles in the scale range of $8 \kpc - 125\mpc$. Thorough investigation shows that there are three scale ranges in the simulations: $\ge 8\mpc$, weakly non-linear regime, where perturbation theory applies with utmost precision, $1\mpc - 8\mpc$, the non-linear plateau, and finally $\le 1\mpc$, a regime where dynamical discreteness effects dominate the higher order statistics. In the physically relevant range of $1\mpc-125\mpc$ the results i) confirm the validity of perturbation theory in the weakly non-linear regime, ii) establish the existence of a plateau in the highly non-linear regime similar to the one observed in scale free simulations iii) show extended perturbation theory to be an excellent approximation for the non-linear regime iv) find the time dependence of the $S_N$'s to be negligible in both regimes v) in comparison with similar measurements in the EDSGC survey, strongly support $\Omega < 1$ with no biasing vi) show that the formulae of Szapudi & Colombi (1996) provide a good approximation for errors on higher order statistics measured in N-body simulations.

We examine the properties of dark matter halos within a rich galaxy cluster using a high resolution simulation that captures the cosmological context of a cold dark matter universe. The mass and force resolution permit the resolution of 150 halos with circular velocities larger than 80 kms within the cluster's virial radius of 2 Mpc. This enables an unprecedented study of the statistical properties of a large sample of dark matter halos evolving in a dense environment. The cumulative fraction of mass attached to these halos varies from 0% at 200 kpc, to 13% at the virial radius. Even at this resolution the overmerging problem persists; halos that pass within 200 kpc of the cluster center are tidally disrupted. Additional substructure is lost at earlier epochs within the massive progenitor halos. The median ratio of apocentric to pericentric radii is 6:1; the orbital distribution is close to isotropic, circular orbits are rare, radial orbits are common. The orbits of halos are unbiased with respect to both position within the cluster and with the orbits of the smooth dark matter background and no velocity bias is detected. The tidal radii of surviving halos are generally well-fit using the simple analytic prediction applied to their orbital pericenters. Halos within clusters have higher concentrations than those in the field. Within the cluster, halo density profiles can be modified by tidal forces and individual encounters with other halos that cause significant mass loss - ``galaxy harassment''. Mergers between halos do not occur inside the clusters virial radius.

We examine the effects of mass resolution and force softening on the density profiles of cold dark matter halos that form within cosmological N-body simulations. As we increase the mass and force resolution, we resolve progenitor halos that collapse at higher redshifts and have very high densities. At our highest resolution we have nearly 3 million particles within the virial radius, several orders of magnitude more than previously used and we can resolve more than one thousand surviving dark matter halos within this single virialised system. The halo profiles become steeper in the central regions and we may not have achieved convergence to a unique slope within the inner 10% of the virialised region. Results from two very high resolution halo simulations yield steep inner density profiles, $\rho(r)\sim r^{-1.4}$. The abundance and properties of arcs formed within this potential will be different from calculations based on lower resolution simulations. The kinematics of disks within such a steep potential may prove problematic for the CDM model when compared with the observed properties of halos on galactic scales.

We investigate the internal structure of cold dark matter halos using high resolution N-body simulations. As the mass and force resolution are increased, halo density profiles become steeper, asymptoting to a slope of $\sim -4/3$ in the central regions and may not have converged to a unique result. At our highest resolution we have nearly 3 million particles within the virial radius, $R_{200}$, and force softening that is $\sim 0.2% R_{200}$. This resolution has also allowed us to resolve a large part of the overmerging problem - we find over 1000 surviving dark halos orbiting within a single cluster potential. These data have given us unprecedented insights into the dynamics and structure of ``halos within halos'', allowing us for the first time, to compare the distribution of dark matter with observations of galaxies in clusters.

We have simulated disk galaxies undergoing continual bombardment by other galaxies in a rich cluster. "Galaxy harassment" leads to dramatic evolution of smaller disk galaxies and provides an extremely effective mechanism to fuel a central quasar. Within a few billion years after a small disk galaxy enters the cluster environment, up to 90% of its gas can be driven into the inner 500 pc. Up to half of the mass can be transferred in a burst lasting just 100-200 Myr. This transport of gas to the center of galaxy is far more efficient than any mechanism proposed before. Galaxy harassment was first proposed to explain the disturbed blue galaxies in clusters seen in clusters at ($z \gsim 0.3$), the "Butcher-Oemler effect". Quasars at the same reshifts lie in more clustered environments than those at lower redshift. Recent HST observations find that roughly half of all observed quasar host galaxiess are fainter than \l*, with many of these less luminous hosts occuring at redshifts $z \gsim 0.3$. We examine 5 quasars that are claimed to have low luminosity hosts and find that 3 are in rich clusters of galaxies, the fourth may be in a cluster but the evidence for this is marginal. The environment of the fifth has not been studied.

Galaxy morphologies in clusters have undergone a remarkable transition over the past several billion years. Distant clusters at $z \sim 0.4$ are filled with small spiral galaxies, many of which are disturbed and show evidence of multiple bursts of star-formation. This population is absent from nearby clusters where spheroidals comprise the faint end of the luminosity function. Our numerical simulations follow the evolution of disk galaxies in a rich cluster owing to encounters with brighter galaxies and the cluster's tidal field---galaxy harassment. After a bursting transient phase, they undergo a complete morphological transformation from "disks" to"spheroidals". We examine the remnants and find support for our theory in detailed comparisons of the photometry and kinematics of the spheroidal galaxies in clusters. Our model naturally accounts for the intermediate age stellar population seen in these spheroidals as well as the trend in dwarf to giant ratio with cluster richness. The final shapes are typically prolate and are flattened primarily by velocity anisotropy. Their mass to light ratios are in the range 3---8 in good agreement with observations.

We present the results from a two-day study in which we discussed various implementations of Smooth Particle Hydrodynamics (SPH), one of the leading methods used across a variety of areas of large-scale astrophysical simulations. In particular, we evaluated the suitability of designing special hardware extensions, to further boost the performance of the high-end general purpose computers currently used for those simulations. We considered a range of hybrid architectures, consisting of a mix of custom LSI and reconfigurable logic, combining the extremely high throughput of Special-Purpose Devices (SPDs) with the flexibility of reconfigurable structures, based on Field Programmable Gate Arrays (FPGAs). The main findings of our workshop consist of a clarification of the decomposition of the computational requirements, together with specific estimates for cost/performance improvements that can be obtained at each stage in this decomposition, by using enabling hardware technology to accelerate the performance of general purpose computers.

Using cosmological N-body simulations of critical (SCDM) and open (Omega_0=0.3, OCDM) cold dark matter models we have identified dark matter halos which are associated with the progenitors of present day bright early type galaxies. By following their merging history, we show how early type galaxies that formed within massive halos at redshift around 3 are now preferentially residing in clusters and groups. On the other hand, those that formed through major merging events between redshift 1 and the present have not yet been accreted into larger, virialized structures. This result is in agreement with analytical predictions in hierarchical clustering models. CDM models are able to explain both the ancient and uniform population of ellipticals that dominates in clusters together with the more recent and heterogeneous population of field ellipticals. Predictions for the comoving number density of bright early-type galaxies are given, and are shown to be consistent with the observed luminosity function. We predict that the number density of interacting bright binary galaxies, from which the field population of ellipticals may have originated, is proportional to (1+z)**4.2+-0.28 and (1+z)**2.5+-0.42 in SCDM and OCDM respectively. This result is consistent with previous analytical estimates and is discussed together with recent observational constraints.

Leapfrog integration has been the method of choice in N-body simulations owing to its low computational cost for a symplectic integrator with second order accuracy. We introduce a new leapfrog integrator that allows for variable timesteps for each particle in large N-body simulations. Tests with single particles in fixed potentials show that it behaves as a symplectic integrator. We then examine the results of both standard leapfrog and our temporally adaptive leapfrog on full N-body integrations of clusters and large scale structure establishing accuracy criteria for both methods. The adaptive method shows significant speed-ups over single step integrations---but the integrator no longer appears to be symplectic or, in the case of large scale structure simulations, accurate. This loss of accuracy appears to be caused by the way that the timestep is chosen, not by the integrator itself. We present a related integration technique that does retain sufficient accuracy. Although it is not symplectic, it is apparently better than previous implementations and is our current integrator of choice for large astrophysical simulations. We also note that the standard leapfrog difference equations used in cosmological N-body integrations in comoving coordinates are not symplectic. We derive an implementation of leapfrog that is in comoving canonical coordinates to correct for this deficiency.

The dynamics of the Local Group and its environment provide a unique challenge to cosmological models. The velocity field within 5h-1 Mpc of the Local Group (LG) is extremely ``cold''. The deviation from a pure Hubble flow, characterized by the observed radial peculiar velocity dispersion, is measured to be about 60km/s. We compare the local velocity field with similarly defined regions extracted from N-body simulations of Universes dominated by cold dark matter (CDM). This test is able to strongly discriminate between models that have different mean mass densities. We find that neither the Omega=1 (SCDM) nor Omega=0.3 (OCDM) cold dark matter models can produce a single candidate Local Group that is embedded in a region with such small peculiar velocities. For these models, we measure velocity dispersions between 500-700km/s and 150-300km/s respectively, more than twice the observed value. Although both CDM models fail to produce environments similar to those of our Local Group on a scale of a few Mpc, they can give rise to many binary systems that have similar orbital properties as the Milky Way--Andromeda system. The local, gravitationally induced bias of halos in the CDM ``Local Group'' environment, if defined within a sphere of 10 Mpc around each Local Group is about 1.5, independent of Omega. No biasing scheme could reconcile the measured velocity dispersions around Local Groups with the observed one. Identification of binary systems using a halo finder (named Skid (http://www-hpcc.astro.washington.edu/tools/DENMAX for a public version)) based on local density maxima instead of a simple linking algorithm, gives a much more complete sample. We show that a standard ``friend of friends'' algorithm would miss 40% of the LG candidates present in the simulations.

Disturbed spiral galaxies with high rates of star formation pervaded clusters of galaxies just a few billion years ago, but nearby clusters exclude spirals in favor of ellipticals. ``Galaxy harassment" (frequent high speed galaxy encounters) drives the morphological transformation of galaxies in clusters, provides fuel for quasars in subluminous hosts and leaves detectable debris arcs. Simulated images of harassed galaxies are strikingly similar to the distorted spirals in clusters at $z \sim 0.4$ observed by the Hubble Space Telescope.

N-body simulations that follow only a collisionless dark matter component have failed to produce galaxy halos or substructure within dense environments. We investigate the `over-merging' problem analytically and with numerical simulations, by calculating dissolution timescales of halos due to physical and artificial dynamical effects. The numerical resolution that has recently been attained is such that mass-loss from two-body relaxation is negligible. We demonstrate that substructure is destroyed in present simulations as a result of large force softening combined with the heating sources of tides and encounters with dissolving substructure. In the limit of infinite numerical resolution, whether or not individual halos or substructure can survive depends sensitively on their inner density profiles. Singular isothermal halos will always survive at some level, however, if halos form with large core radii then the over-merging problem will always exist within dissipationless N-body simulations. In this latter case a dissipational component can increase the halos central density enabling galaxies to survive.

If the dark matter in galactic halos is made up of compact, macroscopic objects (MO), such as black holes with $M_\MO >>M_{stars}$, gravitational scattering will lead to kinematic heating of the stars. Observational constraints on the amount of heating in the disk of the Milky Way put upper limits on $M_\MO \ltorder 10^{6.3}\msun $. We find limits that are three orders of magnitude more stringent by examining the heating limits in low mass stellar systems, where higher densities of dark matter and lower relative velocities disperse the stars in less than a billion years. Specifically, the structure and kinematics of the Local Group member GR8 yield a limit of $M_\MO \ltorder 6\times 10^{3}\msun$. If the properties of the dark matter are universal these results preclude the dominance of dark matter constituents in the cosmologically interesting mass range $\sim 10^6\msun$ and limit them to $M_\MO \ltorder 10^{3.7}\msun$. These results also rule out massive compact halo objects as significant contributors to the kinematic heating of the Galactic disk.