The radial acceleration relation (RAR) locally relates the `observed' acceleration inferred from the dynamics of a system to the acceleration implied by its baryonic matter distribution. The relation as traced by galaxy rotation curves is one-to-one with remarkably little scatter, implying that the dynamics of a system can be predicted simply by measuring its density profile as traced by e.g. stellar light or gas emission lines. Extending the relation to accelerations below those usually probed by practically observable kinematic tracers is challenging, especially once accounting for faintly emitting baryons, such as the putative warm-hot intergalactic medium, becomes important. We show that in the low-acceleration regime, the (inverted) RAR predicts an unphysical, declining enclosed baryonic mass profile for systems with `observed' acceleration profiles steeper than $g_{\rm obs}\propto r^{-1}$ (corresponding to density profiles steeper than isothermal - $\rho(r)\propto r^{-2}$). If the RAR is tantamount to a natural law, such acceleration profiles cannot exist. We apply this argument to test the compatibility of an extrapolation of the rotation curve-derived RAR to low accelerations with data from galaxy-galaxy weak lensing, dwarf spheroidal galaxy stellar kinematic, and outer Milky~Way dynamical measurements, fully independent of the uncertainties inherent in direct measurements of the baryonic matter distribution. In all cases we find that the data weakly favour a break to a steeper low-acceleration slope. Improvements in measurements and modelling of the outer Milky~Way, and weak lensing, seem like the most promising path toward stronger constraints on the low-acceleration behaviour of the RAR.

We address the issue of numerical convergence in cosmological smoothed particle hydrodynamics simulations using a suite of runs drawn from the EAGLE project. Our simulations adopt subgrid models that produce realistic galaxy populations at a fiducial mass and force resolution, but systematically vary the latter in order to study their impact on galaxy properties. We provide several analytic criteria that help guide the selection of gravitational softening for hydrodynamical simulations, and present results from runs that both adhere to and deviate from them. Unlike dark matter-only simulations, hydrodynamical simulations exhibit a strong sensitivity to gravitational softening, and care must be taken when selecting numerical parameters. Our results--which focus mainly on star formation histories, galaxy stellar mass functions and sizes--illuminate three main considerations. First, softening imposes a minimum resolved escape speed, $v_\epsilon$, due to the binding energy between gas particles. Runs that adopt such small softening lengths that $v_\epsilon \gt 10\,{\rm km s^{-1}}$ (the sound speed in ionised $\sim 10^4\,{\rm K}$ gas) suffer from reduced effects of photo-heating. Second, feedback from stars or active galactic nuclei may suffer from numerical over-cooling if the gravitational softening length is chosen below a critical value, $\epsilon_{\rm eFB}$. Third, we note that small softening lengths exacerbate the segregation of stars and dark matter particles in halo centres, often leading to the counter-intuitive result that galaxy sizes {\em increase} as softening is reduced. The structure of dark matter haloes in hydrodynamical runs respond to softening in a way that reflects the sensitivity of their galaxy populations to numerical parameters.

We use a compilation of disc galaxy rotation curves to assess the role of the luminous component ("baryons") in the rotation curve diversity problem. As in earlier work, we find that rotation curve shape correlates with baryonic surface density: high surface density galaxies have rapidly-rising rotation curves consistent with cuspy cold dark matter halos; slowly-rising rotation curves (characteristic of galaxies with inner mass deficits or "cores") occur only in low surface density galaxies. The correlation, however, seems too weak to be the main driver of the diversity. In addition, dwarf galaxies exhibit a clear trend, from "cuspy" systems where baryons are unimportant in the inner mass budget to "cored" galaxies where baryons actually dominate. This trend constrains the various scenarios proposed to explain the diversity, such as (i) baryonic inflows and outflows during galaxy formation; (ii) dark matter self-interactions; (iii) variations in the baryonic mass structure coupled to rotation velocities through the "mass discrepancy-acceleration relation" (MDAR); or (iv) non-circular motions in gaseous discs. Together with analytical modeling and cosmological hydrodynamical simulations, our analysis shows that each of these scenarios has promising features, but none seems to fully account for the observed diversity. The MDAR, in particular, is inconsistent with the observed trend between rotation curve shape and baryonic importance; either the trend is caused by systematic errors in the data or the MDAR does not apply. The origin of the dwarf galaxy rotation curve diversity and its relation to the structure of cold dark matter halos remains an open issue.

Modelling the molecular gas that is routinely detected through CO observations of high-redshift galaxies constitutes a major challenge for ab initio simulations of galaxy formation. We carry out a suite of cosmological hydrodynamic simulations in order to compare three approximate methods that have been used in the literature to track the formation and evolution of the simplest and most abundant molecule, H$_2$. Namely, we consider: i) a semi-empirical procedure that associates H$_2$ to dark-matter haloes based on a series of scaling relations inferred from observations; ii) a model that assumes chemical equilibrium between the H$_2$ formation and destruction rates; and iii) a model that fully solves the out-of-equilibrium rate equations and accounts for the unresolved structure of molecular clouds. We study the impact of finite spatial resolution and show that robust H$_2$ masses at redshift $z\approx 4$ can only be obtained for galaxies that are sufficiently metal enriched in which H$_2$ formation is fast. This corresponds to H$_2$ reservoirs with masses $M_{\mathrm{H_2}}\gtrsim 6\times 10^9 \mathrm{M}_\odot$. In this range, equilibrium and non-equilibrium models predict similar molecular masses (but different galaxy morphologies) while the semi-empirical method produces less H$_2$. The star-formation rates as well as the stellar and H$_2$ masses of the simulated galaxies are in line with those observed in actual galaxies at similar redshifts that are not massive starbursts. The H$_2$ mass functions extracted from the simulations at $z\approx 4$ agree well with recent observations that only sample the high-mass end. However, our results indicate that most molecular material at high-$z$ lies yet undetected in reservoirs with $10^9<M_{\mathrm H_2}<10^{10} \mathrm{M}_\odot$.

The nature of the dark matter can affect the collapse time of dark matter haloes, and can therefore be imprinted in observables such as the stellar population ages and star formation histories of dwarf galaxies. In this paper we use high resolution hydrodynamical simulations of Local Group-analogue (LG) volumes in cold dark matter (CDM), sterile neutrino warm dark matter (WDM) and self-interacting dark matter (SIDM) models with the EAGLE galaxy formation code to study how galaxy formation times change with dark matter model. We are able to identify the same haloes in different simulations, since they share the same initial density field phases. We show that the stellar mass varies systematically with resolution by over a factor of two, in a manner that depends on the final stellar mass. The evolution of the stellar populations in SIDM is largely identical to that of CDM, but in WDM early star formation is instead suppressed. The time at which LG haloes can begin to form stars through atomic cooling is delayed by $\sim$200~Myr in WDM models compared to CDM. 70~per~cent of WDM haloes of mass $>10^{8}M_{\odot}$ collapse early enough to form stars before $z=6$, compared to 90~per~cent of CDM and SIDM galaxies. It will be necessary to measure stellar ages of old populations to a precision of better than 100~Myr, and to address degeneracies with the redshift of reionization, in order to use these observables to distinguish between dark matter models.

Linking the properties of galaxies to the assembly history of their dark matter haloes is a central aim of galaxy evolution theory. This paper introduces a dimensionless parameter $s\in[0,1]$, the "tree entropy", to parametrise the geometry of a halo's entire mass assembly hierarchy, building on a generalisation of Shannon's information entropy. By construction, the minimum entropy ($s=0$) corresponds to smoothly assembled haloes without any mergers. In contrast, the highest entropy ($s=1$) represents haloes grown purely by equal-mass binary mergers. Using simulated merger trees extracted from the cosmological $N$-body simulation SURFS, we compute the natural distribution of $s$, a skewed bell curve peaking near $s=0.4$. This distribution exhibits weak dependences on halo mass $M$ and redshift $z$, which can be reduced to a single dependence on the relative peak height $\delta_{\rm c}/\sigma(M,z)$ in the matter perturbation field. By exploring the correlations between $s$ and global galaxy properties generated by the SHARK semi-analytic model, we find that $s$ contains a significant amount of information on the morphology of galaxies $-$ in fact more information than the spin, concentration and assembly time of the halo. Therefore, the tree entropy provides an information-rich link between galaxies and their dark matter haloes.

A non-parametric higher-order Jeans analysis method, GravSphere, was recently used to constrain the density profiles of Local Group dwarf galaxies, classifying them into those that are more likely to have an inner dark matter cusp and those that are likely to have a core (Read et al.). In this work we test this method using 31 simulated galaxies, comparable to Fornax, selected from the APOSTLE suite of cosmological hydrodynamics simulations, which include CDM and Self-Interacting Dark Matter (SIDM) cosmologies. We find that the mass profiles of the simulated dwarfs are often, but not always, well recovered by GravSphere. The less successful cases may be identified using a chi-squared diagnostic. Although the uncertainties are large in the inner regions, the inferred mass profiles are unbiased and exhibit smaller scatter than comparable Jeans methods. We find that GravSphere recovers the density profiles of simulated dwarfs below the half-light radius and down to the resolution limit of our simulations with better than 10+-30 per cent accuracy, making it a promising Jeans-based approach for modelling dark matter distributions in dwarf galaxies.

We examine the formation of dark matter (DM) cores in dwarf galaxies simulated with the EAGLE model of galaxy formation. As in earlier work, we find that the star formation (SF) gas density threshold ($\rho_{\rm th}$) plays a critical role. At low thresholds (LT), gas is unable to reach densities high enough to dominate the gravitational potential before being dispersed by feedback from supernovae. LT runs show little effect on the inner DM profile, even in systems with extended and bursty SF, two ingredients often cited as critical for core formation. For higher thresholds, gas is able to dominate the gravitational potential before being ejected by feedback. This can lead to a substantial reduction in the inner DM content, but only if the gas is gravitationally important over an extended period of time, allowing the halo to contract before gas removal. Rapid assembly and removal of gas in short SF bursts is less effective at altering the inner DM content. Subsequent gas accretion may draw DM back in and reform a cusp, unless SF is bursty enough to prevent it, preserving the core. Thus, for the EAGLE SF+feedback model, there is no simple relation between core formation and SF history, contrary to recent claims. The dependence of the inner DM content of dwarfs on $\rho_{\rm th}$ hinders robust predictions and the interpretation of observations. A simulation of a $(12 \rm \ Mpc)^3$ volume with high $\rho_{\rm th}$ results in dwarfs with sizeable cores over a limited halo mass range, but with insufficient variety in mass profiles to explain the observed diversity of dwarf galaxy rotation curves.

We study the impact of numerical parameters on the properties of cold dark matter haloes formed in collisionless cosmological simulations. We quantify convergence in the median spherically-averaged circular velocity profiles for haloes of widely varying particle number, as well as in the statistics of their structural scaling relations and mass functions. In agreement with prior work focused on single haloes, our results suggest that cosmological simulations yield robust halo properties for a wide range of gravitational softening parameters, $\epsilon$, provided: 1) $\epsilon$ is not larger than a "convergence radius", $r_{\rm conv}$, which is dictated by 2-body relaxation and determined by particle number, and 2) a sufficient number of timesteps are taken to accurately resolve particle orbits with short dynamical times. Provided these conditions are met, median circular velocity profiles converge to within $\approx 10$ per cent for radii beyond which the local 2-body relaxation timescale exceeds the Hubble time by a factor $\kappa\equiv t_{\rm relax}/t_{\rm H}\gt 0.177$, with better convergence attained for higher $\kappa$. We provide analytic estimates of $r_{\rm conv}$ that build on previous attempts in two ways: first, by highlighting its explicit (but weak) softening-dependence and, second, by providing a simpler criterion in which $r_{\rm conv}$ is determined entirely by the mean inter-particle spacing, $l$; for example, better than $10$ per cent convergence in circular velocity for $r\gt 0.05\,l$. We show how these analytic criteria can be used to assess convergence in structural scaling relations for dark matter haloes as a function of their mass or maximum circular speed.

The impact of 2-body scattering on the innermost density profiles of dark matter haloes is well established. We use a suite of cosmological simulations and idealised numerical experiments to show that 2-body scattering is exacerbated in situations where there are two species of unequal mass. This is a consequence of mass segregation and reflects a flow of kinetic energy from the more to less massive particles. This has important implications for the interpretation of galaxy sizes in cosmological hydrodynamic simulations, which nearly always model stars with less massive particles than are used for the dark matter. We compare idealised models as well as simulations from the EAGLE project that differ only in the mass resolution of the dark matter component, but keep sub-grid physics, baryonic mass resolution and gravitational force softening fixed. If the dark matter particle mass exceeds the mass of stellar particles, then galaxy sizes--quantified by their projected half-mass radii, ${\rm R_{50}}$--increase systematically with time until ${\rm R_{50}}$ exceeds a small fraction of the redshift-dependent mean inter-particle separation, $l$ (${\rm R_{50}}\geq 0.05\times l$). Our conclusions should also apply to simulations that adopt different hydrodynamic solvers, subgrid physics or adaptive softening, but in that case may need quantitative revision. Any simulation employing a stellar-to-dark matter particle mass ratio greater than unity will escalate spurious energy transfer from dark matter to baryons on small scales.

We use two high resolution N-body simulations, one assuming general relativity and the other the Hu-Sawicki form of $f\left(R\right)$ gravity with $\vert \bar{f}_{R} \vert = 10^{-6}$, to investigate the concentration--formation time relation of dark matter haloes. We stack haloes in logarithmically spaced mass bins to fit median density profiles and extract median formation times. At fixed mass, haloes in modified gravity are more concentrated than those in GR, especially at low masses and at low redshift, and do not follow the concentration--formation time relation seen in GR. We assess the sensitivity of the relation to how concentration and formation time are defined, as well as to the segregation of the halo population by the amount of gravitational screening. We find a clear difference between halo concentrations and assembly histories displayed in modified gravity and those in GR. Existing models for the mass--concentration--redshift relation that have gained success in cold and warm dark matter models require revision in $f\left(R\right)$ gravity.

Recently a number of studies have found a similarity between the passive fraction of central and satellite galaxies when controlled for both stellar and halo mass. These results suggest that the quenching processes that affect galaxies are largely agnostic to central/satellite status, which contradicts the traditional picture of increased satellite quenching via environmental processes such as stripping, strangulation and starvation. Here we explore this further using the Galaxy And Mass Assembly (GAMA) survey which extends to ~2dex lower in stellar mass than SDSS, is more complete for closely-separated galaxies (>95% compared to >70%), and identifies lower-halo-mass groups outside of the very local Universe (M$_{\mathrm{halo}}\sim10^{12}$M$_{\odot}$ at $0.1<z<0.2$). As far as possible we aim to replicate the selections, completeness corrections and central/satellite division of one of the previous studies but find clear differences between passive fractions of centrals and satellites. We also find that our passive fractions increase with both halo-to-satellite mass ratio and central-to-second rank mass ratio. This suggests that quenching is more efficient in satellites that are low-mass for their halo ($i.e$ at high halo-to-satellite mass ratio in comparison to low halo-to-satellite mass ratio) and are more likely to be passive in older groups - forming a consistent picture of environmental quenching of satellites. We then discuss potential explanations for the previously observed similarity, such as dependence on the group-finding method.

We apply the model relating halo concentration to formation history proposed by Ludlow et al. to merger trees generated using an algorithm based on excursion set theory. We find that while the model correctly predicts the median relation between halo concentration and mass, it underpredicts the scatter in concentration at fixed mass. Since the same model applied to N-body merger trees predicts the correct scatter, we postulate that the missing scatter is due to the lack of any environmental dependence in merger trees derived from excursion set theory. We show that a simple modification to the merger tree construction algorithm, which makes merger rates dependent on environment, can increase the scatter by the required amount, and simultaneously provide a qualitatively correct correlation between environment and formation epoch in the excursion set merger trees.

Dark matter particles may decay, emitting photons. Drawing on the EAGLE family of hydrodynamic simulations of galaxy formation -- including the APOSTLE and C-EAGLE simulations -- we assess the systematic uncertainties and scatter on the decay flux from different galaxy classes, from Milky Way satellites to galaxy clusters, and compare our results to studies of the 3.55~keV line. We demonstrate that previous detections and non-detections of this line are consistent with a dark matter interpretation. For example, in our simulations the width of the the dark matter decay line for Perseus-analogue galaxy clusters lies in the range 1300-1700~\kms. Therefore, the non-detection of the 3.55~keV line in the centre of the Perseus cluster by the {\it Hitomi} collaboration is consistent with detections by other instruments. We also consider trends with stellar and halo mass and evaluate the scatter in the expected fluxes arising from the anisotropic halo mass distribution and from object-to-object variations. We provide specific predictions for observations with {\it XMM-Newton} and with the planned X-ray telescopes {\it XRISM} and {\it ATHENA}. If future detections of unexplained X-ray lines match our predictions, including line widths, we will have strong evidence that we have discovered the dark matter.

We report on the first earth-scale quantum sensor network based on optical atomic clocks aimed at dark matter (DM) detection. Exploiting differences in the susceptibilities to the fine-structure constant of essential parts of an optical atomic clock, i.e. the cold atoms and the optical reference cavity, we can perform sensitive searches for dark matter signatures without the need of real-time comparisons of the clocks. We report a two orders of magnitude improvement in constraints on transient variations of the fine-structure constant, which considerably improves the detection limit for the standard model (SM) - DM coupling. We use Yb and Sr optical atomic clocks at four laboratories on three continents to search for both topological defect (TD) and massive scalar field candidates. No signal consistent with a dark-matter coupling is identified, leading to significantly improved constraints on the DM-SM couplings.

We study galaxy formation in sterile neutrino dark matter models that differ signifi- cantly from both cold and from `warm thermal relic' models. We use the EAGLE code to carry out hydrodynamic simulations of the evolution of pairs of galaxies chosen to resemble the Local Group, as part of the APOSTLE simulations project. We compare cold dark matter (CDM) with two sterile neutrino models with 7 keV mass: one, the warmest among all models of this mass (LA120) and the other, a relatively cold case (LA10). We show that the lower concentration of sterile neutrino subhalos compared to their CDM counterparts makes the inferred inner dark matter content of galaxies like Fornax (or Magellanic Clouds) less of an outlier in the sterile neutrino cosmologies. In terms of the galaxy number counts the LA10 simulations are indistinguishable from CDM when one takes into account halo-to-halo (or `simulation-to-simulation') scatter. In order for the LA120 model to match the number of Local Group dwarf galaxies, a higher fraction of low mass haloes is required to form galaxies than is predicted by the EAGLE simulations. As the census of the Local Group galaxies nears completion, this population may provide a strong discriminant between cold and warm dark matter models.

We study the vertical structure of polytropic, $P\propto \rho^\Gamma$, centrifugally-supported gaseous discs embedded in cold dark matter (CDM) halos. At fixed radius $R$, the shape of the vertical density profile depends only weakly on whether the disc is self-gravitating (SG) or not (NSG). The disc thickness, set by the midplane sound speed and circular velocity, $(c_s/V_c)R$, in the NSG case, and by the sound speed and surface density, $c_s^2/G\Sigma$, in SG discs, is smaller than either of these scales. SG discs are typically Toomre unstable, NSG discs are stable. Exponential discs in CDM halos with roughly flat circular velocity curves generally "flare" outwards. For the polytropic equation of state of the EAGLE simulations, discs whose mass and size match observational constraints are stable (NSG) for $M_d< 3\times 10^9\, M_\odot$ and unstable (SG) at higher masses, if fully gaseous. We test these analytic results using a set of idealized SPH simulations and find excellent agreement. Our results clarify the role of the gravitational softening on the thickness of simulated discs, and on the onset of radial instabilities. EAGLE low-mass discs are non-self-gravitating so the softening plays no role in their vertical structure. High-mass discs, on the other hand, are expected to be self-gravitating and unstable, and may be artificially thickened and stabilized unless gravity is well resolved. Simulations with spatial resolution high enough to not compromise the vertical structure of a disc also resolve the onset of their instabilities, but the converse is not true: resolving instabilities does not guarantee that the vertical structure is resolved.

We examine the origin of the mass discrepancy--radial acceleration relation (MDAR) of disk galaxies. This is a tight empirical correlation between the disk centripetal acceleration and that expected from the baryonic component. The MDAR holds for most radii probed by disk kinematic tracers, regardless of galaxy mass or surface brightness. The relation has two characteristic accelerations; $a_0$, above which all galaxies are baryon-dominated; and $a_{\rm min}$, an effective minimum aceleration probed by kinematic tracers in isolated galaxies. We use a simple model to show that these trends arise naturally in $\Lambda$CDM. This is because: (i) disk galaxies in $\Lambda$CDM form at the centre of dark matter haloes spanning a relatively narrow range of virial mass; (ii) cold dark matter halo acceleration profiles are self-similar and have a broad maximum at the centre, reaching values bracketed precisely by $a_{\rm min}$ and $a_0$ in that mass range; and (iii) halo mass and galaxy size scale relatively tightly with the baryonic mass of a galaxy in any successful $\Lambda$CDM galaxy formation model. Explaining the MDAR in $\Lambda$CDM does not require modifications to the cuspy inner mass profiles of dark haloes, although these may help to understand the detailed rotation curves of some dwarf galaxies and the origin of extreme outliers from the main relation. The MDAR is just a reflection of the self-similar nature of cold dark matter haloes and of the physical scales introduced by the galaxy formation process.

We simulate neutralino dark matter ($\chi$DM) haloes from their initial collapse, at $\sim$ earth mass, up to a few percent solar. Our results confirm that the density profiles of the first haloes are described by a $\sim r^{-1.5}$ power-law. As haloes grow in mass, their density profiles evolve significantly. In the central regions, they become shallower and reach on average $\sim r^{-1}$, the asymptotic form of an NFW profile. Using non-cosmological controlled simulations, we observe that temporal variations in the gravitational potential caused by major mergers lead to a shallowing of the inner profile. This transformation is more significant for shallower initial profiles and for a higher number of merging systems. Depending on the merger details, the resulting profiles can be shallower or steeper than NFW in their inner regions. Interestingly, mergers have a much weaker effect when the profile is given by a broken power-law with an inner slope of $-1$ (such as NFW or Hernquist profiles). This offers an explanation for the emergence of NFW-like profiles: after their initial collapse, $r^{-1.5}$ $\chi$DM haloes suffer copious major mergers, which progressively shallows the profile. Once an NFW-like profile is established, subsequent merging do not change the profile anymore. This suggests that halo profiles are not universal but rather a combination of (1) the physics of the formation of the microhaloes and (2) their early merger history -- both set by the properties of the dark matter particle -- as well as (3) the resilience of NFW-like profiles to perturbations.

We analyze the total and baryonic acceleration profiles of a set of well-resolved galaxies identified in the EAGLE suite of hydrodynamic simulations. Our runs start from the same initial conditions but adopt different prescriptions for unresolved stellar and AGN feedback, resulting in diverse populations of galaxies by the present day. Some of them reproduce observed galaxy scaling relations, while others do not. However, regardless of the feedback implementation, all of our galaxies follow closely a simple relationship between the total and baryonic acceleration profiles, consistent with recent observations of rotationally supported galaxies. The relation has small scatter: different feedback implementations -- which produce different galaxy populations -- mainly shift galaxies along the relation, rather than perpendicular to it. Furthermore, galaxies exhibit a characteristic acceleration, $g_{\dagger}$, above which baryons dominate the mass budget, as observed. These observations, consistent with simple modified Newtonian dynamics, can be accommodated within the standard cold dark matter paradigm.

We study the mass accretion histories (MAHs) and density profiles of dark matter halos using N-body simulations of self-similar gravitational clustering from scale-free power spectra, $P(k)\propto k^n$. We pay particular attention to the density profile curvature, which we characterize using the shape parameter, $\alpha$, of an Einasto profile. In agreement with previous findings our results suggest that, despite vast differences in their MAHs, the density profiles of virialized halos are remarkably alike. Nonetheless, clear departures from self-similarity are evident: for a given spectral index, $\alpha$ increases slightly but systematically with "peak height", $\nu\equiv\delta_{sc}/\sigma(M,z)$, regardless of mass or redshift. More importantly, however, the "$\alpha-\nu$" relation depends on $n$: the steeper the initial power spectrum, the more gradual the curvature of both the mean MAHs and mean density profiles. These results are consistent with previous findings connecting the shapes of halo mass profiles and MAHs and imply that dark matter halos are not structurally self-similar but, through the merger history, retain a memory of the linear density field from which they form.

We carry out controlled $N$-body simulations that follow the dynamical evolution of binary stars in the dark matter (DM) haloes of ultra-faint dwarf spheroidals (dSphs). We find that wide binaries with semi-major axes $a\gtrsim a_t$ tend to be quickly disrupted by the tidal field of the halo. In smooth potentials the truncation scale, $a_t$, is mainly governed by (i) the mass enclosed within the dwarf half-light radius ($R_h$) and (ii) the slope of the DM halo profile at $R\approx R_h$, and is largely independent of the initial eccentricity distribution of the binary systems and the anisotropy of the stellar orbits about the galactic potential. For the reported velocity dispersion and half-light radius of Segue I, the closest ultra-faint, our models predict $a_t$ values that are a factor 2--3 smaller in cuspy haloes than in any of the cored models considered here. Using mock observations of Segue I we show that measuring the projected two-point correlation function of stellar pairs with sub-arcsecond resolution may provide a useful tool to constrain the amount and distribution of DM in the smallest and most DM-dominated galaxies.

We use a suite of cosmological simulations to study the mass-concentration-redshift relation, $c({\rm M},z)$, of dark matter halos. Our simulations include standard $\Lambda$-cold dark matter (CDM) models, and additional runs with truncated power spectra, consistent with a thermal warm dark matter (WDM) scenario. We find that the mass profiles of CDM and WDM halos are self-similar and well approximated by the Einasto profile. The $c({\rm M},z)$ relation of CDM halos is monotonic: concentrations decrease with increasing virial mass at fixed redshift, and decrease with increasing redshift at fixed mass. The mass accretion histories (MAHs) of CDM halos are also scale-free, and can be used to infer concentrations directly. These results do not apply to WDM halos: their MAHs are not scale-free because of the characteristic scale imposed by the power-spectrum suppression. Further, the WDM $c({\rm M},z)$ relation is non-monotonic: concentrations peak at a mass scale dictated by the truncation scale, and decrease at higher and lower masses. We show that the assembly history of a halo can still be used to infer its concentration, provided that the total mass of its progenitors is considered (the "collapsed mass history"; CMH), rather than just that of its main ancestor. This exploits the scale-free nature of CMHs to derive a simple scaling that reproduces the mass-concentration-redshift relation of both CDM and WDM halos over a vast range of halo masses and redshifts. Our model therefore provides a robust account of the mass, redshift, cosmology and power spectrum dependence of dark matter halo concentrations.

We use high resolution N-Body simulations to study the concentration and spin parameters of dark matter haloes in the mass range $10^8\, {\rm M}_{\odot}\, h^{-1} < {\rm M} < 10^{11}\, {\rm M}_{\odot}\, h^{-1}$ and redshifts $5{<}z{<}10$, corresponding to the haloes of galaxies thought to be responsible for reionization. We build a sub-sample of equilibrium haloes and contrast their properties to the full population that also includes unrelaxed systems. Concentrations are calculated by fitting both NFW and Einasto profiles to the spherically-averaged density profiles of individual haloes. After removing haloes that are out-of-equilibrium, we find a $z{>}5$ concentration$-$mass ($c(M)$) relation that is almost flat and well described by a simple power-law for both NFW and Einasto fits. The intrinsic scatter around the mean relation is $\Delta c_{\rm{vir}}{\sim1}$ (or 20 per cent) at $z=5$. We also find that the analytic model proposed by Ludlow et al. reproduces the mass and redshift-dependence of halo concentrations. Our best-fit Einasto shape parameter, $\alpha$, depends on peak height, $\nu$, in a manner that is accurately described by $\alpha {=}0.0070\nu^2{+}0.1839$. The distribution of the spin parameter, $\lambda$, has a weak dependence on equilibrium state; $\lambda$ peaks at roughly ${\sim}0.033$ for our relaxed sample, and at ${\sim}0.04$ for the full population. The spin--virial mass relation has a mild negative correlation at high redshift.

We introduce a sub-grid model for the non-equilibrium abundance of molecular hydrogen in cosmological simulations of galaxy formation. We improve upon previous work by accounting for the unresolved structure of molecular clouds in a phenomenological way which combines both observational and numerical results on the properties of the turbulent interstellar medium. We apply the model to a cosmological simulation of the formation of a Milky-Way-sized galaxy at z=2, and compare the results to those obtained using other popular prescriptions that compute the equilibrium abundance of H2. In these runs we introduce an explicit link between star formation and the local H2 abundance, and perform an additional simulation in which star formation is linked directly to the density of cold gas. In better agreement with observations, we find that the simulated galaxy produces less stars and harbors a larger gas reservoir when star formation is regulated by molecular hydrogen. In this case, the galaxy is composed of a younger stellar population as early star formation is inhibited in small, metal poor dark-matter haloes which cannot efficiently produce H2. The number of luminous satellites orbiting within the virial radius of the galaxy at z=2 is reduced by 10-30 per cent in models with H2-regulated star formation.

We use two cosmological simulations of structure formation in the LambdaCDM scenario to study the evolutionary histories of dark-matter haloes and to characterize the Lagrangian regions from which they form. We focus on haloes identified at redshift z_id=0 and show that the classic ellipsoidal collapse model systematically overestimates their collapse times. If one imposes that halo collapse takes place at z_id, this model requires starting from a significantly lower linear density contrast than what is measured in the simulations at the locations of halo formation. We attempt to explain this discrepancy by testing two key assumptions of the model. First, we show that the tides felt by collapsing haloes due to the surrounding large-scale structure evolve non-linearly. Although this effect becomes increasingly important for low-mass haloes, accounting for it in the ellipsoidal collapse model only marginally improves the agreement with N-body simulations. Second, we track the time evolution of the physical volume occupied by forming haloes and show that, after turnaround, it generally stabilizes at a well-defined redshift, z_c>z_id, contrary to the basic assumption of extended Press-Schechter theory based on excursion sets. We discuss the implications of this result for understanding the origin of the mass-dependence and scatter in the linear threshold for halo formation. Finally, we show that, when tuned for collapse at z_c, a modified version of the ellipsoidal collapse model that also accounts for the triaxial nature of protohaloes predicts their linear density contrast in an unbiased way.

In the excursion set approach to structure formation initially spherical regions of the linear density field collapse to form haloes of mass $M$ at redshift $z_{\rm id}$ if their linearly extrapolated density contrast, averaged on that scale, exceeds some critical threshold, $\delta_{\rm c}(z_{\rm id})$. The value of $\delta_{\rm c}(z_{\rm id})$ is often calculated from the spherical or ellipsoidal collapse model, which provide well-defined predictions given auxiliary properties of the tidal field at a given location. We use two cosmological simulations of structure growth in a $\Lambda$ cold dark matter scenario to quantify $\delta_{\rm c}(z_{\rm id})$, its dependence on the surrounding tidal field, as well as on the shapes of the Lagrangian regions that collapse to form haloes at $z_{\rm id}$. Our results indicate that the ellipsoidal collapse model provides an accurate description of the mean dependence of $\delta_{\rm c}(z_{\rm id})$ on both the strength of the tidal field and on halo mass. However, for a given $z_{\rm id}$, $\delta_{\rm c}(z_{\rm id})$ depends strongly on the halo's characteristic formation redshift: the earlier a halo forms, the higher its initial density contrast. Surprisingly, the majority of haloes forming $today$ fall below the ellipsoidal collapse barrier, contradicting the model predictions. We trace the origin of this effect to the non-spherical shapes of Lagrangian haloes, which arise naturally due to the asymmetry of the linear tidal field. We show that a modified collapse model, that accounts for the triaxial shape of protohaloes, provides a more accurate description of the measured minimum overdensities of recently collapsed objects.

We use a high-resolution hydrodynamic simulation that tracks the non-equilibrium abundance of molecular hydrogen within a massive high-redshift galaxy to produce mock Atacama Large Millimeter Array (ALMA) maps of the fine-structure lines of atomic carbon, CI 1-0 and CI 2-1. Inspired by recent observational and theoretical work, we assume that CI is thoroughly mixed within giant molecular clouds and demonstrate that its emission is an excellent proxy for H2. Nearly all of the H2 associated with the galaxy can be detected at redshifts z<4 using a compact interferometric configuration with a large synthesized beam (that does not resolve the target galaxy) in less than 4 h of integration time. Low-resolution imaging of the \CI lines (in which the target galaxy is resolved into three to four beams) will detect ~80 per cent of the H2 in less than 12 h of aperture synthesis. In this case, the resulting data cube also provides the crucial information necessary for determining the dynamical state of the galaxy. We conclude that ALMA observations of the CI 1-0 and 2-1 emission are well-suited for extending the interval of cosmic look-back time over which the H2 distributions, the dynamical masses, and the Tully-Fisher relation of galaxies can be robustly probed.

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.

We use the Millennium Simulation series to investigate the mass and redshift dependence of the concentration of equilibrium cold dark matter (CDM) halos. We extend earlier work on the relation between halo mass profiles and assembly histories to show how the latter may be used to predict concentrations for halos of all masses and at any redshift. Our results clarify the link between concentration and the ``collapse redshift'' of a halo as well as why concentration depends on mass and redshift solely through the dimensionless ``peak height'' mass parameter, $\nu(M,z)=\delta_{\rm crit}(z)/\sigma(M,z)$. We combine these results with analytic mass accretion histories to extrapolate the $c(M,z)$ relations to mass regimes difficult to reach through direct simulation. Our model predicts that, at given $z$, $c(M)$ should deviate systematically from a simple power law at high masses, where concentrations approach a constant value, and at low masses, where concentrations are substantially lower than expected from extrapolating published empirical fits. This correction may reduce the expected self-annihilation boost factor from substructure by about one order of magnitude. The model also reproduces the $c(M,z)$ dependence on cosmological parameters reported in earlier work, and thus provides a simple and robust account of the relation between cosmology and the mass-concentration-redshift relation of CDM halos.

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.

We use high resolution N-body simulations to investigate the Lagrangian bias of cold dark matter haloes within the LCDM cosmology. Our analysis focuses on "proto-haloes", which we identify in the simulation initial conditions with the subsets of particles belonging to individual redshift-zero haloes. We then calculate the number-density and velocity-divergence fields of proto-haloes and estimate their auto spectral densities. We also measure the corresponding cross spectral densities with the linear matter distribution. We use our results to test a Lagrangian-bias model presented by Desjacques and Sheth which is based on the assumption that haloes form out of local density maxima of a specific height. Our comparison validates the predicted functional form for the scale-dependence of the bias for both the density and velocity fields. We also show that the bias coefficients are accurately predicted for the velocity divergence. On the contrary, the theoretical values for the density bias parameters do not accurately match the numerical results as a function of halo mass. This is likely due to the simplistic assumptions that relate virialized haloes to density peaks of a given height in the model. We also detect appreciable stochasticity for the Lagrangian density bias, even on very large scales. These are not included in the model at leading order but correspond to higher order corrections.

We use two cosmological simulations of structure formation to study the conditions under which dark matter haloes emerge from the linear density field. Our analysis focuses on matching sites of halo collapse to local density maxima, or "peaks", in the initial conditions of the simulations and provides a crucial test of the central ansatz of the peaks formalism. By identifying peaks on a variety of smoothed, linearly extrapolated density fields we demonstrate that as many as ~70% of well-resolved dark matter haloes form preferentially near peaks whose characteristic masses are similar to that of the halo, with more massive haloes showing a stronger tendency to reside near peaks initially. We identify a small but significant fraction of haloes that appear to evolve from peaks of substantially lower mass than that of the halo itself. We refer to these as "peakless haloes" for convenience. By contrasting directly the properties of these objects with the bulk of the proto-halo population we find two clear differences: 1) their initial shapes are significantly flatter and more elongated than the predominantly triaxial majority, and 2) they are, on average, more strongly compressed by tidal forces associated with their surrounding large scale structure. Using the two-point correlation function we show that peakless haloes tend to emerge from highly clustered regions of the initial density field implying that, at fixed mass, the accretion geometry and mass accretion histories of haloes in highly clustered environments differ significantly from those in the field. This may have important implications for understanding the origin of the halo assembly bias, of galaxy properties in dense environments and how environment affects the morphological transformation of galaxies near groups and rich galaxy clusters.

Cosmological N-body simulations indicate that the spherically-averaged density profiles of cold dark matter halos are accurately described by Einasto profiles, where the logarithmic slope is a power-law of adjustable exponent, \gamma =dln\rho /dlnr ~ r^\alpha $. The pseudo-phase-space density (PPSD) profiles of CDM halos also show remarkable regularity, and are well approximated by simple power laws, Q(r)=\rho /\sigma ^3 ~ r^-\chi . We show that this is expected from dynamical equilibrium considerations, since Jeans' equations predict that the pseudo-phase-space density profiles of Einasto halos should resemble power laws over a wide range of radii. For the values of \alpha typical of CDM halos, the inner Q profiles of equilibrium halos deviate significantly from a power law only very close to the center, and simulations of extremely high-resolution would be needed to detect such deviations unambiguously. We use an ensemble of halos drawn from the Millennium-II simulation to study which of these two alternatives describe best the mass profile of CDM halos. Our analysis indicates that at the resolution of the best available simulations, both Einasto and power-law PPSD profiles (with adjustable exponents \alpha and \chi, respectively) provide equally acceptable fits to the simulations. A full account of the structure of CDM halos requires understanding how the shape parameters that characterize departures from self-similarity, like \alpha or \chi, are determined by evolutionary history, environment or initial conditions.

We use the Aquarius simulation series to study the imprint of assembly history on the structure of Galaxy-mass cold dark matter halos. Our results confirm earlier work regarding the influence of mergers on the mass density profile and the inside-out growth of halos. The inner regions that contain the visible galaxies are stable since early times and are significantly affected only by major mergers. Particles accreted diffusely or in minor mergers are found predominantly in the outskirts of halos. Our analysis reveals trends that run counter to current perceptions of hierarchical halo assembly. For example, major mergers (i.e. those with progenitor mass ratios greater than 1:10) contribute little to the total mass growth of a halo, on average less than 20 per cent for our six Aquarius halos. The bulk is contributed roughly equally by minor mergers and by "diffuse" material which is not resolved into individual objects. This is consistent with modeling based on excursion-set theory which suggests that about half of this diffuse material should not be part of a halo of any scale. Interestingly, the simulations themselves suggest that a significantly fraction is not truly diffuse, since it was ejected from earlier halos by mergers prior to their joining the main system. The Aquarius simulations resolve halos to much lower mass scales than are expected to retain gas or form stars. These results thus confirm that most of the baryons from which visible galaxies form are accreted diffusely, rather than through mergers, and they suggest that only relatively rare major mergers will affect galaxy structure at later times.

We use N-body simulations to investigate the radial dependence of the density and velocity dispersion in cold dark matter (CDM) halos. In particular, we explore how closely Q rho/sigma^3, a surrogate measure of the phase-space density, follows a power-law in radius. Our study extends earlier work by considering, in addition to spherically-averaged profiles, local Q-estimates for individual particles, Q_i; profiles based on the ellipsoidal radius dictated by the triaxial structure of the halo, Q_i(r'); and by carefully removing substructures in order to focus on the profile of the smooth halo, Q^s. The resulting Q_i^s(r') profiles follow closely a power law near the center, but show a clear upturn from this trend near the virial radius, r_{200}. The location and magnitude of the deviations are in excellent agreement with the predictions from Bertschinger's spherical secondary-infall similarity solution. In this model, Q \propto r^{-1.875} in the inner, virialized regions, but departures from a power-law occur near r_{200} because of the proximity of this radius to the location of the first shell crossing - the shock radius in the case of a collisional fluid. Particles there have not yet fully virialized, and so Q departs from the inner power-law profile. Our results imply that the power-law nature of $Q$ profiles only applies to the inner regions and cannot be used to predict accurately the structure of CDM halos beyond their characteristic scale radius.

We study the mass, velocity dispersion, and anisotropy profiles of $\Lambda$CDM halos using a suite of N-body simulations of unprecedented numerical resolution (the {\it Aquarius Project}). Our analysis confirms a number of results claimed by earlier work, and clarifies a few issues where conflicting claims may be found in the recent literature. The spherically-averaged density profile becomes progressively shallower inwards and, at the innermost resolved radius, the logarithmic slope is $\gamma \equiv -$d$\ln\rho/$d$\ln r \simlt 1$. Asymptotic inner slopes as steep as the recently claimed $\rho \propto r^{-1.2}$ are clearly ruled out. The radial dependence of $\gamma$ is well approximated by a power-law, $\gamma \propto r^{\alpha}$ (the Einasto profile). The shape parameter, $\alpha$, varies slightly but significantly from halo to halo, implying that the mass profiles of $\Lambda$CDM halos are not strictly universal: different halos cannot, in general, be rescaled to look identical. Departures from similarity are also seen in velocity dispersion profiles and correlate with those in density profiles so as to preserve a power-law form for the spherically averaged pseudo-phase-space density, $\rho/\sigma^3\propto r^{-1.875}$. Our conclusions are reliable down to radii below 0.4% of the virial radius, providing well-defined predictions for halo structure when baryonic effects are neglected, and thus an instructive theoretical template against which the modifications induced by the baryonic components of real galaxies can be judged.

We study predictions for dark matter phase-space structure near the Sun based on high-resolution simulations of six galaxy halos taken from the Aquarius Project. The local DM density distribution is predicted to be remarkably smooth; the density at the Sun differs from the mean over a best-fit ellipsoidal equidensity contour by less than 15% at the 99.9% confidence level. The local velocity distribution is also very smooth, but it differs systematically from a (multivariate) Gaussian distribution. This is not due to the presence of individual clumps or streams, but to broad features in the velocity modulus and energy distributions that are stable both in space and time and reflect the detailed assembly history of each halo. These features have a significant impact on the signals predicted for WIMP and axion searches. For example, WIMP recoil rates can deviate by ~10% from those expected from the best-fit multivariate Gaussian models. The axion spectra in our simulations typically peak at lower frequencies than in the case of multivariate Gaussian velocity distributions. Also in this case, the spectra show significant imprints of the formation of the halo. This implies that once direct DM detection has become routine, features in the detector signal will allow us to study the dark matter assembly history of the Milky Way. A new field, "dark matter astronomy", will then emerge.

(Abridged) We use cosmological N-body simulations to study the properties of substructure halos in galaxy-sized dark matter halos. We extend prior work on the subject by considering the whole population of subhalos physically associated with the main system. These are defined as subhalos that have at some time in the past been within the virial radius of the halo's main progenitor and that have survived as self-bound entities to $z=0$. We find that this population extends beyond {\it three times} the virial radius, and contains objects on extreme orbits. We trace the origin of these unorthodox orbits to the tidal dissociation of bound groups of subhalos, which results in the ejection of some subhalos along tidal streams. Ejected subhalos are primarily low-mass systems, leading to mass-dependent biases in their spatial distribution and kinematics: the lower the subhalo mass at accretion time, the less centrally concentrated and kinematically hotter their descendant population. The bias is strongest amongst the most massive subhalos, but disappears at the low-mass end. Our findings imply that subhalos identified within the virial radius represent an incomplete census of the substructure physically related to a halo: only about {\it one half} of all associated subhalos are found today within the virial radius of a halo. These results may explain the age dependence of the clustering of low-mass halos, and has implications for (i) the interpretation of the structural parameters and assembly histories of halos neighboring massive systems; (ii) the existence of low-mass dynamical outliers in the Local Group; and (iii) the presence of evidence for evolutionary effects well outside the traditional virial boundary of a galaxy system.

We have performed the largest ever particle simulation of a Milky Way-sized dark matter halo, and present the most comprehensive convergence study for an individual dark matter halo carried out thus far. We have also simulated a sample of 6 ultra-highly resolved Milky-way sized halos, allowing us to estimate the halo-to-halo scatter in substructure statistics. In our largest simulation, we resolve nearly 300,000 gravitationally bound subhalos within the virialized region of the halo. Simulations of the same object differing in mass resolution by factors up to 1800 accurately reproduce the largest subhalos with the same mass, maximum circular velocity and position, and yield good convergence for the abundance and internal properties of dark matter substructures. We detect up to four generations of subhalos within subhalos, but contrary to recent claims, we find less substructure in subhalos than in the main halo when regions of equal mean overdensity are compared. The overall substructure mass fraction is much lower in subhalos than in the main halo. Extrapolating the main halo's subhalo mass spectrum down to an Earth mass, we predict the mass fraction in substructure to be well below 3% within 100 kpc, and to be below 0.1% within the Solar Circle. The inner density profiles of subhalos show no sign of converging to a fixed asymptotic slope and are well fit by gently curving profiles of Einasto form. The mean concentrations of isolated halos are accurately described by the fitting formula of Neto et al. down to maximum circular velocities of 1.5 km/s, an extrapolation over some 5 orders of magnitude in mass. However, at equal maximum circular velocity, subhalos are more concentrated than field halos, with a characteristic density that is typically ~2.6 times larger and increases towards the halo centre.

Dark matter is the dominant form of matter in the universe, but its nature is unknown. It is plausibly an elementary particle, perhaps the lightest supersymmetric partner of known particle species. In this case, annihilation of dark matter in the halo of the Milky Way should produce gamma-rays at a level which may soon be observable. Previous work has argued that the annihilation signal will be dominated by emission from very small clumps (perhaps smaller even than the Earth) which would be most easily detected where they cluster together in the dark matter halos of dwarf satellite galaxies. Here we show, using the largest ever simulation of the formation of a galactic halo, that such small-scale structure will, in fact, have a negligible impact on dark matter detectability. Rather, the dominant and likely most easily detectable signal will be produced by diffuse dark matter in the main halo of the Milky Way. If the main halo is strongly detected, then small dark matter clumps should also be visible, but may well contain no stars, thereby confirming a key prediction of the Cold Dark Matter (CDM) model.

We present spectroscopic observations of the AndXII dwarf spheroidal galaxy using DEIMOS/Keck-II, showing it to be moving rapidly through the Local Group (-556 km/s heliocentric velocity, -281 km/s relative to Andromeda from the MW), falling into the Local Group from ~115 kpc beyond Andromeda's nucleus. AndXII therefore represents a dwarf galaxy plausibly falling into the Local Group for the first time, and never having experienced a dense galactic environment. From Green Bank Telescope observations, a limit on the H{I} gas mass of <3000 Msun suggests that AndXII's gas could have been removed prior to experiencing the tides of the Local Group galaxies. Orbit models suggest the dwarf is close to the escape velocity of M31 for published mass models. AndXII is our best direct evidence for the late infall of satellite galaxies, a prediction of cosmological simulations.

We explore similarities in the luminosity distribution of early type galaxies and the mass profiles of LCDM halos. The spatial structure of these systems may be accurately described by a simple law where the logarithmic slope of the projected density is a power law of radius; the Sersic law. We show that this law provides a significantly better fit than a three-parameter generalization of the NFW profile and derive the best-fitting Sersic parameters for a set of high-resolution LCDM halos spanning a wide range in mass. The mean Sersic n values are 3.0 for dwarf- and galaxy-sized halos and 2.4 for cluster-sized halos, similar to the values that characterize luminous elliptical galaxies. We discuss possible reasons why the same law should describe dark and luminous systems that span a range of over seven decades in mass.