We compare the structure and substructure of dark matter halos in model universes dominated by collisional, strongly self interacting dark matter (SIDM) and collisionless, weakly interacting dark matter (CDM). While SIDM virialised halos are more nearly spherical than CDM halos, they can be rotationally flattened by as much as 20% in their inner regions. Substructure halos suffer ram-pressure truncation and drag which are more rapid and severe than their gravitational counterparts tidal stripping and dynamical friction. Lensing constraints on the size of galactic halos in clusters are a factor of two smaller than predicted by gravitational stripping, and the recent detection of tidal streams of stars escaping from the satellite galaxy Carina suggests that its tidal radius is close to its optical radius of a few hundred parsecs --- an order of magnitude smaller than predicted by CDM models but consistent with SIDM. The orbits of SIDM satellites suffer significant velocity bias $\sigma_{SIDM}/\sigma_{CDM}=0.85$ and are more circular than CDM, $\beta_{SIDM}} \approx 0.5$, in agreement with the inferred orbits of the Galaxy's satellites. In the limit of a short mean free path, SIDM halos have singular isothermal density profiles, thus in its simplest incarnation SIDM is inconsistent with galactic rotation curves.

When only cooling processes are included, smoothed-particle hydrodynamical simulations of galaxy formation in a Cold Dark Matter hierarchical clustering scenario consistently produce collapsed objects that are deficient in angular momentum by a factor of about 25 relative to the disks of observed spiral galaxies. It is widely hoped that proper allowance for star formation feedback effects will resolve this discrepancy. We explore and compare the effects of including two different types of feedback event: uniform reheating of the entire Universe to 5e5 K at a redshift z~6, and gas blow-out from pregalactic gas clouds comparable to present-day dwarf galaxies at a somewhat lower redshift. We find that blow-out is far more successful than early reheating, and that it may even be sufficient to solve the angular momentum problem. We see indications that the remaining angular momentum deficit (a factor of 5 in the blow-out models) is due to current limitations in the numerical method (SPH) used in our and other authors' work. Our most successful models distinguish themselves by the fact that a large fraction of the gas is accreted very gradually, in a cooling flow from a surrounding hot phase, rather than by mergers of massive cold clumps.

We investigate the claim by Burkert and Silk (1997) that the observed rotation curve of the dwarf irregular galaxy DDO154 cannot be reconciled with the universal CDM halo profile of Navarro, Frenk & White (1996,1997) even when allowance is made for the effect of violent gas outflow events on the structure of the galaxy. By means of N-body simulations we show that under certain conditions it is possible to obtain a reasonable fit to the observed rotation curve without invoking Burkert & Silk's proposed spheroidal MACHO component. We are able to best reproduce the observed decline in the rotation curve by postulating additional hidden disc mass, in an amount that is compatible with disc stability requirements. In the process we improve upon the results of Navarro, Eke & Frenk (1996) on the formation of halo cores by mass loss by using actual haloes from Cold Dark Matter simulations instead of Hernquist (1990) distributions.

This article describes a new, fully adaptive Particle-Multiple-Mesh numerical simulation code developed primarily for simulations of small regions (such as a group of galaxies) in a cosmological context. It integrates the equations of motion of a set of particles subject to their mutual gravitational interaction and to an arbitrary external field. The interactions are computed using a hierarchy of nested grids constructed anew at each integration step to enhance the spatial resolution in high-density regions of interest. Significant effort has gone into supporting isolated boundary conditions at the top grid level. This makes our method also applicable to non-cosmological problems, at the cost of some complications which we discuss. We point out the implications of some differences between our approach and those of other authors of similar codes, in particular with respect to the handling of the interface between regions of different spatial resolution. We present a selection of tests performed to verify the correctness and performance of our implementation. The conclusion suggests possible further improvements in the areas of independent time steps and particle softening lengths.

Models of galaxy formation ultimately aim at reproducing the _observed_ properties of galaxies. We report on work in progress to predict luminosities, colours and morphologies of field objects of various masses through N-body + Smoothed Particle Hydrodynamics simulations. We describe our method, illustrate the effects of varying the star formation assumptions, show some preliminary results, draw encouragement from their strengths as well as lessons about what improvements are most needed.