BZPEERPapers, people, evidence and research paths
Menu

Discover

SearchFind papers, people, topics, institutions and opportunities from one place.GraphSee how papers, authors, topics and evidence connect.

Research tools

HomeContinue from the papers, people and topics that matter now.CollectionsTurn saved papers into reading lists, evidence maps and shareable context.ResearchReconstruct the path that led you from search to paper to question.CollaborationMove from a paper, topic or expert request into a focused thread.PublishCreate full-text research articles with metadata, formulas and versions.

Network

ExpertsFind specialists whose work and reviews explain why they can help.CommunitiesJoin research groups organized around papers, topics and shared evidence.

Opportunities

ListingsBrowse roles, calls and collaborations connected to real research context.

Account

Log inReturn to your saved papers, graph views and active threads.Create accountStart saving papers, building a research trail and joining teams.
Log inCreate account

Source author record

Shea Garrison-Kimmel

Shea Garrison-Kimmel appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record
astro-ph.GAastro-ph.COhep-phastro-ph.HEastro-ph.IMMachine Learning

Catalog footprint

What is connected

29works
6topics
4close collaborators

Actions

Connect this record

Log in to claim
Open graphBrowse works

Research graph

See the researcher in context

Open full explorer

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

29 published item(s)

preprint2022arXiv

Chasing Accreted Structures within Gaia DR2 using Deep Learning

In previous work, we developed a deep neural network classifier that only relies on phase-space information to obtain a catalog of accreted stars based on the second data release of Gaia (DR2). In this paper, we apply two clustering algorithms to identify velocity substructure within this catalog. We focus on the subset of stars with line-of-sight velocity measurements that fall in the range of Galactocentric radii $r \in [6.5, 9.5]$ kpc and vertical distances $|z| < 3$ kpc. Known structures such as Gaia Enceladus and the Helmi stream are identified. The largest previously-unknown structure, Nyx, is a vast stream consisting of at least 90 stars in the region of interest. This study displays the power of the machine learning approach by not only successfully identifying known features, but also discovering new kinematic structures that may shed light on the merger history of the Milky Way.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

Evidence for a Vast Prograde Stellar Stream in the Solar Vicinity

Massive dwarf galaxies that merge with the Milky Way on prograde orbits can be dragged into the disk plane before being completely disrupted. Such mergers can contribute to an accreted stellar disk and a dark matter disk. We present evidence for Nyx, a vast new stellar stream in the vicinity of the Sun, that may provide the first indication that such an event occurred in the Milky Way. We identify about 500 stars that have coherent radial and prograde motion in this stream using a catalog of accreted stars built by applying deep learning methods to the second Gaia data release. Nyx is concentrated within $\pm 2$ kpc of the Galactic midplane and spans the full radial range studied (6.5-9.5 kpc). The kinematics of Nyx stars are distinct from those of both the thin and thick disk. In particular, its rotational speed lags the disk by $\sim 80$ km/s and its stars follow more eccentric orbits. A small number of Nyx stars have chemical abundances or inferred ages; from these, we deduce that Nyx stars have a peak metallicity of [Fe/H] $\sim -0.5$ and ages $\sim $10-13 Gyr. Taken together with the kinematic observations, these results strongly favor the interpretation that Nyx is the remnant of a disrupted dwarf galaxy. To further justify this interpretation, we explicitly demonstrate that metal-rich, prograde streams like Nyx can be found in the disk plane of Milky Way-like galaxies using the FIRE hydrodynamic simulations. Future spectroscopic studies will be able to validate whether Nyx stars originate from a single progenitor.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

A dark matter profile to model diverse feedback-induced core sizes of $Λ$CDM haloes

We analyze the cold dark matter density profiles of 54 galaxy halos simulated with FIRE-2 galaxy formation physics, each resolved within $0.5\%$ of the halo virial radius. These halos contain galaxies with masses that range from ultra-faint dwarfs ($M_\star \simeq 10^{4.5} M_{\odot}$) to the largest spirals ($M_\star \simeq 10^{11} M_{\odot}$) and have density profiles that are both cored and cuspy. We characterize our results using a new analytic density profile that extends the standard Einasto form to allow for a pronounced constant-density core in the resolved innermost radius. With one additional core-radius parameter, $r_{c}$, this "core-Einasto" profile is able to characterize the shape and normalization of our feedback-impacted dark matter halos. In order to enable comparisons with observations, we provide fitting functions for $r_{c}$ and other profile parameters as a function of both $M_\star$ and $M_{\star}/M_{\rm halo}$. In agreement with similar studies done in the literature, we find that dark matter core formation is most efficient at the characteristic stellar-mass to halo-mass ratio $M_\star/M_{\rm halo} \simeq 5 \times 10^{-3}$, or $M_{\star} \sim 10^9 \, M_{\odot}$, with cores that are roughly the size of the galaxy half-light radius, $r_{c} \simeq 1-5$ kpc. Furthermore, we find no evidence for core formation at radii $\gtrsim 100\ \rm pc$ in galaxies with $M_{\star}/M_{\rm halo} < 5\times 10^{-4}$ or $M_\star \lesssim 10^6 \, M_{\odot}$. For Milky Way-size galaxies, baryonic contraction often makes halos significantly more concentrated and dense at the stellar half-light radius than dark matter only runs. However, even at the Milky Way scale, FIRE-2 galaxy formation still produces small dark matter cores of $\simeq 0.5-2$ kpc in size. Recent evidence for a ${\sim} 2$ kpc core in the Milky Way's dark matter halo is consistent with this expectation.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

A profile in FIRE: resolving the radial distributions of satellite galaxies in the Local Group with simulations

While many tensions between Local Group (LG) satellite galaxies and LCDM cosmology have been alleviated through recent cosmological simulations, the spatial distribution of satellites remains an important test of physical models and physical versus numerical disruption in simulations. Using the FIRE-2 cosmological zoom-in baryonic simulations, we examine the radial distributions of satellites with Mstar > 10^5 Msun around 8 isolated Milky Way- (MW) mass host galaxies and 4 hosts in LG-like pairs. We demonstrate that these simulations resolve the survival and physical destruction of satellites with Mstar >~ 10^5 Msun. The simulations broadly agree with LG observations, spanning the radial profiles around the MW and M31. This agreement does not depend strongly on satellite mass, even at distances <~ 100 kpc. Host-to-host variation dominates the scatter in satellite counts within 300 kpc of the hosts, while time variation dominates scatter within 50 kpc. More massive host galaxies within our sample have fewer satellites at small distances, likely because of enhanced tidal destruction of satellites via the baryonic disks of host galaxies. Furthermore, we quantify and provide fits to the tidal depletion of subhalos in baryonic relative to dark matter-only simulations as a function of distance. Our simulated profiles imply observational incompleteness in the LG even at Mstar >~ 10^5 Msun: we predict 2-10 such satellites to be discovered around the MW and possibly 6-9 around M31. To provide cosmological context, we compare our results with the radial profiles of satellites around MW analogs in the SAGA survey, finding that our simulations are broadly consistent with most SAGA systems.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

But What About... Cosmic Rays, Magnetic Fields, Conduction, & Viscosity in Galaxy Formation

We present a suite of high-resolution cosmological simulations, using the FIRE-2 feedback physics together with explicit treatment of magnetic fields, anisotropic conduction and viscosity, and cosmic rays (CRs) injected by supernovae (including anisotropic diffusion, streaming, adiabatic, hadronic and Coulomb losses). We survey systems from ultra-faint dwarf ($M_{\ast}\sim 10^{4}\,M_{\odot}$, $M_{\rm halo}\sim 10^{9}\,M_{\odot}$) through Milky Way masses, systematically vary CR parameters (e.g. the diffusion coefficient $κ$ and streaming velocity), and study an ensemble of galaxy properties (masses, star formation histories, mass profiles, phase structure, morphologies). We confirm previous conclusions that magnetic fields, conduction, and viscosity on resolved ($\gtrsim 1\,$pc) scales have small effects on bulk galaxy properties. CRs have relatively weak effects on all galaxy properties studied in dwarfs ($M_{\ast} \ll 10^{10}\,M_{\odot}$, $M_{\rm halo} \lesssim 10^{11}\,M_{\odot}$), or at high redshifts ($z\gtrsim 1-2$), for any physically-reasonable parameters. However at higher masses ($M_{\rm halo} \gtrsim 10^{11}\,M_{\odot}$) and $z\lesssim 1-2$, CRs can suppress star formation by factors $\sim 2-4$, given relatively high effective diffusion coefficients $κ\gtrsim 3\times10^{29}\,{\rm cm^{2}\,s^{-1}}$. At lower $κ$, CRs take too long to escape dense star-forming gas and lose energy to hadronic collisions, producing negligible effects on galaxies and violating empirical constraints from $γ$-ray emission. But around $κ\sim 3\times10^{29}\,{\rm cm^{2}\,s^{-1}}$, CRs escape the galaxy and build up a CR-pressure-dominated halo which supports dense, cool ($T\ll 10^{6}$ K) gas that would otherwise rain onto the galaxy. CR heating (from collisional and streaming losses) is never dominant.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Cataloging Accreted Stars within Gaia DR2 using Deep Learning

The goal of this study is to present the development of a machine learning based approach that utilizes phase space alone to separate the Gaia DR2 stars into two categories: those accreted onto the Milky Way from those that are in situ. Traditional selection methods that have been used to identify accreted stars typically rely on full 3D velocity, metallicity information, or both, which significantly reduces the number of classifiable stars. The approach advocated here is applicable to a much larger portion of Gaia DR2. A method known as "transfer learning" is shown to be effective through extensive testing on a set of mock Gaia catalogs that are based on the FIRE cosmological zoom-in hydrodynamic simulations of Milky Way-mass galaxies. The machine is first trained on simulated data using only 5D kinematics as inputs and is then further trained on a cross-matched Gaia/RAVE data set, which improves sensitivity to properties of the real Milky Way. The result is a catalog that identifies around 767,000 accreted stars within Gaia DR2. This catalog can yield empirical insights into the merger history of the Milky Way and could be used to infer properties of the dark matter distribution.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Synthetic Gaia surveys from the FIRE cosmological simulations of Milky Way-mass galaxies

With Gaia Data Release 2, the astronomical community is entering a new era of multidimensional surveys of the Milky Way. This new phase-space view of our Galaxy demands new tools for comparing observations to simulations of Milky-Way-mass galaxies in a cosmological context, to test the physics of both dark matter and galaxy formation. We present ananke, a framework for generating synthetic phase-space surveys from high-resolution baryonic simulations, and use it to generate a suite of synthetic surveys resembling Gaia DR2 in data structure, magnitude limits, and observational errors. We use three cosmological simulations of Milky-Way-mass galaxies from the Latte suite of the Feedback In Realistic Environments (FIRE) project, which feature self-consistent clustering of star formation in dense molecular clouds and thin stellar/gaseous disks in live cosmological halos with satellite dwarf galaxies and stellar halos. We select three solar viewpoints from each simulation to generate nine synthetic Gaia-like surveys. We sample synthetic stars by assuming each star particle (of mass 7070 $M_{\odot}$) represents a single stellar population. At each viewpoint, we compute dust extinction from the simulated gas metallicity distribution and apply a simple error model to produce a synthetic Gaia-like survey that includes both observational properties and a pointer to the generating star particle. We provide the complete simulation snapshot at $z = 0$ for each simulated galaxy. We describe data access points, the data model, and plans for future upgrades. These synthetic surveys provide a tool for the scientific community to test analysis methods and interpret Gaia data.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

The Edge of the Galaxy

We use cosmological simulations of isolated Milky Way-mass galaxies, as well as Local Group analogues, to define the "edge" -- a caustic manifested in a drop in density or radial velocity -- of Galactic-sized haloes, both in dark matter and in stars. In the dark matter, we typically identify two caustics: the outermost caustic located at ~1.4r_200m corresponding to the "splashback" radius, and a second caustic located at ~0.6r_200m which likely corresponds to the edge of the virialized material which has completed at least two pericentric passages. The splashback radius is ill defined in Local Group type environments where the halos of the two galaxies overlap. However, the second caustic is less affected by the presence of a companion, and is a more useful definition for the boundary of the Milky Way halo. Curiously, the stellar distribution also has a clearly defined caustic, which, in most cases, coincides with the second caustic of the dark matter. This can be identified in both radial density and radial velocity profiles, and should be measurable in future observational programmes. Finally, we show that the second caustic can also be identified in the phase-space distribution of dwarf galaxies in the Local Group. Using the current dwarf galaxy population, we predict the edge of the Milky Way halo to be 292 +/- 61 kpc.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Organized Chaos: Scatter in the relation between stellar mass and halo mass in small galaxies

We use Local Group galaxy counts together with the ELVIS N-body simulations to explore the relationship between the scatter and slope in the stellar mass vs. halo mass relation at low masses, $M_\star \simeq 10^5 - 10^8 M_\odot$. Assuming models with log-normal scatter about a median relation of the form $M_\star \propto M_\mathrm{halo}^α$, the preferred log-slope steepens from $α\simeq 1.8$ in the limit of zero scatter to $α\simeq 2.6$ in the case of $2$ dex of scatter in $M_\star$ at fixed halo mass. We provide fitting functions for the best-fit relations as a function of scatter, including cases where the relation becomes increasingly stochastic with decreasing mass. We show that if the scatter at fixed halo mass is large enough ($\gtrsim 1$ dex) and if the median relation is steep enough ($α\gtrsim 2$), then the "too-big-to-fail" problem seen in the Local Group can be self-consistently eliminated in about $\sim 5-10\%$ of realizations. This scenario requires that the most massive subhalos host unobservable ultra-faint dwarfs fairly often; we discuss potentially observable signatures of these systems. Finally, we compare our derived constraints to recent high-resolution simulations of dwarf galaxy formation in the literature. Though simulation-to-simulation scatter in $M_\star$ at fixed $M_\mathrm{halo}$ is large among separate authors ($\sim 2$ dex), individual codes produce relations with much less scatter and usually give relations that would over-produce local galaxy counts.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Under Pressure: Quenching Star Formation in Low-Mass Satellite Galaxies via Stripping

Recent studies of galaxies in the local Universe, including those in the Local Group, find that the efficiency of environmental (or satellite) quenching increases dramatically at satellite stellar masses below ~ $10^8\ {\rm M}_{\odot}$. This suggests a physical scale where quenching transitions from a slow "starvation" mode to a rapid "stripping" mode at low masses. We investigate the plausibility of this scenario using observed HI surface density profiles for a sample of 66 nearby galaxies as inputs to analytic calculations of ram-pressure and viscous stripping. Across a broad range of host properties, we find that stripping becomes increasingly effective at $M_{*} < 10^{8-9}\ {\rm M}_{\odot}$, reproducing the critical mass scale observed. However, for canonical values of the circumgalactic medium density ($n_{\rm halo} < 10^{-3.5}$ ${\rm cm}^{-3}$), we find that stripping is not fully effective; infalling satellites are, on average, stripped of < 40 - 70% of their cold gas reservoir, which is insufficient to match observations. By including a host halo gas distribution that is clumpy and therefore contains regions of higher density, we are able to reproduce the observed HI gas fractions (and thus the high quenched fraction and short quenching timescale) of Local Group satellites, suggesting that a host halo with clumpy gas may be crucial for quenching low-mass systems in Local Group-like (and more massive) host halos.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

When and where did GW150914 form?

The recent LIGO detection of gravitational waves (GW150914), likely originating from the merger of two $\sim 30 M_\odot$ black holes suggests progenitor stars of low metallicity ($[Z/Z_\odot] \lesssim 0.3$), constraining when and where the progenitor of GW150914 may have formed. We combine estimates of galaxy properties (metallicity, star formation rate and merger rate) across cosmic time to predict the low redshift black hole - black hole merger rate as a function of present day host galaxy mass, $M_\mathrm{gal}$, and the formation redshift of the progenitor system $z_\mathrm{form}$ for different progenitor metallicities $Z_\mathrm{c}$. At $Z_\mathrm{c}=0.1 Z_\odot$, the signal is dominated by binaries in massive galaxies with $z_\mathrm{form}\simeq 2$, with a small contribution from binaries formed around $z_\mathrm{form}\simeq 0.5$ in dwarf galaxies. For $Z_\mathrm{c}=0.01Z_\odot$, fast mergers are possible and very recent star formation in dwarfs likely dominates. Additional gravitational wave detections from merging massive black holes will provide constraints on the mass-metallicity relation and massive star formation at high redshifts.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Properties of Resonantly Produced Sterile Neutrino Dark Matter Subhalos

The anomalous 3.55 keV X-ray line recently detected towards a number of massive dark matter objects may be interpreted as the radiative decays of 7.1 keV mass sterile neutrino dark matter. Depending on its parameters, the sterile neutrino can range from cold to warm dark matter with small-scale suppression that differs in form from commonly-adopted thermal warm dark matter. Here, we numerically investigate the subhalo properties for 7.1 keV sterile neutrino dark matter produced via the resonant Shi-Fuller mechanism. Using accurate matter power spectra, we run cosmological zoom-in simulations of a Milky Way-sized halo and explore the abundance of massive subhalos, their radial distributions, and their internal structure. We also simulate the halo with thermal 2.0 keV warm dark matter for comparison and discuss quantitative differences. We find that the resonantly produced sterile neutrino model for the 3.55 keV line provides a good description of structures in the Local Group, including the number of satellite dwarf galaxies and their radial distribution, and largely mitigates the too-big-to-fail problem. Future searches for satellite galaxies by deep surveys, such as the Dark Energy Survey, Large Synoptic Survey Telescope, and Wide Field Infrared Survey Telescope, will be a strong direct test of warm dark matter scenarios.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Resonant Sterile Neutrino Dark Matter in the Local and High-z Universe

Sterile neutrinos comprise an entire class of dark matter models that, depending on their production mechanism, can be hot, warm, or cold dark matter. We simulate the Local Group and representative volumes of the Universe in a variety of sterile neutrino models, all of which are consistent with the possible existence of a radiative decay line at ~3.5 keV. We compare models of production via resonances in the presence of a lepton asymmetry (suggested by Shi & Fuller 1999) to "thermal" models. We find that properties in the highly nonlinear regime - e.g., counts of satellites and internal properties of halos and subhalos - are insensitive to the precise fall-off in power with wavenumber, indicating that nonlinear evolution essentially washes away differences in the initial (linear) matter power spectrum. In the quasi-linear regime at higher redshifts, however, quantitative differences in the 3D matter power spectra remain, raising the possibility that such models can be tested with future observations of the Lyman-alpha forest. While many of the sterile neutrino models largely eliminate multiple small-scale issues within the Cold Dark Matter (CDM) paradigm, we show that these models may be ruled out in the near future via discoveries of additional dwarf satellites in the Local Group.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Satellite Dwarf Galaxies in a Hierarchical Universe: Infall Histories, Group Preprocessing, and Reionization

In the Local Group (LG), almost all satellite dwarf galaxies that are within the virial radius of the Milky Way (MW) and Andromeda (M31) exhibit strong environmental influence. The orbital histories of these satellites provide the key to understanding the role of the MW/M31 halo, lower-mass groups, and cosmic reionization on the evolution of dwarf galaxies. We examine the virial-infall histories of satellites with M_star = 10^{3-9} M_sun using the ELVIS suite of cosmological zoom-in dissipationless simulations of 48 MW/M31-like halos. Satellites at z = 0 fell into the MW/M31 halos typically 5 - 8 Gyr ago at z = 0.5 - 1. However, they first fell into any host halo typically 7 - 10 Gyr ago at z = 0.7 - 1.5. This difference arises because many satellites experienced "group preprocessing" in another host halo, typically of M_vir ~ 10^{10-12} M_sun, before falling into the MW/M31 halos. Satellites with lower mass and/or those closer to the MW/M31 fell in earlier and are more likely to have experienced group preprocessing; half of all satellites with M_star < 10^6 M_sun were preprocessed in a group. Infalling groups also drive most satellite-satellite mergers within the MW/M31 halos. Finally, none of the surviving satellites at z = 0 were within the virial radius of their MW/M31 halo during reionization (z > 6), and only < 4% were satellites of any other host halo during reionization. Thus, effects of cosmic reionization versus host-halo environment on the formation histories of surviving dwarf galaxies in the LG occurred at distinct epochs, separated typically by 2 - 4 Gyr, so they are separable theoretically and, in principle, observationally.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Satellites of LMC-Mass Dwarfs: Close Friendships Ruined by Milky Way Mass Halos

Motivated by the recent discovery of several dwarfs near the Large Magellanic Cloud (LMC), we study the accretion of massive satellites onto Milky Way (MW)/M31-like halos using the ELVIS suite of N-body simulations. We identify 25 surviving LMC-mass subhalos, and investigate the lower-mass satellites that were associated with these subhalos before they fell into the MW/M31 halos. Typically, 7% of the overall z=0 satellite population of MW/M31 halos were in a surviving LMC-group before falling into the MW/M31 halo. This fraction can vary between 1% and 25%, being higher for groups with higher mass and/or more recent infall times. Groups of satellites disperse rapidly in phase space after infall, and their distances and velocities relative to the group center become statistically similar to the overall satellite population after 4-8 Gyr. We quantify the likelihood that satellites were associated with an LMC-mass group as a function of both distance and velocity relative to the LMC at z=0. The close proximity in distance of the nine Dark Energy Survey candidate dwarf galaxies to the LMC suggest that ~2-4 are likely associated with the LMC. Furthermore, if several of these dwarfs are genuine members, then the LMC-group probably fell into the MW very recently, <2 Gyr ago. If the connection with the LMC is established with follow-up velocity measurements, these "satellites of satellites" represent prime candidates to study the affects of group pre-processing on lower mass dwarfs.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Sweating the small stuff: simulating dwarf galaxies, ultra-faint dwarf galaxies, and their own tiny satellites

We present FIRE/Gizmo hydrodynamic zoom-in simulations of isolated dark matter halos, two each at the mass of classical dwarf galaxies ($M_{\rm vir} \simeq 10^{10} M_{\odot}$) and ultra-faint galaxies ($M_{\rm vir} \simeq 10^9 M_{\odot}$), and with two feedback implementations. The resultant central galaxies lie on an extrapolated abundance matching relation from $M_{\star} \simeq 10^6$ to $10^4 M_{\odot}$ without a break. Every host is filled with subhalos, many of which form stars. Our dwarfs with $M_{\star} \simeq 10^6 M_{\odot}$ each have 1-2 well-resolved satellites with $M_{\star} = 3-200 \times 10^3 M_{\odot}$. Even our isolated ultra-faint galaxies have star-forming subhalos. If this is representative, dwarf galaxies throughout the universe should commonly host tiny satellite galaxies of their own. We combine our results with the ELVIS simulations to show that targeting $\sim 50~ \rm kpc$ regions around nearby isolated dwarfs could increase the chances of discovering ultra-faint galaxies by $\sim 35\%$ compared to random halo pointings, and specifically identify the region around the Phoenix dwarf galaxy as a good potential target. The well-resolved ultra-faint galaxies in our simulations ($M_{\star} \simeq 3 - 30 \times 10^3 M_{\odot}$) form within $M_{\rm peak} \simeq 0.5 - 3 \times 10^9 M_{\odot}$ halos. Each has a uniformly ancient stellar population ($ > 10~ \rm Gyr$) owing to reionization-related quenching. More massive systems, in contrast, all have late-time star formation. Our results suggest that $M_{\rm halo} \simeq 5 \times 10^9 M_{\odot}$ is a probable dividing line between halos hosting reionization "fossils" and those hosting dwarfs that can continue to form stars in isolation after reionization.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Taking Care of Business in a Flash: Constraining the Timescale for Low-Mass Satellite Quenching with ELVIS

The vast majority of dwarf satellites orbiting the Milky Way and M31 are quenched, while comparable galaxies in the field are gas-rich and star-forming. Assuming that this dichotomy is driven by environmental quenching, we use the ELVIS suite of N-body simulations to constrain the characteristic timescale upon which satellites must quench following infall into the virial volumes of their hosts. The high satellite quenched fraction observed in the Local Group demands an extremely short quenching timescale (~ 2 Gyr) for dwarf satellites in the mass range Mstar ~ 10^6-10^8 Msun. This quenching timescale is significantly shorter than that required to explain the quenched fraction of more massive satellites (~ 8 Gyr), both in the Local Group and in more massive host halos, suggesting a dramatic change in the dominant satellite quenching mechanism at Mstar < 10^8 Msun. Combining our work with the results of complementary analyses in the literature, we conclude that the suppression of star formation in massive satellites (Mstar ~ 10^8 - 10^11 Msun) is broadly consistent with being driven by starvation, such that the satellite quenching timescale corresponds to the cold gas depletion time. Below a critical stellar mass scale of ~ 10^8 Msun, however, the required quenching times are much shorter than the expected cold gas depletion times. Instead, quenching must act on a timescale comparable to the dynamical time of the host halo. We posit that ram-pressure stripping can naturally explain this behavior, with the critical mass (of Mstar ~ 10^8 Msun) corresponding to halos with gravitational restoring forces that are too weak to overcome the drag force encountered when moving through an extended, hot circumgalactic medium.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Core Formation in Dwarf Halos with Self Interacting Dark Matter: No Fine-Tuning Necessary

We investigate the effect of self-interacting dark matter (SIDM) on the density profiles of $V_{\rm max} \simeq 40~km~s^{-1}$ isolated dwarf dark matter halos -- the scale of relevance for the too big to fail problem (TBTF) -- using very high-resolution cosmological zoom simulations. Each halo has millions of particles within its virial radius. We find that SIDM models with cross sections per unit mass spanning the range σ/m = $0.5 - 50$ $cm^2~g^{-1}$ alleviate TBTF and produce constant density cores of size 300-1000 pc, comparable to the half-light radii of $M_\star$ ~ $10^{5-7}$ $M_\odot$ dwarfs. The largest, lowest density cores develop for cross sections in the middle of this range, σ/m ~ $5-10~cm^2~g^{-1}$. Our largest SIDM cross section run (σ/m = $50~cm^2~g^{-1}$) develops a slightly denser core owing to mild core-collapse behavior, but it remains less dense than the CDM case and retains a constant density core profile. Our work suggests that SIDM cross sections as large or larger than $50~cm^2~g^{-1}$ remain viable on velocity scales of dwarf galaxies ($v_{\rm rms}$ ~ $40~km~s^{-1}$). The range of SIDM cross sections that alleviate TBTF and the cusp/core problem spans at least two orders of magnitude and therefore need not be particularly fine-tuned.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Cosmological Simulations of Decaying Dark Matter: Implications for Small-scale Structure of Dark Matter Halos

We present a set of N-body simulations of a class of models in which an unstable dark matter particle decays into a stable non-interacting dark matter particle, with decay lifetime comparable to the Hubble time. We study the effects of the kinematic recoil velocity received by the stable dark matter on the structures of dark matter halos ranging from galaxy-cluster to Milky Way mass scales. For Milky Way-mass halos, we use high-resolution, zoom-in simulations to explore the effects of decays on Galactic substructure. In general, halos with circular velocities comparable to the magnitude of kick velocity are most strongly affected by decays. We show that decaying dark matter models with lifetimes comparable to Hubble time and recoil speeds about 20-40 km/s can significantly reduce both the abundance of Galactic subhalos and the internal densities of the subhalos. We also compare subhalo circular velocity profiles with observational constraints on the Milky Way dwarf satellite galaxies. Interestingly, we find that decaying dark matter models that do not violate current astrophysical constraints, can significantly mitigate both the well-documented "missing satellites problem" and the more recent "too big to fail problem" associated with the abundances and densities of Local Group dwarf satellite galaxies. A relatively unique feature of late decaying dark matter models is that they predict significant evolution of halos as a function of time. This is an important consideration because at high redshifts, prior to decays, decaying models exhibit the same sequence of structure formation as cold dark matter. We conclude that models of decaying dark matter make predictions that are relevant for the interpretation of observations of small galaxies in the Local Group and can be tested or constrained by the kinematics of Local Group dwarf galaxies as well as by forthcoming large-scale surveys.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

ELVIS: Exploring the Local Volume in Simulations

We introduce a set of high-resolution dissipationless simulations that model the Local Group (LG) in a cosmological context: Exploring the Local Volume in Simulations (ELVIS). The suite contains 48 Galaxy-size halos, each within high-resolution volumes that span 2-5 Mpc in size, and each resolving thousands of systems with masses below the atomic cooling limit. Half of the ELVIS galaxy halos are in paired configurations similar to the Milky Way (MW) and M31; the other half are isolated, mass-matched analogs. We find no difference in the abundance or kinematics of substructure within the virial radii of isolated versus paired hosts. On Mpc scales, however, LG-like pairs average almost twice as many companions and the velocity field is kinematically hotter and more complex. We present a refined abundance matching relation between stellar mass and halo mass that reproduces the observed satellite stellar mass functions of the MW and M31 down to the regime where incompleteness is an issue, $M_\star \sim 5\times 10^5 \, M_\odot$. Within a larger region spanning approximately 3 Mpc, the same relation predicts that there should be $\sim$ 1000 galaxies with $M_\star > 10^{3}\,M_\odot$ awaiting discovery. We show that up to 50% of halos within 1 Mpc of the MW or M31 could be systems that have previously been within the virial radius of either giant. By associating never-accreted halos with gas-rich dwarfs, we show that there are plausibly 50 undiscovered dwarf galaxies with HI masses $> 10^5\,M_\odot$ within the Local Volume. The radial velocity distribution of these predicted gas-rich dwarfs can be used to inform follow-up searches based on ultra-compact high-velocity clouds found in the ALFALFA survey.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Near-Field Limits on the Role of Faint Galaxies in Cosmic Reionization

Reionizing the Universe with galaxies appears to require significant star formation in low-mass halos at early times, while local dwarf galaxy counts tell us that star formation has been minimal in small halos around us today. Using simple models and the ELVIS simulation suite, we show that reionization scenarios requiring appreciable star formation in halos with $M_{\rm vir} \approx 10^{8}\,M_{\odot}$ at $z=8$ are in serious tension with galaxy counts in the Local Group. This tension originates from the seemingly inescapable conclusion that 30 - 60 halos with $M_{\rm vir} > 10^{8}\,M_{\odot}$ at $z=8$ will survive to be distinct bound satellites of the Milky Way at $z = 0$. Reionization models requiring star formation in such halos will produce dozens of bound galaxies in the Milky Way's virial volume today (and 100 - 200 throughout the Local Group), each with $\gtrsim 10^{5}\,M_{\odot}$ of old stars ($\gtrsim 13$ Gyr). This exceeds the stellar mass function of classical Milky Way satellites today, even without allowing for the (significant) post-reionization star formation observed in these galaxies. One possible implication of these findings is that star formation became sharply inefficient in halos smaller than $\sim 10^9 \,M_{\odot}$ at early times, implying that the high-$z$ luminosity function must break at magnitudes brighter than is often assumed (at ${\rm M_{UV}} \approx -14$). Our results suggest that JWST (and possibly even HST with the Frontier Fields) may realistically detect the faintest galaxies that drive reionization. It remains to be seen how these results can be reconciled with the most sophisticated simulations of early galaxy formation at present, which predict substantial star formation in $M_{\rm vir} \sim 10^8 \, M_{\odot}$ halos during the epoch of reionization.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Running with BICEP2: Implications for Small-Scale Problems in CDM

The BICEP2 results, when interpreted as a gravitational wave signal and combined with other CMB data, suggest a roll-off in power towards small scales in the primordial matter power spectrum. Among the simplest possibilities is a running of the spectral index. Here we show that the preferred level of running alleviates small-scale issues within the $Λ$CDM model, more so even than viable WDM models. We use cosmological zoom-in simulations of a Milky Way-size halo along with full-box simulations to compare predictions among four separate cosmologies: a BICEP2-inspired running index model ($α_s$ = -0.024), two fixed-tilt $Λ$CDM models motivated by Planck, and a 2.6 keV thermal WDM model. We find that the running BICEP2 model reduces the central densities of large dwarf-size halos ($V_\mathrm{max}$ ~ 30 - 80 km s$^{-1}$) and alleviates the too-big-to-fail problem significantly compared to our adopted Planck and WDM cases. Further, the BICEP2 model suppresses the count of small subhalos by ~50% relative to Planck models, and yields a significantly lower "boost" factor for dark matter annihilation signals. Our findings highlight the need to understand the shape of the primordial power spectrum in order to correctly interpret small-scale data.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Satellite Dwarf Galaxies in a Hierarchical Universe: The Prevalence of Dwarf-Dwarf Major Mergers

Mergers are a common phenomenon in hierarchical structure formation, especially for massive galaxies and clusters, but their importance for dwarf galaxies in the Local Group remains poorly understood. We investigate the frequency of major mergers between dwarf galaxies in the Local Group using the ELVIS suite of cosmological zoom-in dissipationless simulations of Milky Way- and M31-like host halos. We find that ~10% of satellite dwarf galaxies with M_star > 10^6 M_sun that are within the host virial radius experienced a major merger of stellar mass ratio closer than 0.1 since z = 1, with a lower fraction for lower mass dwarf galaxies. Recent merger remnants are biased towards larger radial distance and more recent virial infall times, because most recent mergers occurred shortly before crossing within the virial radius of the host halo. Satellite-satellite mergers also occur within the host halo after virial infall, catalyzed by the large fraction of dwarf galaxies that fell in as part of a group. The merger fraction doubles for dwarf galaxies outside of the host virial radius, so the most distant dwarf galaxies in the Local Group are the most likely to have experienced a recent major merger. We discuss the implications of these results on observable dwarf merger remnants, their star formation histories, the gas content of mergers, and massive black holes in dwarf galaxies.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Too Big to Fail in the Local Group

We compare the dynamical masses of dwarf galaxies in the Local Group (LG) to the predicted masses of halos in the ELVIS suite of $Λ$CDM simulations, a sample of 48 Galaxy-size hosts, 24 of which are in paired configuration similar to the LG. We enumerate unaccounted-for dense halos ($V_\mathrm{max} \gtrsim 25$ km s$^{-1}$) in these volumes that at some point in their histories were massive enough to have formed stars in the presence of an ionizing background ($V_\mathrm{peak} > 30$ km s$^{-1}$). Within 300 kpc of the Milky Way, the number of unaccounted-for massive halos ranges from 2 - 25 over our full sample. Moreover, this "too big to fail" count grows as we extend our comparison to the outer regions of the Local Group: within 1.2 Mpc of either giant we find that there are 12-40 unaccounted-for massive halos. This count excludes volumes within 300 kpc of both the MW and M31, and thus should be largely unaffected by any baryonically-induced environmental processes. According to abundance matching -- specifically abundance matching that reproduces the Local Group stellar mass function -- all of these missing massive systems should have been quite bright, with $M_\star > 10^6M_\odot$. Finally, we use the predicted density structure of outer LG dark matter halos together with observed dwarf galaxy masses to derive an $M_\star-V_\mathrm{max}$ relation for LG galaxies that are outside the virial regions of either giant. We find that there is no obvious trend in the relation over three orders of magnitude in stellar mass (a "common mass" relation), from $M_\star \sim 10^8 - 10^5 M_\odot$, in drastic conflict with the tight relation expected for halos that are unaffected by reionization. Solutions to the too big to fail problem that rely on ram pressure stripping, tidal effects, or statistical flukes appear less likely in the face of these results.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Can Feedback Solve the Too Big to Fail Problem?

The observed central densities of Milky Way dwarf spheroidal galaxies (dSphs) are significantly lower than the densities of the largest (Vmax about 35 km/s) subhalos found in dissipationless simulations of Galaxy-size dark matter hosts. One possible explanation is that gas removal from feedback can lower core densities enough to match observations. We model the dynamical effects of supernova feedback through the use of a time-varying central potential in high resolution, idealized numerical simulations and explore the resulting impact on the mass distributions of dwarf dark matter halos. We find that in order to match the observed central masses of M_star about 10^6 M_sun dSphs, the energy equivalent of more than 40,000 supernovae must be delivered with 100% efficiency directly to the dark matter. This energy requirement exceeds the number of supernovae that have ever exploded in most dSphs for typical initial mass functions. We also find that, per unit energy delivered and per cumulative mass removed from the galaxy, single blow-out events are more effective than repeated small bursts in reducing central dark matter densities. We conclude that it is unlikely that supernova feedback alone can solve the "Too Big to Fail" problem for Milky Way subhalos.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

How to Zoom: Bias, Contamination, and Lagrange Volumes in Multimass Cosmological Simulations

We perform a suite of multimass cosmological zoom simulations of individual dark matter halos and explore how to best select Lagrangian regions for resimulation without contaminating the halo of interest with low-resolution particles. Such contamination can lead to significant errors in the gas distribution of hydrodynamical simulations, as we show. For a fixed Lagrange volume, we find that the chance of contamination increases systematically with the level of zoom. In order to avoid contamination, the Lagrangian volume selected for resimulation must increase monotonically with the resolution difference between parent box and the zoom region. We provide a simple formula for selecting Lagrangian regions (in units of the halo virial volume) as a function of the level of zoom required. We also explore the degree to which a halo's Lagrangian volume correlates with other halo properties (concentration, spin, formation time, shape, etc.) and find no significant correlation. There is a mild correlation between Lagrange volume and environment, such that halos living in the most clustered regions have larger Lagrangian volumes. Nevertheless, selecting halos to be isolated is not the best way to ensure inexpensive zoom simulations. We explain how one can safely choose halos with the smallest Lagrangian volumes, which are the least expensive to resimulate, without biasing one's sample.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

On the stark difference in satellite distributions around the Milky Way and Andromeda

We compare spherically-averaged radial number counts of bright (> 10^5 Lsun) dwarf satellite galaxies within 400 kpc of the Milky Way (MW) and M31 and find that the MW satellites are much more centrally concentrated. Remarkably, the two satellite systems are almost identical within the central 100 kpc, while M31 satellites outnumber MW satellites by about a factor of four at deprojected distances spanning 100 - 400 kpc. We compare the observed distributions to those predicted for LCDM suhbalos using a suite of 44 high-resolution ~10^12 halo zoom simulations, 22 of which are in pairs like the MW and M31. We find that the radial distribution of satellites around M31 is fairly typical of those predicted for subhalos, while the Milky Way's distribution is more centrally concentrated that any of our simulated LCDM halos. One possible explanation is that our census is bright (> 10^5 Lsun) MW dwarf galaxies is significantly incomplete beyond ~ 100 kpc of the Sun. If there were ~8 - 20 more bright dwarfs orbiting undetected at 100 - 400 kpc, then the Milky Way's radial distribution would fall within the range expected from subhalo distributions and alos look very much like the known M31 system. We use our simulations to demonstrate that there is enough area left unexplored by the Sloan Digital Sky Survey and its extensions that the discovery of ~10 new bright dwarfs is not implausible given the expected range of angular anisotropy of subhalos in the sky.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Sterile neutrino dark matter bounds from galaxies of the Local Group

We show that the canonical oscillation-based (non-resonant) production of sterile neutrino dark matter is inconsistent at $>99$% confidence with observations of galaxies in the Local Group. We set lower limits on the non-resonant sterile neutrino mass of $2.5$ keV (equivalent to $0.7$ keV thermal mass) using phase-space densities derived for dwarf satellite galaxies of the Milky Way, as well as limits of $8.8$ keV (equivalent to $1.8$ keV thermal mass) based on subhalo counts of $N$-body simulations of M 31 analogues. Combined with improved upper mass limits derived from significantly deeper X-ray data of M 31 with full consideration for background variations, we show that there remains little room for non-resonant production if sterile neutrinos are to explain $100$% of the dark matter abundance. Resonant and non-oscillation sterile neutrino production remain viable mechanisms for generating sufficient dark matter sterile neutrinos.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Cosmological Simulations with Self-Interacting Dark Matter I: Constant Density Cores and Substructure

We use cosmological simulations to study the effects of self-interacting dark matter (SIDM) on the density profiles and substructure counts of dark matter halos from the scales of spiral galaxies to galaxy clusters, focusing explicitly on models with cross sections over dark matter particle mass σ/m = 1 and 0.1 cm^2/g. Our simulations rely on a new SIDM N-body algorithm that is derived self-consistently from the Boltzmann equation and that reproduces analytic expectations in controlled numerical experiments. We find that well-resolved SIDM halos have constant-density cores, with significantly lower central densities than their CDM counterparts. In contrast, the subhalo content of SIDM halos is only modestly reduced compared to CDM, with the suppression greatest for large hosts and small halo-centric distances. Moreover, the large-scale clustering and halo circular velocity functions in SIDM are effectively identical to CDM, meaning that all of the large-scale successes of CDM are equally well matched by SIDM. From our largest cross section runs we are able to extract scaling relations for core sizes and central densities over a range of halo sizes and find a strong correlation between the core radius of an SIDM halo and the NFW scale radius of its CDM counterpart. We construct a simple analytic model, based on CDM scaling relations, that captures all aspects of the scaling relations for SIDM halos. Our results show that halo core densities in σ/m = 1 cm^2/g models are too low to match observations of galaxy clusters, low surface brightness spirals (LSBs), and dwarf spheroidal galaxies. However, SIDM with σ/m ~ 0.1 cm^2/g appears capable of reproducing reported core sizes and central densities of dwarfs, LSBs, and galaxy clusters without the need for velocity dependence. (abridged)

0 citations0 reviewsNo code linkedOpen access