We study the accretion along streams from the cosmic web into high-redshift massive galaxies using three sets of AMR hydro-cosmological simulations. We find that the streams keep a roughly constant accretion rate as they penetrate into the halo centre. The mean accretion rate follows the mass and redshift dependence predicted for haloes by the EPS approximation, dM / dt is proportional to Mvir^{1.25} (1 + z)^{2.5}. The distribution of the accretion rates can well be described by a sum of two Gaussians, the primary corresponding to "smooth inflow" and the secondary to "mergers". The same functional form was already found for the distributions of specific star formation rates in observations. The mass fraction in the smooth component is 60 - 90 %, insensitive to redshift or halo mass. The simulations with strong feedback show clear signs of reaccretion due to recycling of galactic winds. The mean accretion rate for the mergers is a factor 2 - 3 larger than that of the smooth component. The standard deviation of the accretion rate is 0.2 - 0.3 dex, showing no trend with mass or redshift. For the smooth component it is 0.12 - 0.24 dex.

Recent observations have shown that the majority of the Andromeda galaxy's satellites are aligned in a thin plane. On the theoretical side it has been proposed that galaxies acquire their gas via cold streams. In addition, numerical simulations show that the same streams also deliver gas clumps which could potentially develop into satellite galaxies. Assuming that cold streams are a major source of satellite systems around galaxies we calculate the probabilities in two different models to find a certain fraction of satellites within a thin plane around the central galaxy of the host halo with and without having the same sense of rotation. Using simple geometrical considerations and adopting a random orientation of the streams we demonstrate that the vast thin disk of satellites detected around Andromeda can naturally be explained within this framework. In fact, without any satellite scattering, two streams would lead to too many satellites in the thin plane, compared with the observations. Three streams reproduce the observations very well. Natural implications from our model are that all massive galaxies should have a thin plane of satellites and that the satellites should naturally distribute themselves not only into a single plane but into several inclined ones. We estimate the effect of additional satellites accreted from random directions and find it to be of minor relevance for a mild inflow of satellites from random directions.

We study the velocities of the accretion along streams from the cosmic web into massive galaxies at high redshift with the help of three different suites of AMR hydrodynamical cosmological simulations. The results are compared to free-fall velocities and to the sound speeds of the hot ambient medium. The sound speed of the hot ambient medium is calculated using two different methods to determine the medium's temperature. We find that the simulated cold stream velocities are in violent disagreement with the corresponding free-fall profiles. The sound speed is a better albeit not always correct description of the cold flows' velocity. Using these calculations as a first order approximation for the gas inflow velocities vinflow = 0.9 vvir is given. We conclude from the hydrodynamical simulations as our main result that the velocity profiles for the cold streams are constant with radius. These constant inflow velocities seem to have a "parabola-like" dependency on the host halo mass in units of the virial velocity that peaks at Mvir = 10^{12} Msun and we also propose that the best fit functional form for the dependency of the inflow velocity on the redshift is a square root power law relation: vinflow is proportional to sqrt(z + 1) vvir.

Cold gas streaming along the dark-matter filaments of the cosmic web is predicted to be the major provider of resources for disc buildup, violent disk instability and star formation in massive galaxies in the early universe. We study to what extent these cold streams are traceable in the extended circum-galactic environment of galaxies via Ly alpha absorption and selected low ionisation metal absorption lines. We model the expected absorption signatures using high resolution zoom-in AMR cosmological simulations. In the postprocessing, we distinguish between self-shielded gas and unshielded gas. In the self-shielded gas, which is optically thick to Lyman continuum radiation, we assume pure collisional ionisation for species with an ionisation potential greater than 13.6 eV. In the optically thin, unshielded gas these species are also photoionised by the meta-galactic radiation. In addition to absorption of radiation from background quasars, we compute the absorption line profiles of radiation emitted by the galaxy at the centre of the same halo. We predict the strength of the absorption signal for individual galaxies without stacking. We find that the Ly alpha absorption profiles produced by the streams are consistent with observations of absorption and emission Ly alpha profiles in high redshift galaxies. Due to the low metallicities in the streams, and their low covering factors, the metal absorption features are weak and difficult to detect.

We perform a detailed investigation into the disruption of central cusps via the transfer of energy from sinking massive objects. Constant density inner regions form at the radius where the enclosed mass approximately matches the mass of the infalling body. We explore parameter space using numerical simulations and give an empirical relation for the size of the resulting core within structures that have different initial cusp slopes. We find that infalling bodies always stall at the edge of these newly formed cores, experiencing no dynamical friction over many dynamical times. As applications, we consider the resulting decrease in the dark matter annihilation flux due to centrally destroyed cusps, and we present a new theory for the formation of close binary nuclei -- the `stalled binary' model. We focus on one particularly interesting binary nucleus system, the dwarf spheroidal galaxy VCC 128 which is dark matter dominated at all radii. We show that its nuclei would rapidly coalesce within a few million years if it has a central dark matter cusp slope steeper than r^{-1}. However, if its initial dark matter cusp is slightly shallower than a log slope of -0.75 at ~0.1% of the virial radius, then the sinking nuclei naturally create a core equal to their observed separation and stall. This is close to the log slope measured a recent billion particle CDM halo simulation.

We use high-resolution cosmological hydrodynamical AMR simulations to predict the characteristics of La emission from the cold gas streams that fed galaxies in massive haloes at high redshift. The La luminosity in our simulations is powered by the release of gravitational energy as gas flows from the intergalactic medium into the halo potential wells. The UV background contributes only <20% to the gas heating. The La emissivity is due primarily to electron-impact excitation cooling radiation in gas ~2x10^4K. We calculate the La emissivities assuming collisional ionisation equilibrium (CIE) at all gas temperatures. The simulated streams are self-shielded against the UV background, so photoionisation and recombination contribute negligibly to the La line formation. We produce theoretical maps of the La surface brightnesses, assuming that ~85% of the La photons are directly observable. We find that typical haloes of mass Mv~10^12-13 Msun at z~3 emit as La blobs (LABs) with luminosities 10^43-44 erg/s. Most of the La comes from the extended narrow, partly clumpy, inflowing, cold streams that feed the growing galaxies. The predicted LAB morphology is therefore irregular, with dense clumps and elongated extensions. The linewidth is expected to range from 10^2 to more than 10^3 km/s with a large variance. The typical La surface brightness profile is proportional to r^-1.2 where r is the distance from the halo centre. Our simulated LABs are similar in luminosity, morphology and extent to the observed LABs, with distinct kinematic features. The predicted La luminosity function is consistent with observations, and the predicted areas and linewidths roughly recover the observed scaling relations. This mechanism for producing LABs appears inevitable in many high-z galaxies. Some of the LABs may thus be regarded as direct detections of the cold streams that drove galaxy evolution at high z.

The massive galaxies in the young universe, ten billion years ago, formed stars at surprising intensities. Although this is commonly attributed to violent mergers, the properties of many of these galaxies are incompatible with such events, showing gas-rich, clumpy, extended rotating disks not dominated by spheroids (Genzel et al. 2006, 2008). Cosmological simulations and clustering theory are used to explore how these galaxies acquired their gas. Here we report that they are stream-fed galaxies, formed from steady, narrow, cold gas streams that penetrate the shock-heated media of massive dark matter haloes (Dekel & Birnboim 2006; Keres et al. 2005). A comparison with the observed abundance of star-forming galaxies implies that most of the input gas must rapidly convert to stars. One-third of the stream mass is in gas clumps leading to mergers of mass ratio greater than 1:10, and the rest is in smoother flows. With a merger duy cycle of 0.1, three-quarters of the galaxies forming stars at a given rate are fed by smooth streams. The rarer, submillimetre galaxies that form stars even more intensely are largely merger-induced starbursts. Unlike destructive mergers, the streams are likely to keep the rotating disk configuration intact, although turbulent and broken into giant star-forming clumps that merge into a central spheroid (Noguchi 1999; Genzel et al. 2008, Elmegreen, Bournaud & Elmegreen 2008, Dekel, Sari & Ceverino 2009). This stream-driven scenario for the formation of disks and spheroids is an alternative to the merger picture.

Ultra compact dwarf galaxies (UCDs) have similar properties as massive globular clusters or the nuclei of nucleated galaxies. Recent observations suggesting a high dark matter content and a steep spatial distribution within groups and clusters provide new clues as to their origins. We perform high-resolution N-body / smoothed particle hydrodynamics simulations designed to elucidate two possible formation mechanisms for these systems: the merging of globular clusters in the centre of a dark matter halo, or the massively stripped remnant of a nucleated galaxy. Both models produce density profiles as well as the half light radii that can fit the observational constraints. However, we show that the first scenario results to UCDs that are underluminous and contain no dark matter. This is because the sinking process ejects most of the dark matter particles from the halo centre. Stripped nuclei give a more promising explanation, especially if the nuclei form via the sinking of gas, funneled down inner galactic bars, since this process enhances the central dark matter content. Even when the entire disk is tidally stripped away, the nucleus stays intact and can remain dark matter dominated even after severe stripping. Total galaxy disruption beyond the nuclei only occurs on certain orbits and depends on the amount of dissipation during nuclei formation. By comparing the total disruption of CDM subhaloes in a cluster potential we demonstrate that this model also leads to the observed spatial distribution of UCDs which can be tested in more detail with larger data sets.

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.

If the dark matter particle is a neutralino then the first structures to form are cuspy cold dark matter (CDM) haloes collapsing after redshifts z ~ 100 in the mass range 10^{-6} - 10^{-3} Msun. We carry out a detailed study of the survival of these micro-haloes in the Galaxy as they experience tidal encounters with stars, molecular clouds, and other dark matter substructures. We test the validity of analytic impulsive heating calculations using high resolution N-body simulations. A major limitation of analytic estimates is that mean energy inputs are compared to mean binding energies, instead of the actual mass lost from the system. This energy criterion leads to an overestimate of the stripped mass and underestimate of the disruption timescale since CDM haloes are strongly bound in their inner parts. We show that a significant fraction of material from CDM micro-haloes can be unbound by encounters with Galactic substructure and stars, however the cuspy central regions remain relatively intact. Furthermore, the micro-haloes near the solar radius are those which collapse significantly earlier than average and will suffer very little mass loss. Thus we expect a fraction of surviving bound micro-haloes, a smooth component with narrow features in phase space, which may be uncovered by direct detection experiments, as well as numerous surviving cuspy cores with proper motions of arc-minutes per year, which can be detected indirectly via their annihilation into gamma-rays.

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.