Maximum Circular Velocity Research Articles

The Copernicus Complexio: a high-resolution view of the small-scale Universe

We introduce Copernicus Complexio (COCO), a high-resolution cosmological N-body simulation of structure formation in the $\Lambda{\rm CDM}{}$ model. COCO follows an approximately spherical region of radius $\sim 17.4h^{-1}\,{\rm Mpc}$ embedded in a much larger periodic cube that is followed at lower resolution. The high resolution volume has a particle mass of $1.135\times10^5h^{-1}{\rm M}_{\odot}$ (60 times higher than the Millennium-II simulation). COCO gives the dark matter halo mass function over eight orders of magnitude in halo mass; it forms $\sim 60$ haloes of galactic size, each resolved with about 10 million particles. We confirm the power-law character of the subhalo mass function, $\bar{N}(>\mu)\propto\mu^{-s}$, down to a reduced subhalo mass $M_{sub}/M_{200}\equiv\mu=10^{-6}$, with a best-fit power-law index, $s=0.94$, for hosts of mass $\langle M_{200}\rangle=10^{12}h^{-1}{\rm M}_{\odot}$. The concentration-mass relation of COCO haloes deviates from a single power law for masses $M_{200}<\textrm{a few}\times 10^{8}h^{-1}{\rm M}_{\odot}$, where it flattens, in agreement with results by Sanchez-Conde et al. The host mass invariance of the reduced maximum circular velocity function of subhaloes, $\nu\equiv V_{max}/V_{200}$, hinted at in previous simulations, is clearly demonstrated over five orders of magnitude in host mass. Similarly, we find that the average, normalised radial distribution of subhaloes is approximately universal (i.e. independent of subhalo mass), as previously suggested by the Aquarius simulations of individual haloes. Finally, we find that at fixed physical subhalo size, subhaloes in lower mass hosts typically have lower central densities than those in higher mass hosts.

The APOSTLE simulations: solutions to the Local Group's cosmic puzzles

The Local Group of galaxies offer some of the most discriminating tests of models of cosmic structure formation. For example, observations of the Milky Way (MW) and Andromeda satellite populations appear to be in disagreement with N-body simulations of the "Lambda Cold Dark Matter" ({\Lambda}CDM) model: there are far fewer satellite galaxies than substructures in cold dark matter halos (the "missing satellites" problem); dwarf galaxies seem to avoid the most massive substructures (the "too-big-to-fail" problem); and the brightest satellites appear to orbit their host galaxies on a thin plane (the "planes of satellites" problem). Here we present results from APOSTLE (A Project Of Simulating The Local Environment), a suite of cosmological hydrodynamic simulations of twelve volumes selected to match the kinematics of the Local Group (LG) members. Applying the Eagle code to the LG environment, we find that our simulations match the observed abundance of LG galaxies, including the satellite galaxies of the MW and Andromeda. Due to changes to the structure of halos and the evolution in the LG environment, the simulations reproduce the observed relation between stellar mass and velocity dispersion of individual dwarf spheroidal galaxies without necessitating the formation of cores in their dark matter profiles. Satellite systems form with a range of spatial anisotropies, including one similar to that of the MW, confirming that such a configuration is not unexpected in {\Lambda}CDM. Finally, based on the observed velocity dispersion, size, and stellar mass, we provide new estimates of the maximum circular velocity for the halos of nine MW dwarf spheroidals.

Structure formation in a mixed dark matter model with decaying sterile neutrino: the 3.5 keV X-ray line and the Galactic substructure

We perform a set of cosmological simulations of structure formation in a mixed dark matter (MDM) model. Our model is motivated by the recently identified 3.5 keV X-ray line, which can be explained by the decay of non-resonantly produced sterile neutrinos accounting for 20–60% of the dark matter in the Universe. These non-resonantly produced sterile neutrinos have a sizable free-streaming length and hence behave effectively as warm dark matter (WDM). Assuming the rest of dark matter is composed of some cold dark matter (CDM) particles, we follow the coevolution of a mixed WDM plus CDM cosmology. Specifically, we consider the models with the warm component fraction of rwarm=0.25 and 0.50. Our MDM models predict that the comoving Jeans length at the matter-radiation equality is close to that of the thermally produced warm dark matter model with particle mass mWDM=2.4 keV, but the suppression in the fluctuation power spectrum is weaker. We perform large N-body simulations to study the structure of non-linear dark halos in the MDM models. The abundance of substructure is significantly reduced in the MDM models, and hence the so-called small-scale crisis is mitigated. The cumulative maximum circular velocity function (CVF) of at least one halo in the MDM models is in good agreement with the CVFs of the observed satellites in the Milky Way and the Andromeda galaxy. We argue that the MDM models open an interesting possibility to reconcile the reported 3.5 keV line and the internal structure of galaxies.

The chosen few: the low-mass haloes that host faint galaxies

Since reionization prevents star formation in most halos below 3 x 10^9 solar masses, dwarf galaxies only populate a fraction of existing dark matter halos. We use hydrodynamic cosmological simulations of the Local Group to study the discriminating factors for galaxy formation in the early Universe and connect them to the present-day properties of galaxies and halos. A combination of selection effects related to reionization, and the subsequent evolution of halos in different environments, introduces strong biases between the population of halos that host dwarf galaxies, and the total halo population. Halos that host galaxies formed earlier and are more concentrated. In addition, halos more affected by tidal stripping are more likely to host a galaxy for a given mass or maximum circular velocity, vmax, today. Consequently, satellite halos are populated more frequently than field halos, and satellite halos of 10^8 - 10^9 solar masses or vmax of 12 - 20 km/s, similar to the Local Group dwarf spheroidals, have experienced a greater than average reduction in both mass and vmax after infall. They are on closer, more radial orbits with higher infall velocities and earlier infall times. Together, these effects make dwarf galaxies highly biased tracers of the underlying dark matter distribution.

A matter of measurement: rotation velocities and the velocity function of dwarf galaxies

The velocity function derived from large scale surveys can be compared with the predictions of LCDM cosmology, by matching the measured rotation velocities Vrot of galaxies to the maximum circular velocity of dark matter (DM) halos Vmax. For Vrot<50km/s, a major discrepancy arises between the observed and LCDM velocity functions. However, the manner in which different observational measures of Vrot are associated with Vmax is not straight forward in dwarf galaxies. We instead relate galaxies to DM halos using the empirical baryon- mass to halo-mass relation, and show that different observational measures of Vrot result in very different velocity functions. We show how the W50 velocity function, i.e. using the HI profile line width at 50% of peak HI flux to measure Vrot, can be reconciled with a LCDM cosmology. Our semi-empirical methodology allows us to determine the region of rotation curves that are probed by HI measurements (RHI), and shows that the Vrot of dwarfs are generally measured at a fraction of Rmax, explaining their tendency to have rising rotation curves. We provide fitting formulae for relating RHI and Reff (the effective radius) to the virial radius of DM halos. To continue to use velocity functions as a probe of LCDM cosmology, it is necessary to be precise about how the different measures of rotation velocity are probing the mass of the DM halos, dropping the assumption that any measure of rotational velocity can be equally used as a proxy for Vmax.

On the segregation of dark matter substructure

We present the first comprehensive analysis of the segregation of dark matter subhaloes in their host haloes. Using numerical simulations, we examine the segregation of twelve different subhalo properties with respect to both orbital energy and halo-centric radius (in real space as well as in projection). Subhaloes are strongly segregated by accretion redshift, which is an outcome of the inside-out assembly of their host haloes. Since subhaloes that were accreted earlier have experienced more tidal stripping, subhaloes that have lost a larger fraction of their mass at infall are on more bound orbits. Subhaloes are also strongly segregated in their masses and maximum circular velocities at accretion. We demonstrate that part of this segregation is already imprinted in the infall conditions. For massive subhaloes it is subsequently boosted by dynamical friction, but only during their first radial orbit. The impact of these two effects is counterbalanced, though, by the fact that subhaloes with larger accretion masses are accreted later. Because of tidal stripping, subhaloes reveal little to no segregation by present-day mass or maximum circular velocity, while the corresponding torques cause subhaloes on more bound orbits to have smaller spin. There is a weak tendency for subhaloes that formed earlier to be segregated towards the center of their host halo, which is an indirect consequence of the fact that (sub)halo formation time is correlated with other, strongly segregated properties. We discuss the implications of our results for the segregation of satellite galaxies in galaxy groups and clusters.

Comprehensive assessment of the too big to fail problem

We use a semi-analytical model for the substructure of dark matter haloes to assess the too-big-to-fail (TBTF) problem. The model accurately reproduces the average subhalo mass and velocity functions, as well as their halo-to-halo variance, in N-body simulations. We construct thousands of realizations of Milky Way (MW) size host haloes, allowing us to investigate the TBTF problem with unprecedented statistical power. We examine the dependence on host halo mass and cosmology, and explicitly demonstrate that a reliable assessment of TBTF requires large samples of hundreds of host haloes. We argue that previous statistics used to address TBTF suffer from the look-elsewhere effect and/or disregard certain aspects of the data on the MW satellite population. We devise a new statistic that is not hampered by these shortcomings, and, using only data on the 9 known MW satellite galaxies with $V_{\rm max}>15{\rm kms}^{-1}$, demonstrate that $1.4^{+3.3}_{-1.1}\%$ of MW-size host haloes have a subhalo population in statistical agreement with that of the MW. However, when using data on the MW satellite galaxies down to $V_{\rm max}=8{\rm kms}^{-1}$, this MW consistent fraction plummets to $<5\times10^{-4}$ (at 68% CL). Hence, if it turns out that the inventory of MW satellite galaxies is complete down to 8km/s, then the maximum circular velocities of MW satellites are utterly inconsistent with $\Lambda$CDM predictions, unless baryonic effects can drastically increase the spread in $V_{\rm max}$ values of satellite galaxies compared to that of their subhaloes.

A STUDY IN BLUE: THE BARYON CONTENT OF ISOLATED LOW-MASS GALAXIES

We study the baryon content of low mass galaxies selected from the Sloan Digital Sky Survey (SDSS DR8), focusing on galaxies in isolated environments where the complicating physics of galaxy-galaxy interactions are minimized. We measure neutral hydrogen (HI) gas masses and line-widths for 148 isolated galaxies with stellar mass between $10^7$ and $10^{9.5} M_{\odot}$. We compare isolated low mass galaxies to more massive galaxies and galaxies in denser environments by remeasuring HI emission lines from the Arecibo Legacy Fast ALFA (ALFALFA) survey 40% data release. All isolated low mass galaxies either have large atomic gas fractions or large atomic gas fractions cannot be ruled out via their upper limits. We measure a median atomic gas fraction of $f_{\rm gas} = 0.82 \pm 0.13$ for our isolated low mass sample with no systems below 0.30. At all stellar masses, the correlations between galaxy radius, baryonic mass and velocity width are not significantly affected by environment. Finally, we estimate a median baryon to total dynamical mass fraction of $f_{\rm baryon, disk} = 0.15 \pm 0.18$. We also estimate two different median baryon to halo mass fractions using the results of semi-analytic models $(f_{\rm baryon, halo} = 0.04 \pm 0.06)$ and abundance matching $(f_{\rm baryon, halo} = 0.04 \pm 0.02)$. Baryon fractions estimated directly using HI observations appear independent of environment and maximum circular velocity, while baryon fractions estimated using abundance matching show a significant depletion of baryons at low maximum circular velocities.

The unexpected diversity of dwarf galaxy rotation curves

We examine the circular velocity profiles of galaxies in {\Lambda}CDM cosmological hydrodynamical simulations from the EAGLE and LOCAL GROUPS projects and compare them with a compilation of observed rotation curves of galaxies spanning a wide range in mass. The shape of the circular velocity profiles of simulated galaxies varies systematically as a function of galaxy mass, but shows remarkably little variation at fixed maximum circular velocity. This is especially true for low-mass dark matter-dominated systems, reflecting the expected similarity of the underlying cold dark matter haloes. This is at odds with observed dwarf galaxies, which show a large diversity of rotation curve shapes, even at fixed maximum rotation speed. Some dwarfs have rotation curves that agree well with simulations, others do not. The latter are systems where the inferred mass enclosed in the inner regions is much lower than expected for cold dark matter haloes and include many galaxies where previous work claims the presence of a constant density "core". The "cusp vs core" issue is thus better characterized as an "inner mass deficit" problem than as a density slope mismatch. For several galaxies the magnitude of this inner mass deficit is well in excess of that reported in recent simulations where cores result from baryon-induced fluctuations in the gravitational potential. We conclude that one or more of the following statements must be true: (i) the dark matter is more complex than envisaged by any current model; (ii) current simulations fail to reproduce the effects of baryons on the inner regions of dwarf galaxies; and/or (iii) the mass profiles of "inner mass deficit" galaxies inferred from kinematic data are incorrect.

A universality of dark-halo surface density for the Milky Way and Andromeda dwarf satellites as a probe of the coldness of dark matter

AbstractWe propose a new astrophysical test on the nature of dark matter based on the properties of dark halos associated with dwarf spheroidal galaxies. The method adopts a mean surface density of a dark halo defined within a radius of maximum circular velocity, which is derivable for a wide variety of galaxies with any dark-matter density profiles. We find that even though dark halo density profiles are derived based on the different assumptions for each galaxy sample, this surface density is generally constant across a wide mass range of galaxy. We find that at higher halo-mass scales, this constancy for real galaxies can be naturally reproduced by both cold and warm dark matter (CDM and WDM) models. However, at low-mass scales, for which we have estimated from the Milky Way and Andromeda dwarf satellites, the mean surface density derived from WDM models largely deviates from the observed constancy, whereas CDM models are in reasonable agreement with observations.

A robust determination of Milky Way satellite properties using hierarchical mass modelling

We introduce a new methodology to robustly determine the mass profile, as well as the overall distribution, of Local Group satellite galaxies. Specifically we employ a statistical multilevel modelling technique, Bayesian hierarchical modelling, to simultaneously constrain the properties of individual Local Group Milky Way satellite galaxies and the characteristics of the Milky Way satellite population. We show that this methodology reduces the uncertainty in individual dwarf galaxy mass measurements up to a factor of a few for the faintest galaxies. We find that the distribution of Milky Way satellites inferred by this analysis, with the exception of the apparent lack of high-mass haloes, is consistent with the Lambda cold dark matter (Lambda-CDM) paradigm. In particular we find that both the measured relationship between the maximum circular velocity and the radius at this velocity, as well as the inferred relationship between the mass within 300 pc and luminosity, match the values predicted by Lambda-CDM simulations for halos with maximum circular velocities below 20 km/sec. Perhaps more striking is that this analysis seems to suggest a more cusped "average" halo shape that is shared by these galaxies. While this study reconciles many of the observed properties of the Milky Way satellite distribution with that of Lambda-CDM simulations, we find that there is still a deficit of satellites with maximum circular velocities of 20-40 km/sec.

A COMMON SURFACE-DENSITY SCALE FOR THE MILKY WAY AND ANDROMEDA DWARF SATELLITES AS A CONSTRAINT ON DARK MATTER MODELS

In an attempt to place an explicit constraint on dark matter models, we define and estimate a mean surface density of a dark halo within a radius of maximum circular velocity, which is derivable for various galaxies with any dark-matter density profiles. We find that this surface density is generally constant across a wide range of maximum circular velocities of $\sim$ 10 to 400 km s$^{-1}$, irrespective of different density distribution in each of the galaxies. This common surface density at high halo-mass scales is found to be naturally reproduced by both cold and warm dark matter (CDM and WDM) models, even without employing any fitting procedures. However, the common surface density at dwarf-galaxy scales, for which we have derived from the Milky Way and Andromeda dwarf satellites, is reproduced only in a massive range of WDM particle masses, whereas CDM provides a reasonable agreement with the observed constancy. This is due to the striking difference between mass-concentration relations for CDM and WDM halos at low halo-mass scales. In order to explain the universal surface density of dwarf-galaxy scales in WDM models, we suggest that WDM particles need to be heavier than 3 keV.

THE BLACK HOLE–DARK MATTER HALO CONNECTION

We explore the connection between the central supermassive blackholes (SMBH) in galaxies and the dark matter halo through the relation between the masses of the SMBHs and the maximum circular velocities of the host galaxies, as well as the relationship between stellar velocity dispersion of the spheroidal component and the circular velocity. Our assumption here is that the circular velocity is a proxy for the mass of the dark matter halo. We rely on a heterogeneous sample containing galaxies of all types. The only requirement is that the galaxy has a direct measurement of the mass of its SMBH and a direct measurement of its circular velocity and its velocity dispersion. Previous studies have analyzed the connection between the SMBH and dark matter halo through the relationship between the circular velocity and the bulge velocity dispersion, with the assumption that the bulge velocity dispersion stands in for the mass of the SMBH, via the well{}-established SMBH mass{}-bulge velocity dispersion relation. Using intermediate relations may be misleading when one is studying them to decipher the active ingredients of galaxy formation and evolution. We believe that our approach will provide a more direct probe of the SMBH and the dark matter halo connection. We find that the correlation between the mass of supermassive blackholes and the circular velocities of the host galaxies is extremely weak, leading us to state the dark matter halo may not play a major role in regulating the blackhole growth in the present Universe.

Subhalo statistics of galactic haloes: beyond the resolution limit

We study the substructure population of Milky Way (MW)-mass halos in the $\Lambda$CDM cosmology using a novel procedure to extrapolate subhalo number statistics beyond the resolution limit of N-body simulations. The technique recovers the mean and the variance of the subhalo abundance, but not its spatial distribution. It extends the dynamic range over which precise statistical predictions can be made by the equivalent of performing a simulation with 50 times higher resolution, at no additional computational cost. We apply this technique to MW-mass halos, but it can easily be applied to halos of any mass. We find up to $20\%$ more substructures in MW-mass halos than found in previous studies. Our analysis lowers the mass of the MW halo required to accommodate the observation that the MW has only three satellites with a maximum circular velocity $V_{max}\ge30 km/s$ in the $\Lambda$CDM cosmology. The probability of having a subhalo population similar to that in the MW is $20\%$ for a virial mass, $M_{200}=1\times10^{12} M_\odot$ and practically zero for halos more massive than $M_{200}=2\times10^{12} M_\odot$.

Milky Way mass constraints from the Galactic satellite gap

We use the distribution of maximum circular velocities, $V_{max}$, of satellites in the Milky Way (MW) to constrain the virial mass, $M_{200}$, of the Galactic halo under an assumed prior of a $\Lambda$CDM universe. This is done by analysing the subhalo populations of a large sample of halos found in the Millennium II cosmological simulation. The observation that the MW has at most three subhalos with $V_{max}\ge30 km/s$ requires a halo mass $M_{200}\le1.4\times10^{12} M_\odot$, while the existence of the Magellanic Clouds (assumed to have $V_{max}\ge60 km/s$) requires $M_{200}\ge1.0\times10^{12} M_\odot$. The first of these conditions is necessary to avoid the "too-big-to-fail" problem highlighted by Boylan-Kolchin et al., while the second stems from the observation that massive satellites like the Magellanic Clouds are rare. When combining both requirements, we find that the MW halo mass must lie in the range $0.25 \le M_{200}/(10^{12} M_\odot) \le 1.4$ at $90\%$ confidence. The gap in the abundance of Galactic satellites between $30 km/s\le V_{max} \le 60 km/s$ places our galaxy in the tail of the expected satellite distribution.

Coming of age in the dark sector: how dark matter haloes grow their gravitational potential wells

We present a detailed study of how dark matter haloes assemble their mass and grow their (central) potential well. We characterize these via their mass accretion histories (MAHs) and potential well growth histories (PWGHs), which we extract from the Bolshoi simulation and from semi-analytical merger trees supplemented with a method to compute the maximum circular velocity, Vmax, of progenitor haloes. The results of both methods are in excellent agreement, both in terms of the average and the scatter. We show that the MAH and PWGH are tightly correlated, and that growth of the central potential precedes the assembly of mass; the maximum circular velocity is already half the present-day value by the time the halo has accreted only 2 per cent of its final mass. Finally, we demonstrate that MAHs have a universal form, which we use to develop a new and improved universal model that can be used to compute the average or median MAH and PWGH for a halo of any mass in any Λ cold dark matter cosmology, without having to run a numerical simulation or a set of halo merger trees.

THE EFFECTS OF VARYING COSMOLOGICAL PARAMETERS ON HALO SUBSTRUCTURE

We investigate how different cosmological parameters, such as those delivered by the WMAP and Planck missions, affect the nature and evolution of dark matter halo substructure. We use a series of flat $\Lambda$ cold dark matter ($\Lambda$CDM) cosmological $N$-body simulations of structure formation, each with a different power spectrum but the same initial white noise field. Our fiducial simulation is based on parameters from the WMAP 7th year cosmology. We then systematically vary the spectral index, $n_s$, matter density, $\Omega_M$, and normalization of the power spectrum, $\sigma_8$, for 7 unique simulations. Across these, we study variations in the subhalo mass function, mass fraction, maximum circular velocity function, spatial distribution, concentration, formation times, accretion times, and peak mass. We eliminate dependence of subhalo properties on host halo mass and average over many hosts to reduce variance. While the "same" subhalos from identical initial overdensity peaks in higher $\sigma_8, n_s$, and $\Omega_m$ simulations accrete earlier and end up less massive and closer to the halo center at $z=0$, the process of continuous subhalo accretion and destruction leads to a steady state distribution of these properties across all subhalos in a given host. This steady state mechanism eliminates cosmological dependence on all properties listed above except subhalo concentration and $V_{max}$, which remain greater for higher $\sigma_8, n_s$ and $\Omega_m$ simulations, and subhalo formation time, which remains earlier. We also find that the numerical technique for computing scale radius and the halo finder used can significantly affect the concentration-mass relationship computed for a simulation.

THE MASSES OF LOCAL GROUP DWARF SPHEROIDAL GALAXIES: THE DEATH OF THE UNIVERSAL MASS PROFILE

We investigate the claim that all dwarf spheroidal galaxies (dSphs) reside within halos that share a common, universal mass profile as has been derived for dSphs of the Galaxy. By folding in kinematic information for 25 Andromeda dSphs, more than doubling the previous sample size, we find that a singular mass profile can not be found to fit all the observations well. Further, the best-fit dark matter density profile measured for solely the Milky Way dSphs is marginally discrepant (at just beyond the 1 sigma level) with that of the Andromeda dSphs, where a profile with lower maximum circular velocity, and hence mass, is preferred. The agreement is significantly better when three extreme Andromeda outliers, And XIX, XXI and XXV, all of which have large half-light radii (>600pc) and low velocity dispersions (sigma_v < 5km/s) are omitted from the sample. We argue that the unusual properties of these outliers are likely caused by tidal interactions with the host galaxy.

A model for halo formation with axion mixed dark matter

There are several issues to do with dwarf galaxy predictions in the standard $\Lambda$CDM cosmology that have suscitated much recent debate about the possible modification of the nature of dark matter as providing a solution. We explore a novel solution involving ultra-light axions that can potentially resolve the missing satellites problem, the cusp-core problem, and the `too big to fail' problem. We discuss approximations to non-linear structure formation in dark matter models containing a component of ultra-light axions across four orders of magnitude in mass, $10^{-24}\text{ eV}\lesssim m_a \lesssim 10^{-20}\text{ eV}$, a range too heavy to be well constrained by linear cosmological probes such as the CMB and matter power spectrum, and too light/non-interacting for other astrophysical or terrestrial axion searches. We find that an axion of mass $m_a\approx 10^{-21}\text{ eV}$ contributing approximately $85\%$ of the total dark matter can introduce a significant kpc scale core in a typical Milky Way satellite galaxy in sharp contrast to a thermal relic with a transfer function cut off at the same scale, while still allowing such galaxies to form in significant number. Therefore ultra-light axions do not suffer from the \emph{Catch 22} that applies to using a warm dark matter as a solution to the small scale problems of cold dark matter. Our model simultaneously allows formation of enough high redshift galaxies to allow reconciliation with observational constraints, and also reduces the maximum circular velocities of massive dwarfs so that baryonic feedback may more plausibly resolve the predicted overproduction of massive MWG dwarf satellites.

The dark side of galaxy colour

We present age distribution matching, a theoretical formalism for predicting how galaxies of luminosity L and color C occupy dark matter halos. Our model supposes that there are just two fundamental properties of a halo that determine the color and brightness of the galaxy it hosts: the maximum circular velocity Vmax, and the redshift z_starve that correlates with the epoch at which the star formation in the galaxy ceases. The halo property z_starve is intended to encompass physical characteristics of halo mass assembly that may deprive the galaxy of its cold gas supply and, ultimately, quench its star formation. The new, defining feature of the model is that, at fixed luminosity, galaxy color is in monotonic correspondence with z_starve, with the larger values of z_starve being assigned redder colors. We populate an N- body simulation with a mock galaxy catalog based on age distribution matching, and show that the resulting mock galaxy distribution accurately describes a variety of galaxy statistics. Our model suggests that halo and galaxy assembly are indeed correlated. We make publicly available our low-redshift, SDSS M_r <-19 mock galaxy catalog, and main progenitor histories of all z=0 halos, at http://logrus.uchicago.edu/~aphearin