Source author record

Alyson M. Brooks

Alyson M. Brooks appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record

Catalog footprint

What is connected

27works
7topics
4close collaborators

Actions

Connect this record

Log in to claim

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

27 published item(s)

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.

preprint2022arXiv

A hidden population of massive black holes in simulated dwarf galaxies

We explore the characteristics of actively accreting MBHs within dwarf galaxies in the \textsc{Romulus25} cosmological hydrodynamic simulation. We examine the MBH occupation fraction, x-ray active fractions, and AGN scaling relations within dwarf galaxies of stellar mass $10^{8} < M_{\rm star} < 10^{10} M_\odot$ out to redshift $z=2$. In the local universe, the MBH occupation fraction is consistent with observed constraints, dropping below unity at $M_{\rm star} < 3\times10^{10} M_{\odot}$, $M_{\rm 200} < 3\times10^{11} M_\odot$. Local dwarf AGN in \textsc{Romulus25} follow observed scaling relations between AGN x-ray luminosity, stellar mass, and star formation rate, though they exhibit slightly higher active fractions and number densities than comparable x-ray observations. Since $z=2$, the MBH occupation fraction has decreased, the population of dwarf AGN has become overall less luminous, and as a result, the overall number density of dwarf AGN has diminished. We predict the existence of a large population of MBHs in the local universe with low x-ray luminosities and high contamination from x-ray binaries and the hot interstellar medium that are undetectable by current x-ray surveys. These hidden MBHs make up $76\%$ of all MBHs in local dwarf galaxies, and include many MBHs that are undermassive relative to their host galaxy's stellar mass. Their detection relies not only on greater instrument sensitivity but on better modeling of x-ray contaminants or multi-wavelength surveys. Our results indicate dwarf AGN were substantially more active in the past despite being low-luminosity today, and indicate future deep x-ray surveys may uncover many hidden MBHs in dwarf galaxies out to at least $z=2$.

preprint2021arXiv

The Formation of Isolated Ultra-Diffuse Galaxies in Romulus25

We use the \textsc{Romulus25} cosmological simulation volume to identify the largest-ever simulated sample of {\it field} ultra-diffuse galaxies (UDGs). At $z=0$, we find that isolated UDGs have average star formation rates, colors, and virial masses for their stellar masses and environment. UDGs have moderately elevated HI masses, being 70\% (300\%) more HI-rich than typical isolated dwarf galaxies at luminosities brighter (fainter) than M$_\mathrm{B}$=-14. However, UDGs are consistent with the general isolated dwarf galaxy population and make up $\sim$20\% of all field galaxies with 10$^7$<M$_\star$/M$_\odot$<10$^{9}$. The HI masses, effective radii, and overall appearances of our UDGs are consistent with existing observations of field UDGs, but we predict that many isolated UDGs have been missed by current surveys. Despite their isolation at $z=0$, the UDGs in our sample are the products of major mergers. Mergers are no more common in UDG than non-UDG progenitors, but mergers that create UDGs tend to happen earlier - almost never occurring after $z=1$, produce a temporary boost in spin, and cause star formation to be redistributed to the outskirts of galaxies, resulting in lower central star formation rates. The centers of the galaxies fade as their central stellar populations age, but their global star formation rates are maintained through bursts of star formation at larger radii, producing steeper negative g-r color gradients. This formation channel is unique relative to other proposals for UDG formation in isolated galaxies, demonstrating that UDGs can potentially be formed through multiple mechanisms.

preprint2020arXiv

Roll of the Dice: A Stochastically Sampled IMF Alters the Stellar Content of Simulated Dwarf Galaxies

Cosmological simulations are reaching the resolution necessary to study ultra-faint dwarf galaxies. Observations indicate that in small populations, the stellar initial mass function (IMF) is not fully populated; rather, stars are sampled in a way that can be approximated as coming from an underlying probability density function. To ensure the accuracy of cosmological simulations in the ultra-faint regime, we present an improved treatment of the IMF. We implement a self-consistent, stochastically populated IMF in cosmological hydrodynamic simulations. We test our method using high-resolution simulations of a Milky Way halo, run to $z=6$, yielding a sample of nearly 100 galaxies. We also use an isolated dwarf galaxy to investigate the resulting systematic differences in galaxy properties. We find that a stochastic IMF in simulations makes feedback burstier, strengthening feedback, and quenching star formation earlier in small dwarf galaxies. For galaxies in halos with mass $\lesssim10^{8.5}$ M$_\odot$, a stochastic IMF typically leads to lower stellar mass compared to a continuous IMF, sometimes by more than an order of magnitude. We show that existing methods of ensuring discrete supernovae incorrectly determine the mass of the star particle and its associated feedback. This leads to overcooling of surrounding gas, with at least ${\sim}$10 per cent higher star formation and ${\sim}$30 per cent higher cold gas content. Going forward, to accurately model dwarf galaxies and compare to observations, it will be necessary to incorporate a stochastically populated IMF that samples the full spectrum of stellar masses.

preprint2020arXiv

The Diversity and Variability of Star Formation Histories in Models of Galaxy Evolution

Understanding the variability of galaxy star formation histories (SFHs) across a range of timescales provides insight into the underlying physical processes that regulate star formation within galaxies. We compile the SFHs of galaxies at $z=0$ from an extensive set of models, ranging from cosmological hydrodynamical simulations (Illustris, IllustrisTNG, Mufasa, Simba, EAGLE), zoom simulations (FIRE-2, g14, and Marvel/Justice League), semi-analytic models (Santa Cruz SAM) and empirical models (UniverseMachine), and quantify the variability of these SFHs on different timescales using the power spectral density (PSD) formalism. We find that the PSDs are well described by broken power-laws, and variability on long timescales ($\gtrsim1$ Gyr) accounts for most of the power in galaxy SFHs. Most hydrodynamical models show increased variability on shorter timescales ($\lesssim300$ Myr) with decreasing stellar mass. Quenching can induce $\sim0.4-1$ dex of additional power on timescales $>1$ Gyr. The dark matter accretion histories of galaxies have remarkably self-similar PSDs and are coherent with the in-situ star formation on timescales $>3$ Gyr. There is considerable diversity among the different models in their (i) power due to SFR variability at a given timescale, (ii) amount of correlation with adjacent timescales (PSD slope), (iii) evolution of median PSDs with stellar mass, and (iv) presence and locations of breaks in the PSDs. The PSD framework is a useful space to study the SFHs of galaxies since model predictions vary widely. Observational constraints in this space will help constrain the relative strengths of the physical processes responsible for this variability.

preprint2020arXiv

The Formation of Ultra-Diffuse Galaxies in the RomulusC Galaxy Cluster Simulation

We study the origins of 122 ultra-diffuse galaxies (UDGs) in the {\sc RomulusC} zoom-in cosmological simulation of a galaxy cluster (M$_{200} = 1.15\times10^{14}$ M$_{\odot}$), one of the only such simulations capable of resolving the evolution and structure of dwarf galaxies (M$_{\star} < 10^9$ M$_{\odot}$). We find broad agreement with observed cluster UDGs and predict that they are not separate from the overall cluster dwarf population. UDGs in cluster environments form primarily from dwarf galaxies that experienced early cluster in-fall and subsequent quenching due to ram pressure. The ensuing dimming of these dwarf galaxies due to passive stellar evolution results in a population of very low surface brightness galaxies that are otherwise typical dwarfs. UDGs and non-UDGs alike are affected by tidal interactions with the cluster potential. Tidal stripping of dark matter, as well as mass loss from stellar evolution, results in the adiabatic expansion of stars, particularly in the lowest mass dwarfs. High mass dwarf galaxies show signatures of tidal heating while low mass dwarfs that survive until $z=0$ typically have not experienced such impulsive interactions. There is little difference between UDGs and non-UDGs in terms of their dark matter halos, stellar morphology, colors, and location within the cluster. In most respects cluster UDG and non-UDGs alike are similar to isolated dwarf galaxies, except for the fact that they are typically quenched.

preprint2020arXiv

Ultra-faint dwarfs in a Milky Way context: Introducing the Mint Condition DC Justice League Simulations

We present results from the "Mint" resolution DC Justice League suite of Milky Way-like zoom-in cosmological simulations, which extend our study of nearby galaxies down into the ultra-faint dwarf (UFD) regime for the first time. The mass resolution of these simulations is the highest ever published for cosmological Milky Way zoom-in simulations run to $z=0$, with initial star (dark matter) particle masses of 994 (17900) M$_\odot$, and a force resolution of 87 pc. We study the surrounding dwarfs and UFDs, and find the simulations match the observed dynamical properties of galaxies with $-3 < M_V < -19$, and reproduce the scatter seen in the size-luminosity plane for $r_h\gtrsim200$ pc. We predict the vast majority of nearby galaxies will be observable by the Vera Rubin Observatory's co-added Legacy Survey of Space and Time (LSST). We additionally show that faint dwarfs with velocity dispersions $\lesssim5$ km/s result from severe tidal stripping of the host halo. We investigate the quenching of UFDs in a hydrodynamical Milky Way context, and find that the majority of UFDs are quenched prior to interactions with the Milky Way, though some of the quenched UFDs retain their gas until infall. Additionally these simulations yield some unique dwarfs that are the first of their kind to be simulated, e.g., an HI-rich field UFD, a late-forming UFD that has structural properties similar to Crater 2, as well as a compact dwarf satellite that has no dark matter at $z=0$.

preprint2018arXiv

Understanding Dwarf Galaxies in order to Understand Dark Matter

Much progress has been made in recent years by the galaxy simulation community in making realistic galaxies, mostly by more accurately capturing the effects of baryons on the structural evolution of dark matter halos at high resolutions. This progress has altered theoretical expectations for galaxy evolution within a Cold Dark Matter (CDM) model, reconciling many earlier discrepancies between theory and observations. Despite this reconciliation, CDM may not be an accurate model for our Universe. Much more work must be done to understand the predictions for galaxy formation within alternative dark matter models.

preprint2017arXiv

Going, going, gone dark: Quantifying the scatter in the faintest dwarf galaxies

We predict the stellar mass-halo mass (SMHM) relationship for dwarf galaxies and their satellites residing in halos down to M$_{halo} =$ 10$^7$ M$_{\odot}$ with 10$^4$ M$_{\odot} <$ M$_{star}$($z=0$) $< 10^8$ M$_{\odot}$, and quantify the predicted scatter in the relation at the low mass end, using cosmological simulations. The galaxies were drawn from a cosmological simulation of dwarf galaxies, run with the N-body + SPH code, ChaNGA, at a high resolution of 60 pc. For M$_{halo} > 10^9$ M$_{\odot}$, the simulated SMHM relationship agrees with literature determinations, including exhibiting a small scatter. However, the scatter in the SMHM relation increases dramatically for lower-mass halos. We find that some of this scatter is due to {\em dark dwarfs}, halos devoid of stars. However, even when only considering well-resolved halos that contain a stellar population, the scatter in stellar mass reaches nearly 1 dex for M$_{halo}$($z=0$) 10$^7$ M$_{\odot}$. Much of this scatter is due to including satellites of the dwarf galaxies that have had their halo masses reduced through tidal stripping. The fraction of dark dwarfs (those that contain no stars) increases substantially with decreasing halo mass. When these dark halos are considered, the true scatter in the SMHM at low masses is even larger. At the faintest end of the SMHM relation probed by our simulations, a galaxy cannot be assigned a unique halo mass based solely on its luminosity. We provide a formula to stochastically populate low-mass halos following our results. Our predicted large scatter at low halo masses increases the slope of the resulting stellar mass function on the ultra-faint dwarf galaxy scales currently being probed by such surveys as the Dark Energy Survey or the Hyper-Suprime Cam Subaru Strategic Program, and in the future by the Large Synoptic Survey Telescope.

preprint2016arXiv

Assessing Astrophysical Uncertainties in Direct Detection with Galaxy Simulations

We study the local dark matter velocity distribution in simulated Milky Way-mass galaxies, generated at high resolution with both dark matter and baryons. We find that the dark matter in the Solar neighborhood is influenced appreciably by the inclusion of baryons, increasing the speed of dark matter particles compared to dark matter-only simulations. The gravitational potential due to the presence of a baryonic disk increases the amount of high velocity dark matter, resulting in velocity distributions which are more similar to the Maxwellian Standard Halo Model than predicted from dark matter-only simulations. Further, the velocity structures present in baryonic simulations possess a greater diversity than expected from dark matter-only simulation. We show the impact on the direct detection experiments LUX, DAMA/Libra, and CoGeNT using our simulated velocity distributions, and explore how resolution and halo mass within the Milky Way's estimated mass range impact the results. A Maxwellian fit to the velocity distribution tends to overpredict the amount of dark matter in the high velocity tail, even with baryons, and thus leads to overly optimistic direct detection bounds on models which are dependent on this region of phase space for an experimental signal. Our work further demonstrates that it is critical to transform simulated velocity distributions to the lab frame of reference, due to the fact that velocity structure in the Solar neighborhood appears when baryons are included. There is more velocity structure present when baryons are included than in dark matter-only simulations. Even when baryons are included, the importance of the velocity structure is not as apparent in the Galactic frame of reference as in the Earth frame.

preprint2016arXiv

Search for Gamma-ray Emission from Dark Matter Annihilation in the Small Magellanic Cloud with the Fermi Large Area Telescope

The Small Magellanic Cloud (SMC) is the second-largest satellite galaxy of the Milky Way and is only 60 kpc away. As a nearby, massive, and dense object with relatively low astrophysical backgrounds, it is a natural target for dark matter indirect detection searches. In this work, we use six years of Pass 8 data from the Fermi Large Area Telescope to search for gamma-ray signals of dark matter annihilation in the SMC. Using data-driven fits to the gamma-ray backgrounds, and a combination of N-body simulations and direct measurements of rotation curves to estimate the SMC DM density profile, we found that the SMC was well described by standard astrophysical sources, and no signal from dark matter annihilation was detected. We set conservative upper limits on the dark matter annihilation cross section. These constraints are in agreement with stronger constraints set by searches in the Large Magellanic Cloud and approach the canonical thermal relic cross section at dark matter masses lower than 10 GeV in the $b\bar{b}$ and $τ^+τ^-$ channels.

preprint2015arXiv

Bulge Formation via Mergers in Cosmological Simulations

The latest generation of cosmological simulations are on the verge of being able to resolve the structure of bulges for the first time. Hence, we review the current state of bulge formation in cosmological simulations, and discuss open questions that can be addressed in the near future by simulators, with a particular focus on merger-driven bulge growth. Galaxy mergers have long been assumed to produce classical bulges in disk galaxies. Under this bulge-formation model, though, the high rates of mergers in Cold Dark Matter galaxy formation theory predict many more classical bulges than are observed. Furthermore, simulations of galaxy formation continue to generally produce too massive of bulges. Feedback offers a promising avenue for reducing merger-driven bulge growth by maintaining high gas fractions in galaxies and ejecting low-angular momentum gas driven to the centers of galaxies. After reviewing the results of relevant research that has been published to date, we use cosmological simulations to explore the ability of feedback to reduce or even prevent bulge growth during mergers. In dwarf galaxies, mergers actually reduce the central concentration of galaxies as the induced burst of star formation drives out low angular momentum material. This result shows the potential for feedback to reduce central mass growth. However, we also demonstrate that it is very difficult for current stellar feedback models to reproduce the small bulges observed in more massive disk galaxies like the Milky Way. We argue that feedback models need to be improved, or an additional source of feedback such as AGN is necessary to generate the required outflows.

preprint2015arXiv

Consequences of bursty star formation on galaxy observables at high redshifts

The star formation histories (SFHs) of dwarf galaxies are thought to be \emph{bursty}, with large -- order of magnitude -- changes in the star formation rate on timescales similar to O-star lifetimes. As a result, the standard interpretations of many galaxy observables (which assume a slowly varying SFH) are often incorrect. Here, we use the SFHs from hydro-dynamical simulations to investigate the effects of bursty SFHs on sample selection and interpretation of observables and make predictions to confirm such SFHs in future surveys. First, because dwarf galaxies' star formation rates change rapidly, the mass-to-light ratio is also changing rapidly in both the ionizing continuum and, to a lesser extent, the non-ionizing UV continuum. Therefore, flux limited surveys are highly biased toward selecting galaxies in the \emph{burst} phase and very deep observations are required to detect all dwarf galaxies at a given stellar mass. Second, we show that a $\log_{10}[νL_ν(1500{\rm Å})/L_{{\rm H}α}]>2.5$ implies a very recent quenching of star formation and can be used as evidence of stellar feedback regulating star formation. Third, we show that the ionizing continuum can be significantly higher than when assuming a constant SFH, which can affect the interpretation of nebular emission line equivalent widths and direct ionizing continuum detections. Finally, we show that a star formation rate estimate based on continuum measurements only (and not on nebular tracers such as the hydrogen Balmer lines) will not trace the rapid changes in star formation and will give the false impression of a star-forming main sequence with low dispersion.

preprint2015arXiv

Particle tagging and its implications for stellar population dynamics

We establish a controlled comparison between the properties of galactic stellar halos obtained with hydrodynamical simulations and with `particle tagging'. Tagging is a fast way to obtain stellar population dynamics: instead of tracking gas and star formation, it `paints' stars directly onto a suitably defined subset of dark matter particles in a collisionless, dark-matter-only simulation.Our study shows that there are conditions under which particle tagging generates good fits to the hydrodynamical stellar density profiles of a central Milky-Way-like galaxy and its most prominent substructure. Phase-space diffusion processes are crucial to reshaping the distribution of stars in infalling spheroidal systems and hence the final stellar halo. We conclude that the success of any particular tagging scheme hinges on this diffusion being taken into account, at a minimum by making use of `live' tagging schemes, in which particles are regularly tagged throughout the evolution of a galaxy.

preprint2015arXiv

Search for Gamma-ray Emission from Dark Matter Annihilation in the Large Magellanic Cloud with the Fermi Large Area Telescope

At a distance of 50 kpc and with a dark matter mass of $\sim10^{10}$ M$_{\odot}$, the Large Magellanic Cloud (LMC) is a natural target for indirect dark matter searches. We use five years of data from the Fermi Large Area Telescope (LAT) and updated models of the gamma-ray emission from standard astrophysical components to search for a dark matter annihilation signal from the LMC. We perform a rotation curve analysis to determine the dark matter distribution, setting a robust minimum on the amount of dark matter in the LMC, which we use to set conservative bounds on the annihilation cross section. The LMC emission is generally very well described by the standard astrophysical sources, with at most a $1-2σ$ excess identified near the kinematic center of the LMC once systematic uncertainties are taken into account. We place competitive bounds on the dark matter annihilation cross section as a function of dark matter particle mass and annihilation channel.

preprint2014arXiv

Delayed Star Formation in Isolated Dwarf Galaxies: HST Star Formation History of the Aquarius Dwarf Irregular

We have obtained deep images of the highly isolated (d = 1 Mpc) Aquarius dwarf irregular galaxy (DDO 210) with the Hubble Space Telescope (HST) Advanced Camera for Surveys (ACS). The resulting color-magnitude diagram (CMD) reaches more than a magnitude below the oldest main-sequence turnoff, allowing us to derive the star formation history (SFH) over the entire lifetime of the galaxy with a timing precision of ~10% of the lookback time. Using a maximum likelihood fit to the CMD we find that only ~10% of all star formation in Aquarius took place more than 10 Gyr ago (lookback time equivalent to redshift z ~2). The star formation rate increased dramatically ~6-8 Gyr ago (z ~ 0.7-1.1) and then declined until the present time. The only known galaxy with a more extreme confirmed delay in star formation is Leo A, a galaxy of similar M(HI)/M(stellar), dynamical mass, mean metallicity, and degree of isolation. The delayed stellar mass growth in these galaxies does not track the mean dark matter accretion rate from CDM simulations. The similarities between Leo A and Aquarius suggest that if gas is not removed from dwarf galaxies by interactions or feedback, it can linger for several gigayears without cooling in sufficient quantity to form stars efficiently. We discuss possible causes for the delay in star formation including suppression by reionization and late-time mergers. We find reasonable agreement between our measured SFHs and select cosmological simulations of isolated dwarfs. Because star formation and merger processes are both stochastic in nature, delayed star formation in various degees is predicted to be a characteristic (but not a universal) feature of isolated small galaxies.

preprint2014arXiv

Faint dwarfs as a test of DM models: WDM vs. CDM

We use high resolution Hydro$+$N-Body cosmological simulations to compare the assembly and evolution of a small field dwarf (stellar mass ~ 10$^{6-7}$ M$\odot$, total mass 10$^{10}$ M$\odot$ in $Λ$ dominated CDM and 2keV WDM cosmologies. We find that star formation (SF) in the WDM model is reduced and delayed by 1-2 Gyr relative to the CDM model, independently of the details of SF and feedback. Independent of the DM model, but proportionally to the SF efficiency, gas outflows lower the central mass density through `dynamical heating', such that all realizations have circular velocities $<$ 20kms at 500$~$pc, in agreement with local kinematic constraints. As a result of dynamical heating, older stars are less centrally concentrated than younger stars, similar to stellar population gradients observed in nearby dwarf galaxies. Introducing an important diagnostic of SF and feedback models, we translate our simulations into artificial color-magnitude diagrams and star formation histories in order to directly compare to available observations. The simulated galaxies formed most of their stars in many $\sim$10 Myr long bursts. The CDM galaxy has a global SFH, HI abundance and Fe/H and alpha-elements distribution well matched to current observations of dwarf galaxies. These results highlight the importance of directly including `baryon physics' in simulations when 1) comparing predictions of galaxy formation models with the kinematics and number density of local dwarf galaxies and 2) differentiating between CDM and non-standard models with different DM or power spectra.

preprint2014arXiv

The effect of models of the interstellar media on the central mass distribution of galaxies

We compare the central mass distribution of galaxies simulated with three different models of the interstellar medium (ISM) with increasing complexity: primordial (H+He) cooling down to 10^4K, additional cooling via metal lines and to lower temperatures, and molecular hydrogen (H_2) with shielding of atomic and molecular hydrogen, in addition to metal line cooling. In order to analyze the effect of these models, we follow the evolution of four field galaxies with V_peak < 120 km/s to a redshift of zero using high-resolution Smoothed Particle Hydrodynamic simulations in a fully cosmological LCDM context. The spiral galaxies produced in simulations with either primordial cooling or H_2 physics have realistic, rising rotation curves. In contrast, the simulations with metal line cooling and otherwise similar feedback and star formation produced galaxies with the peaked rotation curves typical of most previous LCDM simulations of spiral galaxies. The less-massive bulges and non-peaked rotation curves in the galaxies simulated with primordial cooling or H_2 are linked to changes in the angular momentum distribution of the baryons. These galaxies had smaller amounts of low-angular momentum baryons because of increased gas loss from stellar feedback. When there is only primordial cooling, the star forming gas is hotter and the feedback-heated gas cools more slowly than when metal line cooling is included and so requires less energy to be expelled. When H_2 is included, the accompanying shielding produces large amounts of clumpy, cold gas where H_2 forms. Star formation in clumpy gas results in more concentrated supernova feedback and greater efficiency of mass loss. The higher feedback efficiency causes a decrease of low-angular momentum material. (abridged)

preprint2014arXiv

The Milky Way Tomography with SDSS. V. Mapping the Dark Matter Halo

We present robust constraints from the Sloan Digital Sky Survey (SDSS) on the shape and distribution of the dark matter halo within the Milky Way (MW). Using the number density distribution and kinematics of SDSS halo stars, we probe the dark matter distribution to heliocentric distances exceeding 10 kpc and galactocentric distances exceeding 20 kpc. Our analysis utilizes Jeans equations to generate two-dimensional acceleration maps throughout the volume; this approach is thoroughly tested on a cosmologically derived N-body+SPH simulation of a MW-like galaxy. We show that the known accelerations (gradients of the gravitational potential) can be successfully recovered in such a realistic system. Leveraging the baryonic gravitational potential derived by Bovy & Rix (2013), we show that the gravitational potential implied by the SDSS observations cannot be explained, assuming Newtonian gravity, by visible matter alone: the gravitational force experienced by stars at galactocentric distances of 20 kpc is as much as three times stronger than what can be attributed to purely visible matter. We also show that the SDSS data provide a strong constraint on the shape of the dark matter halo potential. Within galactocentric distances of 20 kpc, the dark matter halo potential is well described as an oblate halo with axis ratio qDM=0.7+/-0.1; this corresponds to an axis ratio qDM=0.4+/-0.1 for the dark matter density distribution. Because of our precise two-dimensional measurements of the acceleration of the halo stars, we can reject several MOND models as an explanation of the observed behavior.

preprint2014arXiv

Why Baryons Matter: The Kinematics of Dwarf Spheroidal Satellites

We use high resolution cosmological simulations of Milky Way-mass galaxies that include both baryons and dark matter to show that baryonic physics (energetic feedback from supernovae and subsequent tidal stripping) significantly reduces the dark matter mass in the central regions of luminous satellite galaxies. The reduced central masses of the simulated satellites reproduce the observed internal dynamics of Milky Way and M31 satellites as a function of luminosity. We use these realistic satellites to update predictions for the observed velocity and luminosity functions of satellites around Milky Way-mass galaxies when baryonic effects are accounted for. We also predict that field dwarf galaxies in the same luminosity range as the Milky Way classical satellites should not exhibit velocities as low as the satellites, since the field dwarfs do not experience tidal stripping. Additionally, the early formation times of the satellites compared to field galaxies at the same luminosity may be apparent in the star formation histories of the two populations. Including baryonic physics in Cold Dark Matter models naturally explains the observed low dark matter densities in the Milky Way's dwarf spheroidal population. Our simulations therefore resolve the tension between kinematics predicted in Cold Dark Matter theory and observations of satellites, without invoking alternative forms of dark matter.

preprint2013arXiv

A Baryonic Solution to the Missing Satellites Problem

It has been demonstrated that the inclusion of baryonic physics can alter the dark matter densities in the centers of low-mass galaxies, making the central dark matter slope more shallow than predicted in pure cold dark matter simulations. This flattening of the dark matter profile can occur in the most luminous subhalos around Milky Way-mass galaxies. Zolotov et al. (2012) have suggested a correction to be applied to the central masses of dark matter-only satellites in order to mimic the affect of (1) the flattening of the dark matter cusp due to supernova feedback in luminous satellites, and (2) enhanced tidal stripping due to the presence of a baryonic disk. In this paper, we apply this correction to the z=0 subhalo masses from the high resolution, dark matter-only Via Lactea II (VL2) simulation, and find that the number of massive subhalos is dramatically reduced. After adopting a stellar mass to halo mass relationship for the VL2 halos, and identifying subhalos that are (1) likely to be destroyed by stripping and (2) likely to have star formation suppressed by photo-heating, we find that the number of massive, luminous satellites around a Milky Way-mass galaxy is in agreement with the number of observed satellites around the Milky Way or M31. We conclude that baryonic processes have the potential to solve the missing satellites problem.

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.

preprint2013arXiv

The Comparative Chemical Evolution of an Isolated Dwarf Galaxy: A VLT and Keck Spectroscopic Survey of WLM

Building on our previous spectroscopic and photometric analysis of the isolated Local Group dwarf irregular (dIrr) galaxy WLM, we present a comparison of the metallicities of its RGB stars with respect to the well studied Local Group dwarf spheroidal galaxies (dSphs) and Magellanic Clouds. We calculate a mean metallicity of [Fe/H]$ = -1.28 \pm 0.02$, and intrinsic spread in metallicity of $σ= 0.38 \pm 0.04$ dex, similar to the mean and spread observed in the massive dSph Fornax and the Small Magellanic Cloud. Thus, despite its isolated environment the global metallicity still follows expectations for WLM's mass and its global chemical evolution is similar to other nearby luminous dwarf galaxies (gas-rich or gas-poor). The data also show a radial gradient in [Fe/H] of $\rm{d[Fe/H]/dr_{c}} = -0.04 \pm 0.04$ dex $\rm{r_{c}^{-1}}$, which is flatter than that seen in the unbiased and spatially extended surveys of dSphs. Comparison of the spatial distribution of [Fe/H] in WLM, the Magellanic Clouds, and a sample of Local Group dSphs, shows an apparent dichotomy in the sense that the dIrrs have statistically flatter radial [Fe/H] gradients than the low angular momentum dSphs. The correlation between angular momentum and radial metallicity gradient is further supported when considering the Local Group dEs. This chemodynamic relationship offers a new and useful constraint for environment driven dwarf galaxy evolution models in the Local Group.

preprint2012arXiv

Baryons Matter: Why Luminous Satellite Galaxies Have Reduced Central Masses

Using high resolution cosmological hydrodynamical simulations of Milky Way-massed disk galaxies, we demonstrate that supernovae feedback and tidal stripping lower the central masses of bright (-15 < M_V < -8) satellite galaxies. These simulations resolve high density regions, comparable to giant molecular clouds, where stars form. This resolution allows us to adopt a prescription for H_2 formation and destruction that ties star formation to the presence of shielded, molecular gas. Before infall, supernova feedback from the clumpy, bursty star formation captured by this physically motivated model leads to reduced dark matter (DM) densities and shallower inner density profiles in the massive satellite progenitors (Mvir > 10^9 Msun, Mstar > 10^7 Msun) compared to DM-only simulations. The progenitors of the lower mass satellites are unable to maintain bursty star formation histories, due to both heating at reionization and gas loss from initial star forming events, preserving the steep inner density profile predicted by DM-only simulations. After infall, tidal stripping acts to further reduce the central densities of the luminous satellites, particularly those that enter with cored dark matter halos, increasing the discrepancy in the central masses predicted by baryon+DM and DM-only simulations. We show that DM-only simulations, which neglect the baryonic effects described in this work, produce denser satellites with larger central velocities. We provide a simple correction to the central DM mass predicted for satellites by DM-only simulations. We conclude that DM-only simulations should be used with great caution when interpreting kinematic observations of the Milky Way's dwarf satellites.

preprint2012arXiv

Constraints on the Shape of the Milky Way Dark Matter Halo from Jeans Equations Applied to SDSS Data

We search for evidence of dark matter in the Milky Way by utilizing the stellar number density distribution and kinematics measured by the Sloan Digital Sky Survey (SDSS) to heliocentric distances exceeding ~10 kpc. We employ the cylindrically symmetric form of Jeans equations and focus on the morphology of the resulting acceleration maps, rather than the normalization of the total mass as done in previous, mostly local, studies. Jeans equations are first applied to a mock catalog based on a cosmologically derived N-body + SPH simulation, and the known acceleration (gradient of gravitational potential) is successfully recovered. The same simulation is also used to quantify the impact of dark matter on the total acceleration. We use Galfast, a code designed to quantitatively reproduce SDSS measurements and selection effects, to generate a synthetic stellar catalog. We apply Jeans equations to this catalog and produce two-dimensional maps of stellar acceleration. These maps reveal that in a Newtonian framework, the implied gravitational potential cannot be explained by visible matter alone. The acceleration experienced by stars at galactocentric distances of ~20 kpc is three times larger than what can be explained by purely visible matter. The application of an analytic method for estimating the dark matter halo axis ratio to SDSS data implies an oblate halo with q_DM = 0.47 +/- 0.14 within the same distance range. These techniques can be used to map the dark matter halo to much larger distances from the Galactic center using upcoming deep optical surveys, such as LSST.

preprint2012arXiv

The Resolved Structure and Dynamics of an Isolated Dwarf Galaxy: A VLT and Keck Spectroscopic Survey of WLM

We present spectroscopic data for 180 red giant branch stars in the isolated dwarf irregular galaxy WLM. Observations of the Calcium II triplet lines in spectra of RGB stars covering the entire galaxy were obtained with FORS2 at the VLT and DEIMOS on Keck II allowing us to derive velocities, metallicities, and ages for the stars. With accompanying photometric and radio data we have measured the structural parameters of the stellar and gaseous populations over the full galaxy. The stellar populations show an intrinsically thick configuration with $0.39 \leq q_{0} \leq 0.57$. The stellar rotation in WLM is measured to be $17 \pm 1$ km s$^{-1}$, however the ratio of rotation to pressure support for the stars is $V/σ\sim 1$, in contrast to the gas whose ratio is seven times larger. This, along with the structural data and alignment of the kinematic and photometric axes, suggests we are viewing WLM as a highly inclined oblate spheroid. Stellar rotation curves, corrected for asymmetric drift, are used to compute a dynamical mass of $4.3\pm 0.3\times10^{8} $M$_{\odot}$ at the half light radius ($r_{h} = 1656 \pm 49$ pc). The stellar velocity dispersion increases with stellar age in a manner consistent with giant molecular cloud and substructure interactions producing the heating in WLM. Coupled with WLM's isolation, this suggests that the extended vertical structure of its stellar and gaseous components and increase in stellar velocity dispersion with age are due to internal feedback, rather than tidally driven evolution. These represent some of the first observational results from an isolated Local Group dwarf galaxy which can offer important constraints on how strongly internal feedback and secular processes modulate SF and dynamical evolution in low mass isolated objects.

preprint2010arXiv

Misaligned Angular Momentum in Hydrodynamic Cosmological Simulations: Warps, Outer Discs, and Thick Discs

We present a detailed analysis of a disc galaxy forming in a high-resolution fully cosmological simulation to investigate the nature of the outer regions of discs and their relevance for the disc formation process. Specifically, we focus on the phenomenon of misaligned disc components and find that the outer disc warp is a consequence of the misalignment between the inner disc and the surrounding hot gaseous halo. As the infalling cold gas sinks toward the centre of the galaxy, it is strongly torqued by the hot gas halo. By the time the fresh gas reaches the central disc-forming region its angular momentum is completely aligned with the spin of the hot gas halo. If the spin of the hot gas halo, in turn, is not aligned with that of the inner disc, a misaligned outer disc forms comprised of newly accreted material. The inner and outer components are misaligned with each other because they respond differently to infalling substructure and accretion. The warped disc feeds the main gas disc due to viscous angular momentum losses, but small amounts of star formation in the warp itself form a low-metallicity thick disc. We show that observations of resolved stellar populations in warped galaxies in the local universe could provide evidence for the presence of these processes and therefore indirectly reveal ongoing gas accretion and the existence of hot gas halos.