We use the Millennium Simulation to measure the cross-correlation between halo centres and mass (or equivalently the average density profiles of dark haloes) in a LCDM cosmology. We present results for radii in the range 10 kpc/h < r < 30 Mpc/h for halo masses in the range 4e10 Msol/h < M200 < 4e14 Msol/h. Both at z=0 and at z=0.76 these cross-correlations are surprisingly well fit by approximating the inner region by a density profile of NFW or Einasto form, the outer region by a biased version of the linear mass autocorrelation function, and by adopting the maximum of the two where they are comparable. We use a simulation of the formation of galaxies within the Millennium Simulation to explore how these results are reflected in cross-correlations between galaxies and mass. These are directly observable through galaxy-galaxy lensing. Here also we find that simple models can represent the simulation results remarkably well, typically to < 10%. Such models can be used to extend our results to other redshifts, to cosmologies with other parameters, and to other assumptions about how galaxies populate dark haloes. The characteristic features predicted in the galaxy-galaxy lensing signal should provide a strong test of the LCDM cosmology as well as a route to understanding how galaxies form within it.

We use a set of cosmological N-body simulations to investigate the structural shape of galaxy-sized cold dark matter (CDM) halos. Unlike most previous work on the subject - which dealt with shapes as measured by the inertia tensor - we focus here on the shape of the gravitational potential, a quantity more directly relevant to comparison with observational probes. A further advantage is that the potential is less sensitive to the effects of substructure and, as a consequence, the isopotential surfaces are typically smooth and well approximated by concentric ellipsoids. Our main result is that the asphericity of the potential increases rapidly towards the center of the halo. The radial trend is more pronounced than expected from constant flattening in the mass distribution, and reflects a strong tendency for dark matter halos to become increasingly aspherical inwards. Near the center the halo potential is approximately prolate ((c/a)_0=0.72 +/- 0.04, (b/a)_0=0.78 +/- 0.08), but it becomes increasingly spherical in the outer regions. The principal axes of the isopotential surfaces remain well aligned, and in most halos the angular momentum tends to be parallel to the minor axis and perpendicular to the major axis. This suggests that galactic disks may form in a plane where the potential is elliptical and where its ellipticity varies rapidly with radius. This can result in significant deviations from circular motion in systems such as low surface brightness galaxies (LSBs), even for relatively minor deviations from circular symmetry. Simulated long-slit rotation curves can appear similar to those of LSBs often cited as evidence for constant density "cores". This suggests that taking into account the 3D shape of the dark mass distribution might help to reconcile such evidence with the cuspy mass profile of CDM halos.

We study the kinematics of gaseous disks in triaxial dark matter halos using the closed-loop orbit solutions in non-axisymmetric potentials. The orbits are in general non-circular and, for given triaxiality, their ellipticity depends on the ratio of escape to circular velocities, V_esc^2/V_c^2. This ratio increases steeply towards the center for cold dark matter (CDM) halo density profiles, implying that even minor deviations from spherical symmetry may induce large deviations from circular orbits in the velocity field of a gaseous disk, especially near the center. This result suggests that caution should be exercised when interpreting constraints on the presence of density cusps in the dark halo derived from the innermost velocity profile. Simulated long-slit rotation curves vary greatly in shape, depending primarily on the viewing angle of the disk and on its orientation relative to the principal axes of the potential. "Solid-body" rotation curves - typically interpreted as a signature of a constant density core in the dark matter distribution - are often obtained when the slit samples velocities near the major axis of the closed loop orbits. Triaxial potentials imprint specific symmetries in 2D velocity fields, generally inducing "twists" in the isovelocity contours and anti-symmetric patterns in opposite quadrants. We suggest that triaxial halos may be responsible for the variety of shapes of long-slit rotation curves of low surface brightness (LSB) galaxies, as well as for the complex central kinematics of LSBs, which are sometimes ascribed to the presence of "radial motions" in the gas. We argue that LSB rotation curves might be reconciled with the structure of CDM halos once the effects of halo triaxiality on the dynamics of gaseous disks are properly taken into account.

The galaxy cluster 1E 0657-56 has a bullet-like subcluster that is moving away from the centre of the main cluster at high speed. Markevitch et al. (2004) recently estimated a relative velocity of V_bullet = 4500 +1100/-800 km/s, based on observations of the bow shock in front of the subcluster. The weak lensing analysis of Clowe et al. (2004) indicates that a substantial secondary mass peak is associated with this subcluster. We estimate the likelihood of such a configuration by examining the distribution of subhalo velocities for clusters in the Millennium Run, a large LCDM cosmological simulation. We find that the most massive subhalo has a velocity as high as that of the bullet subcluster in only about 1 out of every 100 cluster-sized halos. This estimate is strongly dependent on the precise velocity adopted for the bullet. One of the ten most massive subhalos has such a high velocity about 40% of the time. We conclude that the velocity of the bullet subcluster is not exceptionally high for a cluster substructure, and can be accommodated within the currently favoured LCDM comogony.

Subramanian, Pujari and Becker (2004) claim that the correct expression for the angular momentum transport in an accretion disc, which is proportional to dOmega/dR, can be derived on the basis of the analysis of the epicyclic motion of gas parcels in adjacent eddies in the disc. We study their argument and show that their derivation contains several fundamental errors: 1) the biased choice of the desired formula from an infinite number of formulae; 2) the biased choice of parcel trajectories; and 3) confusion regarding the reference frames. Following 1) we could derive, for example, a (invalid) formula in which the angular momentum transport is proportional to dv_phi/dR, and from 2) we could even prove that the angular momentum transport is either inward or null. We present the correct approach to the problem of angular momentum transport in an accretion disc in terms of mean free path theory.

Clarke and Pringle (2004) derived a proper viscosity formula in a rotating gas by applying mean free path theory. We study their argument in detail and show that their result can be derived with a much simpler calculational procedure and physically clearer picture.

We use N-body hydrodynamical simulations to study the structure of disks in triaxial potentials resembling CDM halos. Our analysis focuses on the accuracy of the dark mass distribution inferred from rotation curves derived from simulated long-slit spectra. We consider a massless disk embedded in a halo with axis ratios of 0.5:0.6:1.0 and with its rotation axis aligned with the minor axis of the halo. Closed orbits for the gaseous particles deviate from coplanar circular symmetry, resulting in a variety of long-slit rotation curve shapes, depending on the orientation of the disk relative to the line of sight. Rotation curves may thus differ significantly from the spherically-averaged circular velocity profile of the dark matter halo. "Solid-body" rotation curves--typically interpreted as a signature of a constant density core in the dark matter distribution--are obtained about 25% of the time for random orientations although the dark matter follows the cuspy density profile proposed by Navarro, Frenk & White (NFW). We conclude that the discrepancies reported between the shape of the rotation curve of low surface brightness galaxies and the structure of CDM halos may be resolved once the complex effects of halo triaxiality on the dynamics of the gas component is properly taken into account.

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 investigate the mass profile of LambdaCDM halos using a suite of numerical simulations spanning five decades in halo mass, from dwarf galaxies to rich galaxy clusters. Our analysis confirms the proposal of Navarro, Frenk & White (NFW) that the shape of LambdaCDM halo mass profiles differs strongly from a power law and depends little on mass. The logarithmic slope of the spherically-averaged density profile, as measured by beta=-dln(rho)/dln(r), decreases monotonically towards the center and becomes shallower than isothermal (beta<2) inside a characteristic radius, r_{-2}. Although the fitting formula proposed by NFW provides a reasonably good approximation to the density and circular velocity profiles of individual halos, systematic deviations from the best NFW fits are also noticeable. Inside r_{-2}, the profile of simulated halos gets shallower with radius more gradually than predicted and, as a result, NFW fits tend to underestimate the dark matter density in these regions. This discrepancy has been interpreted as indicating a steeply divergent cusp, but our results suggest a different interpretation. We use the density and enclosed mass at our innermost resolved radii to place strong constraints on beta_{0}: density cusps as steep as r^{-1.5} are inconsistent with most of our simulations, although beta_{0}=1 is still consistent with our data. Our density profiles show no sign of converging to a well-defined asymptotic inner power law. We propose a simple formula that reproduces the radial dependence of the slope better than the NFW profile, and so may minimize errors when extrapolating our results inward to radii not yet reliably probed by numerical simulations.

We use a set of high-resolution cosmological N-body simulations to investigate the inner mass profile of galaxy-sized cold dark matter (CDM) halos. These simulations extend the thorough numerical convergence study presented in Paper I of this series (Power et al. 2003), and demonstrate that the mass profile of CDM halos can be robustly estimated beyond a minimum converged radius of order r_conv ~ 1 kpc/h in our highest resolution runs. The density profiles of simulated halos become progressively shallow from the virial radius inwards, and show no sign of approaching a well-defined power-law behaviour near the centre. At r_conv, the logarithmic slope of the density profile is steeper than the asymptotic \rho \propto r^-1 expected from the formula proposed by Navarro, Frenk, and White (1996), but significantly shallower than the steeply divergent \rho \propto r^-1.5 cusp proposed by Moore et al. (1999). We perform a direct comparison of the spherically-averaged dark matter circular velocity (V_c) profiles with rotation curves of low surface brightness (LSB) galaxies from the samples of de Blok et al. (2001), de Blok and Bosma (2002), and Swaters et al. (2003). Most (about two-thirds) LSB galaxies in this dataset are roughly consistent with CDM halo V_c profiles. However, about one third of LSBs in these samples feature a sharp transition between the rising and flat part of the rotation curve that is not seen in the V_c profiles of CDM halos. This discrepancy has been interpreted as excluding the presence of cusps, but we argue that it might simply reflect the difference between circular velocity and gas rotation speed likely to arise in gaseous disks embedded within realistic, triaxial CDM halos.

We investigate the evolution of substructure in cold dark matter halos using N-body simulations of tidal stripping of substructure halos (subhalos) within a static host potential. We find that halos modeled following the Navarro, Frenk & White (NFW) mass profile lose mass continuously due to tides from the massive host, leading to the total disruption of satellite halos with small tidal radii. The structure of stripped NFW halos depends mainly on the fraction of mass lost, and can be expressed in terms of a simple correction to the original NFW profile. We apply these results to substructure in the Milky Way, and conclude that the dark matter halos surrounding its dwarf spheroidal (dSph) satellites have circular velocity curves that peak well beyond the luminous radius at velocities significantly higher than expected from the stellar velocity dispersion. Our modeling suggests that the true tidal radii of dSphs lie well beyond the putative tidal cutoff observed in the surface brightness profile, suggesting that the latter are not really tidal in origin but rather features in the light profile of limited dynamical relevance. For Draco, in particular, our modeling implies that its tidal radius is much larger than derived by Irwin & Hatzidimitriou (1995), lending support to the interpretation of recent Sloan survey data by Odenkirchen et al. (2001). Similarly, our model suggests that Carina's halo has a peak circular velocity of ~55 km/s, which may help explain how this small galaxy has managed to retain enough gas to undergo several bursts of star formation. Our results imply a close correspondence between the most massive subhalos expected in a CDM universe and the known satellites of the Milky Way, and suggest that only subhalos with peak circular velocities below 35 km/s lack readily detectable luminous counterparts.

The derivation of the viscosity formula of accretion disks appearing in a number of textbooks is based on a mean free path theory of gas particles. However, this procedure, when followed precisely, leads to the incorrect conclusion that the angular momentum of an accretion disk flows from the outer part to the inner part.

The history of hydrodynamic numerical simulations for accretion disks in close binary systems is reviewed, in which emphasis is placed, in particular, on the facts that spiral shock waves were numerically found in 1986 by researchers including one of the present authors and that spiral structure was discovered in IP Pegasi in 1997 by Steeghs et al. The results of our two and three-dimensional numerical simulations in recent years are then summarized, with comparison being made with observations.