Constraining ultra-compact dwarf galaxy formation with galaxy clusters in the local universe

We compare the predictions of a semi-analytic model for ultra-compact dwarf galaxy (UCD) formation by tidal stripping to the observed properties of globular clusters (GCs) and UCDs in the Fornax and Virgo clusters. For Fornax we find the predicted number of stripped nuclei agrees very well with the excess number of GCs$+$UCDs above the GC luminosity function. GCs$+$UCDs with masses $>10^{7.3}$ M$_\odot$ are consistent with being entirely formed by tidal stripping. Stripped nuclei can also account for Virgo UCDs with masses $>10^{7.3}$ M$_\odot$ where numbers are complete by mass. For both Fornax and Virgo, the predicted velocity dispersions and radial distributions of stripped nuclei are consistent with that of UCDs within $\sim$50-100 kpc but disagree at larger distances where dispersions are too high and radial distributions too extended. Stripped nuclei are predicted to have radially biased anisotropies at all radii, agreeing with Virgo UCDs at clustercentric distances larger than 50 kpc. However, ongoing disruption is not included in our model which would cause orbits to become tangentially biased at small radii. We find the predicted metallicities and central black hole masses of stripped nuclei agree well with the metallicities and implied black hole masses of UCDs for masses $>10^{6.5}$ M$_\odot$. The predicted black hole masses also agree well with that of M60-UCD1, the first UCD with a confirmed central black hole. These results suggest that observed GC$+$UCD populations are a combination of genuine GCs and stripped nuclei, with the contribution of stripped nuclei increasing toward the high-mass end.Comment: 18 pages, 12 figures, accepted for publication in MNRA

Contribution of stripped nuclear clusters to globular cluster and ultra-compact dwarf galaxy populations

We use the Millennium II cosmological simulation combined with the semi-analytic galaxy formation model of Guo et al. (2011) to predict the contribution of galactic nuclei formed by the tidal stripping of nucleated dwarf galaxies to globular cluster (GC) and ultra-compact dwarf galaxy (UCD) populations of galaxies. We follow the merger trees of galaxies in clusters back in time and determine the absolute number and stellar masses of disrupted galaxies. We assume that at all times nuclei have a distribution in nucleus-to-galaxy mass and nucleation fraction of galaxies similar to that observed in the present day universe. Our results show stripped nuclei follow a mass function $N(M) \sim M^{-1.5}$ in the mass range $10^6 < M/M_\odot < 10^8$, significantly flatter than found for globular clusters. The contribution of stripped nuclei will therefore be most important among high-mass GCs and UCDs. For the Milky Way we predict between 1 and 3 star clusters more massive than $10^5 M_\odot$ come from tidally disrupted dwarf galaxies, with the most massive cluster formed having a typical mass of a few times $10^6 M_\odot$, like omega Centauri. For a galaxy cluster with a mass $7 \times 10^{13} M_\odot$, similar to Fornax, we predict $\sim$19 UCDs more massive than $2\times10^6 M_\odot$ and $\sim$9 UCDs more massive than $10^7 M_\odot$ within a projected distance of 300 kpc come from tidally stripped dwarf galaxies. The observed number of UCDs are $\sim$200 and 23, respectively. We conclude that most UCDs in galaxy clusters are probably simply the high mass end of the GC mass function.Comment: 15 pages, 9 figures, accepted for publication in MNRA

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.Comment: 15 pages, 15 figures, Accepted to ApJ on March 15, 201

Predicting the locations of possible long-lived low-mass first stars: Importance of satellite dwarf galaxies

The search for metal-free stars has so far been unsuccessful, proving that if there are surviving stars from the first generation, they are rare, they have been polluted, or we have been looking in the wrong place. To predict the likely location of Population~III (Pop~III) survivors, we semi-analytically model early star formation in progenitors of Milky Way-like galaxies and their environments. We base our model on merger trees from the high-resolution dark matter only simulation suite \textit{Caterpillar}. Radiative and chemical feedback are taken into account self-consistently, based on the spatial distribution of the haloes. Our results are consistent with the non-detection of Pop III survivors in the Milky Way today. We find that possible surviving Population III stars are more common in Milky Way satellites than in the main Galaxy. In particular, low mass Milky Way satellites contain a much larger fraction of Pop~III stars than the Milky Way. Such nearby, low mass Milky Way satellites are promising targets for future attempts to find Pop~III survivors, especially for high-resolution, high signal-to-noise spectroscopic observations. We provide the probabilities for finding a Pop~III survivor in the red giant branch phase for all known Milky Way satellites to guide future observations.Comment: 17 pages, 12 figures, 1 table, submitted to MNRA

The Caterpillar Project: A Large Suite of Milky Way Sized Halos

We present the largest number of Milky Way sized dark matter halos simulated at very high mass ($\sim$$10^4$ M$_\odot$/particle) and temporal resolution ($\sim$5 Myrs/snapshot) done to date, quadrupling what is currently available in the literature. This initial suite consists of the first 24 halos of the $Caterpillar$ $Project$ (www.caterpillarproject.org) whose project goal of 60 - 70 halos will be made public when complete. We resolve $\sim$20,000 gravitationally bound subhalos within the virial radius of each host halo. Over the ranges set by our spatial resolution our convergence is excellent and improvements were made upon current state-of-the-art halo finders to better identify substructure at such high resolutions (e.g., on average we recover $\sim$4 subhalos in each host halo above 10$^8$ M$_\odot$ which would have otherwise not been found using conventional methods). For our relaxed halos, the inner profiles are reasonably fit by Einasto profiles ($\alpha$ = 0.169 $\pm$ 0.023) though this depends on the relaxed nature and assembly history of a given halo. Averaging over all halos, the substructure mass fraction is $f_{m,subs} = 0.121 \pm 0.041$, and mass function slope is d$N$/d$M\propto M^{-1.88 \pm 0.10}$ though we find scatter in the normalizations for fixed halo mass due to more concentrated hosts having less subhalos at fixed subhalo mass. There are no biases stemming from Lagrangian volume selection as all Lagrangian volume types are included in our sample. Our detailed contamination study of 264 low resolution halos has resulted in obtaining very large and unprecedented, high-resolution regions around our host halos for our target resolution (sphere of radius $\sim$$1.4 \pm 0.4$ Mpc) allowing for accurate studies of low mass dwarf galaxies at large galactocentric radii and the very first stellar systems at high redshift ($z \geq$ 10).Comment: 19 pages; 14 figures; 6 tables; Received September 3, 2015; Accepted November 15, 2015; Published February 2, 201

An analysis of the evolving comoving number density of galaxies in hydrodynamical simulations

The cumulative comoving number-density of galaxies as a function of stellar mass or central velocity dispersion is commonly used to link galaxy populations across different epochs. By assuming that galaxies preserve their number-density in time, one can infer the evolution of their properties, such as masses, sizes, and morphologies. However, this assumption does not hold in the presence of galaxy mergers or when rank ordering is broken owing to variable stellar growth rates. We present an analysis of the evolving comoving number density of galaxy populations found in the Illustris cosmological hydrodynamical simulation focused on the redshift range $0\leq z \leq 3$. Our primary results are as follows: 1) The inferred average stellar mass evolution obtained via a constant comoving number density assumption is systematically biased compared to the merger tree results at the factor of $\sim$2(4) level when tracking galaxies from redshift $z=0$ out to redshift $z=2(3)$; 2) The median number density evolution for galaxy populations tracked forward in time is shallower than for galaxy populations tracked backward in time; 3) A similar evolution in the median number density of tracked galaxy populations is found regardless of whether number density is assigned via stellar mass, stellar velocity dispersion, or dark matter halo mass; 4) Explicit tracking reveals a large diversity in galaxies' assembly histories that cannot be captured by constant number-density analyses; 5) The significant scatter in galaxy linking methods is only marginally reduced by considering a number of additional physical and observable galaxy properties as realized in our simulation. We provide fits for the forward and backward median evolution in stellar mass and number density and discuss implications of our analysis for interpreting multi-epoch galaxy property observations.Comment: 18 pages, 11 figures, submitted to MNRAS, comments welcom

Biliary dyskinesia: a potentially unrecognized cause of abdominal pain in children

Biliary dyskinesia is defined as symptomatic biliary colic without cholelithiasis, and is diagnosed during cholescintigraphy by assessing gallbladder emptying with cholecystokinin (CCK) stimulation. Unfortunately, gallbladder emptying is not routinely assessed during cholescintigraphy in pediatric patients. The purpose of this review is to assess the effectiveness of cholecystectomy in patients with chronic abdominal pain and delayed gallbladder emptying and to assess whether these findings correlate with the histologic evidence of chronic cholecystitis. We retrospectively reviewed the medical records of all patients ( n =16) at our institution from October 1997 to August 2001 who underwent quantitative cholescintigraphy with CCK stimulation that demonstrated delayed gallbladder emptying (<35% at 60 min) and who subsequently underwent cholecystectomy. Laparoscopic cholecystectomy was performed in 16 patients with chronic abdominal pain. All 16 patients had delayed gallbladder emptying (mean ejection fraction : 15±8%, range: 3–32%). The mean age was 12±2 years (range: 8–17 years). Presenting symptoms included abdominal pain (86%), fatty food intolerance (27%), emesis (13%), and diarrhea (13%). Mean duration of abdominal pain before operation was 11±19 months (range: 2 weeks–6 years). One patient’s symptoms persisted postoperatively , but abdominal pain resolved in all other patient s. Histologic evidence of chronic cholecystitis was demonstrated in 86% of surgical specimens. Five patients underwent concurrent appendectomy , and all had normal appendiceal histology. Our experience suggests that children with chronic abdominal pain and delayed gallbladder emptying on CCK-stimulated cholescintigraphy are likely to benefit from cholecystectomy and to have histologic evidence of chronic cholecystitis.Peer Reviewedhttp://deepblue.lib.umich.edu/bitstream/2027.42/47166/1/383_2004_Article_1234.pd

An observer's guide to the (Local Group) dwarf galaxies: predictions for their own dwarf satellite populations

A recent surge in the discovery of new ultrafaint dwarf satellites of the Milky Way has inspired the idea of searching for faint satellites, 10³M[subscript ⊙] 99 per cent chance that at least one satellite with stellar mass M* > 10⁵M[subscript ⊙] exists around the combined five Local Group field dwarf galaxies with the largest stellar mass. When considering satellites with M* > 10⁴M[subscript ⊙], we predict a combined 5-25 satellites for the five largest field dwarfs, and 10-50 for the whole Local Group field dwarf population. Because of the relatively small number of predicted dwarfs, and their extended spatial distribution, a large fraction each Local Group dwarf's virial volume will need to be surveyed to guarantee discoveries. We compute the predicted number of satellites in a given field of view of specific Local Group galaxies, as a function of minimum satellite luminosity, and explicitly obtain such values for the Solitary Local dwarfs survey. Uncertainties in abundance-matching and reionization models are large, implying that comprehensive searches could lead to refinements of both models.National Science Foundation (U.S.) (Grant 1122374

Tracing the first stars and galaxies of the Milky Way

We use 30 high-resolution dark matter haloes of the Caterpillar simulation suite to probe the first stars and galaxies of MilkyWay-mass systems. We quantify the environment of the high-z progenitors of theMilkyWay and connect them to the properties of the host and satellites today. We identify the formation sites of the first generation of Population III (Pop III) stars (z ~ 25) and first galaxies (z ~ 22) with several different models based on a minimum halo mass. This includes a simple model for radiative feedback, the primary limitation of the model. Through this method we find approximately 23 000 ± 5000 Pop III potentially star-forming sites per MilkyWay-mass host, though this number is drastically reduced to ~550 star-forming sites if feedback is included. The majority of these haloes identified form in isolation (96 per cent at z = 15) and are not subject to external enrichment by neighbouring haloes (median separation ~1 kpc at z = 15), though half merge with a system larger than themselves within 1.5 Gyr. Using particle tagging, we additionally trace the Pop III remnant population to z = 0 and find an order of magnitude scatter in their number density at small (i.e. r 50 kpc) galactocentric radii. We provide fitting functions for determining the number of progenitor minihalo and atomic cooling halo systems that present-day satellite galaxies might have accreted since their formation. We determine that observed dwarf galaxies with stellar masses below 10[superscript 4.6]M[subscript ⊙] are unlikely to have merged with any other star-forming systems. Keywords: Galaxy: halo; galaxies: evolution; galaxies: formation; dark ages; reionization; first stars; dark matter; cosmology: theoryNational Science Foundation (U.S.) (Grant 1122374)National Science Foundation (U.S.) (Grant PHY-0941373)National Science Foundation (U.S.) (Grant PHY 08-22648)National Science Foundation (U.S.) (Grant PHY-1430152)United States. National Aeronautics and Space Administration (Grant NNX12AC98G)United States. National Aeronautics and Space Administration (Grant NNX15AP39G)United States. National Aeronautics and Space Administration (Grant HST-AR-13261.01-A)United States. National Aeronautics and Space Administration (Grant HST-AR-14315.001-A