It is common practice to describe formal size and mass scales of dark matter halos as spherical overdensities with respect to an evolving density threshold. Here, we critically investigate the evolutionary effects of several such commonly used definitions and compare them to the halo evolution within fixed physical scales as well as to the evolution of other intrinsic physical properties of dark matter halos. It is shown that, in general, the traditional way of characterizing sizes and masses of halos dramatically overpredicts the degree of evolution in the last 10 Gyr, especially for low-mass halos. This pseudo-evolution leads to the illusion of growth even though there are no major changes within fixed physical scales. Such formal size definitions also serve as proxies for the virialized region of a halo in the literature. In general, those spherical overdensity scales do not coincide with the virialized region. A physically more precise nomenclature would be to simply characterize them by their very definition instead of calling such formal size and mass definitions 'virial'. In general, we find a discrepancy between the evolution of the underlying physical structure of dark matter halos seen in cosmological structure formation simulations and pseudo-evolving formal virial quantities. We question the importance of the role of formal virial quantities currently ubiquitously used in descriptions, models and relations that involve properties of dark matter structures. Concepts and relations based on pseudo-evolving formal virial quantities do not properly reflect the actual evolution of dark matter halos and lead to an inaccurate picture of the physical evolution of our universe.

We revisit the formation and evolution of the first galaxies using new hydrodynamic cosmological simulations with the ART code. Our simulations feature a recently developed model for H2 formation and dissociation, and a star formation recipe that is based on molecular rather than atomic gas. Here, we develop and implement a recipe for the formation of metal-free Population III stars in galaxy-scale simulations that resolve primordial clouds with sufficiently high density. We base our recipe on the results of prior zoom-in simulations that resolved the protostellar collapse in pre-galactic objects. We find the epoch during which Pop III stars dominated the energy and metal budget of the first galaxies to be short-lived. Galaxies which host Pop III stars do not retain dynamical signatures of their thermal and radiative feedback for more than 10^8 yr after the lives of the stars end in pair-instability supernovae, even when we consider the maximum reasonable efficiency of the feedback. Though metals ejected by the supernovae can travel well beyond the virial radius of the host galaxy, they typically begin to fall back quickly, and do not enrich a large fraction of the intergalactic medium. Galaxies with total mass in excess of 3 x 10^6 Msun re-accrete most of their baryons and transition to metal-enriched Pop II star formation.

We investigate the transition from primordial Population III (Pop III) star formation to normal Pop II star formation in the first galaxies using new cosmological hydrodynamic simulations. We find that while the first stars seed their host galaxies with metals, they cannot sustain significant outflows to enrich the intergalactic medium, even assuming a top-heavy initial mass function. This means that Pop III star formation could potentially continue until z~6 in different unenriched regions of the universe, before being ultimately shut off by cosmic reionization. Within an individual galaxy, the metal production and stellar feedback from Pop II stars overtake Pop III stars in 20-200 Myr, depending on galaxy mass.

The ever increasing size and complexity of data coming from simulations of cosmic structure formation demands equally sophisticated tools for their analysis. During the past decade, the art of object finding in these simulations has hence developed into an important discipline itself. A multitude of codes based upon a huge variety of methods and techniques have been spawned yet the question remained as to whether or not they will provide the same (physical) information about the structures of interest. Here we summarize and extent previous work of the "halo finder comparison project": we investigate in detail the (possible) origin of any deviations across finders. To this extent we decipher and discuss differences in halo finding methods, clearly separating them from the disparity in definitions of halo properties. We observe that different codes not only find different numbers of objects leading to a scatter of up to 20 per cent in the halo mass and Vmax function, but also that the particulars of those objects that are identified by all finders differ. The strength of the variation, however, depends on the property studied, e.g. the scatter in position, bulk velocity, mass, and the peak value of the rotation curve is practically below a few per cent, whereas derived quantities such as spin and shape show larger deviations. Our study indicates that the prime contribution to differences in halo properties across codes stems from the distinct particle collection methods and -- to a minor extent -- the particular aspects of how the procedure for removing unbound particles is implemented. We close with a discussion of the relevance and implications of the scatter across different codes for other fields such as semi-analytical galaxy formation models, gravitational lensing, and observables in general.

We study the detailed structure of galaxies at redshifts z > 2 using cosmological simulations with improved modeling of the interstellar medium and star formation. The simulations follow the formation and dissociation of molecular hydrogen, and include star formation only in cold molecular gas. The molecular gas is more concentrated towards the center of galaxies than the atomic gas, and as a consequence, the resulting stellar distribution is very compact. For halos with total mass above 10^{11} Mo, the median half-mass radius of the stellar disks is 0.8 kpc at z = 3. The vertical structure of the molecular disk is much thinner than that of the atomic neutral gas. Relative to the non-radiative run, the inner regions of the dark matter halo change shape from prolate to mildly oblate and align with the stellar disk. However, we do not find evidence for a significant dark disk of dark matter around the stellar disk. The outer halo regions retain the orientation acquired during accretion and mergers, and are significantly misaligned with the inner regions. The radial profile of the dark matter halo contracts in response to baryon dissipation, establishing an approximately isothermal profile throughout most of the halo. This effect can be accurately described by a modified model of halo contraction. The angular momentum of a fixed amount of inner dark matter is approximately conserved over time, while in the dissipationless case most of it is transferred outward during mergers. The conservation of the dark matter angular momentum provides supporting evidence for the validity of the halo contraction model in a hierarchical galaxy formation process.

A basic property of objects, like galaxies and halos that form in cosmological structure formation simulations, is their shape. Here, we critically investigate shape determination methods that are commonly used in the literature. It is found that using an enclosed integration volume and weight factors r^{-2} and r_{ell}^{-2} (elliptical radius) for the contribution of each particle or volume element in the shape tensor leads to biased axis ratios and smoothing of details when calculating the local shape as a function of distance from the center. To determine the local shape of matter distributions as a function of distance for well resolved objects (typically more than O(10^4) particles), we advocate a method that (1) uses an ellipsoidal shell (homoeoid) as an integration volume without any weight factors in the shape tensor and (2) removes subhalos.

[abridged] We present a detailed comparison of fundamental dark matter halo properties retrieved by a substantial number of different halo finders. These codes span a wide range of techniques including friends-of-friends (FOF), spherical-overdensity (SO) and phase-space based algorithms. We further introduce a robust (and publicly available) suite of test scenarios that allows halo finder developers to compare the performance of their codes against those presented here. This set includes mock haloes containing various levels and distributions of substructure at a range of resolutions as well as a cosmological simulation of the large-scale structure of the universe. All the halo finding codes tested could successfully recover the spatial location of our mock haloes. They further returned lists of particles (potentially) belonging to the object that led to coinciding values for the maximum of the circular velocity profile and the radius where it is reached. All the finders based in configuration space struggled to recover substructure that was located close to the centre of the host halo and the radial dependence of the mass recovered varies from finder to finder. Those finders based in phase space could resolve central substructure although they found difficulties in accurately recovering its properties. Via a resolution study we found that most of the finders could not reliably recover substructure containing fewer than 30-40 particles. However, also here the phase space finders excelled by resolving substructure down to 10-20 particles. By comparing the halo finders using a high resolution cosmological volume we found that they agree remarkably well on fundamental properties of astrophysical significance (e.g. mass, position, velocity, and peak of the rotation curve).

Merging compact binaries are the most viable and best studied candidates for gravitational wave (GW) detection by the fully operational network of ground-based observatories. In anticipation of the first detections, the expected distribution of GW sources in the local universe is of considerable interest. Here we investigate the full phase space distribution of coalescing compact binaries at $z = 0$ using dark matter simulations of structure formation. The fact that these binary systems acquire large barycentric velocities at birth ("kicks") results in merger site distributions that are more diffusely distributed with respect to their putative hosts, with mergers occurring out to distances of a few Mpc from the host halo. Redshift estimates based solely on the nearest galaxy in projection can, as a result, be inaccurate. On the other hand, large offsets from the host galaxy could aid the detection of faint optical counterparts and should be considered when designing strategies for follow-up observations. The degree of isotropy in the projected sky distributions of GW sources is found to be augmented with increasing kick velocity and to be severely enhanced if progenitor systems possess large kicks as inferred from the known population of pulsars and double compact binaries. Even in the absence of observed electromagnetic counterparts, the differences in sky distributions of binaries produced by disparate kick-velocity models could be discerned by GW observatories, within the expected accuracies and detection rates of advanced LIGO--in particular with the addition of more interferometers.

The velocity distribution function of dark matter particles is expected to show significant departures from a Maxwell-Boltzmann distribution. This can have profound effects on the predicted dark matter - nucleon scattering rates in direct detection experiments, especially for dark matter models in which the scattering is sensitive to the high velocity tail of the distribution, such as inelastic dark matter (iDM) or light (few GeV) dark matter (LDM), and for experiments that require high energy recoil events, such as many directionally sensitive experiments. Here we determine the velocity distribution functions from two of the highest resolution numerical simulations of Galactic dark matter structure (Via Lactea II and GHALO), and study the effects for these scenarios. For directional detection, we find that the observed departures from Maxwell-Boltzmann increase the contrast of the signal and change the typical direction of incoming DM particles. For iDM, the expected signals at direct detection experiments are changed dramatically: the annual modulation can be enhanced by more than a factor two, and the relative rates of DAMA compared to CDMS can change by an order of magnitude, while those compared to CRESST can change by a factor of two. The spectrum of the signal can also change dramatically, with many features arising due to substructure. For LDM the spectral effects are smaller, but changes do arise that improve the compatibility with existing experiments. We find that the phase of the modulation can depend upon energy, which would help discriminate against background should it be found.

Merging compact binaries are the one source of gravitational radiation so far identified. Because short-period systems which will merge in less than a Hubble time have already been observed as binary pulsars, they are important both as gravitational wave sources for observatories such as LIGO but also as progenitors for short gamma-ray bursts (SGRBs). The fact that these systems must have large systemic velocities implies that by the time they merge, they will be far from their formation site. The locations of merging sites depend sensitively on the gravitational potential of the galaxy host, which until now has been assumed to be static. Here we refine such calculations to incorporate the temporal evolution of the host's gravitational potential as well as that of its nearby neighbors using cosmological simulations of structure formation. This results in merger site distributions that are more diffusively distributed with respect to their putative hosts, with locations extending out to distances of a few Mpc for lighter halos. The degree of mixing between neighboring compact binary populations computed in this way is severely enhanced in environments with a high number density of galaxies. We find that SGRB redshift estimates based solely on the nearest galaxy in projection can be very inaccurate, if progenitor systems inhere large systematic kicks at birth.

The coalescence of a massive black hole (MBH) binary leads to the gravitational-wave recoil of the system and its ejection from the galaxy core. We have carried out N-body simulations of the motion of a MBH = 3.7x10^6 Msun MBH remnant in the Via Lactea I simulation, a Milky Way sized dark matter halo. The black hole receives a recoil velocity of Vkick = 80, 120, 200, 300, and 400 km/s at redshift 1.5, and its orbit is followed for over 1 Gyr within a live host halo, subject only to gravity and dynamical friction against the dark matter background. We show that, owing to asphericities in the dark matter potential, the orbit of the MBH is hightly non-radial, resulting in a significantly increased decay timescale compared to a spherical halo. The simulations are used to construct a semi-analytic model of the motion of the MBH in a time-varying triaxial Navarro-Frenk-White dark matter halo plus a spherical stellar bulge, where the dynamical friction force is calculated directly from the velocity dispersion tensor. Such a model should offer a realistic picture of the dynamics of kicked MBHs in situations where gas drag, friction by disk stars, and the flattening of the central cusp by the returning black hole are all negligible effects. We find that MBHs ejected with initial recoil velocities Vkick > 500 km/s do not return to the host center within Hubble time. In a Milky Way-sized galaxy, a recoiling hole carrying a gaseous disk of initial mass ~MBH may shine as a quasar for a substantial fraction of its wandering phase. The long decay timescales of kicked MBHs predicted by this study may thus be favorable to the detection of off-nuclear quasar activity.

We review results from recent high resolution cosmological structure formation simulations, namely the Via Lactea I & II and GHALO projects. These simulations study the formation of Milky Way sized objects within a cosmological framework. We discuss the general properties of cold dark matter halos at redshift z = 0 and focus on new insights into the structure of halos we got due to the unprecedented high resolution in these simulations.

The hierarchical theory of galaxy formation rests on the idea that smaller galactic structures merge to form the galaxies that we see today. The past decade has provided remarkable observational support for this scenario, driven in part by advances in spectroscopic instrumentation. Multi-object spectroscopy enabled the discovery of kinematically cold substructures around the Milky Way and M31 that are likely the debris of disrupting satellites. Improvements in high-resolution spectroscopy have produced key evidence that the abundance patterns of the Milky Way halo and its dwarf satellites can be explained by Galactic chemical evolution models based on hierarchical assembly. These breakthroughs have depended almost entirely on observations of nearby stars in the Milky Way and luminous red giant stars in M31 and Local Group dwarf satellites. In the next decade, extremely large telescopes will allow observations far down the luminosity function in the known dwarf galaxies, and they will enable observations of individual stars far out in the Galactic halo. The chemical abundance census now available for the Milky Way will become possible for our nearest neighbor, M31. Velocity dispersion measurements now available in M31 will become possible for systems beyond the Local Group such as Sculptor and M81 Group galaxies. Detailed studies of a greater number of individual stars in a greater number of spiral galaxies and their satellites will test hierarchical assembly in new ways because dynamical and chemical evolution models predict different outcomes for halos of different masses in different environments.

The final inspiral and coalescence of a black hole binary can produce highly beamed gravitational wave radiation. To conserve linear momentum, the black hole remnant can recoil with "kick" velocity as high as 4000 km/s. We present two sets of full N-body simulations of recoiling massive black holes (MBH) in high-resolution, non-axisymmetric potentials. The host to the first set of simulations is the main halo of the Via Lactea I simulation (Diemand et al. 2007). The nature of the resulting orbits is investigated through a numerical model where orbits are integrated assuming an evolving, triaxial NFW potential, and dynamical friction is calculated directly from the velocity dispersion along the major axes of the main halo of Via Lactea I. By comparing the triaxial case to a spherical model, we find that the wandering time spent by the MBH is significantly increased due to the asphericity of the halo. For kicks larger than 200 km/s, the remnant MBH does not return to the inner 200 pc within 1 Gyr, a timescale an order of magnitude larger than the upper limit of the estimated QSO lifetime. The second set of simulations is run using the outcome of a high-resolution gas-rich merger (Mayer et al. 2007) as host potential. In this case, a recoil velocity of 500 km/s cannot remove the MBH from the nuclear region.

We use the recently completed one billion particle Via Lactea II LambdaCDM simulation to investigate local properties like density, mean velocity, velocity dispersion, anisotropy, orientation and shape of the velocity dispersion ellipsoid, as well as structure in velocity space of dark matter haloes. We show that at the same radial distance from the halo centre, these properties can deviate by orders of magnitude from the canonical, spherically averaged values, a variation that can only be partly explained by triaxiality and the presence of subhaloes. The mass density appears smooth in the central relaxed regions but spans four orders of magnitude in the outskirts, both because of the presence of subhaloes as well as of underdense regions and holes in the matter distribution. In the inner regions the local velocity dispersion ellipsoid is aligned with the shape ellipsoid of the halo. This is not true in the outer parts where the orientation becomes more isotropic. The clumpy structure in local velocity space of the outer halo can not be well described by a smooth multivariate normal distribution. Via Lactea II also shows the presence of cold streams made visible by their high 6D phase space density. Generally, the structure of dark matter haloes shows a high degree of graininess in phase space that cannot be described by a smooth distribution function.

It is a clear unique prediction of the cold dark matter paradigm of cosmological structure formation that galaxies form hierarchically and are embedded in massive, extended dark halos teeming with self-bound substructure or "subhalos". The amount and spatial distribution of subhalos around their host provide unique information and clues on the galaxy assembly process and the nature of the dark matter. Here we present results from the Via Lactea INCITE simulation, a one billion particle, one million cpu-hour simulation of the formation and evolution of a Galactic dark matter halo and its substructure population.

Our recently completed one billion particle Via Lactea II simulation of a Milky Way-sized dark matter halo resolves over 50,000 gravitationally bound clumps orbiting today within the virialized region of the main host. About 2,300 of these subhalos have one or more "progenitors" above 1e6 Msun at redshift 11, i.e. massive enough for their gas to have cooled via excitation of H2 and fragmented prior to the epoch of cosmic reionization. We count 4,500 such progenitors: if these were able to convert a fraction of their gas content into very metal-poor stars with a Salpeter initial mass function (IMF), they would be shining today with a visual magnitude M_V=6.7 per solar mass in stars. Assuming a universal baryon fraction, we show that mean star formation efficiencies as low as 0.1% in progenitors below 1e8 Msun would overproduce the abundance of the faint Galatic dwarf spheroidals observed by the Sloan Digital Sky Survey. Star formation at first light must have occurred either with an IMF lacking stars below 0.9 Msun, or was intrinsically very inefficient in small dark matter halos. If the latter, our results may be viewed as another hint of a minimum scale in galaxy formation.

We perform a series of simulations of a Galactic mass dark matter halo at different resolutions, our largest uses over three billion particles and has a mass resolution of 1000 M_sun. We quantify the structural properties of the inner dark matter distribution and study how they depend on numerical resolution. We can measure the density profile to a distance of 120 pc (0.05% of R_vir) where the logarithmic slope is -0.8 and -1.4 at (0.5% of R_vir). We propose a new two parameter fitting function that has a linearly varying logarithmic density gradient which fits the GHALO and VL2 density profiles extremely well. Convergence in the density profile and the halo shape scales as N^(-1/3), but the shape converges at a radius three times larger at which point the halo becomes more spherical due to numerical resolution. The six dimensional phase-space profile is dominated by the presence of the substructures and does not follow a power law, except in the smooth under-resolved inner few kpc.

In cold dark matter cosmological models, structures form and grow by merging of smaller units. Numerical simulations have shown that such merging is incomplete; the inner cores of halos survive and orbit as "subhalos" within their hosts. Here we report a simulation that resolves such substructure even in the very inner regions of the Galactic halo. We find hundreds of very concentrated dark matter clumps surviving near the solar circle, as well as numerous cold streams. The simulation reveals the fractal nature of dark matter clustering: Isolated halos and subhalos contain the same relative amount of substructure and both have cuspy inner density profiles. The inner mass and phase-space densities of subhalos match those of recently discovered faint, dark matter-dominated dwarf satellite galaxies and the overall amount of substructure can explain the anomalous flux ratios seen in strong gravitational lenses. Subhalos boost gamma-ray production from dark matter annihilation, by factors of 4-15, relative to smooth galactic models. Local cosmic ray production is also enhanced, typically by a factor 1.4, but by more than a factor of ten in one percent of locations lying sufficiently close to a large subhalo. These estimates assume that gravitational effects of baryons on dark matter substructure are small.

We present a simple and efficient method to set up spherical structure models for N-body simulations with a multimass technique. This technique reduces by a substantial factor the computer run time needed in order to resolve a given scale as compared to single-mass models. It therefore allows to resolve smaller scales in N-body simulations for a given computer run time. Here, we present several models with an effective resolution of up to 1.68 x 10^9 particles within their virial radius which are stable over cosmologically relevant time-scales. As an application, we confirm the theoretical prediction by Dehnen (2005) that in mergers of collisonless structures like dark matter haloes always the cusp of the steepest progenitor is preserved. We model each merger progenitor with an effective number of particles of approximately 10^8 particles. We also find that in a core-core merger the central density approximately doubles whereas in the cusp-cusp case the central density only increases by approximately 50%. This may suggest that the central region of flat structures are better protected and get less energy input through the merger process.

We present a new time-stepping criterion for N-body simulations that is based on the true dynamical time of a particle. This allows us to follow the orbits of particles correctly in all environments since it has better adaptivity than previous time-stepping criteria used in N-body simulations. Furthermore, it requires far fewer force evaluations in low density regions of the simulation and has no dependence on artificial parameters such as, for example, the softening length. This can be orders of magnitude faster than conventional ad-hoc methods that employ combinations of acceleration and softening and is ideally suited for hard problems, such as obtaining the correct dynamics in the very central regions of dark matter haloes. We also derive an eccentricity correction for a general leapfrog integration scheme that can follow gravitational scattering events for orbits with eccentricity e -> 1 with high precision. These new approaches allow us to study a range of problems in collisionless and collisional dynamics from few body problems to cosmological structure formation. We present tests of the time-stepping scheme in N-body simulations of 2-body orbits with eccentricity e -> 1 (elliptic and hyperbolic), equilibrium haloes and a hierarchical cosmological structure formation run.

The dark matter dominated Fornax dwarf spheroidal has five globular clusters orbiting at ~1 kpc from its centre. In a cuspy CDM halo the globulars would sink to the centre from their current positions within a few Gyrs, presenting a puzzle as to why they survive undigested at the present epoch. We show that a solution to this timing problem is to adopt a cored dark matter halo. We use numerical simulations and analytic calculations to show that, under these conditions, the sinking time becomes many Hubble times; the globulars effectively stall at the dark matter core radius. We conclude that the Fornax dwarf spheroidal has a shallow inner density profile with a core radius constrained by the observed positions of its globular clusters. If the phase space density of the core is primordial then it implies a warm dark matter particle and gives an upper limit to its mass of ~0.5 keV, consistent with that required to significantly alleviate the substructure problem.

The Milky Way contains several distinct old stellar components that provide a fossil record of its formation. We can understand their spatial distribution and kinematics in a hierarchical formation scenario by associating the proto-galactic fragments envisaged by Searle and Zinn (1978) with the rare peaks able to cool gas in the cold dark matter density field collapsing at redshift z>10. We use hierarchical structure formation simulations to explore the kinematics and spatial distribution of these early star-forming structures in galaxy haloes today. Most of the proto-galaxies rapidly merge, their stellar contents and dark matter becoming smoothly distributed and forming the inner Galactic halo. The metal-poor globular clusters and old halo stars become tracers of this early evolutionary phase, centrally biased and naturally reproducing the observed steep fall off with radius. The most outlying peaks fall in late and survive to the present day as satellite galaxies. The observed radial velocity dispersion profile and the local radial velocity anisotropy of Milky Way halo stars are successfully reproduced in this model. If this epoch of structure formation coincides with a suppression of further cooling into lower sigma peaks then we can reproduce the rarity, kinematics and spatial distribution of satellite galaxies as suggested by Bullock et al. (2000). Reionisation at z=12+/-2 provides a natural solution to the missing satellites problem. Measuring the distribution of globular clusters and halo light on scales from galaxies to clusters could be used to constrain global versus local reionisation models. If reionisation occurs contemporary, our model predicts a constant frequency of blue globulars relative to the host halo mass, except for dwarf galaxies where the average relative frequencies become smaller.

We resolve the inner region of a massive cluster forming in a cosmological LCDM simulation with a mass resolution of 2*10^6 Msun and before z=4.4 even 3*10^5 Msun. This is a billion times less than the clusters final virial mass and a substantial increase over current LCDM simulations. We achieve this resolution using a new multi-mass refinement procedure and are now able to probe a dark matter halo density profile down to 0.1 percent of the virial radius. The inner density profile of this cluster halo is well fitted by a power-law rho ~ r^-gamma down to the smallest resolved scale. An inner region with roughly constant logarithmic slope is now resolved, which suggests that cuspy profiles describe the inner profile better than recently proposed profiles with a core. The cluster studied here is one out of a sample of six high resolution cluster simulations of Diemand et al. (2004) and it's inner slope of gamma = 1.2 lies close to the sample average.

Several general trends have been identified for equilibrated, self-gravitating collisionless systems, such as density or anisotropy profiles. These are integrated quantities which naturally depend on the underlying velocity distribution function (VDF) of the system. We study this VDF through a set of numerical simulations, which allow us to extract both the radial and the tangential VDF. We find that the shape of the VDF is universal, in the sense that it depends only on two things namely the dispersion (radial or tangential) and the local slope of the density. Both the radial and the tangential VDF's are universal for a collection of simulations, including controlled collisions with very different initial conditions, radial infall simulation, and structures formed in cosmological simulations.