Globular Cluster Systems and the Missing Satellite Problem: Implications for Cold Dark Matter Models

We analyze the metallicity distributions of globular clusters belonging to 28 early-type galaxies in the survey of Kundu & Whitmore (2001). A Monte Carlo algorithm which simulates the chemical evolution of galaxies that grow hierarchically via dissipationless mergers is used to determine the most probable protogalactic mass function for each galaxy. Contrary to the claims of Kundu & Whitmore, we find that the observed metallicity distributions are in close agreement with the predictions of such hierarchical formation models. The mass spectrum of protogalactic fragments for the galaxies in our sample has a power-law behavior, with an exponent of roughly -2. This spectrum is indistinguishable from the mass spectrum of dark matter halos predicted by cold dark matter models for structure formation. We argue that these protogalactic fragments, the likely sites of globular cluster formation in the early universe, are the disrupted remains of the "missing" satellite galaxies predicted by cold dark matter models. Our findings suggest that the solution to the missing satellite problem is through the suppression of gas accretion in low-mass halos after reionization, or via self-interacting dark matter, and argue against models with suppressed small-scale power or warm dark matter.Comment: 28 pages, 19 postscript figures. Accepted for publication in the Astrophysical Journa

Three-Dimensional Geometric Morphometric Analysis of Fossil Canid Mandibles and Skulls

Acknowledgements We thank C.P. Klingenberg for critical discussion of methodology. A. Drake and R. Losey were supported by a grant from the Social Sciences and Humanities Research Council of Canada grant (#SSHRC IG 435-2014-0075) and a European Research Council Grant to D. Anderson (#295458). M. Sablin acknowledges participation of ZIN RAS (state assignment № АААА-А17-117022810195-3) to this research. Supplementary information accompanies this paper at doi:10.1038/s41598-017-10232-1Peer reviewedPublisher PD

The Clustering of Massive Halos

The clustering properties of dark matter halos are a firm prediction of modern theories of structure formation. We use two large volume, high-resolution N-body simulations to study how the correlation function of massive dark matter halos depends upon their mass and formation history. We find that halos with the lowest concentrations are presently more clustered than those of higher concentration, the size of the effect increasing with halo mass; this agrees with trends found in studies of lower mass halos. The clustering dependence on other characterizations of the full mass accretion history appears weaker than the effect with concentration. Using the integrated correlation function, marked correlation functions, and a power-law fit to the correlation function, we find evidence that halos which have recently undergone a major merger or a large mass gain have slightly enhanced clustering relative to a randomly chosen population with the same mass distribution.Comment: 10 pages, 8 figures; text improved, references and one figure added; accepted for publication in Ap

Robustness of Cosmological Simulations I: Large Scale Structure

The gravitationally-driven evolution of cold dark matter dominates the formation of structure in the Universe over a wide range of length scales. While the longest scales can be treated by perturbation theory, a fully quantitative understanding of nonlinear effects requires the application of large-scale particle simulation methods. Additionally, precision predictions for next-generation observations, such as weak gravitational lensing, can only be obtained from numerical simulations. In this paper, we compare results from several N-body codes using test problems and a diverse set of diagnostics, focusing on a medium resolution regime appropriate for studying many observationally relevant aspects of structure formation. Our conclusions are that -- despite the use of different algorithms and error-control methodologies -- overall, the codes yield consistent results. The agreement over a wide range of scales for the cosmological tests is test-dependent. In the best cases, it is at the 5% level or better, however, for other cases it can be significantly larger than 10%. These include the halo mass function at low masses and the mass power spectrum at small scales. While there exist explanations for most of the discrepancies, our results point to the need for significant improvement in N-body errors and their understanding to match the precision of near-future observations. The simulation results, including halo catalogs, and initial conditions used, are publicly available.Comment: 32 pages, 53 figures, data from the simulations is available at http://t8web.lanl.gov/people/heitmann/arxiv, accepted for publication in ApJS, several minor revisions, reference added, main conclusions unchange

Precision Determination of the Mass Function of Dark Matter Halos

The predicted mass function of dark matter halos is essential in connecting observed galaxy cluster counts and models of galaxy clustering to the properties of the primordial density field. We determine the mass function in the concordance $\Lambda$CDM cosmology, as well as its uncertainty, using sixteen $1024^3$-particle nested-volume dark-matter simulations, spanning a mass range of over five orders of magnitude. Using the nested volumes and single-halo tests, we find and correct for a systematic error in the friends-of-friends halo-finding algorithm. We find a fitting form and full error covariance for the mass function that successfully describes the simulations' mass function and is well-behaved outside the simulations' resolutions. Estimated forecasts of uncertainty in cosmological parameters from future cluster count surveys have negligible contribution from remaining statistical uncertainties in the central cosmology multiplicity function. There exists a potentially non-negligible cosmological dependence (non-universality) of the halo multiplicity function.Comment: 4 pages, 3 figures, submitted to ApJ

Close Pairs as Proxies for Galaxy Cluster Mergers

Galaxy cluster merger statistics are an important component in understanding the formation of large-scale structure. Unfortunately, it is difficult to study merger properties and evolution directly because the identification of cluster mergers in observations is problematic. We use large N-body simulations to study the statistical properties of massive halo mergers, specifically investigating the utility of close halo pairs as proxies for mergers. We examine the relationship between pairs and mergers for a wide range of merger timescales, halo masses, and redshifts (0<z<1). We also quantify the utility of pairs in measuring merger bias. While pairs at very small separations will reliably merge, these constitute a small fraction of the total merger population. Thus, pairs do not provide a reliable direct proxy to the total merger population. We do find an intriguing universality in the relation between close pairs and mergers, which in principle could allow for an estimate of the statistical merger rate from the pair fraction within a scaled separation, but including the effects of redshift space distortions strongly degrades this relation. We find similar behavior for galaxy-mass halos, making our results applicable to field galaxy mergers at high redshift. We investigate how the halo merger rate can be statistically described by the halo mass function via the merger kernel (coagulation), finding an interesting environmental dependence of merging: halos within the mass resolution of our simulations merge less efficiently in overdense environments. Specifically, halo pairs with separations less than a few Mpc/h are more likely to merge in underdense environments; at larger separations, pairs are more likely to merge in overdense environments.Comment: 12 pages, 9 figures; Accepted for publication in ApJ. Significant additions to text and two figures changed. Added new findings on the universality of pair mergers and added analysis of the effect of FoF linking length on halo merger

The Environmental Dependence of the Infrared Luminosity and Stellar Mass Functions

We investigate the dependence of the galaxy infrared luminosity function (LF) and the associated stellar mass function (SMF) on environment and spectral type using photometry from the Two Micron All Sky Survey and redshifts from the Las Campanas Redshift Survey for galaxies brighter than M_J<-19+5 log h. In the field environment, galaxies with emission lines have LFs with much steeper faint end slopes (alpha_J=-1.39) than galaxies without emission lines (alpha_J=-0.59). In the cluster environment, however, even the non-emission line galaxies have a steep faint-end LF (alpha_J=-1.22). There is also a significant (95%) difference between the overall cluster and field LFs, $\Delta \alpha_J=-0.34, \Delta M_J^\ast=-0.54$. All of these variations are more pronounced in the SMFs, which we compute by relating the strength of the 4000 A break in the optical spectra to a mass-to-light ratio.Comment: 14 pages, 13 figures, emulateapj style ApJ, post-referee. Very minor changes, mostly typographical in natur

The Use of Sequential Surveys to Shorten Implementation Time for Healthcare System-Level Enhanced Recovery After Surgery (ERAS) Pathways

Background Enhanced Recovery After Surgery (ERAS) pathways improve healthcare quality, safety, and cost-effectiveness. We hypothesized that the RAND Method (a hybrid Delphi approach), involving anonymous sequential surveys and face-to-face meetings, would allow for more rapid agreement and initiation of new ERAS pathways. Methods Using the ERAS Society guidelines for cesarean section as a baseline, our institution’s ERAS Leadership Team (ELT) compiled published literature and institutional practices to design a 32-component survey that was sent to obstetricians, nurse midwives, anesthesiologists, pharmacists, and nurses. Components that did not reach 90% consensus were included in a second survey the following week, and meetings were held to review results. At the conclusion of this process, time to agreement was retrospectively compared to the colorectal ERAS pathway process at this institution. Results ERAS pathway components were compiled and reviewed by 121 stakeholders at 7 hospitals using iterative surveys with review meetings over a 13-week period. Survey response rates were 61% and 50% in the initial and follow-up surveys, respectively. There was agreement on 28/32 and 32/32 items on the initial and follow-up surveys. Using the RAND Method, time to agreement decreased by 54.1% (24 vs 13 weeks) compared to prior system-wide efforts to standardize the colorectal surgery ERAS pathway. Discussion With rapidly expanding healthcare systems, effective methods to gain consensus and adopt ERAS pathways are critical to implementation of ERAS guidelines. We demonstrate that the RAND Method allows for a transparent and efficient means of agreement across a diverse group of clinicians practicing in several settings

Using Cluster Abundances and Peculiar Velocities to Test the Gaussianity of the Cosmological Density Field

(Abridged) By comparing the frequency of typical events with that of unusual events, one can test whether the cosmological density distribution function is consistent with the normally made assumption of Gaussianity. To this end, we compare the consistency of the tail-inferred (from clusters) and measured values (from large-scale flows) of the rms level of mass fluctuations for two distribution functions: a Gaussian, and a texture (positively-skewed) PDF. Averaging the recent large-scale flow measurements, we find that observations of the rms and the tail at the 10 h^-1 Mpc scale disfavor a texture PDF at ~1.5 sigma in all cases. However, taking only the most recent measurement of the rms, that from Willick et al. (1997b), the comparison disfavors textures for low Omega_0=0.3, and disfavors Gaussian models if Omega_0=1 (again at ~1.5 sigma). Predictions for evolution of high temperature clusters can also be made for the models considered, and strongly disfavor Omega_0=1 in Gaussian models and marginally disfavor Omega_0=1 in texture models. Only Omega_0=0.3 Gaussian models are consistent with all the data considered.Comment: 34 pg incl. 8 embedded figures, LaTeX, aaspp4.sty, submitted to Ap