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

James S. Bullock

James S. Bullock 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.COastro-phhep-phastro-ph.HEastro-ph.IMastro-ph.SRastro-ph.EPgr-qc

Catalog footprint

What is connected

91works
9topics
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

91 published item(s)

preprint2024arXiv

An analytic surface density profile for $Λ$CDM halos and gravitational lensing studies

We introduce an analytic surface density profile for dark matter halos that accurately reproduces the structure of simulated halos of mass $M_{\rm vir} = 10^{7-11}\ M_\odot$, making it useful for modeling line-of-sight perturbers in strong gravitational lensing models. The two-parameter function has an analytic deflection potential and is more accurate than the projected Navarro, Frenk & White (NFW) profile commonly adopted at this mass scale for perturbers, especially at the small radii of most relevant for lensing perturbations. Using a characteristic radius, $R_{-1}$, where the log slope of surface density is equal to $-1$, and an associated surface density, $Σ_{-1}$, we can represent the expected lensing signal from line-of-sight halos statistically, for an ensemble of halo orientations, using a distribution of {\em projected concentration} parameters, $\mathcal{C}_{\rm vir} := r_{\rm vir}/ R_{-1}$. Though an individual halo can have a projected concentration that varies with orientation with respect to the observer, the range of projected concentrations correlates with the usual three-dimensional halo concentration in a way that enables ease of use.

0 citations0 reviewsNo code linkedOpen access
preprint2023arXiv

The JWST Resolved Stellar Populations Early Release Science Program II. Survey Overview

We present the JWST Resolved Stellar Populations Early Release Science (ERS) science program. We obtained 27.5 hours of NIRCam and NIRISS imaging of three targets in the Local Group (Milky Way globular cluster M92, ultra-faint dwarf galaxy Draco II, star-forming dwarf galaxy WLM), which span factors of $\sim10^5$ in luminosity, $\sim10^4$ in distance, and $\sim10^5$ in surface brightness. We describe the survey strategy, scientific and technical goals, implementation details, present select NIRCam color-magnitude diagrams (CMDs), and validate the NIRCam exposure time calculator (ETC). Our CMDs are among the deepest in existence for each class of target. They touch the theoretical hydrogen burning limit in M92 ($<0.08$ $M_{\odot}$; SNR $\sim5$ at $m_{F090W}\sim28.2$; $M_{F090W}\sim+13.6$), include the lowest-mass stars observed outside the Milky Way in Draco II (0.09 $M_{\odot}$; SNR $=10$ at $m_{F090W}\sim29$; $M_{F090W}\sim+12.1$), and reach $\sim1.5$ magnitudes below the oldest main sequence turnoff in WLM (SNR $=10$ at $m_{F090W}\sim29.5$; $M_{F090W}\sim+4.6$). The PARSEC stellar models provide a good qualitative match to the NIRCam CMDs, though are $\sim0.05$ mag too blue compared to M92 F090W$-$F150W data. The NIRCam ETC (v2.0) matches the SNRs based on photon noise from DOLPHOT stellar photometry in uncrowded fields, but the ETC may not be accurate in more crowded fields, similar to what is known for HST. We release beta versions of DOLPHOT NIRCam and NIRISS modules to the community. Results from this ERS program will establish JWST as the premier instrument for resolved stellar populations studies for decades to come.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

Amplified J-factors in the Galactic Center for velocity-dependent darkmatter annihilation in FIRE simulations

We use FIRE-2 zoom cosmological simulations of Milky Way size galaxy halos to calculate astrophysical J-factors for dark matter annihilation and indirect detection studies. In addition to velocity-independent (s-wave) annihilation cross sections $σ_v$, we also calculate effective J-factors for velocity-dependent models, where the annihilation cross section is either either p-wave ($\propto v^2/c^2$) or d-wave ($\propto v^4/c^4$). We use 12 pairs of simulations, each run with dark-matter-only (DMO) physics and FIRE-2 physics. We observe FIRE runs produce central dark matter velocity dispersions that are systematically larger than in DMO runs by factors of $\sim 2.5-4$. They also have a larger range of central ($\sim 400$ pc) dark matter densities than the DMO runs ($ρ_{\rm FIRE}/ρ_{\rm DMO} \simeq 0.5 - 3$) owing to the competing effects of baryonic contraction and feedback. At 3 degrees from the Galactic Center, FIRE J-factors are $5-50$ (p-wave) and $15-500$ (d-wave) times higher than in the DMO runs. The change in s-wave signal at 3 degrees is more modest and can be higher or lower ($\sim 0.3-6$), though the shape of the emission profile is flatter (less peaked towards the Galactic Center) and more circular on the sky in FIRE runs. Our results for s-wave are broadly consistent with the range of assumptions in most indirect detection studies. We observe p-wave J-factors that are significantly enhanced compared to most past estimates. We find that thermal models with p-wave annihilation may be within range of detection in the near future.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

Comparing Implementations of Self-Interacting Dark Matter in the Gizmo and Arepo Codes

Self-interacting dark matter (SIDM) models have received great attention over the past decade as solutions to the small-scale puzzles of astrophysics. Though there are different implementations of dark matter (DM) self-interactions in N-body codes of structure formation, there has not been a systematic study to compare the predictions of these different implementations. We investigate the implementation of dark matter self-interactions in two simulation codes: Gizmo and Arepo. We begin with identical initial conditions for an isolated $10^{10}$ M$_\odot$ dark matter halo and investigate the evolution of the density and velocity dispersion profiles in Gizmo and Arepo for SIDM cross-section over mass of 1, 5, and 50 $\rm cm^2 g^{-1}$. Our tests are restricted to the core expansion phase where the core density decreases and core radius increases with time. We find better than 30% agreement between the codes for the density profile in this phase of evolution, with the agreement improving at higher resolution. We find that varying code-specific SIDM parameters changes the central halo density by less than 10% outside of the convergence radius. We argue that SIDM core formation is robust across the two different schemes and conclude that these codes can reliably differentiate between cross-sections of 1, 5, and 50 $\rm cm^2 g^{-1}$ but finer distinctions would require further investigation.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

Galaxies lacking dark matter produced by close encounters in a cosmological simulation

The standard cold dark matter plus cosmological constant model predicts that galaxies form within dark-matter haloes, and that low-mass galaxies are more dark-matter dominated than massive ones. The unexpected discovery of two low-mass galaxies lacking dark matter immediately provoked concerns about the standard cosmology and ignited explorations of alternatives, including self-interacting dark matter and modified gravity. Apprehension grew after several cosmological simulations using the conventional model failed to form adequate numerical analogues with comparable internal characteristics (stellar masses, sizes, velocity dispersions and morphologies). Here we show that the standard paradigm naturally produces galaxies lacking dark matter with internal characteristics in agreement with observations. Using a state-of-the-art cosmological simulation and a meticulous galaxy-identification technique, we find that extreme close encounters with massive neighbours can be responsible for this. We predict that approximately 30 percent of massive central galaxies (with at least 1e11 solar masses in stars) harbour at least one dark-matter-deficient satellite (with 1e8 - 1e9 solar masses in stars). This distinctive class of galaxies provides an additional layer in our understanding of the role of interactions in shaping galactic properties. Future observations surveying galaxies in the aforementioned regime will provide a crucial test of this scenario.

0 citations0 reviewsNo code linkedOpen access
preprint2021arXiv

Out of sight, out of mind? The impact of correlated clustering in substructure lensing

A promising route for revealing the existence of dark matter structures on mass scales smaller than the faintest galaxies is through their effect on strong gravitational lenses. We examine the role of local, lens-proximate clustering in boosting the lensing probability relative to contributions from substructure and unclustered line-of-sight (LOS) halos. Using two cosmological simulations that can resolve halo masses of $M_{\rm halo} \simeq 10^{9}\ M_{\odot}$ (in a simulation box of length $L_{\rm box}{\sim}100\,{\rm Mpc}$) and $10^{7}\ M_{\odot}$ ($L_{\rm box}\sim20\,{\rm Mpc}$), we demonstrate that clustering in the vicinity of the lens host produces a clear enhancement relative to an assumption of unclustered halos that persists to $> 20\,R_{\rm vir}$. This enhancement exceeds estimates that use a two-halo term to account for clustering, particularly within $2-5\,R_{\rm vir}$. We provide an analytic expression for this excess, clustered contribution. We find that local clustering boosts the expected count of $10^9 \ M_\odot$ perturbing halos by ${\sim}35\%$ compared to substructure alone, a result that will significantly enhance expected signals for low-redshift ($z_l \simeq 0.2$) lenses, where substructure contributes substantially compared to LOS halos. We also find that the orientation of the lens with respect to the line of sight (e.g., whether the line of sight passes through the major axis of the lens) can also have a significant effect on the lensing signal, boosting counts by an additional $\sim 50\%$ compared to a random orientations. This could be important if discovered lenses are biased to be oriented along their principal axis.

0 citations0 reviewsNo code linkedOpen access
preprint2021arXiv

The effects of LMC-mass environments on their dwarf satellite galaxies in the FIRE simulations

Characterizing the predicted environments of dwarf galaxies like the Large Magellanic Cloud (LMC) is becoming increasingly important as next generation surveys push sensitivity limits into this low-mass regime at cosmological distances. We study the environmental effects of LMC-mass halos ($M_{200m} \sim 10^{11}$ M$_\odot$) on their populations of satellites ($M_\star \geq 10^4$ M$_\odot$) using a suite of zoom-in simulations from the Feedback In Realistic Environments (FIRE) project. Our simulations predict significant hot coronas with $T\sim10^5$ K and $M_\text{gas}\sim10^{9.5}$ M$_\odot$. We identify signatures of environmental quenching in dwarf satellite galaxies, particularly for satellites with intermediate mass ($M_\star = 10^{6-7}$ M$_\odot$). The gas content of such objects indicates ram-pressure as the likely quenching mechanism, sometimes aided by star formation feedback. Satellites of LMC-mass hosts replicate the stellar mass dependence of the quiescent fraction found in satellites of MW mass hosts (i.e. that the quiescent fraction increases as stellar mass decreases). Satellites of LMC-mass hosts have a wider variety of quenching times when compared to the strongly bi-modal distribution of quenching times of nearby centrals. Finally, we identify significant tidal stellar structures around four of our six LMC-analogs, suggesting that stellar streams may be common. These tidal features originated from satellites on close orbits, extend to $\sim$80 kpc from the central galaxy, and contain $\sim10^{6-7}$ M$_\odot$~of stars.

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

A Relationship Between Stellar Metallicity Gradients and Galaxy Age in Dwarf Galaxies

We explore the origin of stellar metallicity gradients in simulated and observed dwarf galaxies. We use FIRE-2 cosmological baryonic zoom-in simulations of 26 isolated galaxies as well as existing observational data for 10 Local Group dwarf galaxies. Our simulated galaxies have stellar masses between $10^{5.5}$ and $10^{8.6} \msun$. Whilst gas-phase metallicty gradients are generally weak in our simulated galaxies, we find that stellar metallicity gradients are common, with central regions tending to be more metal-rich than the outer parts. The strength of the gradient is correlated with galaxy-wide median stellar age, such that galaxies with younger stellar populations have flatter gradients. Stellar metallicty gradients are set by two competing processes: (1) the steady "puffing" of old, metal-poor stars by feedback-driven potential fluctuations, and (2) the accretion of extended, metal-rich gas at late times, which fuels late-time metal-rich star formation. If recent star formation dominates, then extended, metal-rich star formation washes out pre-existing gradients from the "puffing" process. We use published results from ten Local Group dwarf galaxies to show that a similar relationship between age and stellar metallicity-gradient strength exists among real dwarfs. This suggests that observed stellar metallicity gradients may be driven largely by the baryon/feedback cycle rather than by external environmental effects.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Accurate mass estimates from the proper motions of dispersion-supported galaxies

Starting with the spherical Jeans equation, we show that there exists a radius where the mass enclosed depends only on the projected tangential velocity dispersion, assuming that the anisotropy profile slowly varies. This is well-approximated at the radius where the log-slope of the stellar tracer profile is $-2$: $r_{-2}$. The associated mass is $M(r_{-2}) = 2 G^{-1} \langle σ_{\mathcal{T}}^{2}\rangle^{*} r_{-2}$ and the circular velocity is $V^{2}({r_{-2}}) = 2\langle σ_{\mathcal{T}}^{2}\rangle^{*}$. For a Plummer profile $r_{-2} \simeq 4R_e/5$. Importantly, $r_{-2}$ is smaller than the characteristic radius for line-of-sight velocities derived by Wolf et al. 2010. Together, the two estimators can constrain the mass profiles of dispersion-supported galaxies. We illustrate its applicability using published proper motion measurements of dwarf galaxies Draco and Sculptor, and find that they are consistent with inhabiting cuspy NFW subhalos of the kind predicted in CDM but we cannot rule out a core. We test our combined mass estimators against previously-published, non-spherical cosmological dwarf galaxy simulations done in both CDM and SIDM. For CDM, the estimates for the dynamic rotation curves are found to be accurate to $10\%$ while SIDM are accurate to $15\%$. Unfortunately, this level of accuracy is not good enough to measure slopes at the level required to distinguish between cusps and cores of the type predicted in viable SIDM models without stronger priors. However, we find that this provides good enough accuracy to distinguish between the normalization differences predicted at small radii ($r \simeq r_{-2} < r_{\rm core}$) for interesting SIDM models. As the number of galaxies with internal proper motions increases, mass estimators of this kind will enable valuable constraints on SIDM and CDM models.

0 citations0 reviewsNo code linkedOpen access
preprint2019arXiv

Stars made in outflows may populate the stellar halo of the Milky Way

We study stellar-halo formation using six Milky Way-mass galaxies in FIRE-2 cosmological zoom simulations. We find that $5-40\%$ of the outer ($50-300$ kpc) stellar halo in each system consists of $\textit{in-situ}$ stars that were born in outflows from the main galaxy. Outflow stars originate from gas accelerated by super-bubble winds, which can be compressed, cool, and form co-moving stars. The majority of these stars remain bound to the halo and fall back with orbital properties similar to the rest of the stellar halo at $z=0$.In the outer halo, outflow stars are more spatially homogeneous, metal rich, and alpha-element-enhanced than the accreted stellar halo. At the solar location, up to $\sim 10 \%$ of our kinematically-identified halo stars were born in outflows; the fraction rises to as high as $\sim 40\%$ for the most metal-rich local halo stars ([Fe/H] $> -0.5$). We conclude that the Milky Way stellar halo could contain local counterparts to stars that are observed to form in molecular outflows in distant galaxies. Searches for such a population may provide a new, near-field approach to constraining feedback and outflow physics. A stellar halo contribution from outflows is a phase-reversal of the classic halo formation scenario of Eggen, Lynden-Bell $\&$ Sandange, who suggested that halo stars formed in rapidly $\textit{infalling}$ gas clouds. Stellar outflows may be observable in direct imaging of external galaxies and could provide a source for metal-rich, extreme velocity stars in the Milky Way.

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

The Local Group: The Ultimate Deep Field

Near-field cosmology -- using detailed observations of the Local Group and its environs to study wide-ranging questions in galaxy formation and dark matter physics -- has become a mature and rich field over the past decade. There are lingering concerns, however, that the relatively small size of the present-day Local Group ($\sim 2$ Mpc diameter) imposes insurmountable sample-variance uncertainties, limiting its broader utility. We consider the region spanned by the Local Group's progenitors at earlier times and show that it reaches $3' \approx 7$ co-moving Mpc in linear size (a volume of $\approx 350\,{\rm Mpc}^3$) at $z=7$. This size at early cosmic epochs is large enough to be representative in terms of the matter density and counts of dark matter halos with $M_{\rm vir}(z=7) \lesssim 2\times 10^{9}\,M_{\odot}$. The Local Group's stellar fossil record traces the cosmic evolution of galaxies with $10^{3} \lesssim M_{\star}(z=0) / M_{\odot} \lesssim 10^{9}$ (reaching $M_{1500} > -9$ at $z\sim7$) over a region that is comparable to or larger than the Hubble Ultra-Deep Field (HUDF) for the entire history of the Universe. It is highly complementary to the HUDF, as it probes much fainter galaxies but does not contain the intrinsically rarer, brighter sources that are detectable in the HUDF. Archaeological studies in the Local Group also provide the ability to trace the evolution of individual galaxies across time as opposed to evaluating statistical connections between temporally distinct populations. In the JWST era, resolved stellar populations will probe regions larger than the HUDF and any deep JWST fields, further enhancing the value of near-field cosmology.

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
preprint2015arXiv

Are rotating planes of satellite galaxies ubiquitous?

We compare the dynamics of satellite galaxies in the Sloan Digital Sky Survey to simple models in order to test the hypothesis that a large fraction of satellites co-rotate in coherent planes. We confirm the previously-reported excess of co-rotating satellite pairs located near diametric opposition with respect to the host, but show that this signal is unlikely to be due to rotating discs (or planes) of satellites. In particular, no overabundance of co-rotating satellites pairs is observed within $\sim 20^{\circ}-50^{\circ}$ of direct opposition, as would be expected for planar distributions inclined relative to the line-of-sight. Instead, the excess co-rotation for satellite pairs within $\sim 10^{\circ}$ of opposition is consistent with random noise associated with undersampling of an underlying isotropic velocity distribution. We conclude that at most $10\%$ of the hosts in our sample harbor co-rotating satellite planes (as traced by the luminous satellite population).

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Beacons In the Dark: Using Novae and Supernovae to Detect Dwarf Galaxies in the Local Universe

We propose that luminous transients, including novae and supernovae, can be used to detect the faintest galaxies in the universe. Beyond a few Mpc, dwarf galaxies with stellar masses $<10^6 M_{\odot}$ will likely be too faint and/or too low in surface brightness to be directly detected in upcoming large area ground-based photometric surveys. However, single epoch LSST photometry will be able to detect novae to distances of $\sim30$ Mpc and SNe to Gpc-scale distances. Depending on the form of the stellar mass-halo mass relation and the underlying star formation histories of low mass dwarfs, the expected nova rates will be a few to $\sim100$ yr$^{-1}$ and the expected SN rates (including both type Ia and core-collapse) will be $\sim10^2-10^4$ within the observable ($4π$ sr) volume. The transient rate associated with intrahalo stars will be comparably large, but these transients will be located close to bright galaxies, in contrast to the dwarfs, which should trace the underlying large scale structure of the cosmic web. Aggressive follow-up of hostless transients has the potential to uncover the predicted enormous population of low mass field dwarf galaxies.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Forged in FIRE: cusps, cores, and baryons in low-mass dwarf galaxies

We present ultra-high resolution cosmological hydrodynamic simulations of $M_*\simeq10^{4-6.3}M_{\odot}$ dwarf galaxies that form within $M_{v}=10^{9.5-10}M_{\odot}$ dark matter halos. Our simulations rely on the FIRE implementation of star formation feedback and were run with high enough force and mass resolution to directly resolve stellar and dark matter structure on the ~200 pc scales of interest for classical and ultra-faint dwarfs in the Local Group. The resultant galaxies sit on the $M_*$ vs. $M_{v}$ relation required to match the Local Group stellar mass function. They have bursty star formation histories and also form with half-light radii and metallicities that broadly match those observed for local dwarfs at the same stellar mass. We demonstrate that it is possible to create a large (~1 kpc) dark matter core in a cosmological simulation of an $M_*\simeq10^{6.5}M_{\odot}$ dwarf galaxy that resides within an $M_{v}=10^{10}M_{\odot}$ halo -- precisely the scale of interest for resolving the Too Big to Fail problem. However, these large cores are not ubiquitous and appear to correlate closely with the star formation histories of the dwarfs: dark matter cores are largest in systems that form their stars late ($z\lesssim2$), after the early epoch of cusp building mergers has ended. Our $M_*\simeq10^4M_{\odot}$ dwarf retains a cuspy dark matter halo density profile that matches almost identically that of a dark-matter only run of the same system. Despite forming in a field environment, this very low mass dwarf has observable properties that match closely to those of ultra-faint satellite galaxies of the Milky Way, including a uniformly old stellar population (>10 Gyr). Though ancient, most of the stars in our ultra-faint form after reionization; the UV field acts mainly to suppress fresh gas accretion, not to boil away gas that is already present in the proto-dwarf.

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

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
preprint2015arXiv

The Local Group as a time machine: studying the high-redshift Universe with nearby galaxies

We infer the UV luminosities of Local Group galaxies at early cosmic times ($z \sim 2$ and $z \sim 7$) by combining stellar population synthesis modeling with star formation histories derived from deep color-magnitude diagrams constructed from Hubble Space Telescope (HST) observations. Our analysis provides a basis for understanding high-$z$ galaxies - including those that may be unobservable even with the James Webb Space Telescope (JWST) - in the context of familiar, well-studied objects in the very low-$z$ Universe. We find that, at the epoch of reionization, all Local Group dwarfs were less luminous than the faintest galaxies detectable in deep HST observations of blank fields. We predict that JWST will observe $z \sim 7$ progenitors of galaxies similar to the Large Magellanic Cloud today; however, the HST Frontier Fields initiative may already be observing such galaxies, highlighting the power of gravitational lensing. Consensus reionization models require an extrapolation of the observed blank-field luminosity function at $z \approx 7$ by at least two orders of magnitude in order to maintain reionization. This scenario requires the progenitors of the Fornax and Sagittarius dwarf spheroidal galaxies to be contributors to the ionizing background at $z \sim 7$. Combined with numerical simulations, our results argue for a break in the UV luminosity function from a faint-end slope of $α\sim -2$ at $M_{\rm UV} < -13$ to $α\sim -1.2$ at lower luminosities. Applied to photometric samples at lower redshifts, our analysis suggests that HST observations in lensing fields at $z \sim 2$ are capable of probing galaxies with luminosities comparable to the expected progenitor of Fornax.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

The MOSFIRE Deep Evolution Field (MOSDEF) Survey: Rest-Frame Optical Spectroscopy for ~1500 H-Selected Galaxies at 1.37 < z < 3.8

In this paper we present the MOSFIRE Deep Evolution Field (MOSDEF) survey. The MOSDEF survey aims to obtain moderate-resolution (R=3000-3650) rest-frame optical spectra (~3700-7000 Angstrom) for ~1500 galaxies at 1.37<z<3.80 in three well-studied CANDELS fields: AEGIS, COSMOS, and GOODS-N. Targets are selected in three redshift intervals: 1.37<z<1.70, 2.09<z<2.61, and 2.95<z<3.80, down to fixed H_AB (F160W) magnitudes of 24.0, 24.5 and 25.0, respectively, using the photometric and spectroscopic catalogs from the 3D-HST survey. We target both strong nebular emission lines (e.g., [OII], Hbeta, [OIII], 5008, Halpha, [NII], and [SII]) and stellar continuum and absorption features (e.g., Balmer lines, Ca-II H and K, Mgb, 4000 Angstrom break). Here we present an overview of our survey, the observational strategy, the data reduction and analysis, and the sample characteristics based on spectra obtained during the first 24 nights. To date, we have completed 21 masks, obtaining spectra for 591 galaxies. For ~80% of the targets we derive a robust redshift from either emission or absorption lines. In addition, we confirm 55 additional galaxies, which were serendipitously detected. The MOSDEF galaxy sample includes unobscured star-forming, dusty star-forming, and quiescent galaxies and spans a wide range in stellar mass (~10^9-10^11.5 Msol) and star formation rate (~10^0-10^3 Msol/yr). The spectroscopically confirmed sample is roughly representative of an H-band limited galaxy sample at these redshifts. With its large sample size, broad diversity in galaxy properties, and wealth of available ancillary data, MOSDEF will transform our understanding of the stellar, gaseous, metal, dust, and black hole content of galaxies during the time when the universe was most active.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

XMM-Newton Survey of Local O VII Absorption Lines in the Spectra of Active Galactic Nuclei

Highly ionized, z=0 metal absorption lines detected in the X-ray spectra of background active galactic nuclei (AGNs) provide an effective method to probe the hot ($T\sim10^6$ K) gas and its metal content in and around the Milky Way. We present an all-sky survey of the $K_α$ transition of the local O VII absorption lines obtained by Voigt-profile fitting archival XMM-Newton observations. A total of 43 AGNs were selected, among which 12 are BL Lac-type AGNs, and the rest are Seyfert 1 galaxies. At above the $3σ$ level the local O VII absorption lines were detected in 21 AGNs, among which 7 were newly discovered in this work. The sky covering fraction, defined as the ratio between the number of detections and the sample size, increases from at about 40% for all targets to 100% for the brightest targets, suggesting a uniform distribution of the O VII absorbers. We correlate the line equivalent width with the Galactic coordinates and do not find any strong correlations between these quantities. Some AGNs have warm absorbers that may complicate the analysis of the local X-ray absorber since the recession velocity can be compensated by the outflow velocity, especially for the nearby targets. We discuss the potential impact of the warm absorbers on our analysis. A comprehensive theoretical modelling of the X-ray absorbers will be presented in a later paper.

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

Galaxies on FIRE (Feedback In Realistic Environments): Stellar Feedback Explains Cosmologically Inefficient Star Formation

We present a series of high-resolution cosmological simulations of galaxy formation to z=0, spanning halo masses ~10^8-10^13 M_sun, and stellar masses ~10^4-10^11. Our simulations include fully explicit treatment of both the multi-phase ISM (molecular through hot) and stellar feedback. The stellar feedback inputs (energy, momentum, mass, and metal fluxes) are taken directly from stellar population models. These sources of stellar feedback, with zero adjusted parameters, reproduce the observed relation between stellar and halo mass up to M_halo~10^12 M_sun (including dwarfs, satellites, MW-mass disks, and small groups). By extension, this leads to reasonable agreement with the stellar mass function for M_star<10^11 M_sun. We predict weak redshift evolution in the M_star-M_halo relation, consistent with current constraints to z>6. We find that the M_star-M_halo relation is insensitive to numerical details, but is sensitive to the feedback physics. Simulations with only supernova feedback fail to reproduce the observed stellar masses, particularly in dwarf and high-redshift galaxies: radiative feedback (photo-heating and radiation pressure) is necessary to disrupt GMCs and enable efficient coupling of later supernovae to the gas. Star formation rates agree well with the observed Kennicutt relation at all redshifts. The galaxy-averaged Kennicutt relation is very different from the numerically imposed law for converting gas into stars in the simulation, and is instead determined by self-regulation via stellar feedback. Feedback reduces star formation rates considerably and produces a reservoir of gas that leads to rising late-time star formation histories significantly different from the halo accretion history. Feedback also produces large short-timescale variability in galactic SFRs, especially in dwarfs. Many of these properties are not captured by common 'sub-grid' galactic wind models.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Global Properties of M31's Stellar Halo from the SPLASH Survey: II. Metallicity Profile

We present the metallicity distribution of red giant branch (RGB) stars in M31's stellar halo, derived from photometric metallicity estimates for over 1500 spectroscopically confirmed RGB halo stars. The stellar sample comes from 38 halo fields observed with the Keck/DEIMOS spectrograph, ranging from 9 to 175 kpc in projected distance from M31's center, and includes 52 confirmed M31 halo stars beyond 100 kpc. While a wide range of metallicities is seen throughout the halo, the metal-rich peak of the metallicity distribution function becomes significantly less prominent with increasing radius. The metallicity profile of M31's stellar halo shows a continuous gradient from 9 to ~100 kpc, with a magnitude of -0.01 dex/kpc. The stellar velocity distributions in each field are used to identify stars that are likely associated with tidal debris features. The removal of tidal debris features does not significantly alter the metallicity gradient in M31's halo: a gradient is maintained in fields spanning 10 to 90 kpc. We analyze the halo metallicity profile, as well as the relative metallicities of stars associated with tidal debris features and the underlying halo population, in the context of current simulations of stellar halo formation. We argue that the large scale gradient in M31's halo implies M31 accreted at least one relatively massive progenitor in the past, while the field to field variation seen in the metallicity profile indicates that multiple smaller progenitors are likely to have contributed substantially to M31's outer halo.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

M31 Satellite Masses Compared to LCDM Subhaloes

We compare the kinematics of M31's satellite galaxies to the mass profiles of the subhaloes they are expected to inhabit in LCDM. We consider the most massive subhaloes of an approximately M31-sized halo, following the assumption of a monotonic galaxy luminosity-to-subhalo mass mapping. While this abundance matching relation is consistent with the kinematic data for galaxies down to the luminosity of the bright satellites of the Milky Way and M31, it is not consistent with kinematic data for fainter dwarf galaxies (those with L <~ 10^8 Lsun). Comparing the kinematics of M31's dwarf Spheroidal (dSph) satellites to the subhaloes reveals that M31's dSph satellites are too low density to be consistent with the subhaloes' mass profiles. A similar discrepancy has been reported between Milky Way dSphs and their predicted subhaloes, the "too big to fail" problem. By contrast, total mass profiles of the dwarf Elliptical (and similarly bright) satellites are consistent with the subhaloes. However, they suffer from large systematic uncertainties in their dark matter content because of substantial (and potentially dominant) contributions from baryons within their half-light radii.

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

The dynamics of isolated Local Group galaxies

We measured velocities of 862 individual red giant stars in seven isolated dwarf galaxies in the Local Group: NGC 6822, IC 1613, VV 124 (UGC 4879), the Pegasus dwarf irregular galaxy (DDO 216), Leo A, Cetus, and Aquarius (DDO 210). We also computed velocity dispersions, taking into account the measurement uncertainties on individual stars. None of the isolated galaxies is denser than the densest Local Group satellite galaxy. Furthermore, the isolated dwarf galaxies have no obvious distinction in the velocity dispersion--half-light radius plane from the satellite galaxies of the Milky Way and M31. The similarity of the isolated and satellite galaxies' dynamics and structural parameters imposes limitations on environmental solutions to the too-big-to-fail problem, wherein there are fewer dense dwarf satellite galaxies than would be expected from cold dark matter simulations. This data set also has many other applications for dwarf galaxy evolution, including the transformation of dwarf irregular into dwarf spheroidal galaxies. We intend to explore these issues in future work.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

The High-$z$ Universe Confronts Warm Dark Matter: Galaxy Counts, Reionization and the Nature of Dark Matter

We use $N$-body simulations to show that high-redshift galaxy counts provide an interesting constraint on the nature of dark matter, specifically Warm Dark Matter (WDM), owing to the lack of early structure formation these models. Our simulations include three WDM models with thermal-production masses of 0.8 keV, 1.3 keV, and 2.6 keV, as well as CDM. Assuming a relationship between dark halo mass and galaxy luminosity that is set by the observed luminosity function at bright magnitudes, we find that 0.8 keV WDM is disfavored by direct galaxy counts in the Hubble Ultra Deep Field at $>\!\!10σ$. Similarly, 1.3 keV WDM is statistically inconsistent at $2.2σ$. Future observations with JWST (and possibly HST via the Frontier Fields) could rule out $1.3$ keV WDM at high significance, and may be sensitive to WDM masses greater than 2.6 keV. We also examine the ability of galaxies in these WDM models to reionize the universe, and find that 0.8 keV and 1.3 keV WDM produce optical depths to the Cosmic Microwave Background (CMB) that are inconsistent at 68% C.L. with current Planck results, even with extremely high ionizing radiation escape fractions, and 2.6 keV WDM requires an optimistic escape fraction to yield an optical depth consistent with Planck data. Although CMB optical depth calculations are model dependent, we find a strong challenge for stellar processes alone to reionize the universe in a 0.8 keV and 1.3 keV WDM cosmology.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

The surprising inefficiency of dwarf satellite quenching

We study dwarf satellite galaxy quenching using observations from the Geha et al. (2012) NSA/SDSS catalog together with LCDM cosmological simulations to facilitate selection and interpretation. We show that fewer than 30% of dwarfs (M* ~ 10^8.5-10^9.5 Msun) identified as satellites within massive host halos (Mhost ~ 10^12.5-10^14 Msun) are quenched, in spite of the expectation from simulations that half of them should have been accreted more than 6 Gyr ago. We conclude that whatever the action triggering environmental quenching of dwarf satellites, the process must be highly inefficient. We investigate a series of simple, one-parameter quenching models in order to understand what is required to explain the low quenched fraction and conclude that either the quenching timescale is very long (> 9.5 Gyr, a "slow starvation" scenario) or that the environmental trigger is not well matched to accretion within the virial volume. We discuss these results in light of the fact that most of the low mass dwarf satellites in the Local Group are quenched, a seeming contradiction that could point to a characteristic mass scale for satellite quenching.

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

A Dichotomy in Satellite Quenching Around L* Galaxies

We examine the star formation properties of bright (~0.1 L*) satellites around isolated ~L* hosts in the local Universe using spectroscopically confirmed systems in the Sloan Digital Sky Survey DR7. Our selection method is carefully designed with the aid of N-body simulations to avoid groups and clusters. We find that satellites are significantly more likely to be quenched than a stellar mass-matched sample of isolated galaxies. Remarkably, this quenching occurs only for satellites of hosts that are themselves quenched: while star formation is unaffected in the satellites of star-forming hosts, satellites around quiescent hosts are more than twice as likely to be quenched than stellar-mass matched field samples. One implication of this is that whatever shuts down star formation in isolated, passive L* galaxies also plays at least an indirect role in quenching star formation in their bright satellites. The previously-reported tendency for "galactic conformity" in color/morphology may be a by-product of this host-specific quenching dichotomy. The Sérsic indices of quenched satellites are statistically identical to those of field galaxies with the same specific star formation rates, suggesting that environmental and secular quenching give rise to the same morphological structure. By studying the distribution of pairwise velocities between the hosts and satellites, we find dynamical evidence that passive host galaxies reside in dark matter halos that are ~45% more massive than those of star-forming host galaxies of the same stellar mass. We emphasize that even around passive hosts, the mere fact that galaxies become satellites does not typically result in star formation quenching: we find that only ~30% of ~0.1 L* galaxies that fall in from the field are quenched around passive hosts, compared with ~0% around star forming hosts.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Angular Momentum Acquisition in Galaxy Halos

We use high-resolution cosmological hydrodynamic simulations to study the angular momentum acquisition of gaseous halos around Milky Way sized galaxies. We find that cold mode accreted gas enters a galaxy halo with ~70% more specific angular momentum than dark matter averaged over cosmic time (though with a very large dispersion). In fact, we find that all matter has a higher spin parameter when measured at accretion than when averaged over the entire halo lifetime, and is well characterized by λ~0.1, at accretion. Combined with the fact that cold flow gas spends a relatively short time (1-2 dynamical times) in the halo before sinking to the center, this naturally explains why cold flow halo gas has a specific angular momentum much higher than that of the halo and often forms "cold flow disks". We demonstrate that the higher angular momentum of cold flow gas is related to the fact that it tends to be accreted along filaments.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Cold dark matter: controversies on small scales

The cold dark matter (CDM) cosmological model has been remarkably successful in explaining cosmic structure over an enormous span of redshift, but it has faced persistent challenges from observations that probe the innermost regions of dark matter halos and the properties of the Milky Way's dwarf galaxy satellites. We review the current observational and theoretical status of these "small scale controversies." Cosmological simulations that incorporate only gravity and collisionless CDM predict halos with abundant substructure and central densities that are too high to match constraints from galaxy dynamics. The solution could lie in baryonic physics: recent numerical simulations and analytic models suggest that gravitational potential fluctuations tied to efficient supernova feedback can flatten the central cusps of halos in massive galaxies, and a combination of feedback and low star-formation efficiency could explain why most of the dark matter subhalos orbiting the Milky Way do not host visible galaxies. However, it is not clear that this solution can work in the lowest mass galaxies where discrepancies are observed. Alternatively, the small-scale conflicts could be evidence of more complex physics in the dark sector itself. For example, elastic scattering from strong dark matter self-interactions can alter predicted halo mass profiles, leading to good agreement with observations across a wide range of galaxy mass. Gravitational lensing and dynamical perturbations of tidal streams in the stellar halo provide evidence for an abundant population of low mass subhalos in accord with CDM predictions. These observational approaches will get more powerful over the next few years.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Groups of two galaxies in SDSS: implications of colours on star formation quenching time-scales

We have devised a method to select galaxies that are isolated in their dark matter halo (N=1 systems) and galaxies that reside in a group of exactly two (N=2 systems). Our N=2 systems are widely-separated (up to $\sim$\,200\,$h^{-1}$\,kpc), where close galaxy-galaxy interactions are not dominant. We apply our selection criteria to two volume-limited samples of galaxies from SDSS DR6 with $M_{r}-5 \log_{10} h \leq$ -19 and -20 to study the effects of the environment of very sparse groups on galaxy colour. For satellite galaxies in a group of two, we find a red excess attributed to star formation quenching of 0.15\,$\pm$\,0.01 and 0.14\,$\pm$\,0.01 for the -19 and -20 samples, respectively, relative to isolated galaxies of the same stellar mass. Assuming N=1 systems are the progenitors of N=2 systems, an immediate-rapid star formation quenching scenario is inconsistent with these observations. A delayed-then-rapid star formation quenching scenario with a delay time of 3.3 and 3.7\,Gyr for the -19 and -20 samples, respectively, yields a red excess prediction in agreement with the observations. The observations also reveal that central galaxies in a group of two have a slight blue excess of 0.06\,$\pm$\,0.02 and 0.02\,$\pm$\,0.01 for the -19 and -20 samples, respectively, relative to N=1 populations of the same stellar mass. Our results demonstrate that even the environment of very sparse groups of luminous galaxies influence galaxy evolution and in-depth studies of these simple systems are an essential step towards understanding galaxy evolution in general.

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

Segue 2: The Least Massive Galaxy

Segue 2, discovered by Belokurov et al. (2009), is a galaxy with a luminosity of only 900 L_sun. We present Keck/DEIMOS spectroscopy of 25 members of Segue 2--a threefold increase in spectroscopic sample size. The velocity dispersion is too small to be measured with our data. The upper limit with 90% (95%) confidence is sigma_v < 2.2 (2.6) km/s, the most stringent limit for any galaxy. The corresponding limit on the mass within the 3-D half-light radius (46 pc) is M_1/2 < 1.5 (2.1) x 10^5 M_sun. Segue 2 is the least massive galaxy known. We identify Segue 2 as a galaxy rather than a star cluster based the wide dispersion in [Fe/H] (from -2.85 to -1.33) among the member stars. The stars' [alpha/Fe] ratios decline with increasing [Fe/H], indicating that Segue 2 retained Type Ia supernova ejecta despite its presently small mass and that star formation lasted for at least 100 Myr. The mean metallicity, <[Fe/H]> = -2.22 +/- 0.13 (about the same as the Ursa Minor galaxy, 330 times more luminous than Segue 2), is higher than expected from the luminosity-metallicity relation defined by more luminous dwarf galaxy satellites of the Milky Way. Segue 2 may be the barest remnant of a tidally stripped, Ursa Minor-sized galaxy. If so, it is the best example of an ultra-faint dwarf galaxy that came to be ultra-faint through tidal stripping. Alternatively, Segue 2 could have been born in a very low-mass dark matter subhalo (v_max < 10 km/s), below the atomic hydrogen cooling limit.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

The Outer Limits of the M31 System: Kinematics of the Dwarf Galaxy Satellites And XXVIII and And XXIX

We present Keck/DEIMOS spectroscopy of resolved stars in the M31 satellites And XXVIII & And XXIX. We show that these are likely self-bound galaxies based on 18 and 24 members in And XXVIII & And XXIX, respectively. And XXVIII has a systemic velocity of -331.1 +/- 1.8 km/s and velocity dispersion of 4.9+/-1.6 km/s, implying a mass-to-light ratio (within r_1/2) of ~ 44 +/- 41. And XXIX has a systemic velocity of -194.4 +/- 1.5 km/s and velocity dispersion of 5.7+/-1.2 km/s, implying a mass-to-light ratio (within r_1/2) of ~ 124 +/- 72. The internal kinematics and implied masses of And XXVIII & And XXIX are similar to dwarf spheroidals (dSphs) of comparable luminosities, implying that these objects are dark matter-dominated dwarf galaxies. Despite the large projected distances from their host (380 and 188 kpc), the kinematics of these dSph suggest that they are bound M31 satellites.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

The Space Motion of Leo I: Hubble Space Telescope Proper Motion and Implied Orbit

We present the first absolute proper motion measurement of Leo I, based on two epochs of HST ACS/WFC images separated by ~5 years. The average shift of Leo I stars with respect to ~100 background galaxies implies a proper motion of (mu_W, mu_N) = (0.1140 +/- 0.0295, -0.1256 +/- 0.0293) mas/yr. The implied Galactocentric velocity vector, corrected for the reflex motion of the Sun, has radial and tangential components V_rad = 167.9 +/- 2.8 km/s and V_tan = 101.0 +/- 34.4 km/s, respectively. We study the detailed orbital history of Leo I by solving its equations of motion backward in time for a range of plausible mass models for the Milky Way and its surrounding galaxies. Leo I entered the Milky Way virial radius 2.33 +/- 0.21 Gyr ago, most likely on its first infall. It had a pericentric approach 1.05 +/- 0.09 Gyr ago at a Galactocentric distance of 91 +/- 36 kpc. We associate these time scales with characteristic time scales in Leo I's star formation history, which shows an enhanced star formation activity ~2 Gyr ago and quenching ~1 Gyr ago. There is no indication from our calculations that other galaxies have significantly influenced Leo I's orbit, although there is a small probability that it may have interacted with either Ursa Minor or Leo II within the last ~1 Gyr. For most plausible Milky Way masses, the observed velocity implies that Leo I is bound to the Milky Way. However, it may not be appropriate to include it in models of the Milky Way satellite population that assume dynamical equilibrium, given its recent infall. Solution of the complete (non-radial) timing equations for the Leo I orbit implies a Milky Way mass M_MW,vir = 3.15 (-1.36, +1.58) x 10^12 Msun, with the large uncertainty dominated by cosmic scatter. In a companion paper, we compare the new observations to the properties of Leo I subhalo analogs extracted from cosmological simulations.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

The Space Motion of Leo I: The Mass of the Milky Way's Dark Matter Halo

We combine our Hubble Space Telescope measurement of the proper motion of the Leo I dwarf spheroidal galaxy (presented in a companion paper) with the highest resolution numerical simulations of Galaxy-size dark matter halos in existence to constrain the mass of the Milky Way's dark matter halo (M_MW). Despite Leo I's large Galacto-centric space velocity (200 km/s) and distance (261 kpc), we show that it is extremely unlikely to be unbound if Galactic satellites are associated with dark matter substructure, as 99.9% of subhalos in the simulations are bound to their host. The observed position and velocity of Leo I strongly disfavor a low mass Milky Way: if we assume that Leo I is the least bound of the Milky Way's classical satellites, then we find that M_MW > 10^{12} M_sun at 95% confidence for a variety of Bayesian priors on M_MW. In lower mass halos, it is vanishingly rare to find subhalos at 261 kpc moving as fast as Leo I. Should an additional classical satellite be found to be less bound than Leo I, this lower limit on M_MW would increase by 30%. Imposing a mass weighted LCDM prior, we find a median Milky Way virial mass of M_MW=1.6 x 10^{12} M_sun, with a 90% confidence interval of [1.0-2.4] x 10^{12} M_sun. We also confirm a strong correlation between subhalo infall time and orbital energy in the simulations and show that proper motions can aid significantly in interpreting the infall times and orbital histories of satellites.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

The Universal Stellar Mass-Stellar Metallicity Relation for Dwarf Galaxies

We present spectroscopic metallicities of individual stars in seven gas-rich dwarf irregular galaxies (dIrrs), and we show that dIrrs obey the same mass-metallicity relation as the dwarf spheroidal (dSph) satellites of both the Milky Way and M31: Z_* ~ M_*^(0.30 +/- 0.02). The uniformity of the relation is in contradiction to previous estimates of metallicity based on photometry. This relationship is roughly continuous with the stellar mass-stellar metallicity relation for galaxies as massive as M_* = 10^12 M_sun. Although the average metallicities of dwarf galaxies depend only on stellar mass, the shapes of their metallicity distributions depend on galaxy type. The metallicity distributions of dIrrs resemble simple, leaky box chemical evolution models, whereas dSphs require an additional parameter, such as gas accretion, to explain the shapes of their metallicity distributions. Furthermore, the metallicity distributions of the more luminous dSphs have sharp, metal-rich cut-offs that are consistent with the sudden truncation of star formation due to ram pressure stripping.

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
preprint2012arXiv

Cosmological Simulations with Self-Interacting Dark Matter II: Halo Shapes vs. Observations

If dark matter has a large self-interaction scattering cross section, then interactions among dark-matter particles will drive galaxy and cluster halos to become spherical in their centers. Work in the past has used this effect to rule out velocity-independent, elastic cross sections larger than sigma/m ~ 0.02 cm^2/g based on comparisons to the shapes of galaxy cluster lensing potentials and X-ray isophotes. In this paper, we use cosmological simulations to show that these constraints were off by more than an order of magnitude because (a) they did not properly account for the fact that the observed ellipticity gets contributions from the triaxial mass distribution outside the core set by scatterings, (b) the scatter in axis ratios is large and (c) the core region retains more of its triaxial nature than estimated before. Including these effects properly shows that the same observations now allow dark matter self-interaction cross sections at least as large as sigma/m = 0.1 cm^2/g. We show that constraints on self-interacting dark matter from strong-lensing clusters are likely to improve significantly in the near future, but possibly more via central densities and core sizes than halo shapes.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Dark matter concentrations and a search for cores in Milky Way dwarf satellites

We investigate the mass distributions within eight classical Milky Way dwarf spheroidal galaxies (MW dSphs) using an equilibrium Jeans analysis and we compare our results to the mass distributions predicted for subhalos in dissipationless ΛCDM simulations. In order to match the dark matter density concentrations predicted, the stars in these galaxies must have a fairly significant tangential velocity dispersion anisotropy (β~-1.5). For the limiting case of an isotropic velocity dispersion (β=0), the classical MW dSphs predominantly prefer to live in halos that are less concentrated than ΛCDM predictions. We also investigate whether the dSphs prefer to live in halos with constant density cores in the limit of isotropic velocity dispersion. Interestingly, even in this limit, not all of the dSphs prefer large constant-density cores: the Sculptor dSph prefers a cusp while Carina, Draco and Leo I prefer cores. The other four dSphs do not show a statistically significant preference for either cuspy or cored profiles. Finally, we re-examine the hypothesis that the density profiles of these eight MW dSphs can be quantified by a common dark matter halo.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

On the Hot Gas Content of the Milky Way Halo

The Milky Way appears to be missing baryons, as the observed mass in stars and gas is well below the cosmic mean. One possibility is that a substantial fraction of the Galaxy's baryons are embedded within an extended, million-degree hot halo, an idea supported indirectly by observations of warm gas clouds in the halo and gas-free dwarf spheroidal satellites. X-ray observations have established that hot gas does exist in our Galaxy beyond the local hot bubble; however, it may be distributed in a hot disk configuration. Moreover, recent investigations into the X-ray constraints have suggested that any Galactic corona must be insignificant. Here we re-examine the observational data, particularly in the X-ray and radio bands, in order to determine whether it is possible for a substantial fraction of the Galaxy's baryons to exist in ~ 10^6 K gas. In agreement with past studies, we find that a baryonically closed halo is clearly ruled out if one assumes that the hot corona is distributed with a cuspy NFW profile. However, if the hot corona of the galaxy is in an extended, low-density distribution with a large central core, as expected for an adiabatic gas in hydrostatic equilibrium, then it may contain up to 10^11 M_sun of material, possibly accounting for all of the missing Galactic baryons. We briefly discuss some potential avenues for discriminating between a massive, extended hot halo and a local hot disk.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Signatures of minor mergers in the Milky Way disc I: The SEGUE stellar sample

It is now known that minor mergers are capable of creating structure in the phase-space distribution of their host galaxy's disc. In order to search for such imprints in the Milky Way, we analyse the SEGUE F/G-dwarf and the Schuster et al. (2006) stellar samples. We find similar features in these two completely independent stellar samples, consistent with the predictions of a Milky Way minor-merger event. We next apply the same analyses to high-resolution, idealised N-body simulations of the interaction between the Sagittarius dwarf galaxy and the Milky Way. The energy distributions of stellar particle samples in small spatial regions in the host disc reveal strong variations of structure with position. We find good matches to the observations for models with a mass of Sagittarius' dark matter halo progenitor $\lessapprox 10^{11}$ M$_{\odot}$. Thus, we show that this kind of analysis could be used to provide unprecedentedly tight constraints on Sagittarius' orbital parameters, as well as place a lower limit on its mass.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Stellar Kinematics of the Andromeda II Dwarf Spheroidal Galaxy

We present kinematical profiles and metallicity for the M31 dwarf spheroidal (dSph) satellite galaxy Andromeda II (And II) based on Keck DEIMOS spectroscopy of 531 red giant branch stars. Our kinematical sample is among the largest for any M31 satellite and extends out to two effective radii (r_eff = 5.3' = 1.1 kpc). We find a mean systemic velocity of -192.4+-0.5 km/s and an average velocity dispersion of sigma_v = 7.8+-1.1 km/s. While the rotation velocity along the major axis of And II is nearly zero (<1 km/s), the rotation along the minor axis is significant with a maximum rotational velocity of v_max=8.6+-1.8 km/s. We find a kinematical major axis, with a maximum rotational velocity of v_max=10.9+-2.4 km/s, misaligned by 67 degrees to the isophotal major axis. And II is thus the first dwarf galaxy with evidence for nearly prolate rotation with a v_max/sigma_v = 1.1, although given its ellipticity of epsilon = 0.10, this object may be triaxial. We measured metallicities for a subsample of our data, finding a mean metallicity of [Fe/H] = -1.39+- 0.03 dex and an internal metallicity dispersion of 0.72+-0.03 dex. We find a radial metallicity gradient with metal-rich stars more centrally concentrated, but do not observe a significant difference in the dynamics of two metallicity populations. And II is the only known dwarf galaxy to show minor axis rotation making it a unique system whose existence offers important clues on the processes responsible for the formation of dSphs.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

The Milky Way's bright satellites as an apparent failure of LCDM

We use the Aquarius simulations to show that the most massive subhalos in galaxy-mass dark matter halos in LCDM are grossly inconsistent with the dynamics of the brightest Milky Way dwarf spheroidal galaxies. While the best-fitting hosts of the dwarf spheroidals all have 12 < Vmax < 25 km/s, LCDM simulations predict at least ten subhalos with Vmax > 25 km/s. These subhalos are also among the most massive at earlier times, and significantly exceed the UV suppression mass back to z ~ 10. No LCDM-based model of the satellite population of the Milky Way explains this result. The problem lies in the satellites' densities: it is straightforward to match the observed Milky Way luminosity function, but doing so requires the dwarf spheroidals to have dark matter halos that are a factor of ~5 more massive than is observed. Independent of the difficulty in explaining the absence of these dense, massive subhalos, there is a basic tension between the derived properties of the bright Milky Way dwarf spheroidals and LCDM expectations. The inferred infall masses of these galaxies are all approximately equal and are much lower than standard LCDM predictions for systems with their luminosities. Consequently, their implied star formation efficiencies span over two orders of magnitude, from 0.2% to 20% of baryons converted into stars, in stark contrast with expectations gleaned from more massive galaxies. We explore possible solutions to these problems within the context of LCDM and find them to be unconvincing. In particular, we use controlled simulations to demonstrate that the small stellar masses of the bright dwarf spheroidals make supernova feedback an unlikely explanation for their low inferred densities.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

The SPLASH Survey: Kinematics of Andromeda's Inner Spheroid

The combination of large size, high stellar density, high metallicity, and Sersic surface brightness profile of the spheroidal component of the Andromeda galaxy (M31) within R_proj ~ 20 kpc suggest that it is unlike any subcomponent of the Milky Way. In this work we capitalize on our proximity to and external view of M31 to probe the kinematical properties of this "inner spheroid." We employ a Markov chain Monte Carlo (MCMC) analysis of resolved stellar kinematics from Keck/DEIMOS spectra of 5651 red giant branch stars to disentangle M31's inner spheroid from its stellar disk. We measure the mean velocity and dispersion of the spheroid in each of five spatial bins after accounting for a locally cold stellar disk as well as the Giant Southern Stream and associated tidal debris. For the first time, we detect significant spheroid rotation (v_rot ~ 50 km/s) beyond R_proj ~ 5 kpc. The velocity dispersion decreases from about 140 km/s at R_proj = 7 kpc to 120 km/s at R_proj = 14 kpc, consistent to 2 sigma with existing measurements and models. We calculate the probability that a given star is a member of the spheroid and find that the spheroid has a significant presence throughout the spatial extent of our sample. Lastly, we show that the flattening of the spheroid is due to velocity anisotropy in addition to rotation. Though this suggests that the inner spheroid of M31 more closely resembles an elliptical galaxy than a typical spiral galaxy bulge, it should be cautioned that our measurements are much farther out (2 - 14 r_eff) than for the comparison samples.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

The SPLASH Survey: Spectroscopy of 15 M31 Dwarf Spheroidal Satellite Galaxies

We present a resolved-star spectroscopic survey of 15 dwarf spheroidal (dSph) satellites of the Andromeda Galaxy (M31) as part of the Spectroscopic and Photometric Landscape of Andromeda's Stellar Halo (SPLASH) project. We filter foreground contamination from Milky Way (MW) stars, noting that MW substructure is evident in this contaminant sample. We also filter M31 halo field giant stars, and identify the remainder as probable dSph members. We then use these members to determine the kinematical properties of the dSphs. For the first time, we confirm that And XVIII, XXI, and XXII show kinematics consistent with bound, dark matter-dominated galaxies. From the velocity dispersions for the full sample of dSphs we determine masses, which we combine with the size and luminosity of the galaxies to produce mass-size-luminosity scaling relations. With these scalings we determine that the M31 dSphs are fully consistent with the MW dSphs, suggesting that the well-studied MW satellite population provides a fair sample for broader conclusions. We also estimate dark matter halo masses of the satellites, and find that there is no sign that the luminosity of these galaxies depends on their dark halo mass, a result consistent with what is seen for MW dwarfs. Two of the M31 dSphs (And XV, XVI) have estimated maximum circular velocities smaller than 12 km/s (to 1sigma), which likely places them within the lowest mass dark matter halos known to host stars (along with Bootes I of the MW). Finally, we use the systemic velocities of the M31 satellites to estimate the total mass of the M31 halo: within 139 kpc, M31's mass is 8 +4.1 - 3.7 x 10^11 corresponding to a virial mass for M31's dark matter halo of 1.2 +0.9 -0.7 x 10^12, consistent with previous results.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Vertical density waves in the Milky Way disc induced by the Sagittarius Dwarf Galaxy

Recently, Widrow and collaborators announced the discovery of vertical density waves in the Milky Way disk. Here we investigate a scenario where these waves were induced by the Sagittarius dwarf galaxy as it plunged through the Galaxy. Using numerical simulations, we find that the Sagittarius impact produces North-South asymmetries and vertical wave-like behavior that qualitatively agrees with what is observed. The extent to which vertical modes can radially penetrate into the disc, as well as their amplitudes, depend on the mass of the perturbing satellite. We show that the mean height of the disc is expected to vary more rapidly in the radial than in the azimuthal direction. If the observed vertical density asymmetry is indeed caused by vertical oscillations, we predict radial and azimuthal variations of the mean vertical velocity, correlating with the spatial structure. These variations can have amplitudes as large as 8 km/s.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

A Complete Spectroscopic Survey of the Milky Way Satellite Segue 1: The Darkest Galaxy

We present the results of a comprehensive Keck/DEIMOS spectroscopic survey of the ultra-faint Milky Way satellite galaxy Segue 1. We have obtained velocity measurements for 98.2% of the stars within 67 pc (10 arcmin, or 2.3 half-light radii) of the center of Segue 1 that have colors and magnitudes consistent with membership, down to a magnitude limit of r=21.7. Based on photometric, kinematic, and metallicity information, we identify 71 stars as probable Segue 1 members, including some as far out as 87 pc. After correcting for the influence of binary stars using repeated velocity measurements, we determine a velocity dispersion of 3.7^{+1.4}_{-1.1} km/s, with a corresponding mass within the half-light radius of 5.8^{+8.2}_{-3.1} x 10^5 Msun. The stellar kinematics of Segue 1 require very high mass-to-light ratios unless the system is far from dynamical equilibrium, even if the period distribution of unresolved binary stars is skewed toward implausibly short periods. With a total luminosity less than that of a single bright red giant and a V-band mass-to-light ratio of 3400 Msun/Lsun, Segue 1 is the darkest galaxy currently known. We critically re-examine recent claims that Segue 1 is a tidally disrupting star cluster and that kinematic samples are contaminated by the Sagittarius stream. The extremely low metallicities ([Fe/H] < -3) of two Segue 1 stars and the large metallicity spread among the members demonstrate conclusively that Segue 1 is a dwarf galaxy, and we find no evidence in favor of tidal effects. We also show that contamination by the Sagittarius stream has been overestimated. Segue 1 has the highest measured dark matter density of any known galaxy and will therefore be a prime testing ground for dark matter physics and galaxy formation on small scales.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Exploring the links between star formation and minor companions around isolated galaxies

Previous studies have shown that galaxies with minor companions exhibit an elevated star formation rate. We reverse this inquiry, constructing a volume-limited sample of \simL\star (Mr \leq -19.5 + 5 log h) galaxies from the Sloan Digital Sky Survey that are isolated with respect to other luminous galaxies. Cosmological simulations suggest that 99.8% of these galaxies are alone in their dark matter haloes with respect to other luminous galaxies. We search the area around these galaxies for photometric companions. Matching strongly star forming (EW(Hα)\geq 35 Å) and quiescent (EW(Hα)< 35 Å) samples for stellar mass and redshift using a Monte Carlo resampling technique, we demonstrate that rapidly star-forming galaxies are more likely to have photometric companions than other galaxies. The effect is relatively small; about 11% of quiescent, isolated galaxies have minor photometric companions at radii \leq 60 kpc h$^{-1}$ kpc while about 16% of strongly star-forming ones do. Though small, the cumulative difference in satellite counts between strongly star-forming and quiescent galaxies is highly statistically significant (PKS = 1.350 \times10$^{-3}$) out to to radii of \sim 100 h$^{-1}$ kpc. We discuss explanations for this excess, including the possibility that \sim 5% of strongly star-forming galaxies have star formation that is causally related to the presence of a minor companion.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Infall Times for Milky Way Satellites From Their Present-Day Kinematics

We analyze subhalos in the Via Lactea II (VL2) cosmological simulation to look for correlations among their infall times and z = 0 dynamical properties. We find that the present day orbital energy is tightly correlated with the time at which subhalos last crossed into the virial radius. This energy-infall correlation provides a means to infer infall times for Milky Way satellite galaxies. Assuming that the Milky Way's assembly can be modeled by VL2, we show that the infall times of some satellites are well constrained given only their Galactocentric positions and line-of-sight velocities. The constraints sharpen for satellites with proper motion measurements. We find that Carina, Ursa Minor, and Sculptor were all accreted early, more than 8 Gyr ago. Five other dwarfs, including Sextans and Segue 1, are also probable early accreters, though with larger uncertainties. On the other extreme, Leo T is just falling into the Milky Way for the first time while Leo I fell in \sim 2 Gyr ago and is now climbing out of the Milky Way's potential after its first perigalacticon. The energies of several other dwarfs, including Fornax and Hercules, point to intermediate infall times, 2 - 8 Gyr ago. We compare our infall time estimates to published star formation histories and find hints of a dichotomy between ultrafaint and classical dwarfs. The classical dwarfs appear to have quenched star formation after infall but the ultrafaint dwarfs tend to be quenched long before infall, at least for the cases in which our uncertainties allow us to discern differences. Our analysis suggests that the Large Magellanic Cloud crossed inside the Milky Way virial radius recently, within the last \sim 4 billion years.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Observing the End of Cold Flow Accretion using Halo Absorption Systems

We use cosmological SPH simulations to study the cool, accreted gas in two Milky Way-size galaxies through cosmic time to z=0. We find that gas from mergers and cold flow accretion results in significant amounts of cool gas in galaxy halos. This cool circum-galactic component drops precipitously once the galaxies cross the critical mass to form stable shocks, Mvir = Msh ~ 10^12 Msun. Before reaching Msh, the galaxies experience cold mode accretion (T<10^5 K) and show moderately high covering fractions in accreted gas: f_c ~ 30-50% for R<50 co-moving kpc and N_HI>10^16 cm^-2. These values are considerably lower than observed covering fractions, suggesting that outflowing gas (not included here) is important in simulating galaxies with realistic gaseous halos. Within ~500 Myr of crossing the Msh threshold, each galaxy transitions to hot mode gas accretion, and f_c drops to ~5%. The sharp transition in covering fraction is primarily a function of halo mass, not redshift. This signature should be detectable in absorption system studies that target galaxies of varying host mass, and may provide a direct observational tracer of the transition from cold flow accretion to hot mode accretion in galaxies.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Orbiting Circum-galactic Gas as a Signature of Cosmological Accretion

We use cosmological SPH simulations to study the kinematic signatures of cool gas accretion onto a pair of well-resolved galaxy halos. Cold-flow streams and gas-rich mergers produce a circum-galactic component of cool gas that generally orbits with high angular momentum about the galaxy halo before falling in to build the disk. This signature of cosmological accretion should be observable using background-object absorption line studies as features that are offset from the galaxy's systemic velocity by ~100 km/s. Accreted gas typically co-rotates with the central disk in the form of a warped, extended cold flow disk, such that the observed velocity offset is in the same direction as galaxy rotation, appearing in sight lines that avoid the galactic poles. This prediction provides a means to observationally distinguish accreted gas from outflow gas: the accreted gas will show large one-sided velocity offsets in absorption line studies while radial/bi-conical outflows will not (except possibly in special polar projections). This rotation signature has already been seen in studies of intermediate redshift galaxy-absorber pairs; we suggest that these observations may be among the first to provide indirect observational evidence for cold accretion onto galactic halos. Cold mode halo gas typically has ~3-5 times more specific angular momentum than the dark matter. The associated cold mode disk configurations are likely related to extended HI/XUV disks seen around galaxies in the local universe. The fraction of galaxies with extended cold flow disks and associated offset absorption-line gas should decrease around bright galaxies at low redshift, as cold mode accretion dies out.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Quantifying Kinematic Substructure in the Milky Way's Stellar Halo

We present and analyze the positions, distances, and radial velocities for over 4000 blue horizontal-branch (BHB) stars in the Milky Way's halo, drawn from SDSS DR8. We search for position-velocity substructure in these data, a signature of the hierarchical assembly of the stellar halo. Using a cumulative "close pair distribution" (CPD) as a statistic in the 4-dimensional space of sky position, distance, and velocity, we quantify the presence of position-velocity substructure at high statistical significance among the BHB stars: pairs of BHB stars that are close in position on the sky tend to have more similar distances and radial velocities compared to a random sampling of these overall distributions. We make analogous mock-observations of 11 numerical halo formation simulations, in which the stellar halo is entirely composed of disrupted satellite debris, and find a level of substructure comparable to that seen in the actually observed BHB star sample. This result quantitatively confirms the hierarchical build-up of the stellar halo through a signature in phase (position-velocity) space. In detail, the structure present in the BHB stars is somewhat less prominent than that seen in most simulated halos, quite possibly because BHB stars represent an older sub-population. BHB stars located beyond 20 kpc from the Galactic center exhibit stronger substructure than at $\rm r_{gc} < 20$ kpc.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Small-Scale Structure in the SDSS and LCDM: Isolated L* Galaxies with Bright Satellites

We use a volume-limited spectroscopic sample of isolated galaxies in the Sloan Digital Sky Survey (SDSS) to investigate the frequency and radial distribution of luminous (M_r <~ -18.3) satellites like the Large Magellanic Cloud (LMC) around ~L* Milky Way analogs and compare our results object-by-object to LCDM predictions based on abundance matching in simulations. We show that 12% of Milky Way-like galaxies host an LMC-like satellite within 75 kpc (projected), and 42 % within 250 kpc (projected). This implies ~10% have a satellite within the distance of the LMC, and ~40% of L* galaxies host a bright satellite within the virialized extent of their dark matter halos. Remarkably, the simulation reproduces the observed frequency, radial dependence, velocity distribution, and luminosity function of observed secondaries exceptionally well, suggesting that LCDM provides an accurate reproduction of the observed Universe to galaxies as faint as L~10^9 Lsun on ~50 kpc scales. When stacked, the observed projected pairwise velocity dispersion of these satellites is sigma~160 km/s, in agreement with abundance-matching expectations for their host halo masses. Finally, bright satellites around L* primaries are significantly redder than typical galaxies in their luminosity range, indicating that environmental quenching is operating within galaxy-size dark matter halos that typically contain only a single bright satellite. This redness trend is in stark contrast to the Milky Way's LMC, which is unusually blue even for a field galaxy. We suggest that the LMC's discrepant color might be further evidence that it is undergoing a triggered star-formation event upon first infall.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

The Sagittarius impact as an architect of spirality and outer rings in the Milky Way

Like many galaxies of its size, the Milky Way is a disk with prominent spiral arms rooted in a central bar, although our knowledge of its structure and origin is incomplete. Traditional attempts to understand the Galaxy's morphology assume that it has been unperturbed by major external forces. Here we report simulations of the response of the Milky Way to the infall of the Sagittarius dwarf galaxy (Sgr), which results in the formation of spiral arms, influences the central bar and produces a flared outer disk. Two ring-like wrappings emerge towards the Galactic anti-Center in our model that are reminiscent of the low- latitude arcs observed in the same area of the Milky Way. Previous models have focused on Sgr itself to reproduce the dwarf's orbital history and place associated constraints on the shape of the Milky Way gravitational potential, treating the Sgr impact event as a trivial influence on the Galactic disk. Our results show that the Milky Way's morphology is not purely secular in origin and that low-mass minor mergers predicted to be common throughout the Universe probably have a similarly important role in shaping galactic structure.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Too big to fail? The puzzling darkness of massive Milky Way subhaloes

We show that dissipationless LCDM simulations predict that the majority of the most massive subhaloes of the Milky Way are too dense to host any of its bright satellites (L_V > 10^5 L_sun). These dark subhaloes have circular velocities at infall of 30-70 km/s and infall masses of [0.2-4] x 10^10 M_sun. Unless the Milky Way is a statistical anomaly, this implies that galaxy formation becomes effectively stochastic at these masses. This is in marked contrast to the well-established monotonic relation between galaxy luminosity and halo circular velocity (or halo mass) for more massive haloes. We show that at least two (and typically four) of these massive dark subhaloes are expected to produce a larger dark matter annihilation flux than Draco. It may be possible to circumvent these conclusions if baryonic feedback in dwarf satellites or different dark matter physics can reduce the central densities of massive subhaloes by order unity on a scale of 0.3 - 1 kpc.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Accurate masses for dispersion-supported galaxies

We derive an accurate mass estimator for dispersion-supported stellar systems and demonstrate its validity by analyzing resolved line-of-sight velocity data for globular clusters, dwarf galaxies, and elliptical galaxies. Specifically, by manipulating the spherical Jeans equation we show that the dynamical mass enclosed within the 3D deprojected half-light radius r_1/2 can be determined with only mild assumptions about the spatial variation of the stellar velocity dispersion anisotropy. We find M_1/2 = 3 σ_los^2 r_1/2 / G ~ 4 σ_los^2 R_eff / G, where σ_los^2 is the luminosity-weighted square of the line-of-sight velocity dispersion and R_eff is the 2D projected half-light radius. While deceptively familiar in form, this formula is not the virial theorem, which cannot be used to determine accurate masses unless the radial profile of the total mass is known a priori. We utilize this finding to show that all of the Milky Way dwarf spheroidal galaxies (MW dSphs) are consistent with having formed within a halo of mass approximately 3 x 10^9 M_sun in Lambda CDM cosmology. The faintest MW dSphs seem to have formed in dark matter halos that are at least as massive as those of the brightest MW dSphs, despite the almost five orders of magnitude spread in luminosity. We expand our analysis to the full range of observed dispersion-supported stellar systems and examine their I-band mass-to-light ratios (M/L). The M/L vs. M_1/2 relation for dispersion-supported galaxies follows a U-shape, with a broad minimum near M/L ~ 3 that spans dwarf elliptical galaxies to normal ellipticals, a steep rise to M/L ~ 3,200 for ultra-faint dSphs, and a more shallow rise to M/L ~ 800 for galaxy cluster spheroids.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Counts-in-Cylinders in the Sloan Digital Sky Survey with Comparisons to N-body Simulations

Environmental statistics provide a necessary means of comparing the properties of galaxies in different environments and a vital test of models of galaxy formation within the prevailing, hierarchical cosmological model. We explore counts-in-cylinders, a common statistic defined as the number of companions of a particular galaxy found within a given projected radius and redshift interval. Galaxy distributions with the same two-point correlation functions do not necessarily have the same companion count distributions. We use this statistic to examine the environments of galaxies in the Sloan Digital Sky Survey, Data Release 4. We also make preliminary comparisons to four models for the spatial distributions of galaxies, based on N-body simulations, and data from SDSS DR4 to study the utility of the counts-in-cylinders statistic. There is a very large scatter between the number of companions a galaxy has and the mass of its parent dark matter halo and the halo occupation, limiting the utility of this statistic for certain kinds of environmental studies. We also show that prevalent, empirical models of galaxy clustering that match observed two- and three-point clustering statistics well fail to reproduce some aspects of the observed distribution of counts-in-cylinders on 1, 3 and 6-Mpc/h scales. All models that we explore underpredict the fraction of galaxies with few or no companions in 3 and 6-Mpc/h cylinders. Roughly 7% of galaxies in the real universe are significantly more isolated within a 6 Mpc/h cylinder than the galaxies in any of the models we use. Simple, phenomenological models that map galaxies to dark matter halos fail to reproduce high-order clustering statistics in low-density environments.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

From Galaxy Clusters to Ultra-Faint Dwarf Spheroidals: A Fundamental Curve Connecting Dispersion-supported Galaxies to Their Dark Matter Halos

We examine scaling relations of dispersion-supported galaxies over more than eight orders of magnitude in luminosity by transforming standard fundamental plane parameters into a space of mass (M1/2), radius (r1/2), and luminosity (L1/2). We find that from ultra-faint dwarf spheroidals to giant cluster spheroids, dispersion-supported galaxies scatter about a one-dimensional "fundamental curve" through this MRL space. The weakness of the M1/2-L1/2 slope on the faint end may imply that potential well depth limits galaxy formation in small galaxies, while the stronger dependence on L1/2 on the bright end suggests that baryonic physics limits galaxy formation in massive galaxies. The mass-radius projection of this curve can be compared to median dark matter halo mass profiles of LCDM halos in order to construct a virial mass-luminosity relationship (Mvir-L) for galaxies that spans seven orders of magnitude in Mvir. Independent of any global abundance or clustering information, we find that (spheroidal) galaxy formation needs to be most efficient in halos of Mvir ~ 10^12 Msun and to become inefficient above and below this scale. Moreover, this profile matching technique is most accurate at the high and low luminosity extremes (where dark matter fractions are highest) and is therefore quite complementary to statistical approaches that rely on having a well-sampled luminosity function. We also consider the significance and utility of the scatter about this relation, and find that in the dSph regime observational errors are almost at the point where we can explore the intrinsic scatter in the luminosity-virial mass relation. Finally, we note that purely stellar systems like Globular Clusters and Ultra Compact Dwarfs do not follow the fundamental curve relation. This allows them to be easily distinguished from dark-matter dominated dSph galaxies in MRL space. (abridged)

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Heated Disc Stars in the Stellar Halo

Minor accretion events with mass ratio M_sat : M_host ~ 1:10 are common in the context of LCDM cosmology. We use high-resolution simulations of Galaxy-analogue systems to show that these mergers can dynamically eject disk stars into a diffuse light component that resembles a stellar halo both spatially and kinematically. For a variety of orbital configurations, we find that ~3-5e8 M_sun of primary stellar disk material is ejected to a distance larger than 5 kpc above the galactic plane. This ejected contribution is similar to the mass contributed by the tidal disruption of the satellite galaxy itself, though it is less extended. If we restrict our analysis to the approximate solar neighborhood in the disk plane, we find that ~1% of the initial disk stars in that region would be classified kinematically as halo stars. Our results suggest that the inner parts of galactic stellar halos contain ancient disk stars and that these stars may have been liberated in the very same events that delivered material to the outer stellar halo.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Notes on the Missing Satellites Problem

The Missing Satellites Problem (MSP) broadly refers to the overabundance of predicted Cold Dark Matter (CDM) subhalos compared to satellite galaxies known to exist in the Local Group. The most popular interpretation of the MSP is that the smallest dark matter halos in the universe are extremely inefficient at forming stars. The question from that standpoint is to identify the feedback source that makes small halos dark and to identify any obvious mass scale where the truncation in the efficiency of galaxy formation occurs. Among the most exciting developments in near-field cosmology in recent years is the discovery of a new population satellite galaxies orbiting the Milky Way and M31. Wide field, resolved star surveys have more than doubled the dwarf satellite count in less than a decade, revealing a population of ultrafaint galaxies that are less luminous that some star clusters. For the first time, there are empirical reasons to believe that there really are missing satellite galaxies in the Local Group, lurking just beyond our ability to detect them, or simply inhabiting a region of the sky that has yet to have been surveyed. Both kinematic studies and completeness-correction studies seem to point to a characteristic potential well depth for satellite subhalos that is quite close to the mass scale where photoionization and atomic cooling should limit galaxy formation. Among the more pressing problems associated with this interpretation is to understand the selection biases that limit our ability to detect the lowest mass galaxies. The least massive satellite halos are likely to host stealth galaxies with very-low surface brightness and this may be an important limitation in the hunt for low-mass fossils from the epoch of reionization.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Stealth Galaxies in the Halo of the Milky Way

We predict that there is a population of low-luminosity dwarf galaxies orbiting within the halo of the Milky Way that have surface brightnesses low enough to have escaped detection in star-count surveys. The overall count of stealth galaxies is sensitive to the presence (or lack) of a low-mass threshold in galaxy formation. These systems have luminosities and stellar velocity dispersions that are similar to those of known ultrafaint dwarf galaxies but they have more extended stellar distributions (half light radii greater than about 100 pc) because they inhabit dark subhalos that are slightly less massive than their higher surface brightness counterparts. As a result, the typical peak surface brightness is fainter than 30 mag per square arcsec. One implication is that the inferred common mass scale for Milky Way dwarfs may be an artifact of selection bias. If there is no sharp threshold in galaxy formation at low halo mass, then ultrafaint galaxies like Segue 1 represent the high-mass, early forming tail of a much larger population of objects that could number in the hundreds and have typical peak circular velocities of about 8 km/s and masses within 300 pc of about 5 million solar masses. Alternatively, if we impose a low-mass threshold in galaxy formation in order to explain the unexpectedly high densities of the ultrafaint dwarfs, then we expect only a handful of stealth galaxies in the halo of the Milky Way. A complete census of these objects will require deeper sky surveys, 30m-class follow-up telescopes, and more refined methods to identify extended, self-bound groupings of stars in the halo.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

The Case Against Warm or Self-Interacting Dark Matter as Explanations for Cores in Low Surface Brightness Galaxies

Warm dark matter (WDM) and self-interacting dark matter (SIDM) are often motivated by the inferred cores in the dark matter halos of low surface brightness (LSB) galaxies. We test thermal WDM, non-thermal WDM, and SIDM using high-resolution rotation curves of nine LSB galaxies. We fit these dark matter models to the data and determine the halo core radii and central densities. While the minimum core size in WDM models is predicted to decrease with halo mass, we find that the inferred core radii increase with halo mass and also cannot be explained with a single value of the primordial phase space density. Moreover, if the core size is set by WDM particle properties, then even the smallest cores we infer would require primordial phase space density values that are orders of magnitude smaller than lower limits obtained from the Lyman alpha forest power spectra. We also find that the dark matter halo core densities vary by a factor of about 30 from system to system while showing no systematic trend with the maximum rotation velocity of the galaxy. This strongly argues against the core size being directly set by large self-interactions (scattering or annihilation) of dark matter. We therefore conclude that the inferred cores do not provide motivation to prefer WDM or SIDM over other dark matter models.

0 citations0 reviewsNo code linkedOpen access
preprint2009arXiv

LSST Science Book, Version 2.0

A survey that can cover the sky in optical bands over wide fields to faint magnitudes with a fast cadence will enable many of the exciting science opportunities of the next decade. The Large Synoptic Survey Telescope (LSST) will have an effective aperture of 6.7 meters and an imaging camera with field of view of 9.6 deg^2, and will be devoted to a ten-year imaging survey over 20,000 deg^2 south of +15 deg. Each pointing will be imaged 2000 times with fifteen second exposures in six broad bands from 0.35 to 1.1 microns, to a total point-source depth of r~27.5. The LSST Science Book describes the basic parameters of the LSST hardware, software, and observing plans. The book discusses educational and outreach opportunities, then goes on to describe a broad range of science that LSST will revolutionize: mapping the inner and outer Solar System, stellar populations in the Milky Way and nearby galaxies, the structure of the Milky Way disk and halo and other objects in the Local Volume, transient and variable objects both at low and high redshift, and the properties of normal and active galaxies at low and high redshift. It then turns to far-field cosmological topics, exploring properties of supernovae to z~1, strong and weak lensing, the large-scale distribution of galaxies and baryon oscillations, and how these different probes may be combined to constrain cosmological models and the physics of dark energy.

0 citations0 reviewsNo code linkedOpen access
preprint2009arXiv

Lyman break galaxy close and interacting pairs at z ~ 3

To date, the identification of interactions at z ~ 3 and above has relied on morphological analysis. Here, we present five serendipitous spectroscopic z ~ 3 Lyman break galaxy (LBG) pairs with projected proper separations < 15 h^-1 kpc in our survey of nine separate Keck fields. The data consist of 140 of our highest signal-to-noise ratio LBG spectra and ~500 of our most confident colour-selected LBGs. We show that the pairs are composed of two distinct close and/or interacting LBGs from a detailed analysis of the rest-frame ultraviolet spectra and images. In addition, we show that the pair number and separation distribution is expected from (1) the angular correlation function when applied to our survey and ~2500 colour-selected LBGs from the literature and (2) an analysis of a carefully matched high-resolution hybrid numerical and analytical cosmological simulation. Because the spectroscopic slitlets have random orientations with respect to the close pairs on the sky, the serendipitous pairs provide an unbiased sampling of the underlying close pair fraction. Finally, we discover two Ly-a emitters (LAEs) in our slitlets and find that they reside within 50 projected h^-1 kpc of LBGs. In this work, we uncover a strong relationship between Ly-a emission and pair separation. All confirmed and all candidate LBG pairs with separations of < 15 projected h^-1 kpc exhibit Ly-a in emission and we find an overabundance of Ly-a emission in pairs with < 50 projected h^-1 kpc separations. This relationship suggests a picture in which a measurable fraction of the Ly-a emission of LBGs, and potentially LAEs, is generated via interaction mechanisms. As a result, serendipitous spectroscopic close pairs provide a unique means to help identify and study high-redshift galaxy interactions using conventional ground-based optical data.

0 citations0 reviewsNo code linkedOpen access
preprint2009arXiv

Redistributing hot gas around galaxies: do cool clouds signal a solution to the overcooling problem?

We present a pair of high-resolution smoothed particle hydrodynamics (SPH) simulations that explore the evolution and cooling behavior of hot gas around Milky-Way size galaxies. The simulations contain the same total baryonic mass and are identical other than their initial gas density distributions. The first is initialised with a low entropy hot gas halo that traces the cuspy profile of the dark matter, and the second is initialised with a high-entropy hot halo with a cored density profile as might be expected in models with pre-heating feedback. Galaxy formation proceeds in dramatically different fashion depending on the initial setup. While the low-entropy halo cools rapidly, primarily from the central region, the high-entropy halo is quasi-stable for ~4 Gyr and eventually cools via the fragmentation and infall of clouds from ~100 kpc distances. The low-entropy halo's X-ray surface brightness is ~100 times brighter than current limits and the resultant disc galaxy contains more than half of the system's baryons. The high-entropy halo has an X-ray brightness that is in line with observations, an extended distribution of pressure-confined clouds reminiscent of observed populations, and a final disc galaxy that has half the mass and ~50% more specific angular momentum than the disc formed in the low-entropy simulation. The final high-entropy system retains the majority of its baryons in a low-density hot halo. The hot halo harbours a trace population of cool, mostly ionised, pressure-confined clouds that contain ~10% of the halo's baryons after 10 Gyr of cooling. The covering fraction for HI and MgII absorption clouds in the high-entropy halo is ~0.4 and ~0.6, respectively, although most of the mass that fuels disc growth is ionised, and hence would be under counted in HI surveys.

0 citations0 reviewsNo code linkedOpen access
preprint2009arXiv

The Dark Disk of the Milky Way

Massive satellite accretions onto early galactic disks can lead to the deposition of dark matter in disk-like configurations that co-rotate with the galaxy. This phenomenon has potentially dramatic consequences for dark matter detection experiments. We utilize focused, high-resolution simulations of accretion events onto disks designed to be Galaxy analogues, and compare the resultant disks to the morphological and kinematic properties of the Milky Way's thick disk in order to bracket the range of co-rotating accreted dark matter. We find that the Milky Way's merger history must have been unusually quiescent compared to median LCDM expectations and therefore its dark disk must be relatively small: the fraction of accreted dark disk material near the Sun is about 20% of the host halo density or smaller and the co-rotating dark matter fraction near the Sun, defined as particles moving with a rotational velocity lag less than 50 km/s, is enhanced by about 30% or less compared to a standard halo model. Such a dark disk could contribute dominantly to the low energy (of order keV for a dark matter particle with mass 100 GeV) nuclear recoil event rate of direct dectection experiments, but it will not change the likelihood of detection significantly. These dark disks provide testable predictions of weakly-interacting massive particle dark matter models and should be considered in detailed comparisons to experimental data. Our findings suggest that the dark disk of the Milky Way may provide a detectable signal for indirect detection experiments, contributing up to about 25% of the dark matter self-annihilation signal in the direction of the center of the Galaxy, lending the signal a noticeably oblate morphology.

0 citations0 reviewsNo code linkedOpen access
preprint2008arXiv

Cold Dark Matter Substructure and Galactic Disks

We perform a set of high-resolution, dissipationless N-body simulations to investigate the influence of cold dark matter (CDM) substructure on the dynamical evolution of thin galactic disks. Our method combines cosmological simulations of galaxy-sized CDM halos to derive the properties of substructure populations and controlled numerical experiments of consecutive subhalo impacts onto initially-thin, fully-formed disk galaxies. We demonstrate that close encounters between massive subhalos and galactic disks since z~1 should be common occurrences in LCDM models. In contrast, extremely few satellites in present-day CDM halos are likely to have a significant impact on the disk structure. One typical host halo merger history is used to seed controlled N-body experiments of subhalo-disk encounters. As a result of these accretion events, the disk thickens considerably at all radii with the disk scale height increasing in excess of a factor of 2 in the solar neighborhood. We show that interactions with the subhalo population produce a wealth of distinctive morphological signatures in the disk stars including: conspicuous flares; bars; low-lived, ring-like features in the outskirts; and low-density, filamentary structures above the disk plane. We compare a resulting dynamically-cold, ring-like feature in our simulations to the Monoceros ring stellar structure in the MW. The comparison shows quantitative agreement in both spatial distribution and kinematics, suggesting that such observed complex stellar components may arise naturally as disk stars are excited by encounters with subhalos. These findings highlight the significant role of CDM substructure in setting the structure of disk galaxies and driving galaxy evolution.

0 citations0 reviewsNo code linkedOpen access
preprint2008arXiv

Cold Dark Matter Substructure and Galactic Disks I: Morphological Signatures of Hierarchical Satellite Accretion

(Abridged) We conduct a series of high-resolution, dissipationless N-body simulations to investigate the cumulative effect of substructure mergers onto thin disk galaxies in the context of the LCDM paradigm of structure formation. Our simulation campaign is based on a hybrid approach. Substructure properties are culled directly from cosmological simulations of galaxy-sized cold dark matter (CDM) halos. In contrast to what can be inferred from statistics of the present-day substructure populations, accretions of massive subhalos onto the central regions of host halos, where the galactic disk resides, since z~1 should be common occurrences. One host halo merger history is subsequently used to seed controlled numerical experiments of repeated satellite impacts on an initially-thin Milky Way-type disk galaxy. We show that these accretion events produce several distinctive observational signatures in the stellar disk including: a ring-like feature in the outskirts; a significant flare; a central bar; and faint filamentary structures that (spuriously) resemble tidal streams. The final distribution of disk stars exhibits a complex vertical structure that is well-described by a standard ``thin-thick'' disk decomposition. We conclude that satellite-disk encounters of the kind expected in LCDM models can induce morphological features in galactic disks that are similar to those being discovered in the Milky Way, M31, and in other disk galaxies. These results highlight the significant role of CDM substructure in setting the structure of disk galaxies and driving galaxy evolution. Upcoming galactic structure surveys and astrometric satellites may be able to distinguish between competing cosmological models by testing whether the detailed structure of galactic disks is as excited as predicted by the CDM paradigm.

0 citations0 reviewsNo code linkedOpen access
preprint2008arXiv

Merger Histories of Galaxy Halos and Implications for Disk Survival

We study the merger histories of galaxy dark matter halos using a high resolution LCDM N-body simulation. Our merger trees follow ~17,000 halos with masses M_0 = (10^11--10^13) Msun at z=0 and track accretion events involving objects as small as m = 10^10 Msun. We find that mass assembly is remarkably self-similar in m/M_0, and dominated by mergers that are ~10% of the final halo mass. While very large mergers, m > 0.4 M_0, are quite rare, sizeable accretion events, m ~ 0.1 M_0, are common. Over the last 10 Gyr, an overwhelming majority (~95%) of Milky Way-sized halos with M_0 = 10^12 Msun have accreted at least one object with greater total mass than the Milky Way disk (m > 5x10^10 Msun), and approximately 70% have accreted an object with more than twice that mass (m > 10^11 Msun). Our results raise serious concerns about the survival of thin-disk dominated galaxies within the current paradigm for galaxy formation in a CDM universe. In order to achieve a ~70% disk-dominated fraction in Milky Way-sized CDM halos, mergers involving m ~ 2x10^11 Msun objects must not destroy disks. Considering that most thick disks and bulges contain old stellar populations, the situation is even more restrictive: these mergers must not heat disks or drive gas into their centers to create young bulges.

0 citations0 reviewsNo code linkedOpen access
preprint2008arXiv

Mergers and Disk Survival in LCDM

Disk galaxies are common in our universe and this is a source of concern for hierarchical formation models like LCDM. Here we investigate this issue as motivated by raw merger statistics derived for galaxy-size dark matter halos from LCDM simulations. Our analysis shows that a majority (~70%) of galaxy halos with M = 10^12 M_sun at z=0 should have accreted at least one object with mass m > 10^11 M_sun ~ 3 M_disk over the last 10 Gyr. Mergers involving larger objects should have been rare, and this pinpoints m/M = 0.1 mass-ratio mergers as the most worrying ones for the survival of thin galactic disks. Motivated by these results, we use use high-resolution, dissipationless N-body simulations to study the response of stellar Milky-Way type disks to these common mergers and show that thin disks do not survive the bombardment. The remnant galaxies are roughly three times as thick and twice as kinematically hot as the observed thin disk of the Milky Way. Finally, we evaluate the suggestion that disks may be preserved if the mergers involve gas-rich progenitors. Using empirical measures to assign stellar masses and gas masses to dark matter halos as a function of redshift, we show that the vast majority of large mergers experienced by 10^12 M_sun halos should be gas-rich (f_gas > 0.5), suggesting that this is a potentially viable solution to the disk formation conundrum. Moreover, gas-rich mergers should become increasingly rare in more massive halos > 10^12.5 M_sun, and this suggest that merger gas fractions may play an important role in establishing morphological trends with galaxy luminosity.

0 citations0 reviewsNo code linkedOpen access
preprint2008arXiv

The Metallicity of Diffuse Intrahalo Light

We make predictions for the metallicity of diffuse stellar components in systems ranging from small spiral galaxies to rich galaxy clusters. We extend the formalism of Purcell et al. (2007), in which diffuse stellar mass is produced via galaxy disruption, and we convolve this result with the observed mass-metallicity relation for galaxies in order to analyze the chemical abundance of intrahalo light (IHL) in host halos with virial mass 10^10.5 M_sun < M_host < 10^15 M_sun. We predict a steep rise of roughly two dex in IHL metallicity from the scales of small to large spiral galaxies. In terms of the total dynamical mass M_host of the host systems under consideration, we predict diffuse light metallicities ranging from Z_IHL < -2.5 for M_host ~ 10^11 M_sun, to Z_IHL ~ -1.0 for M_host ~ 10^12 M_sun. In larger systems, we predict a gradual flattening of this trend with Z_IHL ~ -0.4 for M_host ~ 10^13 M_sun, increasing to Z_IHL ~ 0.1 for M_host ~ 10^15 M_sun. This behavior is coincident with a narrowing of the intrahalo metallicity distribution as host mass increases. The observable distinction in surface brightness between old, metal-poor IHL stars and more metal-rich, dynamically-younger tidal streams is of crucial importance when estimating the chemical abundance of an intrahalo population with multiple origins.

0 citations0 reviewsNo code linkedOpen access
preprint2007arXiv

Isolating Triggered Star Formation

Galaxy pairs provide a potentially powerful means of studying triggered star formation from galaxy interactions. We use a large cosmological N-body simulation coupled with a well-tested semi-analytic substructure model to demonstrate that the majority of galaxies in close pairs reside within cluster or group-size halos and therefore represent a biased population, poorly suited for direct comparison to ``field'' galaxies. Thus, the frequent observation that some types of galaxies in pairs have redder colors than ``field'' galaxies is primarily a selection effect. We select galaxy pairs that are isolated in their dark matter halos with respect to other massive subhalos (N=2 halos) and a control sample of isolated galaxies (N=1 halos) for comparison. We then apply these selection criteria to a volume-limited subset of the 2dF Galaxy Redshift Survey with M_Bj <= -19 and obtain the first clean measure of the typical fraction of galaxies affected by triggered star formation and the average elevation in the star formation rate. We find that 24% (30.5%) of these L^\star and sub-L^{\star} galaxies in isolated 50 (30) kpc/h pairs exhibit star formation that is boosted by a factor of >~ 5 above their average past value, while only 10% of isolated galaxies in the control sample show this level of enhancement. Thus, 14% (20 %) of the galaxies in these close pairs show clear triggered star formation. The isolation criteria we develop provide a means to constrain star formation and feedback prescriptions in hydrodynamic simulations and a very general method of understanding the importance of triggered star formation in a cosmological context. (Abridged.)

0 citations0 reviewsNo code linkedOpen access
preprint2007arXiv

Redefining the Missing Satellites Problem

Numerical simulations of Milky-Way size Cold Dark Matter (CDM) halos predict a steeply rising mass function of small dark matter subhalos and a substructure count that greatly outnumbers the observed satellites of the Milky Way. Several proposed explanations exist, but detailed comparison between theory and observation in terms of the maximum circular velocity (Vmax) of the subhalos is hampered by the fact that Vmax for satellite halos is poorly constrained. We present comprehensive mass models for the well-known Milky Way dwarf satellites, and derive likelihood functions to show that their masses within 0.6 kpc (M_0.6) are strongly constrained by the present data. We show that the M_0.6 mass function of luminous satellite halos is flat between ~ 10^7 and 10^8 M_\odot. We use the ``Via Lactea'' N-body simulation to show that the M_0.6 mass function of CDM subhalos is steeply rising over this range. We rule out the hypothesis that the 11 well-known satellites of the Milky Way are hosted by the 11 most massive subhalos. We show that models where the brightest satellites correspond to the earliest forming subhalos or the most massive accreted objects both reproduce the observed mass function. A similar analysis with the newly-discovered dwarf satellites will further test these scenarios and provide powerful constraints on the CDM small-scale power spectrum and warm dark matter models.

0 citations0 reviewsNo code linkedOpen access
preprint2007arXiv

Shredded Galaxies as the Source of Diffuse Intrahalo Light On Varying Scales

We make predictions for diffuse stellar mass fractions in dark matter halos from the scales of small spiral galaxies to those of large galaxy clusters. We use an extensively-tested analytic model for subhalo infall and evolution and empirical constraints from galaxy survey data to set the stellar mass in each accreted subhalo to model diffuse light. We add stellar mass to the diffuse light as subhalos become disrupted due to interactions within their host halos. We predict that the stellar mass fraction in diffuse, intrahalo light should rise on average from ~0.5% to approximately 20% from small galaxy halos to poor groups. The trend with mass flattens considerably beyond the group scale, increasing weakly from a fraction of ~20% in poor galaxy clusters (~10^14 M_sun) to roughly ~30% in massive clusters (~10^15 M_sun). The mass-dependent diffuse light fraction is governed primarily by the empirical fact that the mass-to-light ratio in galaxy halos must vary as a function of halo mass. Galaxy halos have little diffuse light because they accrete most of their mass in small subhalos that themselves have high mass-to-light ratios; stellar halos around galaxies are built primarily from disrupted dwarf-irregular-type galaxies with M*~10^8.5 M_sun. The diffuse light in group and cluster halos is built from satellite galaxies that form stars efficiently and have correspondingly low mass-to-light ratios; intracluster light is dominated by material liberated from massive galaxies with M*~10^11 M_sun. Our results are consistent with existing observations spanning the galaxy, group, and cluster scale; however, they can be tested more rigorously in future deep surveys for faint diffuse light.

0 citations0 reviewsNo code linkedOpen access
preprint2006arXiv

A Large Dark Matter Core in the Fornax Dwarf Spheroidal Galaxy?

We use measurements of the stellar velocity dispersion profile of the Fornax dwarf spheroidal galaxy to derive constraints on its dark matter distribution. Though the data are unable to distinguish between models with small cores and those with cusps, we show that a large > 1 kpc dark matter core in Fornax is highly implausible. Irrespective of the origin of the core, reasonable dynamical limits on the mass of the Fornax halo constrain its core radius to be no larger than 700 pc. We derive an upper limit core radius of 300 pc by demanding that the central phase space density of Fornax not exceed that directly inferred from the rotation curves of low-mass spiral galaxies. Further, if the halo is composed of warm dark matter then phase-space constraints force the core to be quite small in order to avoid conservative limits from the Ly alpha forest power spectrum, implying a core radius < 85 pc. We discuss our results in the context of the idea that the extended globular cluster distribution in Fornax can be explained by the presence of a large 1.5 kpc core. A self-consistent core of this size would be drastically inconsistent with the expectations of standard warm or cold dark matter models, and would also require an unreasonably massive dark matter halo, with a maximum circular velocity of 200 km/s.

0 citations0 reviewsNo code linkedOpen access
preprint2006arXiv

Close Galaxy Counts as a Probe of Hierarchical Structure Formation

Standard LCDM predicts that the major merger rate of galaxy-size dark matter halos rises rapidly with redshift. The average number of close companions per galaxy, Nc, is often used to infer the galaxy merger rate, however, recent observational studies suggest that Nc evolves very little with redshift. Here we use a "hybrid" N- body simulation plus analytic substructure model to predict Nc directly. We identify dark matter subhalos with galaxies and show that the observed lack of close pair count evolution arises because the high merger rate per halo at early times is counteracted by a decrease in the number of halos massive enough to host a galaxy pair. We compare our results to data compiled from the DEEP2, SSRS2, and the UZC redshift surveys. Observed pair counts match our predictions if we assume a monotonic mapping between galaxy luminosity and the maximum circular velocity that each subhalo had when it was first accreted onto its host halo. This suggests that satellite galaxies are significantly more resilient to destruction than are dissipationless dark matter subhalos. We argue that while Nc does not provide a direct measure of the halo merger rate, it offers a powerful means to constrain the Halo Occupation Distribution and the spatial distribution of galaxies within halos. Interpreted in this way, close pair counts provide a useful test of galaxy formation processes on < 100 kpc scales.

0 citations0 reviewsNo code linkedOpen access
preprint2004arXiv

Multi-Phase Galaxy Formation: High Velocity Clouds and the Missing Baryon Problem

The standard treatment of cooling in Cold Dark Matter halos assumes that all of the gas within a ``cooling radius'' cools and contracts monolithically to fuel galaxy formation. Here we take into account the expectation that the hot gas in galactic halos is thermally unstable and prone to fragmentation during cooling and show that the implications are more far-reaching than previously expected: allowing multi-phase cooling fundamentally alters expectations about gas infall in halos and naturally explains the bright-end cutoff in the galaxy luminosity function. We argue that cooling should proceed via the formation of high-density, 10^4 K clouds, pressure-confined within a hot gas background. The background medium has a low density, and can survive as a stable corona with a long cooling time. The fraction of baryons contained in the residual hot core grows with halo mass because the cooling density increases, and this leads to an upper-mass limit in quiescent, non-merged galaxies of ~10^11 Msun. In this scenario, galaxy formation is fueled by the infall of pressure-supported clouds. For Milky-Way-size systems, clouds of mass ~ 5x10^6 Msun that formed or merged within the last several Gyrs should still exist as a residual population in the halo, with a total mass in clouds of ~ 2 x 10^10 Msun. The mass of the Milky Way galaxy is explained naturally in this model, and is a factor of two smaller than would result in the standard treatment without feedback. We expect clouds in galactic halos to be ~ 1 kpc in size and to extend ~150 kpc from galactic centers. The predicted properties of clouds match well the observed radial velocities, angular sizes, column densities, and velocity widths of High Velocity Clouds around our Galaxy. The clouds also explain high-ion absorption systems at z<1.

0 citations0 reviewsNo code linkedOpen access
preprint2000arXiv

A universal angular momentum profile for galactic halos

[Abridged] We study the angular-momentum profiles of a statistical sample of halos drawn from a high-resolution N-body simulation of the LCDM cosmology. We find that the cumulative mass distribution of specific angular momentum, j, in a halo of mass Mv is well fit by a universal function, M(<j) = Mv μj/(j_0+j). This profile is defined by one shape parameter (μor j_0) in addition to the global spin parameter λ. It follows a power-law over most of the mass, and flattens at large j, with the flattening more pronounced for small values of μ. Compared to a uniform sphere in solid-body rotation, most halos have a higher fraction of their mass in the low- and high-j tails of the distribution. The spatial distribution of angular momentum in halos tends to be cylindrical and is well-aligned within each halo for ~80% of the halos. We investigate two ideas for the origin of this profile. The first is based on a revised version of linear tidal-torque theory combined with extended Press-Schechter mass accretion, and the second focuses on j transport in minor mergers. Finally, we briefly explore implications of the M(<j) profile on the formation of galactic disks assuming that j is conserved during an adiabatic baryonic infall. The implied gas density profile deviates from an exponential disk, with a higher density at small radii and a tail extending to large radii. The steep central density profiles may imply disk scale lengths that are smaller than observed. This is reminiscent of the "angular-momentum problem" seen in hydrodynamic simulations, even though we have assumed perfect j conservation. A possible solution is to associate the central excesses with bulge components and the outer regions with extended gaseous disks.

0 citations0 reviewsNo code linkedOpen access
preprint2000arXiv

Probing Galaxy Formation with High Energy Gamma-Rays

We discuss how measurements of the absorption of $γ$-rays from GeV to TeV energies via pair production on the extragalactic background light (EBL) can probe important issues in galaxy formation. We use semi-analytic models (SAMs) of galaxy formation, set within the hierarchical structure formation scenario, to obtain predictions of the EBL for 0.1-1000$μ$m. SAMs incorporate simplified physical treatments of the key processes of galaxy formation --- including gravitational collapse and merging of dark matter halos, gas cooling and dissipation, star formation, supernova feedback and metal production --- and have been shown to reproduce key observations at low and high redshift. Here we also introduce improved modelling of the spectral energy distributions in the mid-to-far-IR arising from emission by dust grains. Assuming a flat \lcdm cosmology with $Ω_m=0.3$ and Hubble parameter $h=0.65$, we investigate the consequences of variations in input assumptions such as the stellar initial mass function (IMF) and the efficiency of converting cold gas into stars. We conclude that observational studies of the absorption of $γ$-rays with energies from 10s of Gev to 10s of TeV will help to determine the EBL, and also help to explain its origin by constraining some of the most uncertain features of galaxy formation theory, including the IMF, the history of star formation, and the reprocessing of light by dust.

0 citations0 reviewsNo code linkedOpen access
preprint2000arXiv

Profiles of dark haloes: evolution, scatter, and environment

We study dark-matter halo density profiles in a high-resolution N-body simulation of an LCDM cosmology. Our statistical sample contains ~5000 haloes in the range 10^{11}-10^{14} M_sun. The profiles are parameterized by an NFW form with two parameters, an inner radius r_s and a virial radius r_v, and we define the halo concentration c_v = r_v/r_s. First, we find that, for a given halo mass, the redshift dependence of the median concentration is c_v ~ 1/(1+z), corresponding to a roughly constant r_s with redshift. We present an improved analytic treatment of halo formation that fits the measured relations between halo parameters and their redshift dependence. The implications are that high-redshift galaxies are predicted to be more extended and dimmer than expected before. Second, we find that the scatter in log(c_v) is ~0.18, corresponding to a scatter in maximum rotation velocities of dV/V ~ 0.12. We discuss implications for modelling the Tully-Fisher relation, which has a smaller reported intrinsic scatter. Third, haloes in dense environments tend to be more concentrated than isolated haloes. These results suggest that c_v is an essential parameter for the theory of galaxy modelling, and we briefly discuss implications for the universality of the Tully-Fisher relation, the formation of low surface brightness galaxies, and the origin of the Hubble sequence.

0 citations0 reviewsNo code linkedOpen access
preprint1996arXiv

Non-Gaussian Fluctuations and Primordial Black Holes from Inflation

We explore the role of non-Gaussian fluctuations in primordial black hole (PBH) formation and show that the standard Gaussian assumption, used in all PBH formation papers to date, is not justified. Since large spikes in power are usually associated with flat regions of the inflaton potential, quantum fluctuations become more important in the field dynamics, leading to mode-mode coupling and non-Gaussian statistics. Moreover, PBH production requires several sigma (rare) fluctuations in order to prevent premature matter dominance of the universe, so we are necessarily concerned with distribution tails, where any intrinsic skewness will be especially important. We quantify this argument by using the stochastic slow-roll equation and a relatively simple analytic method to obtain the final distribution of fluctuations. We work out several examples with toy models that produce PBH's, and show that the naive Gaussian assumption can result in errors of many orders of magnitude. For models with spikes in power, our calculations give sharp cut-offs in the probability of large positive fluctuations, meaning that Gaussian distributions would vastly over-produce PBH's. The standard results that link inflation-produced power spectra and PBH number densities must then be reconsidered, since they rely quite heavily on the Gaussian assumption. We point out that since the probability distributions depend on the nature of the potential, it is impossible to obtain results for general models. However, calculating the distribution of fluctuations for any specific model seems to be relatively straightforward, at least in the single inflaton case.

0 citations0 reviewsNo code linkedOpen access