We expand our previous study on the relationship between changes in the orientation of the angular momentum vector of dark matter haloes ("spin flips") and changes in their mass (Bett & Frenk 2012), to cover the full range of halo masses in a simulation cube of length 100 $h^{-1}$ Mpc. Since strong disturbances to a halo (such as might be indicated by a large change in the spin direction) are likely also to disturb the galaxy evolving within, spin flips could be a mechanism for galaxy morphological transformation without involving major mergers. We find that 35% of haloes have, at some point in their lifetimes, had a spin flip of at least $45\deg$ that does not coincide with a major merger. Over 75% of large spin flips coincide with non-major mergers; only a quarter coincide with major mergers. We find a similar picture for changes to the inner-halo spin orientation, although here there is an increased likelihood of a flip occurring. Changes in halo angular momentum orientation, and other such measures of halo perturbation, are therefore very important quantities to consider, in addition to halo mergers, when modelling the formation and evolution of galaxies and confronting such models with observations.

We present weak lensing constraints on the ellipticity of galaxy-scale matter haloes and the galaxy-halo misalignment. Using data from the Canada-France-Hawaii Telescope Lensing Survey (CFHTLenS), we measure the weighted-average ratio of the aligned projected ellipticity components of galaxy matter haloes and their embedded galaxies, $f_\mathrm{h}$, split by galaxy type. We then compare our observations to measurements taken from the Millennium Simulation, assuming different models of galaxy-halo misalignment. Using the Millennium Simulation we verify that the statistical estimator used removes contamination from cosmic shear. We also detect an additional signal in the simulation, which we interpret as the impact of intrinsic shape-shear alignments between the lenses and their large-scale structure environment. These alignments are likely to have caused some of the previous observational constraints on $f_\mathrm{h}$ to be biased high. From CFHTLenS we find $f_\mathrm{h}=-0.04 \pm 0.25$ for early-type galaxies, which is consistent with current models for the galaxy-halo misalignment predicting $f_\mathrm{h}\simeq 0.20$. For late-type galaxies we measure $f_\mathrm{h}=0.69_{-0.36}^{+0.37}$ from CFHTLenS. This can be compared to the simulated results which yield $f_\mathrm{h}\simeq 0.02$ for misaligned late-type models.

Intrinsic galaxy alignments constitute the major astrophysical systematic of forthcoming weak gravitational lensing surveys but also yield unique insights into galaxy formation and evolution. We build analytic models for the distribution of galaxy shapes based on halo properties extracted from the Millennium Simulation, differentiating between early- and late-type galaxies as well as central galaxies and satellites. The resulting ellipticity correlations are investigated for their physical properties and compared to a suite of current observations. The best-faring model is then used to predict the intrinsic alignment contamination of planned weak lensing surveys. We find that late-type galaxy models generally have weak intrinsic ellipticity correlations, marginally increasing towards smaller galaxy separation and higher redshift. The signal for early-type models at fixed halo mass strongly increases by three orders of magnitude over two decades in galaxy separation, and by one order of magnitude from z=0 to z=2. The intrinsic alignment strength also depends strongly on halo mass, but not on galaxy luminosity at fixed mass, or galaxy number density in the environment. We identify models that are in good agreement with all observational data, except that all models over-predict alignments of faint early-type galaxies. The best model yields an intrinsic alignment contamination of a Euclid-like survey between 0.5-10% at z>0.6 and on angular scales larger than a few arcminutes. Cutting 20% of red foreground galaxies using observer-frame colours can suppress this contamination by up to a factor of two.

We use the Millennium Simulation series to study the relation between the accretion history (MAH) and mass profile of cold dark matter halos. We find that the mean density within the scale radius, r_{-2} (where the halo density profile has isothermal slope), is directly proportional to the critical density of the Universe at the time when the main progenitor's virial mass equals the mass enclosed within r_{-2}. Scaled to these characteristic values of mass and density, the mean MAH, expressed in terms of the critical density of the Universe, M(\rho_{crit}(z)), resembles that of the enclosed density profile, M(<\rho >), at z=0. Both follow closely the NFW profile, suggesting that the similarity of halo mass profiles originates from the mass-independence of halo MAHs. Support for this interpretation is provided by outlier halos whose accretion histories deviate from the NFW shape; their mass profiles show correlated deviations from NFW and are better approximated by Einasto profiles. Fitting both M(<\rho >) and M(\rho_{crit}) with either NFW or Einasto profiles yield concentration and shape parameters that are correlated, confirming and extending earlier work linking the concentration of a halo with its accretion history. These correlations also confirm that halo structure is insensitive to initial conditions: only halos whose accretion histories differ greatly from the NFW shape show noticeable deviations from NFW in their mass profiles. As a result, the NFW profile provides acceptable fits to hot dark matter halos, which do not form hierarchically, and for fluctuation power spectra other than CDM. Our findings, however, predict a subtle but systematic dependence of mass profile shape on accretion history which, if confirmed, would provide strong support for the link between accretion history and halo structure we propose here.

The statistical properties of the ellipticities of galaxy images depend on how galaxies form and evolve, and therefore constrain models of galaxy morphology, which are key to the removal of the intrinsic alignment contamination of cosmological weak lensing surveys, as well as to the calibration of weak lensing shape measurements. We construct such models based on the halo properties of the Millennium Simulation and confront them with a sample of 90,000 galaxies from the COSMOS Survey, covering three decades in luminosity and redshifts out to z=2. The ellipticity measurements are corrected for effects of point spread function smearing, spurious image distortions, and measurement noise. Dividing galaxies into early, late, and irregular types, we find that early-type galaxies have up to a factor of two lower intrinsic ellipticity dispersion than late-type galaxies. None of the samples shows evidence for redshift evolution, while the ellipticity dispersion for late-type galaxies scales strongly with absolute magnitude at the bright end. The simulation-based models reproduce the main characteristics of the intrinsic ellipticity distributions although which model fares best depends on the selection criteria of the galaxy sample. We observe fewer close-to-circular late-type galaxy images in COSMOS than expected for a sample of randomly oriented circular thick disks and discuss possible explanations for this deficit.

We use the Millennium Simulation series to study how the dynamical state of dark matter halos affects the relation between mass and concentration. We find that a large fraction of massive systems are identified when they are substantially out of equilibrium and in a particular phase of their dynamical evolution: the more massive the halo, the more likely it is found at a transient stage of high concentration. This state reflects the recent assembly of massive halos and corresponds to the first pericentric passage of recently-accreted material when, before virialization, the kinetic and potential energies reach maximum and minimum values, respectively. This result explains the puzzling upturn in the mass-concentration relation reported in recent work for massive halos; indeed, the upturn disappears when only dynamically-relaxed systems are considered in the analysis. Our results warn against applying simple equilibrium models to describe the structure of rare, massive galaxy clusters and urges caution when extrapolating scaling laws calibrated on lower-mass systems, where such deviations from equilibrium are less common. The evolving dynamical state of galaxy clusters ought to be carefully taken into account if cluster studies are to provide precise cosmological constraints.

During the growth of a cold dark matter halo, the direction of its spin can undergo rapid changes. These could disrupt or even destroy a stellar disc forming in the halo, possibly resulting in the generation of a bulge or spheroid. We investigate the frequency of significant changes in the orientation of the angular momentum vector of dark matter haloes ("spin flips"), and their degree of correlation with mergers. We focus on haloes of mass similar to that of the Milky Way (MW) halo at redshift z=0 (log M/(Msol/h) = 12.0 to 12.5) and consider flips in the spin of the whole halo or just its inner parts. We find that a greater fraction of major mergers are associated with large spin flips than minor mergers. However, since major mergers are rare, the vast majority (93%) of large whole-halo spin flips (angle change >= 45 degrees) coincide with small mass changes, not major mergers. The spin vector of the inner halo experiences much more frequent flips than the halo as a whole. Over their entire lifetimes (i.e. after a halo acquires half of its final mass), over 10% of halos experience a flip of at least 45 degrees in the spin of the entire halo and nearly 60 percent experience a flip this large in the inner halo. These numbers are reduced to 9 percent for the whole halo and 47 percent for the inner halo when we consider only haloes with no major mergers after formation. Our analysis suggests that spin flips (whose effects are not currently included in galaxy formation models) could be an important factor in the morphological transformation of disc galaxies.

We analyse the impact of galaxy--halo misalignment on the ability of weak lensing studies to constrain the shape of dark matter haloes, using a combination of the Millennium dark matter N-body simulation and different semi-analytic galaxy formation models, as well as simpler Monte Carlo tests. Since the distribution of galaxy--halo alignments is not known in detail, we test various alignment models, together with different methods of determining the halo shape. In addition to alignment, we examine the interplay of halo mass and shape, and galaxy colour and morphology with the resulting stacked projected halo shape. We find that only in the case where significant numbers of galaxy and halo minor axes are parallel does the stacked, projected halo axis ratio fall below 0.95. When using broader misalignment distributions, such as those found in recent simulations of galaxy formation, the halo ellipticity signal is washed out and would be extremely difficult to measure observationally. It is important to note that the spread in stacked halo axis ratio due to theoretical unknowns (differences between semi-analytic models, and between alignment models) are much bigger than any statistical uncertainty: It is naive to assume that, simply because LCDM predicts aspherical haloes, the stacked projected shape will be elliptical. In fact, there is no robust LCDM prediction yet for this procedure, and the interpretation of any such elliptical halo signal from lensing in terms of physical halo properties will be extremely difficult.

We investigate the magnitude and internal alignment of the angular momentum of cold dark matter haloes in simulations with and without baryons. We analyse the cumulative angular momentum profiles of hundreds of thousands of haloes in the Millennium simulation and in a smaller, but higher resolution, simulation, in total spanning 5 orders of magnitude in mass. For haloes of a given mass, the median specific angular momentum increases as j(<r) proportional to r. The direction of the vector varies considerably with radius: the median angle between the inner (< 0.25 Rvir) and total (< Rvir) angular momentum vectors is about 25degr. To investigate how baryons affect halo spin, we use another high resolution simulation, which includes gas cooling, star formation and feedback. This simulation produces a sample of galaxies with a realistic distribution of disc-to-total ratios. The formation of the galaxy spins up the dark matter within 0.1 Rvir such that the specific halo angular momentum increases by approx 50% in the median. The dark matter angular momentum becomes better aligned, but there remains a broad distribution of (mis-)alignments between the halo and the central galaxy, with a median angle between their angular momenta of ~ 30degr. Galaxies have a range of orientations relative to the shape of the halo: half of them have their minor axes misaligned by more than 45degr. Finally, we align a sample of haloes according to the orientation of their galaxies and stack the projected mass distributions. Although the individual haloes are aspherical, galaxy--halo misalignments produce a stacked mass distribution that is indistinguishable from circular. If the misalignment found in our simulations is realistic, it will be extremely difficult for weak lensing studies to measure the shape of dark haloes using this technique.

Semi-analytic models of galaxy formation typically form the spheroidal components of galaxies ("bulges"), solely through galactic major mergers. However, it is possible that non-merger events (e.g. a "fly-by" by a smaller halo) can perturb a galaxy--halo system sufficiently to form a bulge. We present a preliminary investigation into the frequency of major changes in halo and galaxy spin direction, which could be signatures of such events.

We use the Millennium Simulation (MS) to study the statistics of LCDM halo concentrations at z = 0. Our results confirm that the average halo concentration declines monotonically with mass; a power-law fits well the concentration-mass relation for over 3 decades in mass, up to the most massive objects to form in a LCDM universe (~ 10^15 h^-1 Msol). This is in clear disagreement with the predictions of the model proposed by Bullock et al. for these rare objects, and agrees better with the original predictions of Navarro, Frenk, & White. The large volume surveyed, together with the unprecedented numerical resolution of the MS, allow us to estimate with confidence the distribution of concentrations and, consequently, the abundance of systems with unusual properties. About one in a hundred cluster haloes (M200 >~ 3x10^14 h^-1 Msol) have concentrations exceeding c200 = 7.5, a result that may be used to interpret the likelihood of unusually strong massive gravitational lenses, such as Abell 1689, in the LCDM cosmogony. A similar fraction (1 in 100) of galaxy-sized haloes (M200 ~ 10^12 h^-1 Msol) have c200 < 4.5, an important constraint on models that attempt to reconcile the rotation curves of low surface-brightness galaxies by appealing to haloes of unexpectedly low concentration. We find that halo concentrations are independent of spin once haloes manifestly out of equilibrium are removed from the sample. Compared to their relaxed brethren, the concentrations of out-of-equilibrium haloes tend to be lower and to have more scatter, while their spins tend to be higher... (continue)

We investigate the spins and shapes of over a million dark matter haloes identified at z=0 in the Millennium simulation. Our sample spans halo masses ranging from dwarf galaxies to rich galaxy clusters. The very large dynamic range of this simulation enables the distribution of spins and shapes and their variation with halo mass and environment to be characterised with unprecedented precision. We compare results for haloes identified using different algorithms, including a novel one based on the branches of the halo merger trees. We investigate (and remove) biases in the estimate of angular momentum introduced by both the algorithms themselves and by numerical effects. We find that for this many objects, the traditional lognormal function is no longer an adequate description of the distribution of the dimensionless spin parameter lambda, and we provide a different function that gives a better fit. The variation of spin with halo mass is weak but detectable, although the trend depends strongly on the halo definition used. The haloes exhibit a range of shapes, with a preference for prolateness over oblateness. More massive haloes tend to be less spherical and more prolate. We find that the more spherical haloes have spin in the median, and those closest to spherical have a spin independent of mass. The most massive have a spin independent of shape. The majority of haloes have their angular momentum vector aligned with their minor axis and perpendicular to their major axis. We find that higher tend to be more clustered, with a stronger effect for more massive haloes. (abridged)