Search

Galaxy cluster abundance measurements provide a classic test of cosmology. They are most sensitive to the evolved amplitude of fluctuations, usually expressed as $S_8 = \sigma_8\sqrt{\Omega_m/0.3}$. Thus, abundance constraints exhibit a strong degeneracy between $\sigma_8$ and $\Omega_{\rm m}$, as do other similar low-redshift tests such as cosmic shear. The mass distribution in the infall region around galaxy clusters, where material is being accreted from the surrounding field, also exhibits a cosmological dependence, but in this case it is nearly orthogonal to the $S_8$ direction in the $\Omega_m$--$\sigma_8$ plane, making it highly complementary to halo abundance or cosmic shear studies. We explore how weak lensing measurements of the infall region might be used to complement abundance studies, considering three different tests. The splashback radius is a prominent feature of the infall region; we show that detection of this feature in lensing data from the Euclid survey could independently constrain $\Omega_{\rm m}$ and $\sigma_8$ to $\pm 0.05$. Another feature, the depletion radius where the bias reaches a minimum, also shows cosmological dependence, though it is challenging to observe in practice. The strongest constraints come from direct measurements of the shear profile in the infall region at $2$--$4\,r_{200{\rm c}}$. Combining the latter with abundance constraints such as those reported from SRG$/$eROSITA should reduce the area of the error contours by an estimated factor of $1.2$ using a sample of clusters observed by the UNIONS survey, or a factor of $3$ using clusters observed by the Euclid Wide survey over a broader range of redshift., Comment: 14 pages, 10 figures. Accepted for publication in MNRAS

Cosmological parameters such as $\Omega_{\rm{M}}$ and $\sigma_{8}$ can be measured indirectly using various methods, including galaxy cluster abundance and cosmic shear. These measurements constrain the composite parameter $S_{8}$, leading to degeneracy between $\Omega_{\rm{M}}$ and $\sigma_{8}$. However, some structural properties of galaxy clusters also correlate with cosmological parameters, due to their dependence on a cluster's accretion history. In this work, we focus on the splashback radius, an observable cluster feature that represents a boundary between a cluster and the surrounding Universe. Using a suite of cosmological simulations with a range of values for $\Omega_{\rm{M}}$ and $\sigma_{8}$, we show that the position of the splashback radius around cluster-mass halos is greater in cosmologies with smaller values of $\Omega_{\rm{M}}$ or larger values of $\sigma_{8}$. This variation breaks the degeneracy between $\Omega_{\rm{M}}$ and $\sigma_{8}$ that comes from measurements of the $S_{8}$ parameter. We also show that this variation is, in principle, measurable in observations. As the splashback radius can be determined from the same weak lensing analysis already used to estimate $S_{8}$, this new approach can tighten low-redshift constraints on cosmological parameters, either using existing data, or using upcoming data such as that from Euclid and LSST., Comment: 12 pages, 7 figures, 1 table, accepted for publication in ApJ

Numerous metrics exist to quantify the dynamical state of galaxy clusters, both observationally and within simulations. Many of these correlate strongly with one another, but it is not clear whether all of these measures probe the same intrinsic properties. In this work, we use two different statistical approaches -- principal component analysis (PCA) and uniform manifold approximation and projection (UMAP) -- to investigate which dynamical properties of a cluster are in fact the best descriptors of its dynamical state. We use measurements taken directly from The Three Hundred suite of galaxy cluster simulations, as well as morphological properties calculated using mock X-ray and SZ maps of the same simulated clusters. We find that four descriptions of dynamical state naturally arise, and although correlations exist between these, a given cluster can be "dynamically relaxed" according to all, none, or some of these four descriptions. These results demonstrate that it is highly important for future observational and theoretical studies to consider in which sense clusters are dynamically relaxed. Cluster dynamical states are complex and multi-dimensional, and so it is not meaningful to classify them simply as "relaxed" and "unrelaxed" based on a single linear scale., Comment: 18 pages, 7 figures, 3 tables, accepted for publication in MNRAS

Dark matter haloes grow at a rate that depends on the value of the cosmological parameters $\sigma_8$ and $\Omega_{\rm m}$ through the initial power spectrum and the linear growth factor. While halo abundance is routinely used to constrain these parameters, through cluster abundance studies, the halo growth rate is not. In recent work, we proposed constraining the cosmological parameters using observational estimates of the overall dynamical "age" of clusters, expressed, for instance, by their half-mass assembly redshift $z_{50}$. Here we explore the prospects for using the instantaneous growth rate, as estimated from the halo merger rate, from the average growth rate over the last dynamical time, or from the fraction of systems with recent episodes of major growth. We show that the merger rate is mainly sensitive to the amplitude of fluctuations $\sigma_8$, while the rates of recent growth provide constraints in the $\Omega_{\rm m}$-$\sigma_8$ plane that are almost orthogonal to those provided by abundance studies. Data collected for forthcoming cluster abundance studies, or studies of the galaxy merger rate in current and future galaxy surveys, may thus provide additional constraints on the cosmological parameters complementary to those already derived from halo abundance., Comment: 16 pages, 12 figures, 1 table and two figures in appendix

Dwarf galaxies are important tracers of small-scale cosmological structure, yet much of our knowledge about these systems comes from the limited sample of dwarf galaxies within the Local Group. To make a comprehensive inventory of dwarf populations in the local Universe, we require effective methods for deriving distance estimates for large numbers of faint, low surface brightness objects. Here we test the surface brightness fluctuation (SBF) method, traditionally applied to brighter early-type galaxies, on a sample of 20 nearby dwarf galaxies detected in the COSMOS field. These objects are partially resolved in HST ACS images, and have confirmed redshift distances in the range 17-130 Mpc. We discuss the many model choices required in applying the SBF method, and explore how these affect the final distance estimates. Amongst other variations on the method, when applying the SBF method, we alter the standard equation to include a term accounting for the power spectrum of the background, greatly improving our results. For the most robust modelling choices, we find a roughly Gaussian SBF signal that correlates linearly with distance out to distances of 50-100 Mpc, but with only a fraction of the power expected. At larger distances, there is excess power relative to that predicted, probably from undetected point sources. Overall, obtaining accurate SBF distances to faint, irregular galaxies remains challenging, but may yet prove possible with the inclusion of more information about galaxy properties and point source populations, and the use of more advanced techniques., Comment: 18 pages, 26 figures, accepted by MNRAS

A growing body of evidence suggests that the central density of cuspy dark matter subhaloes is conserved in minor mergers. However, empirical models of subhalo evolution, calibrated from simulations, often assume a drop in the central density. Since empirical models of subhaloes are used in galaxy-galaxy lensing studies and dark matter annihilation calculations, we explore the consequences of assuming different subhalo models. We find that dark matter annihilation calculations are very sensitive to the assumed subhalo mass profile, and different models can give more than a magnitude difference in the J-factor and boost factor in individual haloes. On the other hand, the shear and convergence profiles used in galaxy-galaxy lensing are sensitive to the initial profile assumed (e.g., NFW versus Einato) but are otherwise well-approximated by a simple model in which the original profile is sharply truncated. We conclude that since the innermost parts of haloes are difficult to resolve in simulations, it is important to have a theoretical understanding of how subhaloes evolve to make accurate predictions of the dark matter annihilation signal., Comment: 15 pages, 13 figures. Submitted to MNRAS

Accurate models of the structural evolution of dark matter subhaloes, as they orbit within larger systems, are fundamental to understanding the detailed distribution of dark matter at the present day. Numerical simulations of subhalo evolution support the idea that the mass loss associated with tidal stripping is most naturally understood in energy space, with the particles that are the least bound being removed first. Starting from this premise, we recently proposed a zero-parameter "energy-truncation model" for subhalo evolution. We tested this model with simulations of tidal stripping of satellites with initial NFW profiles, and showed that the energy-truncation model accurately predicts both the mass loss and density profiles. In this work, we apply the model to a variety of Hernquist, Einasto and King profiles. We show that it matches the simulation results quite closely in all cases, indicating that it may serve as a universal model to describe tidally stripped collisionless systems. A key prediction of the energy-truncation model is that the central density of dark matter subhaloes is conserved as they lose mass; this has important implications for dark matter annihilation calculations, and for other observational tests of dark matter., Comment: 19 pages, 14 figures. Accepted for publication in MNRAS

We present a new method for quantifying the abundance of satellites around field galaxies and in groups. The method is designed to work with samples, such as local photometric redshift catalogues, that do not have full spectroscopic coverage, but for which some redshift or distance information is available. It consists of identifying the galaxies most likely to be centrals, and using the clustering signal around them as a template to iteratively decompose the full population into satellite and central populations. In that sense it is similar to performing crowded-field photometry, after having first used isolated stars to determine the point spread function of the image. The method does not identify individual satellites or centrals conclusively, but assigns a probability to each galaxy of being one or the other. Averaged over a large sample, it provides a statistical estimate of satellite abundance, even in crowded fields with large redshift uncertainties. We test the method using data from the COSMOS field, which includes a large set of local objects with accurate photometric redshifts. We measure satellite abundance as a function of central stellar or halo mass, as well as the satellite luminosity function, and find results consistent with previous studies, but extending over a broader range of central masses. We also consider a number of possible systematic uncertainties in the method, and show that they are generally smaller than our random errors. Having presented the method in this paper, we will use it to study the properties of the satellite populations in a forthcoming one., Comment: 17 pages, 17 figures

The abundance of galaxy clusters in the low-redshift universe provides an important cosmological test, constraining a product of the initial amplitude of fluctuations and the amount by which they have grown since early times. The degeneracy of the test with respect to these two factors remains a limitation of abundance studies. Clusters will have different mean assembly times, however, depending on the relative importance of initial fluctuation amplitude and subsequent growth. Thus, structural probes of cluster age such as concentration, shape or substructure may provide a new cosmological test that breaks the main degeneracy in number counts. We review analytic predictions for how mean assembly time should depend on cosmological parameters, and test these predictions using cosmological simulations. Given the overall sensitivity expected, we estimate the cosmological parameter constraints that could be derived from the cluster catalogues of forthcoming surveys such as Euclid, the Nancy Grace Roman Space Telescope, eROSITA, or CMB-S4. We show that by considering the structural properties of their cluster samples, such surveys could easily achieve errors of $\Delta \sigma_8$ = 0.01 or better., Comment: 18 pages, 17 figures

The Shelby County iZone is a district-led effort to dramatically improve, or "turn around," 23 of the lowest performing schools in Tennessee in the 15th largest district in the country. Despite circumstances that have derailed many past reforms, iZone schools have made statistically significant and educationally meaningful gains in mathematics and English language arts (ELA) (Zimmer, Henry, Kho, & Viano, 2015). This report presents the findings of a multi-year research project that examined the evolution of the iZone as it shifted its strategy from school-level autonomy to one that featured a common curriculum, shared pedagogy, and collegial learning. The analysis delves into the challenges that iZone leaders, principals, and teachers confronted in coping with the needs of a student population mired in intergeneration poverty, rigorous performance standards, and a stringent accountability system. The results illuminate the importance of multi-level system re-design, continuous improvement, and compromise and negotiation among key stakeholder groups.

The orbital parameters of dark matter (DM) subhaloes play an essential role in determining their mass-loss rates and overall spatial distribution within a host halo. Haloes in cosmological simulations grow by a combination of relatively smooth accretion and more violent mergers, and both processes will modify subhalo orbits. To isolate the impact of the smooth growth of the host halo from other relevant mechanisms, we study subhalo orbital evolution using numerical calculations in which subhaloes are modelled as massless particles orbiting in a time-varying spherical potential. We find that the radial action of the subhalo orbit decreases over the first few orbits, indicating that the response to the growth of the host halo is not adiabatic during this phase. The subhalo orbits can shrink by a factor of $\sim$1.5 in this phase. Subsequently, the radial action is well conserved and orbital contraction slows down. We propose a model accurately describing the orbital evolution. Given these results, we consider the spatial distribution of the population of subhaloes identified in high-resolution cosmological simulations. We find that it is consistent with this population having been accreted at z < 3, indicating that any subhaloes accreted earlier are unresolved in the simulations. We also discuss tidal stripping as a formation scenario for NGC1052-DF2, an ultra diffuse galaxy significantly lacking DM, and find that its expected DM mass could be consistent with observational constraints if its progenitor was accreted early enough, z > 1.5, although it should still be a relatively rare object., Comment: 15 pages, 11 figures, 1 table, accepted for publication in MNRAS

The ability to accurately predict the evolution of tidally stripped haloes is important for understanding galaxy formation and testing the properties of dark matter. Most studies of substructure evolution make predictions based on empirical models of tidal mass loss that are calibrated using numerical simulations. This approach can be accurate in the cases considered, but lacks generality and does not provide a physical understanding of the processes involved. Recently, we demonstrated that truncating NFW distribution functions sharply in energy results in density profiles that resemble those of tidally stripped systems, offering a path to constructing physically motivated models of tidal mass loss. In this work, we review calculations of mass loss based on energy truncation alone, and then consider what secondary effects may modulate mass loss beyond this. We find that a combination of dependence on additional orbital parameters and variations in individual particle energies over an orbit results in a less abrupt truncation in energy space as a subhalo loses mass. Combining the energy truncation approach with a simple prediction for the mass-loss rate, we construct a full model of mass loss that can accurately predict the evolution of a subhalo in terms of a single parameter $\eta_{\rm eff}$. This parameter can be fully determined from the initial orbital and halo properties, and does not require calibration with numerical simulations., Comment: 20 pages, 14 figures, accepted to MNRAS

Kinematic studies have produced accurate measurements of the total dark matter mass and mean dark matter density within the optical extent of galaxies, for large samples of objects. Here we consider theoretical predictions for the latter quantity, $\bar{\rho}_{dm}$, measured within the isophotal radius $R_{23.5}$, for isolated halos with universal density profiles. Through a combination of empirical scaling relations, we show that $\bar{\rho}_{dm}$ is expected to depend weakly on halo mass and strongly on redshift. When galaxy halos fall into larger groups or clusters they become tidally stripped, reducing their total dark matter mass, but this process is expected to preserve central density until an object is close to disruption. We confirm this with collisonless simulations of cluster formation, finding that subhalos have values of $\bar{\rho}_{dm}$ close to the analytic predictions. This suggests that $\bar{\rho}_{dm}$ may be a useful indicator of infall redshift onto the cluster. We test this hypothesis with data from the SHIVir survey, which covers a reasonable fraction of the Virgo cluster. We find that galaxies with high $\bar{\rho}_{dm}$ do indeed trace the densest regions of the cluster, with a few notable exceptions. Samples selected by environment have higher densities at a significance of 3.5-4$\sigma$, while samples selected by density are more clustered at 3-3.5$\sigma$ significance. We conclude that halo density can be a powerful tracer of the assembly history of clusters and their member galaxies., Comment: 17 pages, 18 figures; accepted for publication in MNRAS

The structural properties of individual dark matter haloes, including shape, spin, concentration, and substructure, are linked to the halo's growth history, but the exact connection between the two is unclear. One open question, in particular, is the effect of major mergers on halo structure. We have performed a large set of simulations of binary equal-mass mergers between isolated haloes with various density profiles, to map out the relationship between the initial conditions and merger parameters and the structure of the final remnant. In this paper we describe our initial set-up and analysis methods, and report on the results for the size, shape, and spin of the merger remnant. The outcomes of mergers are most easily understood in terms of a scaled dimensionless energy parameter $\kappa$ and an angular momentum (or spin) parameter $\lambda$. We find that the axis ratio $c/a$ scales roughly linearly with energy $\kappa$ while the axis ratio $c/b$ scales linearly with spin $\lambda$. Qualitatively, mergers on radial orbits produce prolate remnants, while mergers on tangential orbits produce oblate remnants. The spin of the remnant can be predicted from angular momentum conservation, while the overall size changes as $\sim \kappa^{-5}$, as expected from self-similar scaling at constant mean density. We discuss potential cosmological applications for these simple patterns., Comment: 18 pages, 20 figures; accepted for publication, MNRAS

Several lines of evidence suggest that as dark matter haloes grow their scale radius increases, and that the density in their central region drops. Major mergers seem an obvious mechanism to explain both these phenomena, and the resulting patterns in the concentration--mass--redshift relation. To test this possibility, we have simulated equal-mass mergers between haloes with a variety of cosmological density profiles, placed on various different orbits. The remnants typically have higher densities than the initial conditions, but differ only slightly from self-similar scaling predictions. They are reasonably well fit by Einasto profiles, but have parameters distinct from those of the initial conditions. The net internal energy available to the merger remnant, relative to the internal energy of the initial conditions, $\kappa$, has the greatest influence on the properties of the final mass distribution. As expected, energetic encounters produce more extended remnants while mergers of strongly bound systems produce compact remnants. Surprisingly, however, the scale radius of the density profile shows the opposite trend, {\it increasing} in the remnants of low-energy encounters relative to energetic ones. Also even in the most energetic encounters, the density within the scale radius decreases only slightly (by 10--20\%), while for very low-energy systems it increases significantly after the merger. We conclude that while major mergers can produce remnants that are more diffuse at large radii, they are relatively ineffective at changing the central densities of haloes, and seem unlikely to explain the mean trends in the concentration--mass--redshift relation., Comment: 19 pages, 16 figures; accepted for publication, MNRAS

Faint dwarf satellite galaxies are important as tracers of small-scale structure, but remain poorly characterized outside the Local Group, due to the difficulty of identifying them consistently at larger distances. We review a recently proposed method for estimating the average satellite population around a given sample of nearby bright galaxies, using a combination of size and magnitude cuts (to select low-redshift dwarf galaxies preferentially) and clustering measurements (to estimate the fraction of true satellites in the cut sample). We test this method using the high-precision photometric redshift catalog of the COSMOS survey, exploring the effect of specific cuts on the clustering signal. The most effective of the size-magnitude cuts considered recover the clustering signal around low-redshift primaries (z < 0.15) with about two-thirds of the signal and 80\%\ of the signal-to-noise ratio obtainable using the full COSMOS photometric redshifts. These cuts are also fairly efficient, with more than one third of the selected objects being clustered satellites. We conclude that structural selection represents a useful tool in characterizing dwarf populations to fainter magnitudes and/or over larger areas than are feasible with spectroscopic surveys.} In reviewing the low-redshift content of the COSMOS field, we also note the existence of several dozen objects that appear resolved or partially resolved in the HST imaging, and are confirmed to be local (at distances of $\sim$250 Mpc or less) by their photometric or spectroscopic redshifts. This underlines the potential for future space-based surveys to reveal local populations of intrinsically faint galaxies through imaging alone., Comment: 22 pages, 21 figures; accepted for publication in MNRAS

Individual extragalactic dark matter halos, such as those associated with nearby galaxies and galaxy clusters, are promising targets for searches for gamma-rays from dark matter annihilation. We review the predictions for the annihilation flux from individual halos, focusing on the effect of current uncertainties in the concentration-mass relation and the contribution from halo substructure, and also estimating the intrinsic halo-to-halo scatter expected. After careful consideration of recent simulation results, we conclude that the concentrations of the smallest halos, while well determined at high redshift, are still uncertain by a factor of 4-6 when extrapolated to low redshift. This in turn produces up to two orders of magnitude uncertainty in the predicted annihilation flux for any halo mass above this scale. Substructure evolution, the small-scale cutoff to the power spectrum, cosmology, and baryonic effects all introduce smaller, though cumulative, uncertainties. We then consider intrinsic variations from halo to halo. These arise from variations in concentration and substructure, leading to a scatter of $\sim$ 2.5 in the predicted annihilation luminosity. Finally, we consider the problem of detecting gamma-rays from annihilation, given the expected contributions from other sources. We estimate the signal-to-noise ratio for gamma-ray detection as a function of halo mass, assuming that cosmic rays from star formation are the main noise source in the detection. This calculation suggests that group-scale halos, individually or in stacks, may be a particularly interesting target for the next generation of annihilation searches., Comment: 30 pages, 9 figures, updated to published version

We present here parameter distributions and fundamental scaling relations for 190 galaxies as part of the Spectroscopy and H-bang Imaging of Virgo cluster galaxies (SHIVir) survey. We find the distribution of galaxy velocities to be bimodal about $V_{\rm circ} \sim 125$ km ${\rm s^{-1}}$, hinting at the existence of dynamically unstable modes in the inner regions of galaxies. An analysis of the Tully-Fisher relation (TFR) of late-type galaxies (LTGs) and fundamental plane (FP) of early-type galaxies (ETGs) is also presented, yielding a compendium of galaxy scaling relations. The slope and zero-point of the Virgo TFR match those of field galaxies, while scatter differences likely reflect distinct evolutionary histories. The velocities minimizing scatter for the TFR and FP are measured at large apertures where the baryonic fraction becomes subdominant. While TFR residuals remain independent of any galaxy parameters, FP residuals (i.e. the FP "tilt") correlate strongly with the dynamical-to-stellar mass ratio, yielding stringent galaxy formation constraints. Furthermore, we construct a stellar-to-total mass relation (STMR) for ETGs and LTGs and find linear but distinct trends over the range $M_{*} = 10^{8-11} M_{\odot}$. Stellar-to-halo mass relations (SHMRs), which probe the extended dark matter halo, can be scaled down to masses estimated within the optical radius, showing a tight match with the Virgo STMR at low masses; however, possibly inadequate halo abundance matching prescriptions and broad radial scalings complicate this comparison at all masses. While ETGs appear to be more compact than LTGs of the same stellar mass in projected space, their mass-size relations in physical space are identical. The trends reported here call for validation through well-resolved numerical simulations., Comment: Accepted to ApJ. 28 pages, 16 figures

We propose a new method for generating equilibrium models of spherical systems of collisionless particles that are finite in extent, but whose central regions resemble dark matter halos from cosmological simulations. This method involves iteratively removing unbound particles from a Navarro-Frenk-White profile truncated sharply at some radius. The resulting models are extremely stable, and thus provide a good starting point for N-body simulations of isolated halos. We provide a code to generate such models for NFW and a variety of other common density profiles. We then develop an analytic approximation to this truncated distribution function. Our method proceeds by analogy with the King model, truncating and shifting the original distribution function of an infinitely extended Navarro-Frenk-White profile in energy space. We show that the density profiles of our models closely resemble the tidally truncated density profiles seen previously in studies of satellite evolution. Pursuing this analogy further with a series of simulations of tidal mass loss, we find that our models provide a good approximation to the full distribution function of tidally stripped systems, thus allowing theoretically motivated phase-space calculations for such systems.

We use deep optical photometry from the Next Generation Virgo Cluster Survey [NGVS] to investigate the color-magnitude diagram for the galaxies inhabiting the core of this cluster. The sensitivity of the NGVS imaging allows us to continuously probe galaxy colors over a factor of $\sim 2 \times 10^5$ in luminosity, from brightest cluster galaxies to scales overlapping classical satellites of the Milky Way [$M_{g^{\prime}}$ $\sim$ $-$9; $M_{*}$ $\sim 10^6$ M$_{\odot}$], within a single environment. Remarkably, we find the first evidence that the RS flattens in all colors at the faint-magnitude end [starting between $-$14 $\le$ $M_{g^{\prime}}$ $\le$ $-$13, around $M_{*}$ $\sim 4 \times 10^7$ M$_{\odot}$], with the slope decreasing to $\sim$60% or less of its value at brighter magnitudes. This could indicate that the stellar populations of faint dwarfs in Virgo's core share similar characteristics [e.g. constant mean age] over $\sim$3 mags in luminosity, suggesting that these galaxies were quenched coevally, likely via pre-processing in smaller hosts. We also compare our results to galaxy formation models, finding that the RS in model clusters have slopes at intermediate magnitudes that are too shallow, and in the case of semi-analytic models, do not reproduce the flattening seen at both extremes [bright/faint] of the Virgo RS. Deficiencies in the chemical evolution of model galaxies likely contribute to the model-data discrepancies at all masses, while overly efficient quenching may also be a factor at dwarf scales. Deep UV and near-IR photometry are required to unambiguously diagnose the cause of the faint-end flattening., Comment: 19 pages, 11 figures, accepted to ApJ

We present measurements of the galaxy luminosity and stellar mass function in a 3.71 deg$^2$ (0.3 Mpc$^2$) area in the core of the Virgo cluster, based on $ugriz$ data from the Next Generation Virgo Cluster Survey (NGVS). The galaxy sample consists of 352 objects brighter than $M_g=-9.13$ mag, the 50% completeness limit of the survey. Using a Bayesian analysis, we find a best-fit faint end slope of $\alpha=-1.33 \pm 0.02$ for the g-band luminosity function; consistent results are found for the stellar mass function as well as the luminosity function in the other four NGVS bandpasses. We discuss the implications for the faint-end slope of adding 92 ultra compact dwarfs galaxies (UCDs) -- previously compiled by the NGVS in this region -- to the galaxy sample, assuming that UCDs are the stripped remnants of nucleated dwarf galaxies. Under this assumption, the slope of the luminosity function (down to the UCD faint magnitude limit, $M_g = -9.6$ mag) increases dramatically, up to $\alpha = -1.60 \pm 0.06$ when correcting for the expected number of disrupted non-nucleated galaxies. We also calculate the total number of UCDs and globular clusters that may have been deposited in the core of Virgo due to the disruption of satellites, both nucleated and non-nucleated. We estimate that ~150 objects with $M_g\lesssim-9.6$ mag and that are currently classified as globular clusters, might, in fact, be the nuclei of disrupted galaxies. We further estimate that as many as 40% of the (mostly blue) globular clusters in the core of Virgo might once have belonged to such satellites; these same disrupted satellites might have contributed ~40% of the total luminosity in galaxies observed in the core region today. Finally, we use an updated Local Group galaxy catalog to provide a new measurement of the luminosity function of Local Group satellites, $\alpha=-1.21\pm0.05$., Comment: ApJ, accepted

The Next Generation Virgo Cluster Survey has recently determined the luminosity function of galaxies in the core of the Virgo cluster down to unprecedented magnitude and surface brightness limits. Comparing simulations of cluster formation to the derived central stellar mass function, we attempt to estimate the stellar-to-halo-mass ratio (SHMR) for dwarf galaxies, as it would have been before they fell into the cluster. This approach ignores several details and complications, e.g., the contribution of ongoing star formation to the present-day stellar mass of cluster members, and the effects of adiabatic contraction and/or violent feedback on the subhalo and cluster potentials. The final results are startlingly simple, however; we find that the trends in the SHMR determined previously for bright galaxies appear to extend down in a scale-invariant way to the faintest objects detected in the survey. These results extend measurements of the formation efficiency of field galaxies by two decades in halo mass, or five decades in stellar mass, down to some of the least massive dwarf galaxies known, with stellar masses of $\sim 10^5 M_\odot$., Comment: 18 pages, 12 figures; published in ApJ July 1st 2015

We present a new measurement of the mass-concentration relation and the stellar-to-halo mass ratio over the halo mass range $5\times 10^{12}$ to $2\times 10^{14}M_{\odot}$. To achieve this, we use weak lensing measurements from the CFHT Stripe 82 Survey (CS82), combined with the central galaxies from the redMaPPer cluster catalogue and the LOWZ/CMASS galaxy sample of the Sloan Digital Sky Survey-III Baryon Oscillation Spectroscopic Survey Tenth Data Release. The stacked lensing signals around these samples are modelled as a sum of contributions from the central galaxy, its dark matter halo, and the neighboring halos, as well as a term for possible centering errors. We measure the mass-concentration relation: $c_{200c}(M)=A(\frac{M_{200c}}{M_0})^{B}$ with $A=5.24\pm1.24, B=-0.13\pm0.10$ for $0.2

We report on a large-scale study of the distribution of globular clusters (GCs) throughout the Virgo cluster, based on photometry from the Next Generation Virgo Cluster Survey, a large imaging survey covering Virgo's primary subclusters to their virial radii. Using the g', (g'-i') color-magnitude diagram of unresolved and marginally-resolved sources, we constructed 2-D maps of the GC distribution. We present the clearest evidence to date showing the difference in concentration between red and blue GCs over the extent of the cluster, where the red (metal-rich) GCs are largely located around the massive early-type galaxies, whilst the blue (metal-poor) GCs have a more extended spatial distribution, with significant populations present beyond 83' (215 kpc) along the major axes of M49 and M87. The GC distribution around M87 and M49 shows remarkable agreement with the shape, ellipticity and boxiness of the diffuse light surrounding both galaxies. We find evidence for spatial enhancements of GCs surrounding M87 that may be indicative of recent interactions or an ongoing merger history. We compare the GC map to the locations of Virgo galaxies and the intracluster X-ray gas, and find good agreement between these baryonic structures. The Virgo cluster contains a total population of 67300$\pm$14400 GCs, of which 35% are located in M87 and M49 alone. We compute a cluster-wide specific frequency S_N,CL=$2.8\pm0.7$, including Virgo's diffuse light. The GC-to-baryonic mass fraction is e_b=$5.7\pm1.1\times10^{-4} $and the GC-to-total cluster mass formation efficiency is e_t=$2.9\pm0.5\times10^{-5}$, values slightly lower than, but consistent with, those derived for individual galactic halos. Our results show that the production of the complex structures in the unrelaxed Virgo cluster core (including the diffuse intracluster light) is an ongoing process.(abridged), Comment: 23 pages, 17 figures. Accepted for publication in the Astrophysical Journal. Figure 1 has reduced resolution. Revised version with updated references, corrected typos -- no changes to results

We measure the cross-correlation of cosmic microwave background lensing convergence maps derived from Atacama Cosmology Telescope data with galaxy lensing convergence maps as measured by the Canada-France-Hawaii Telescope Stripe 82 Survey. The CMB-galaxy lensing cross power spectrum is measured for the first time with a significance of 4.2{\sigma}, which corresponds to a 12% constraint on the amplitude of density fluctuations at redshifts ~ 0.9. With upcoming improved lensing data, this novel type of measurement will become a powerful cosmological probe, providing a precise measurement of the mass distribution at intermediate redshifts and serving as a calibrator for systematic biases in weak lensing measurements., Comment: 11 pages, 7 figures. Submitted to PRD. Updates to lensing maps and analysis pipeline, with revised results, increased significance

We present a weak lensing mass map covering ~124 square degrees of the Canada-France-Hawaii Telescope Stripe 82 Survey (CS82). We study the statistics of rare peaks in the map, including peak abundance, the peak-peak correlation functions and the tangential-shear profiles around peaks. We find that the abundance of peaks detected in CS82 is consistent with predictions from a Lambda-CDM cosmological model, once noise effects are properly included. The correlation functions of peaks with different signal-to-noise ratio (SNR) are well described by power laws, and there is a clear cross-correlation between the Sloan Digital Sky Survey III/Constant Mass galaxies and high SNR peaks. The tangential-shear profiles around peaks increase with peak SNR. We fit analytical models to the tangential-shear profiles, including a projected singular isothermal sphere (SIS) model and a projected Navarro, Frenk & White (NFW) model, plus a two-halo term. For the high SNR peaks, the SIS model is rejected at ~3-sigma. The NFW model plus a two-halo term gives more acceptable fits to the data. Some peaks match the positions of optically detected clusters, while others are relatively dark. Comparing dark and matched peaks, we find a difference in lensing signal of a factor of 2, suggesting that about half of the dark peaks are false detections., Comment: 10 pages, 12 figures, Published in MNRAS

The number and distribution of dwarf satellite galaxies remain a critical test of cold dark matter-dominated structure formation on small scales. Until recently, observational information about galaxy formation on these scales has been limited mainly to the Local Group. We have searched for faint analogues of Local Group dwarfs around nearby bright galaxies, using a spatial clustering analysis of the photometric catalog of the Sloan Digital Sky Survey (SDSS) Data Release 8. Several other recent searches of SDSS have detected clustered satellite populations down to $\Delta m_r \equiv ({m}_{r,\, {\rm sat}} -\, {m}_{r,\, {\rm main}}) \sim 6$-$8$, using photometric redshifts to reduce background contamination. SDSS photometric redshifts are relatively imprecise, however, for faint and nearby galaxies. Instead we use angular size to select potential nearby dwarfs, and consider only the nearest isolated bright galaxies as primaries. As a result, we are able to detect an excess clustering signal from companions down to $\Delta m_r = 12$, four magnitudes fainter than most recent studies. We detect an over-density of objects at separations $< 400$ kpc, corresponding to about $4.6 \pm 0.5$ satellites per central galaxy, consistent with the satellite abundance expected from the Local Group given our selection function. Although the sample of satellites detected is incomplete by construction, since it excludes the least and most compact dwarfs, this detection provides a lower bound on the average satellite luminosity function, down to luminosities corresponding to the faintest "classical" dwarfs of the Local Group., Comment: 12 pages, 7 figures, submitted to ApJ

The scaling between X-ray observables and mass for galaxy clusters and groups is instrumental for cluster based cosmology and an important probe for the thermodynamics of the intracluster gas. We calibrate a scaling relation between the weak lensing mass and X-ray spectroscopic temperature for 10 galaxy groups in the COSMOS field, combined with 55 higher mass clusters from the literature. The COSMOS data includes HST imaging and redshift measurements of 46 source galaxies per square arcmin, enabling us to perform unique weak lensing measurements of low mass systems. Our sample extends the mass range of the lensing calibrated M-T relation an order of magnitude lower than any previous study, resulting in a power-law slope of 1.48$^{+0.13}_{-0.09}$. The slope is consistent with the self-similar model, predictions from simulations, and observations of clusters. However, X-ray observations relying on mass measurements derived under the assumption of hydrostatic equilibrium have indicated that masses at group scales are lower than expected. Both simulations and observations suggest that hydrostatic mass measurements can be biased low. Our external weak lensing masses provides the first observational support for hydrostatic mass bias at group level, showing an increasing bias with decreasing temperature and reaching a level of 30-50% at 1 keV., Comment: 13 pages, 12 figures, submitted to ApJ

Primordial, non-Gaussian perturbations can generate scale-dependent bias in the galaxy distribution. This in turn will modify correlations between galaxy positions and peculiar velocities at late times, since peculiar velocities reflect the underlying matter distribution, whereas galaxies are a biased tracer of the same. We study this effect, and show that non-Gaussianity can be constrained by comparing the observed peculiar velocity field to a model velocity field reconstructed from the galaxy density field assuming linear bias. The amplitude of the spatial correlations in the residual map obtained after subtracting one velocity field from the other is directly proportional to the strength of the primordial non-Gaussianity. We construct the corresponding likelihood function use it to constrain the amplitude of the linear flow $\beta$ and the amplitude of local non-Gaussianity $f^{\rm NL}_{\rm local}$. Applying our method to two observational data sets, the Type-Ia supernovae (A1SN) and Spiral Field \textit{I}-band (SFI++) catalogues, we obtain constraints on the linear flow parameter consistent with the values derived previously assuming Gaussianity. The marginalised 1-D distribution of $|f^{\rm NL}_{\rm local}|$ does not show strong evidence for non-zero $f^{\rm NL}_{\rm local}$, and we set 95% upper limits $|f^{\rm NL}_{\rm local}|<51.4$ from A1SN and $|f^{\rm NL}_{\rm local}|<92.6$ from SFI++. These limits on $f^{\rm NL}_{\rm local}$ are as tight as any set by previous large-scale structure measurements. Our method can be applied to any survey with radial velocities and density field data, and provides an independent check of recent CMB constraints on $f^{\rm NL}_{\rm local}$, extending these to smaller spatial scales., Comment: 10 pages, 11 figures

In the theory of structure formation, galaxies are biased tracers of the underlying matter density field. The statistical relation between galaxy and matter density field is commonly referred as galaxy bias. In this paper, we test the linear bias model with weak-lensing and galaxy clustering measurements in the 2 square degrees COSMOS field (Scoville et al. 2007). We estimate the bias of galaxies between redshifts z=0.2 and z=1, and over correlation scales between R=0.2 h^-1 Mpc and R=15 h^-1 Mpc. We focus on three galaxy samples, selected in flux (simultaneous cuts I_814W < 26.5 and K_s < 24), and in stellar-mass (10^9 < M_* < 10^10 h^-2 Msun and 10^10 < M^*< 10^11 h^-2 Msun). At scales R > 2 h^-1 Mpc, our measurements support a model of bias increasing with redshift. The Tinker et al. (2010) fitting function provides a good fit to the data. We find the best fit mass of the galaxy halos to be log(M_200 h^-1 Msun) = 11.7^+0.6_-1.3 and log(M_200 h^-1 Msun) = 12.4^+0.2_-2.9 respectively for the low and high stellar-mass samples. In the halo model framework, bias is scale-dependent with a change of slope at the transition scale between the one and the two halo terms. We detect a scale-dependence of bias with a turn-down at scale R=2.3\pm1.5 h^-1 Mpc, in agreement with previous galaxy clustering studies. We find no significant amount of stochasticity, suggesting that a linear bias model is sufficient to describe our data. We use N-body simulations to quantify both the amount of cosmic variance and systematic errors in the measurement., Comment: 20 pages, 18 figures, accepted for publication in ApJ

Individual dark matter halos in cosmological simulations vary widely in their detailed structural properties such as shape, rotation, substructure and degree of internal relaxation. Recent non-parametric (principal component) analyses suggest that a few principal components explain a large fraction of the scatter in halo properties. The main principal component is closely linked with concentration, which in turn is known to be related to the mass accretion history of the halo. Here we examine more generally the connection between mass accretion history and structural parameters. The space of mass accretion histories has principal components of its own. We find that the strongest two can be interpreted as the overall age of the halo and the acceleration or deceleration of growth at late times. These two components only account for $\sim70$%\ of the scatter in mass accretions histories however, due to the stochastic effect of major mergers. Relating structural parameters to formation history, we find that concentration correlates strongly with the early history of the halo, while relaxation correlates with the late history. We examine the inferences about formation history that can be drawn by splitting haloes into subsamples, based on observable properties such as concentration and shape at some final time. This approach suggests interesting possibilities, such as the possibility of defining young and old samples of galaxy clusters in a rigorous, quantitative way, or testing the dynamical assumptions of galaxy formation models empirically., Comment: 14 pages, 12 figures

Gravitational lensing can provide pure geometric tests of the structure of space-time, for instance by determining empirically the angular diameter distance-redshift relation. This geometric test has been demonstrated several times using massive clusters which produce a large lensing signal. In this case, matter at a single redshift dominates the lensing signal, so the analysis is straightforward. It is less clear how weaker signals from multiple sources at different redshifts can be stacked to demonstrate the geometric dependence. We introduce a simple measure of relative shear which for flat cosmologies separates the effect of lens and source positions into multiplicative terms, allowing signals from many different source-lens pairs to be combined. Applying this technique to a sample of groups and low-mass clusters in the COSMOS survey, we detect a clear variation of shear with distance behind the lens. This represents the first detection of the geometric effect using weak lensing by multiple, low-mass systems. The variation of distance with redshift is measured with sufficient precision to constrain the equation of state of the universe under the assumption of flatness, equivalent to a detection of a dark energy component Omega_X at greater than 99% confidence for an equation-of-state parameter -2.5 < w < -0.1. For the case w = -1, we find a value for the cosmological constant density parameter Omega_Lambda = 0.85+0.044-0.19 (68% C.L.), and detect cosmic acceleration (q_0 < 0) at the 98% C.L.. We consider the systematic uncertainties associated with this technique and discuss the prospects for applying it in forthcoming weak-lensing surveys., Comment: 13 pages, 7 figures; submitted to ApJ

Using data from the COSMOS survey, we perform the first joint analysis of galaxy-galaxy weak lensing, galaxy spatial clustering, and galaxy number densities. Carefully accounting for sample variance and for scatter between stellar and halo mass, we model all three observables simultaneously using a novel and self-consistent theoretical framework. Our results provide strong constraints on the shape and redshift evolution of the stellar-to-halo mass relation (SHMR) from z=0.2 to z=1. At low stellar mass, we find that halo mass scales as Mh M*^0.46 and that this scaling does not evolve significantly with redshift to z=1. We show that the dark-to-stellar ratio, Mh/M*, varies from low to high masses, reaching a minimum of Mh/M*~27 at M*=4.5x10^10 Msun and Mh=1.2x10^12 Msun. This minimum is important for models of galaxy formation because it marks the mass at which the accumulated stellar growth of the central galaxy has been the most efficient. We describe the SHMR at this minimum in terms of the "pivot stellar mass", M*piv, the "pivot halo mass", Mhpiv, and the "pivot ratio", (Mh/M*)piv. Thanks to a homogeneous analysis of a single data set, we report the first detection of mass downsizing trends for both Mhpiv and M*piv. The pivot stellar mass decreases from M*piv=5.75+-0.13x10^10 Msun at z=0.88 to M*piv=3.55+-0.17x10^10 Msun at z=0.37. Intriguingly, however, the corresponding evolution of Mhpiv leaves the pivot ratio constant with redshift at (Mh/M*)piv~27. We use simple arguments to show how this result raises the possibility that star formation quenching may ultimately depend on Mh/M* and not simply Mh, as is commonly assumed. We show that simple models with such a dependence naturally lead to downsizing in the sites of star formation. Finally, we discuss the implications of our results in the context of popular quenching models, including disk instabilities and AGN feedback., Comment: 31 pages, 18 figures, 5 tables. Submitted to ApJ. Comments welcome

We present the results of Suzaku X-ray follow-up observation of weak-lensing-detected halos in the field around galaxy cluster ZwCl0823.2+0425. We clearly detected X-ray emission associated with most of these halos and determined their detailed physical parameters such as X-ray luminosity, temperature, and metal abundance, for the first time. We find that the X-ray luminosity - temperature relation for these halos agrees with former typical results. With mass determined from the weak gravitational lensing data, the mass-temperature relation for them is also investigated and found to be consistent with the prediction from a simple self-similar model and results of the previous studies with both lensing and X-ray data. We would like to emphasize that the self-similar scaling relation of mass and temperature is shown here for the first time using a weak-lensing selected sample, whereas previous studies of the mass scaling relation used X-ray-selected samples of clusters. Therefore, our study demonstrates importance of X-ray follow-up observations of shear-selected clusters, and shows that a joint X-ray and lensing analysis will be crucial for clusters discovered by the forthcoming weak-lensing surveys, such as the one planned with Subaru/Hyper-Suprime-Cam., Comment: 20 pages, 13 figures, accepted for publication in PASJ

The balance of evidence indicates that individual galaxies and groups or clusters of galaxies are embedded in enormous distributions of cold, weakly interacting dark matter. These dark matter 'halos' provide the scaffolding for all luminous structure in the universe, and their properties comprise an essential part of the current cosmological model. I review the internal properties of dark matter halos, focussing on the simple, universal trends predicted by numerical simulations of structure formation. Simulations indicate that halos should all have roughly the same spherically-averaged density profile and kinematic structure, and predict simple distributions of shape, formation history and substructure in density and kinematics, over an enormous range of halo mass and for all common variants of the concordance cosmology. I describe observational progress towards testing these predictions by measuring masses, shapes, profiles and substructure in real halos, using baryonic tracers or gravitational lensing. An important property of simulated halos (possibly the most important property) is their dynamical 'age', or degree of internal relaxation. The age of a halo may have almost as much effect as its mass in determining the state of its baryonic contents, so halo ages are also worth trying to measure observationally. I review recent gravitational lensing studies of galaxy clusters which should measure substructure and relaxation in a large sample of individual cluster halos, producing quantitative measures of age that are well-matched to theoretical predictions. The age distributions inferred from these studies will lead to second-generation tests of the cosmological model, as well as an improved understanding of cluster assembly and the evolution of galaxies within clusters., Comment: v2: additional references and minor corrections to match the published version

The Chandra COSMOS Survey (C-COSMOS) is a large, 1.8 Ms, Chandra} program that has imaged the central 0.5 sq.deg of the COSMOS field (centered at 10h, +02deg) with an effective exposure of ~160ksec, and an outer 0.4sq.deg. area with an effective exposure of ~80ksec. The limiting source detection depths are 1.9e-16 erg cm(-2) s(-1) in the Soft (0.5-2 keV) band, 7.3e(-16) erg cm^-2 s^-1 in the Hard (2-10 keV) band, and 5.7e(-16) erg cm(-2) s(-1) in the Full (0.5-10 keV) band. Here we describe the strategy, design and execution of the C-COSMOS survey, and present the catalog of 1761 point sources detected at a probability of being spurious of <2e(-5) (1655 in the Full, 1340 in the Soft, and 1017 in the Hard bands). By using a grid of 36 heavily (~50%) overlapping pointing positions with the ACIS-I imager, a remarkably uniform (to 12%) exposure across the inner 0.5 sq.deg field was obtained, leading to a sharply defined lower flux limit. The widely different PSFs obtained in each exposure at each point in the field required a novel source detection method, because of the overlapping tiling strategy, which is described in a companion paper. (Puccetti et al. Paper II). This method produced reliable sources down to a 7-12 counts, as verified by the resulting logN-logS curve, with sub-arcsecond positions, enabling optical and infrared identifications of virtually all sources, as reported in a second companion paper (Civano et al. Paper III). The full catalog is described here in detail, and is available on-line., Comment: Revised to omit egregious bold facing and fix missing ',' in author list

There is a vast menagerie of plausible candidates for the constituents of dark matter, both within and beyond extensions of the Standard Model of particle physics. Each of these candidates may have scattering (and other) cross section properties that are consistent with the dark matter abundance, BBN, and the most scales in the matter power spectrum; but which may have vastly different behavior at sub-galactic "cutoff" scales, below which dark matter density fluctuations are smoothed out. The only way to quantitatively measure the power spectrum behavior at sub-galactic scales at distances beyond the local universe, and indeed over cosmic time, is through probes available in multiply imaged strong gravitational lenses. Gravitational potential perturbations by dark matter substructure encode information in the observed relative magnifications, positions, and time delays in a strong lens. Each of these is sensitive to a different moment of the substructure mass function and to different effective mass ranges of the substructure. The time delay perturbations, in particular, are proving to be largely immune to the degeneracies and systematic uncertainties that have impacted exploitation of strong lenses for such studies. There is great potential for a coordinated theoretical and observational effort to enable a sophisticated exploitation of strong gravitational lenses as direct probes of dark matter properties. This opportunity motivates this white paper, and drives the need for: a) strong support of the theoretical work necessary to understand all astrophysical consequences for different dark matter candidates; and b) tailored observational campaigns, and even a fully dedicated mission, to obtain the requisite data., Comment: Science white paper submitted to the Astro2010 Decadal Cosmology & Fundamental Physics Science Frontier Panel