Showing 1–50 of 50 results for author: Christensen, C

Bursting with Feedback: The Relationship between Feedback Model and Bursty Star Formation Histories in Dwarf Galaxies

Authors: Bianca Azartash-Namin, Anna Engelhardt, Ferah Munshi, B. W. Keller, Alyson M. Brooks, Jordan Van Nest, Charlotte R. Christensen, Tom Quinn, James Wadsley

Abstract: We use high-resolution cosmological simulations to compare the effect of bursty star formation histories on dwarf galaxy structure for two different subgrid supernovae (SNe) feedback models in dwarf galaxies with stellar masses from $5000 <$ M$_*$/M$_\odot$ $< 10^{9}$. Our simulations are run using two distinct supernova feedback models: superbubble and blastwave. We show that both models are capa… ▽ More We use high-resolution cosmological simulations to compare the effect of bursty star formation histories on dwarf galaxy structure for two different subgrid supernovae (SNe) feedback models in dwarf galaxies with stellar masses from $5000 <$ M$_*$/M$_\odot$ $< 10^{9}$. Our simulations are run using two distinct supernova feedback models: superbubble and blastwave. We show that both models are capable of producing galaxies that are cored and reproduce observed scaling relations for metallicity, luminosity, mass, and size. We show that continuous bursty star formation and the resulting stellar feedback are able to sustain dark matter cores in the higher dwarf galaxy mass regime, while the majority of ultra-faint and classical dwarfs retain cuspy central dark matter density profiles. We find that both subgrid SN models are able to create bursty star formation histories. We find that effective core formation peaks at M$_*$/M$_\odot$ $\simeq 5 \times 10^{-3}$ for both feedback models. Galaxies simulated with superbubble feedback peak at lower mean burstiness values relative to blastwave feedback, indicating that core formation in the superbubble sample may be less motivated by the burstiness of star formation. △ Less

Submitted 11 January, 2024; originally announced January 2024.

Comments: 18 pages, 7 figures

Environment Matters: Predicted Differences in the Stellar Mass--Halo Mass Relation and History of Star Formation for Dwarf Galaxies

Authors: Charlotte R. Christensen, Alyson Brooks, Ferah Munshi, Claire Riggs, Jordan Van Nest, Hollis Akins, Thomas R Quinn, Lucas Chamberland

Abstract: We are entering an era in which we will be able to detect and characterize hundreds of dwarf galaxies within the Local Volume. It is already known that a strong dichotomy exists in the gas content and star formation properties of field dwarf galaxies versus satellite dwarfs of larger galaxies. In this work, we study the more subtle differences that may be detectable in galaxies as a function of di… ▽ More We are entering an era in which we will be able to detect and characterize hundreds of dwarf galaxies within the Local Volume. It is already known that a strong dichotomy exists in the gas content and star formation properties of field dwarf galaxies versus satellite dwarfs of larger galaxies. In this work, we study the more subtle differences that may be detectable in galaxies as a function of distance from a massive galaxy, such as the Milky Way. We compare smoothed particle hydrodynamic simulations of dwarf galaxies formed in a Local Volume-like environment (several Mpc away from a massive galaxy) to those formed nearer to Milky Way-mass halos. We find that the impact of environment on dwarf galaxies extends even beyond the immediate region surrounding Milky Way-mass halos. Even before being accreted as satellites, dwarf galaxies near a Milky Way-mass halo tend to have higher stellar masses for their halo mass than more isolated galaxies. Dwarf galaxies in high-density environments also tend to grow faster and form their stars earlier. We show observational predictions that demonstrate how these trends manifest in lower quenching rates, higher HI fractions, and bluer colors for more isolated dwarf galaxies. △ Less

Submitted 8 November, 2023; originally announced November 2023.

Comments: 26 pages, 14 figures, submitted to ApJ

The Scatter Matters: Circumgalactic Metal Content in the Context of the $M-σ$ Relation

Authors: N. Nicole Sanchez, Jessica K. Werk, Charlotte Christensen, O. Grace Telford, Michael Tremmel, Thomas Quinn, Jennifer Mead, Ray Sharma, Alyson Brooks

Abstract: The interaction between supermassive black hole (SMBH) feedback and the circumgalactic medium (CGM) continues to be an open question in galaxy evolution. In our study, we use SPH simulations to explore the impact of SMBH feedback on galactic metal retention and the motion of metals and gas into and through the CGM of L$_{*}$ galaxies. We examine 140 galaxies from the 25 Mpc cosmological volume, Ro… ▽ More The interaction between supermassive black hole (SMBH) feedback and the circumgalactic medium (CGM) continues to be an open question in galaxy evolution. In our study, we use SPH simulations to explore the impact of SMBH feedback on galactic metal retention and the motion of metals and gas into and through the CGM of L$_{*}$ galaxies. We examine 140 galaxies from the 25 Mpc cosmological volume, Romulus25, with stellar masses between 3 $\times$ 10$^{9}$ - 3 $\times$ 10$^{11}$ M$_{\odot}$. We measure the fraction of metals remaining in the ISM and CGM of each galaxy, and calculate the expected mass of its SMBH based on the $M-σ$ relation. The deviation of each SMBH from its expected mass, $ΔM_{BH}$ is compared to the potential of its host via $σ$. We find that SMBHs with accreted mass above the empirical $M-σ$ relation are about 15\% more effective at removing metals from the ISM than under-massive SMBHs in star forming galaxies. Over-massive SMBHs suppress the overall star formation of their host galaxies and more effectively move metals from the ISM into the CGM. However, we see little evidence for the evacuation of gas from their halos, in contrast with other simulations. Finally, we predict that C IV column densities in the CGM of L$_{*}$ galaxies may depend on host galaxy SMBH mass. Our results show that the scatter in the low mass end of $M-σ$ relation may indicate how effective a SMBH is at the local redistribution of mass in its host galaxy. △ Less

Submitted 11 May, 2023; originally announced May 2023.

Comments: 16 pages, 14 figures, Submitted to ApJ

The Role of Mass and Environment on Satellite distributions around Milky Way analogs in the Romulus25 simulation

Authors: Jordan Van Nest, Ferah Munshi, Charlotte Christensen, Alyson M. Brooks, Michael Tremmel, Thomas R. Quinn

Abstract: We study satellite counts and quenched fractions for satellites of Milky Way analogs in Romulus25, a large-volume cosmological hydrodynamic simulation. Depending on the definition of a Milky Way analog, we have between 66 and 97 Milky Way analogs in Romulus25, a 25 Mpc per-side uniform volume simulation. We use these analogs to quantify the effect of environment and host properties on satellite po… ▽ More We study satellite counts and quenched fractions for satellites of Milky Way analogs in Romulus25, a large-volume cosmological hydrodynamic simulation. Depending on the definition of a Milky Way analog, we have between 66 and 97 Milky Way analogs in Romulus25, a 25 Mpc per-side uniform volume simulation. We use these analogs to quantify the effect of environment and host properties on satellite populations. We find that the number of satellites hosted by a Milky Way analog increases predominantly with host stellar mass, while environment, as measured by the distance to a Milky Way-mass or larger halo, may have a notable impact in high isolation. Similarly, we find that the satellite quenched fraction for our analogs also increases with host stellar mass, and potentially in higher-density environments. These results are robust for analogs within 3 Mpc of another Milky Way-mass or larger halo, the environmental parameter space where the bulk of our sample resides. We place these results in the context of observations through comparisons to the Exploration of Local VolumE Satellites and Satellites Around Galactic Analogs surveys. Our results are robust to changes in Milky Way analog selection criteria, including those that mimic observations. Finally, as our samples naturally include Milky Way-Andromeda pairs, we examine quenched fractions in pairs vs isolated systems. We find potential evidence, though not conclusive, that pairs, defined as being within 1 Mpc of another Milky Way-mass or larger halo, may have higher satellite quenched fractions. △ Less

Submitted 26 September, 2023; v1 submitted 4 April, 2023; originally announced April 2023.

Stellar Migration and Chemical Enrichment in the Milky Way Disc: A Hybrid Model

Authors: James W. Johnson, David H. Weinberg, Fiorenzo Vincenzo, Jonathan C. Bird, Sarah R. Loebman, Alyson M. Brooks, Thomas R. Quinn, Charlotte R. Christensen, Emily J. Griffith

Abstract: We develop a hybrid model of galactic chemical evolution that combines a multi-ring computation of chemical enrichment with a prescription for stellar migration and the vertical distribution of stellar populations informed by a cosmological hydrodynamic disc galaxy simulation. Our fiducial model adopts empirically motivated forms of the star formation law and star formation history, with a gradien… ▽ More We develop a hybrid model of galactic chemical evolution that combines a multi-ring computation of chemical enrichment with a prescription for stellar migration and the vertical distribution of stellar populations informed by a cosmological hydrodynamic disc galaxy simulation. Our fiducial model adopts empirically motivated forms of the star formation law and star formation history, with a gradient in outflow mass loading tuned to reproduce the observed metallicity gradient. With this approach, the model reproduces many of the striking qualitative features of the Milky Way disc's abundance structure: (i) the dependence of the [O/Fe]-[Fe/H] distribution on radius $R_\text{gal}$ and midplane distance $|z|$; (ii) the changing shapes of the [O/H] and [Fe/H] distributions with $R_\text{gal}$ and $|z|$; (iii) a broad distribution of [O/Fe] at sub-solar metallicity and changes in the [O/Fe] distribution with $R_\text{gal}$, $|z|$, and [Fe/H]; (iv) a tight correlation between [O/Fe] and stellar age for [O/Fe] $>$ 0.1; (v) a population of young and intermediate-age $α$-enhanced stars caused by migration-induced variability in the Type Ia supernova rate; (vi) non-monotonic age-[O/H] and age-[Fe/H] relations, with large scatter and a median age of $\sim$4 Gyr near solar metallicity. Observationally motivated models with an enhanced star formation rate $\sim$2 Gyr ago improve agreement with the observed age-[Fe/H] and age-[O/H] relations, but worsen agreement with the observed age-[O/Fe] relation. None of our models predict an [O/Fe] distribution with the distinct bimodality seen in the observations, suggesting that more dramatic evolutionary pathways are required. All code and tables used for our models are publicly available through the Versatile Integrator for Chemical Evolution (VICE; https://pypi.org/project/vice). △ Less

Submitted 17 March, 2021; originally announced March 2021.

Comments: 24 pages; 18 figures; to be submitted to MNRAS; comments welcome

The Origins of Off-Centre Massive Black Holes in Dwarf Galaxies

Authors: Jillian M. Bellovary, Sarra Hayoune, Katheryn Chafla, Donovan Vincent, Alyson Brooks, Charlotte Christensen, Ferah Munshi, Michael Tremmel, Thomas R. Quinn, Jordan Van Nest, Serena K. Sligh, Michelle Luzuriaga

Abstract: Massive black holes often exist within dwarf galaxies, and both simulations and observations have shown that a substantial fraction of these may be off-center with respect to their hosts. We trace the evolution of off-center massive black holes (MBHs) in dwarf galaxies using cosmological hydrodynamical simulations, and show that the reason for off-center locations is mainly due to galaxy-galaxy me… ▽ More Massive black holes often exist within dwarf galaxies, and both simulations and observations have shown that a substantial fraction of these may be off-center with respect to their hosts. We trace the evolution of off-center massive black holes (MBHs) in dwarf galaxies using cosmological hydrodynamical simulations, and show that the reason for off-center locations is mainly due to galaxy-galaxy mergers. We calculate dynamical timescales and show that off-center MBHs are unlikely to sink to their galaxys' centers within a Hubble time, due to the shape of the hosts' potential wells and low stellar densities. These wandering MBHs are unlikely to be detected electromagnetically, nor is there a measurable dynamical effect on the galaxy's stellar population. We conclude that off-center MBHs may be common in dwarfs, especially if the mass of the MBH is small or the stellar mass of the host galaxy is large. However detecting them is extremely challenging, because their accretion luminosities are very low and they do not measurably alter the dynamics of their host galaxies. △ Less

Submitted 22 October, 2021; v1 submitted 18 February, 2021; originally announced February 2021.

Comments: updated to reflect MNRAS accepted version. Includes carbon footprint estimate due to supercomputer use

Quantifying scatter in galaxy formation at the lowest masses

Authors: Ferah Munshi, Alyson Brooks, Elaad Applebaum, Charlotte Christensen, Jordan P. Sligh, T. Quinn

Abstract: We predict the stellar mass -- halo mass (SMHM) relationship for dwarf galaxies, using simulated galaxies with peak halo masses of M$_{\rm peak} = 10^{11}$ M$_{\odot}$ down into the ultra-faint dwarf range to M$_{\rm peak} =$ 10$^7$ M$_{\odot}$. Our simulated dwarfs have stellar masses of M$_{\rm star} = $ 790 M$_{\odot}$ to $8.2 \times 10^8$ M$_{\odot}$, with corresponding $V$-band magnitudes fro… ▽ More We predict the stellar mass -- halo mass (SMHM) relationship for dwarf galaxies, using simulated galaxies with peak halo masses of M$_{\rm peak} = 10^{11}$ M$_{\odot}$ down into the ultra-faint dwarf range to M$_{\rm peak} =$ 10$^7$ M$_{\odot}$. Our simulated dwarfs have stellar masses of M$_{\rm star} = $ 790 M$_{\odot}$ to $8.2 \times 10^8$ M$_{\odot}$, with corresponding $V$-band magnitudes from $-2$ to $-18.5$. For M$_{\rm peak} > 10^{10}$ M$_{\odot}$, the simulated SMHM relationship agrees with literature determinations, including exhibiting a small scatter of 0.3 dex. However, the scatter in the SMHM relation increases for lower-mass halos. We first present results for well-resolved halos that contain a simulated stellar population, but recognize that whether a halo hosts a galaxy is inherently mass resolution dependent. We thus adopt a probabilistic model to populate "dark" halos below our resolution limit to predict an "intrinsic" slope and scatter for the SMHM relation. We fit linearly growing log-normal scatter in stellar mass, which grows to more than 1 dex at M$_{\rm peak}$ $=$ 10$^8$ M$_{\odot}$. 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. Instead, we provide a formula to stochastically populate low-mass halos following our results. Finally, we show that our growing log-normal scatter steepens the faint-end slope of the predicted stellar mass function. △ Less

Submitted 14 January, 2021; originally announced January 2021.

Comments: 16 pages, 5 figures; submitted to ApJ. arXiv admin note: text overlap with arXiv:1705.06286

One-Two Quench: A Double Minor Merger Scenario

Authors: N. Nicole Sanchez, Michael Tremmel, Jessica K. Werk, Andrew Pontzen, Charlotte Christensen, Thomas Quinn, Sarah Loebman, Akaxia Cruz

Abstract: Using the N-body+Smoothed particle hydrodynamics code, ChaNGa, we identify two merger-driven processes\textemdash disk disruption and supermassive black hole (SMBH) feedback\textemdash which work together to quench L$^*$ galaxies for over 7 Gyr. Specifically, we examine the cessation of star formation in a simulated Milky Way (MW) analog, driven by an interaction with two minor satellites. Both in… ▽ More Using the N-body+Smoothed particle hydrodynamics code, ChaNGa, we identify two merger-driven processes\textemdash disk disruption and supermassive black hole (SMBH) feedback\textemdash which work together to quench L$^*$ galaxies for over 7 Gyr. Specifically, we examine the cessation of star formation in a simulated Milky Way (MW) analog, driven by an interaction with two minor satellites. Both interactions occur within $\sim$100 Myr of each other, and the satellites both have masses 5 to 20 times smaller than that of their MW-like host galaxy. Using the genetic modification process of \cite{Roth2016}, we generate a set of four zoom-in, MW-mass galaxies all of which exhibit unique star formation histories due to small changes to their assembly histories. In two of these four cases, the galaxy is quenched by $z = 1$. Because these are controlled modifications, we are able to isolate the effects of two closely-spaced minor merger events, the relative timing of which determines whether the MW-mass main galaxy quenches. This one-two punch works to: 1. fuel the primary halo's supermassive black hole (SMBH) at its peak accretion rate; and 2. disrupt the cold, gaseous disk of the host galaxy. The end result is that feedback from the SMBH thoroughly and abruptly ends the galaxy's star formation by $z\approx1$. We search for and find a similar quenching event in {\sc Romulus25}, a hydrodynamical $(25\,\mathrm{Mpc})^3$ volume simulation, demonstrating that the mechanism is common enough to occur even in a small sample of MW-mass quenched galaxies at $z=0$. △ Less

Submitted 25 May, 2021; v1 submitted 11 September, 2020; originally announced September 2020.

Comments: 12 pages, 9 figures, accepted to ApJ

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

Authors: Elaad Applebaum, Alyson M. Brooks, Charlotte R. Christensen, Ferah Munshi, Thomas R. Quinn, Sijing Shen, Michael Tremmel

Abstract: 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) particl… ▽ More 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$. △ Less

Submitted 17 November, 2020; v1 submitted 25 August, 2020; originally announced August 2020.

Comments: 33 pages, 17 figures, ApJ accepted version

Quenching timescales of dwarf satellites around Milky Way-mass hosts

Authors: Hollis B. Akins, Charlotte R. Christensen, Alyson M. Brooks, Ferah Munshi, Elaad Applebaum, Anna Engelhardt, Lucas Chamberland

Abstract: Observations of the low-mass satellites in the Local Group have shown high fractions of gas-poor, quiescent galaxies relative to isolated dwarfs, implying that the host halo environment plays an important role in the quenching of dwarf galaxies. In this work, we present measurements of the quenched fractions and quenching timescales of dwarf satellite galaxies in the DC Justice League suite of 4 h… ▽ More Observations of the low-mass satellites in the Local Group have shown high fractions of gas-poor, quiescent galaxies relative to isolated dwarfs, implying that the host halo environment plays an important role in the quenching of dwarf galaxies. In this work, we present measurements of the quenched fractions and quenching timescales of dwarf satellite galaxies in the DC Justice League suite of 4 high-resolution cosmological zoom-in simulations of Milky Way-mass halos. We show that these simulations accurately reproduce the satellite luminosity functions of observed nearby galaxies, as well as the variation in satellite quenched fractions from $M_* \sim 10^5$ solar masses to $10^{10}$ solar masses. We then trace the histories of satellite galaxies back to $z \sim 15$, and find that many satellites with $M_* \sim 10^{6-8}$ solar masses quench within 2 Gyr of infall into the host halo, while others in the same mass range remain star-forming for as long as 5 Gyr. We show that this scatter can be explained by the satellite's gas mass and the ram pressure it feels at infall. Finally, we identify a characteristic stellar mass scale of $10^8$ solar masses above which infalling satellites are largely resistant to rapid environmental quenching. △ Less

Submitted 2 February, 2021; v1 submitted 6 August, 2020; originally announced August 2020.

Comments: 16 pages, 7 figures, accepted to ApJ

Inside Out and Upside-Down: The Roles of Gas Cooling and Dynamical Heating in Shaping the Stellar Age-Velocity Relation

Authors: Jonathan C. Bird, Sarah R. Loebman, David H. Weinberg, Alyson Brooks, Thomas R. Quinn, Charlotte R. Christensen

Abstract: Kinematic studies of disk galaxies, using individual stars in the Milky Way or statistical studies of global disk kinematics over time, provide insight into how disks form and evolve. We use a high-resolution, cosmological zoom-simulation of a Milky Way-mass disk galaxy h277 to tie together local disk kinematics and the evolution of the disk over time. The present-day stellar age-velocity relation… ▽ More Kinematic studies of disk galaxies, using individual stars in the Milky Way or statistical studies of global disk kinematics over time, provide insight into how disks form and evolve. We use a high-resolution, cosmological zoom-simulation of a Milky Way-mass disk galaxy h277 to tie together local disk kinematics and the evolution of the disk over time. The present-day stellar age-velocity relationship (AVR) of h277 is nearly identical to that of the analogous solar-neighborhood measurement in the Milky Way. A crucial element of this success is the simulation's dynamically cold multi-phase ISM, which allows young stars to form with a low velocity dispersion ($σ_{\mathrm{birth}}$$\sim 6 - 8\ \mathrm{km}\, \mathrm{s}^{-1}$) at late times. Older stars are born kinematically hotter (i.e., the disk settles over time in an "upside-down" formation scenario), and are subsequently heated after birth. The disk also grows "inside-out", and many of the older stars in the solar neighborhood at $z=0$ are present because of radial mixing. We demonstrate that the evolution of $σ_{\mathrm{birth}}$ in h277 can be explained by the same model used to describe the general decrease in velocity dispersion observed in disk galaxies from $z\sim 2-3$ to the present-day, in which the disk evolves in quasi-stable equilibrium and the ISM velocity dispersion decreases over time due to a decreasing gas fraction. Thus, our results tie together local observations of the Milky Way's AVR with observed kinematics of high $z$ disk galaxies. △ Less

Submitted 26 May, 2020; originally announced May 2020.

Comments: Submitting to MNRAS. Comments welcomed. This simulation is publicly available. Contact Alyson Brooks for details

Understanding the circumgalactic medium is critical for understanding galaxy evolution

Authors: Molly S. Peeples, Peter Behroozi, Rongmon Bordoloi, Alyson Brooks, James S. Bullock, Joseph N. Burchett, Hsiao-Wen Chen, John Chisholm, Charlotte Christensen, Alison Coil, Lauren Corlies, Aleksandar Diamond-Stanic, Megan Donahue, Claude-André Faucher-Giguère, Henry Ferguson, Drummond Fielding, Andrew J. Fox, David M. French, Steven R. Furlanetto, Mario Gennaro, Karoline M. Gilbert, Erika Hamden, Nimish Hathi, Matthew Hayes, Alaina Henry , et al. (47 additional authors not shown)

Abstract: Galaxies evolve under the influence of gas flows between their interstellar medium and their surrounding gaseous halos known as the circumgalactic medium (CGM). The CGM is a major reservoir of galactic baryons and metals, and plays a key role in the long cycles of accretion, feedback, and recycling of gas that drive star formation. In order to fully understand the physical processes at work within… ▽ More Galaxies evolve under the influence of gas flows between their interstellar medium and their surrounding gaseous halos known as the circumgalactic medium (CGM). The CGM is a major reservoir of galactic baryons and metals, and plays a key role in the long cycles of accretion, feedback, and recycling of gas that drive star formation. In order to fully understand the physical processes at work within galaxies, it is therefore essential to have a firm understanding of the composition, structure, kinematics, thermodynamics, and evolution of the CGM. In this white paper we outline connections between the CGM and galactic star formation histories, internal kinematics, chemical evolution, quenching, satellite evolution, dark matter halo occupation, and the reionization of the larger-scale intergalactic medium in light of the advances that will be made on these topics in the 2020s. We argue that, in the next decade, fundamental progress on all of these major issues depends critically on improved empirical characterization and theoretical understanding of the CGM. In particular, we discuss how future advances in spatially-resolved CGM observations at high spectral resolution, broader characterization of the CGM across galaxy mass and redshift, and expected breakthroughs in cosmological hydrodynamic simulations will help resolve these major problems in galaxy evolution. △ Less

Submitted 13 March, 2019; originally announced March 2019.

Comments: Astro2020 Decadal Science White Paper

Journal ref: Astro2020: Decadal Survey on Astronomy and Astrophysics, science white papers, no. 368; Bulletin of the American Astronomical Society, Vol. 51, Issue 3, id. 368 (2019)

Implementing Dust Shielding as a Criteria for Star Formation

Authors: Lindsey Byrne, Charlotte Christensen, Marios Tsekitsidis, Alyson Brooks, Tom Quinn

Abstract: Star formation is observed to be strongly correlated to dense regions of molecular gas. Although the exact nature of the link between star formation and molecular hydrogen is still unclear, some have suggested that shielding of dense gas by dust grains is the key factor enabling the presence of both. We present a sub-grid model for use in galaxy formation simulations in which star formation is lin… ▽ More Star formation is observed to be strongly correlated to dense regions of molecular gas. Although the exact nature of the link between star formation and molecular hydrogen is still unclear, some have suggested that shielding of dense gas by dust grains is the key factor enabling the presence of both. We present a sub-grid model for use in galaxy formation simulations in which star formation is linked explicitly to local dust shielding. We developed and tested our shielding and star formation models using smoothed particle hydrodynamic simulations of solar and sub-solar metallicity isolated Milky Way-mass disk galaxies. We compared our dust shielding-based star formation model to two other star formation recipes that used gas temperature and H$_2$ fraction as star formation criteria. We further followed the evolution of a dwarf galaxy within a cosmological context using both the shielding and H$_2$-based star formation models. We find that the shielding-based model allows for star formation at higher temperatures and lower densities than a model in which star formation is tied directly to H$_2$ abundance, as requiring H$_2$ formation leads the gas to undergo additional gravitational collapse before star formation. However, the resulting galaxies are very similar for both the shielding and H$_2$-based star formation models, and both models reproduce the resolved Kennicutt-Schmidt law. Therefore, both star formation models appear viable in the context of galaxy formation simulations. △ Less

Submitted 3 January, 2019; originally announced January 2019.

Comments: 13 pages, 12 figures

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

Authors: Elaad Applebaum, Alyson M. Brooks, Thomas R. Quinn, Charlotte R. Christensen

Abstract: 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-fain… ▽ More 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. △ Less

Submitted 20 February, 2020; v1 submitted 31 October, 2018; originally announced November 2018.

Comments: 14 pages, 11 figures, Updated to reflect accepted version, MNRAS

Journal ref: Monthly Notices of the Royal Astronomical Society, 2020, 492, 8-21

Dancing in the Dark: Uncertainty in ultra-faint dwarf galaxy predictions from cosmological simulations

Authors: Ferah Munshi, Alyson M. Brooks, Charlotte Christensen, Elaad Applebaum, Kelly Holley-Bockelmann, Thomas R. Quinn, James Wadsley

Abstract: The existence of ultra-faint dwarf (UFD) galaxies highlights the need to push our theoretical understanding of galaxies to extremely low mass. We examine the formation of UFDs by twice running a fully cosmological simulations of dwarf galaxies, but varying star formation. One run uses a temperature-density threshold for star formation, while the other uses an H$_{2}$-based sub-grid star formation… ▽ More The existence of ultra-faint dwarf (UFD) galaxies highlights the need to push our theoretical understanding of galaxies to extremely low mass. We examine the formation of UFDs by twice running a fully cosmological simulations of dwarf galaxies, but varying star formation. One run uses a temperature-density threshold for star formation, while the other uses an H$_{2}$-based sub-grid star formation model. The total number of dwarf galaxies that forms is different by a factor of 2 between the two runs, but most of these are satellites, leading to a factor of 5 difference in the number of luminous UFD companions around more massive, isolated dwarfs. The first run yields a 47\% chance of finding a satellite around a M$_{halo}$ $\sim 10^{10}$ M$_{\odot}$ host, while the H$_2$ run predicts only a 16\% chance. Metallicity is the primary physical parameter that creates this difference. As metallicity decreases, the formation of H$_2$ is slowed and relegated to higher-density material. Thus, our H$_2$ run is unable to form many (and often, any) stars before reionization removes gas. These results emphasize that predictions for UFD properties made using hydrodynamic simulations, in particular regarding the frequency of satellites around dwarf galaxies, the slope of the stellar mass function at low masses, as well as the properties of ultra-faint galaxies occupying the smallest halos, are extremely sensitive to the subgrid physics of star formation contained within the simulation. However, upcoming discoveries of ultra-faint dwarfs will provide invaluable constraining power on the physics of the first star formation. △ Less

Submitted 8 February, 2019; v1 submitted 29 October, 2018; originally announced October 2018.

Comments: 18 pages, 7 Figures. Accepted version, ApJ

Not So Heavy Metals: Black Hole Feedback Enriches The Circumgalactic Medium

Authors: N. Nicole Sanchez, Jessica K. Werk, Michael Tremmel, Andrew Pontzen, Charlotte Christensen, Thomas Quinn, Akaxia Cruz

Abstract: We examine the effects of SMBH feedback on the CGM using a cosmological hydrodynamic simulation \citep[{\sc Romulus25};][]{Tremmel2017} and a set of four zoom-in `genetically modified' Milky Way-mass galaxies sampling different evolutionary paths. By tracing the distribution of metals in the circumgalactic medium (CGM), we show that \ion{O}{6} is a sensitive indicator of supermassive black hole (S… ▽ More We examine the effects of SMBH feedback on the CGM using a cosmological hydrodynamic simulation \citep[{\sc Romulus25};][]{Tremmel2017} and a set of four zoom-in `genetically modified' Milky Way-mass galaxies sampling different evolutionary paths. By tracing the distribution of metals in the circumgalactic medium (CGM), we show that \ion{O}{6} is a sensitive indicator of supermassive black hole (SMBH) feedback. First, we calculate the column densities of \ion{O}{6} in simulated Milky Way-mass galaxies and compare them with observations from the COS-Halos Survey. Our simulations show column densities of \ion{O}{6} in the CGM consistent with those of COS-Halos star forming and quenched galaxies. These results contrast with those from previous simulation studies which typically underproduce CGM column densities of \ion{O}{6}. We determine that a galaxy's star formation history and assembly record have little effect on the amount of \ion{O}{6} in its CGM. Instead, column densities of \ion{O}{6} are closely tied to galaxy halo mass and BH growth history. The set of zoom-in, genetically modified Milky Way-mass galaxies indicates that the SMBH drives highly metal-enriched material out into its host galaxy's halo which in turn elevates the column densities of \ion{O}{6} in the CGM. △ Less

Submitted 25 May, 2021; v1 submitted 29 October, 2018; originally announced October 2018.

Comments: 15 pages, 12 figures, Accepted to ApJ; Note: Figures have transparency that may not render in browser

Tracing Outflowing Metals in Simulations of Dwarf and Spiral Galaxies

Authors: Charlotte R. Christensen, Romeel Dave, Alyson Brooks, Thomas Quinn, Sijing Shen

Abstract: We analyze the metal accumulation in dwarf and spiral galaxies by following the history of metal enrichment and outflows in a suite of twenty high-resolution simulated galaxies. These simulations agree with the observed stellar and gas-phase mass-metallicity relation, an agreement that relies on large fractions of the produced metals escaping into the CGM. For instance, in galaxies with Mvir ~ 1e9… ▽ More We analyze the metal accumulation in dwarf and spiral galaxies by following the history of metal enrichment and outflows in a suite of twenty high-resolution simulated galaxies. These simulations agree with the observed stellar and gas-phase mass-metallicity relation, an agreement that relies on large fractions of the produced metals escaping into the CGM. For instance, in galaxies with Mvir ~ 1e9.5 -- 1e10 solar masses, we find that about ~ 85% of the available metals are outside of the galactic disk at z = 0, although the fraction decreases to a little less than half in Milky Way-mass galaxies. In many cases, these metals are spread far beyond the virial radius. We analyze the metal deficit within the ISM and stars in the context of previous work tracking the inflow and outflow of baryons. Outflows are prevalent across the entire mass range, as is reaccretion. We find that between 40 and 80% of all metals removed from the galactic disk are later reaccreted. The outflows themselves are metal enriched relative to the ISM by a factor of 0.2 dex because of the correspondence between sites of metal enrichment and outflows. As a result, the metal mass loading factor scales as eta_metals \propto v_circ^-0.91, a somewhat shallower scaling than the total mass loading factor. We analyze the simulated galaxies within the context of analytic chemical evolution models by determining their net metal expulsion efficiencies, which encapsulate the rates of metal loss and reaccretion. We discuss these results in light of the inflow and outflow properties necessary for reproducing the mass-metallicity relation. △ Less

Submitted 23 August, 2018; originally announced August 2018.

Comments: Under review by ApJ. 21 pages, 15 figures

Multimessenger Signatures of Massive Black Holes in Dwarf Galaxies

Authors: Jillian Bellovary, Colleen Cleary, Ferah Munshi, Michael Tremmel, Charlotte Christensen, Alyson Brooks, Thomas Quinn

Abstract: Recent discoveries of massive black holes (MBHs) in dwarf galaxies suggest that they may have a more common presence than once thought. Systematic searches are revealing more candidates, but this process could be accelerated by predictions from simulations. We perform a study of several high-resolution, cosmological, zoom-in simulations focusing on dwarf galaxies that host massive black holes at z… ▽ More Recent discoveries of massive black holes (MBHs) in dwarf galaxies suggest that they may have a more common presence than once thought. Systematic searches are revealing more candidates, but this process could be accelerated by predictions from simulations. We perform a study of several high-resolution, cosmological, zoom-in simulations focusing on dwarf galaxies that host massive black holes at z = 0, with the aim of determining when the black holes are most observable. Larger dwarf galaxies are more likely to host MBHs than those of lower mass. About 50% of the MBHs in dwarfs are not centrally located, but rather are wandering within a few kpc of the galaxy center. The accretion luminosities of MBHs in dwarfs are low throughout cosmic time, rendering them extremely difficult to detect. However, the merger history of these MBHs is optimal for gravitational wave detection by LISA. △ Less

Submitted 19 October, 2018; v1 submitted 1 June, 2018; originally announced June 2018.

Comments: replaced with accepted version

Sowing black hole seeds: Direct collapse black hole formation with realistic Lyman-Werner radiation in cosmological simulations

Authors: Glenna Dunn, Jillian Bellovary, Kelly Holley-Bockelmann, Charlotte Christensen, Thomas Quinn

Abstract: We study the birth of supermassive black holes from the direct collapse process and characterize the sites where these black hole seeds form. In the pre-reionization epoch, molecular hydrogen (H$_2$) is an efficient coolant, causing gas to fragment and form Population III stars, but Lyman-Werner radiation can suppress H$_2$ formation and allow gas to collapse directly into a massive black hole. Th… ▽ More We study the birth of supermassive black holes from the direct collapse process and characterize the sites where these black hole seeds form. In the pre-reionization epoch, molecular hydrogen (H$_2$) is an efficient coolant, causing gas to fragment and form Population III stars, but Lyman-Werner radiation can suppress H$_2$ formation and allow gas to collapse directly into a massive black hole. The critical flux required to inhibit H$_2$ formation, $J_{\rm crit}$, is hotly debated, largely due to the uncertainties in the source radiation spectrum, H$_2$ self-shielding, and collisional dissociation rates. Here, we test the power of the direct collapse model in a self-consistent, time-dependant, non-uniform Lyman-Werner radiation field -- the first time such has been done in a cosmological volume -- using an updated version of the SPH+N-body tree code Gasoline with H$_2$ non-equilibrium abundance tracking, H$_2$ cooling, and a modern SPH implementation. We vary $J_{\rm crit} $ from $30$ to $10^3$ in units of $J_{21}$ to study how this parameter impacts the number of seed black holes and the type of galaxies which host them. We focus on black hole formation as a function of environment, halo mass, metallicity, and proximity of the Lyman-Werner source. Massive black hole seeds form more abundantly with lower $J_{\rm crit}$ thresholds, but regardless of $J_{\rm crit}$, these seeds typically form in halos that have recently begun star formation. Our results do not confirm the proposed atomic cooling halo pair scenario; rather black hole seeds predominantly form in low-metallicity pockets of halos which already host star formation. △ Less

Submitted 11 June, 2018; v1 submitted 2 March, 2018; originally announced March 2018.

Comments: 10 pages, 8 figures, accepted by ApJ

Reignition of Star Formation in Dwarf Galaxies

Authors: Anna C. Wright, Alyson M. Brooks, Daniel R. Weisz, Charlotte R. Christensen

Abstract: The Local Group hosts a number of star-forming dwarf galaxies that show evidence of periods of little to no star formation. We use a suite of cosmological simulations to study how star formation is reignited in such galaxies. We focus on isolated galaxies at $z=0$ with halo masses between 9.2$\times$10$^8$ M$_\odot$ and 8.4$\times$10$^9$ M$_\odot$, where star formation is typically shut off by rei… ▽ More The Local Group hosts a number of star-forming dwarf galaxies that show evidence of periods of little to no star formation. We use a suite of cosmological simulations to study how star formation is reignited in such galaxies. We focus on isolated galaxies at $z=0$ with halo masses between 9.2$\times$10$^8$ M$_\odot$ and 8.4$\times$10$^9$ M$_\odot$, where star formation is typically shut off by reionization or by supernova feedback. Nearly 20% of these simulated galaxies later restart star formation, due to interactions with streams of gas in the intergalactic medium, indicating that this mechanism is relatively common in this mass range and that many isolated dwarfs at $z=0$ may not have been isolated throughout their histories. The source of this gas is not necessarily cosmic filaments. Rather, the dwarfs interact with gas thrown off by nearby galaxy mergers or streams extending from neighboring galaxies. While high ram pressure interactions of this nature lead to stripping, the encounters that reignite star formation are low density and/or low velocity and thus low ram pressure, resulting in compression of the hot gas in the halos of our dwarfs. The gas mass bound up in hot halos can be substantial -- at least an order of magnitude greater than the mass contained in HI. Consequently, we find that dwarfs that have experienced reignition tend to be more HI-rich and have a higher M$_{HI}$/M$_{*}$ ratio at $z=0$ than galaxies with continuous star formation. Using this fact, we identify galaxies in the Local Volume that might have "gappy" star formation histories, and can be studied by the Hubble Space Telescope or the James Webb Space Telescope. △ Less

Submitted 11 October, 2018; v1 submitted 8 February, 2018; originally announced February 2018.

Comments: 16 pages, 11 figures, 1 table, 2 movies (online), accepted to MNRAS

A unified model for age-velocity dispersion relations in Local Group galaxies: Disentangling ISM turbulence and latent dynamical heating

Authors: Ryan Leaman, J. Trevor Mendel, Emily Wisnioski, Alyson M. Brooks, Michael A. Beasley, Else Starkenburg, Marie Martig, Giuseppina Battaglia, Charlotte Christensen, Andrew A. Cole, T. J. L. de Boer, Drew Wills

Abstract: We analyze age-velocity dispersion relations (AVRs) from kinematics of individual stars in eight Local Group galaxies ranging in mass from Carina ($M_{*} \sim 10^{6}$) to M31 ($M_{*} \sim 10^{11}$). Observationally the $σ$ vs. stellar age trends can be interpreted as dynamical heating of the stars by GMCs, bars/spiral arms, or merging subhalos; alternatively the stars could have simply been born o… ▽ More We analyze age-velocity dispersion relations (AVRs) from kinematics of individual stars in eight Local Group galaxies ranging in mass from Carina ($M_{*} \sim 10^{6}$) to M31 ($M_{*} \sim 10^{11}$). Observationally the $σ$ vs. stellar age trends can be interpreted as dynamical heating of the stars by GMCs, bars/spiral arms, or merging subhalos; alternatively the stars could have simply been born out of a more turbulent ISM at high redshift and retain that larger velocity dispersion till present day - consistent with recent IFU studies. To ascertain the dominant mechanism and better understand the impact of instabilities and feedback, we develop models based on observed SFHs of these Local Group galaxies in order to create an evolutionary formalism which describes the ISM velocity dispersion due to a galaxy's evolving gas fraction. These empirical models relax the common assumption that the stars are born from gas which has constant velocity dispersion at all redshifts. Using only the observed SFHs as input, the ISM velocity dispersion and a mid-plane scattering model fits the observed AVRs of low mass galaxies without fine tuning. Higher mass galaxies above $M_{vir} > 10^{11}$ need a larger contribution from latent dynamical heating processes (for example minor mergers), in excess of the ISM model. Using the SFHs we also find that supernovae feedback does not appear to be a dominant driver of the gas velocity dispersion compared to gravitational instabilities - at least for dispersions $σ\gtrsim 25$ km/s. Together our results point to stars being born with a velocity dispersion close to that of the gas at the time of their formation, with latent dynamical heating operating with a galaxy mass-dependent efficiency. These semi-empirical relations may help constrain the efficiency of feedback and its impact on the physics of disk settling in galaxy formation simulations. △ Less

Submitted 3 October, 2017; originally announced October 2017.

Comments: 21 pages, 10 figures. accepted for publication in MNRAS

Beta Dips in the Gaia Era: Simulation Predictions of the Galactic Velocity Anisotropy Parameter for Stellar Halos

Authors: Sarah R. Loebman, Monica Valluri, Kohei Hattori, Victor P. Debattista, Eric F. Bell, Greg Stinson, Charlotte R. Christensen, Alyson Brooks, Thomas R. Quinn, Fabio Governato

Abstract: The velocity anisotropy parameter, beta, is a measure of the kinematic state of orbits in the stellar halo which holds promise for constraining the merger history of the Milky Way (MW). We determine global trends for beta as a function of radius from three suites of simulations, including accretion only and cosmological hydrodynamic simulations. We find that both types of simulations are consisten… ▽ More The velocity anisotropy parameter, beta, is a measure of the kinematic state of orbits in the stellar halo which holds promise for constraining the merger history of the Milky Way (MW). We determine global trends for beta as a function of radius from three suites of simulations, including accretion only and cosmological hydrodynamic simulations. We find that both types of simulations are consistent and predict strong radial anisotropy (<beta>~0.7) for Galactocentric radii greater than 10 kpc. Previous observations of beta for the MW's stellar halo claim a detection of an isotropic or tangential "dip" at r~20 kpc. Using the N-body+SPH simulations, we investigate the temporal persistence, population origin, and severity of "dips" in beta. We find dips in the in situ stellar halo are long-lived, while dips in the accreted stellar halo are short-lived and tied to the recent accretion of satellite material. We also find that a major merger as early as z~1 can result in a present day low (isotropic to tangential) value of beta over a wide range of radii and angular expanse. While all of these mechanisms are plausible drivers for the beta dip observed in the MW, in the simulations, each mechanism has a unique metallicity signature associated with it, implying that future spectroscopic surveys could distinguish between them. Since an accurate knowledge of beta(r) is required for measuring the mass of the MW halo, we note significant transient dips in beta could cause an overestimate of the halo's mass when using spherical Jeans equation modeling. △ Less

Submitted 3 February, 2018; v1 submitted 20 April, 2017; originally announced April 2017.

Comments: 14 pages, 7 figures, accepted to ApJ, companion paper to Hattori et al. (2017), appendix with Gaia specific predictions

How to Reconcile the Observed Velocity Function of Galaxies with Theory

Authors: Alyson M. Brooks, Emmanouil Papastergis, Charlotte R. Christensen, Fabio Governato, Adrienne Stilp, Thomas R. Quinn, James Wadsley

Abstract: Within a Lambda Cold Dark Matter (LCDM) scenario, we use high resolution cosmological simulations spanning over four orders of magnitude in galaxy mass to understand the deficit of dwarf galaxies in observed velocity functions. We measure velocities in as similar a way as possible to observations, including generating mock HI data cubes for our simulated galaxies. We demonstrate that this apples-t… ▽ More Within a Lambda Cold Dark Matter (LCDM) scenario, we use high resolution cosmological simulations spanning over four orders of magnitude in galaxy mass to understand the deficit of dwarf galaxies in observed velocity functions. We measure velocities in as similar a way as possible to observations, including generating mock HI data cubes for our simulated galaxies. We demonstrate that this apples-to-apples comparison yields an "observed" velocity function in agreement with observations, reconciling the large number of low-mass halos expected in a LCDM cosmological model with the low number of observed dwarfs at a given velocity. We then explore the source of the discrepancy between observations and theory, and conclude that the dearth of observed dwarf galaxies is primarily explained by two effects. The first effect is that galactic rotational velocities derived from the HI linewidth severely underestimate the maximum halo velocity. The second effect is that a large fraction of halos at the lowest masses are too faint to be detected by current galaxy surveys. We find that cored dark matter density profiles can contribute to the lower observed velocity of galaxies, but only for galaxies in which the velocity is measured interior to the size of the core (~3 kpc). △ Less

Submitted 25 October, 2017; v1 submitted 26 January, 2017; originally announced January 2017.

Comments: updated after acceptance for publication in ApJ

The Role of Baryons in Creating Statistically Significant Planes of Satellites around Milky Way-Mass Galaxies

Authors: Sheehan H. Ahmed, Alyson M. Brooks, Charlotte R. Christensen

Abstract: We investigate whether the inclusion of baryonic physics influences the formation of thin, coherently rotating planes of satellites such as those seen around the Milky Way and Andromeda. For four Milky Way-mass simulations, each run both as dark matter-only and with baryons included, we are able to identify a planar configuration that significantly maximizes the number of plane satellite members.… ▽ More We investigate whether the inclusion of baryonic physics influences the formation of thin, coherently rotating planes of satellites such as those seen around the Milky Way and Andromeda. For four Milky Way-mass simulations, each run both as dark matter-only and with baryons included, we are able to identify a planar configuration that significantly maximizes the number of plane satellite members. The maximum plane member satellites are consistently different between the dark matter-only and baryonic versions of the same run due to the fact that satellites are both more likely to be destroyed and to infall later in the baryonic runs. Hence, studying satellite planes in dark matter-only simulations is misleading, because they will be composed of different satellite members than those that would exist if baryons were included. Additionally, the destruction of satellites in the baryonic runs leads to less radially concentrated satellite distributions, a result that is critical to making planes that are statistically significant compared to a random distribution. Since all planes pass through the centre of the galaxy, it is much harder to create a plane from a random distribution if the satellites have a low radial concentration. We identify Andromeda's low radial satellite concentration as a key reason why the plane in Andromeda is highly significant. Despite this, when co-rotation is considered, none of the satellite planes identified for the simulated galaxies are as statistically significant as the observed planes around the Milky Way and Andromeda, even in the baryonic runs. △ Less

Submitted 10 October, 2016; originally announced October 2016.

Comments: 16 pages, 9 figures, submitted to MNRAS

The effects of host galaxy properties on merging compact binaries detectable by LIGO

Authors: Richard O'Shaughnessy, Jillian Bellovary, Alyson Brooks, Sijing Shen, Fabio Governato, Charlotte Christensen

Abstract: Cosmological simulations of galaxy formation can produce present-day galaxies with a large range of assembly and star formation histories. A detailed study of the metallicity evolution and star formation history of such simulations can assist in predicting LIGO-detectable compact object binary mergers. Recent simulations of compact binary evolution suggest the compact object merger rate depends se… ▽ More Cosmological simulations of galaxy formation can produce present-day galaxies with a large range of assembly and star formation histories. A detailed study of the metallicity evolution and star formation history of such simulations can assist in predicting LIGO-detectable compact object binary mergers. Recent simulations of compact binary evolution suggest the compact object merger rate depends sensitively on the progenitor's metallicity. Rare low-metallicity star formation during galaxy assembly can produce more detected compact binaries than typical star formation. Using detailed simulations of galaxy and chemical evolution, we determine how sensitively the compact binary populations of galaxies with similar present-day appearance depend on the details of their assembly. We also demonstrate by concrete example the extent to which dwarf galaxies overabundantly produce compact binary mergers, particularly binary black holes, relative to more massive galaxies. We discuss the implications for transient multimessenger astronomy with compact binary sources. △ Less

Submitted 21 September, 2016; originally announced September 2016.

Comments: 9 pages, re-submitted to MNRAS after responding to referee comments

Bulge Formation via Mergers in Cosmological Simulations

Authors: Alyson M. Brooks, Charlotte R. Christensen

Abstract: 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 pr… ▽ More 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. △ Less

Submitted 12 November, 2015; originally announced November 2015.

Comments: The first 5 sections are review, but Section 6 introduces new analysis. We demonstrate that mergers prevent bulge formation in low mass galaxies. Approximately 15 pages in MNRAS format

Journal ref: Review article written for "Galactic Bulges," published by Springer, 2015, edited by E. Laurikainen, R.F. Peletier, D.A. Gadotti

In-N-Out: the gas cycle from dwarfs to spiral galaxies

Authors: Charlotte R. Christensen, Romeel Davé, Fabio Governato, Andrew Pontzen, Alyson Brooks, Ferah Munshi, Thomas Quinn, James Wadsley

Abstract: We examine the scalings of galactic outflows with halo mass across a suite of twenty high-resolution cosmological zoom galaxy simulations covering halo masses from 10^9.5 - 10^12 M_sun. These simulations self-consistently generate outflows from the available supernova energy in a manner that successfully reproduces key galaxy observables including the stellar mass-halo mass, Tully-Fisher, and mass… ▽ More We examine the scalings of galactic outflows with halo mass across a suite of twenty high-resolution cosmological zoom galaxy simulations covering halo masses from 10^9.5 - 10^12 M_sun. These simulations self-consistently generate outflows from the available supernova energy in a manner that successfully reproduces key galaxy observables including the stellar mass-halo mass, Tully-Fisher, and mass-metallicity relations. We quantify the importance of ejective feedback to setting the stellar mass relative to the efficiency of gas accretion and star formation. Ejective feedback is increasingly important as galaxy mass decreases; we find an effective mass loading factor that scales as v_circ^(-2.2), with an amplitude and shape that is invariant with redshift. These scalings are consistent with analytic models for energy-driven wind, based solely on the halo potential. Recycling is common: about half the outflow mass across all galaxy masses is later re-accreted. The recycling timescale is typically about 1 Gyr, virtually independent of halo mass. Recycled material is re-accreted farther out in the disk and with typically about 2-3 times more angular momentum. These results elucidate and quantify how the baryon cycle plausibly regulates star formation and alters the angular momentum distribution of disk material across the halo mass range where most of cosmic star formation occurs. △ Less

Submitted 30 May, 2018; v1 submitted 31 July, 2015; originally announced August 2015.

Comments: 21 pages, 17 figures, accepted to ApJ

Consequences of bursty star formation on galaxy observables at high redshifts

Authors: Alberto Domínguez, Brian Siana, Alyson M. Brooks, Charlotte R. Christensen, Gustavo Bruzual, Daniel P. Stark, Anahita Alavi

Abstract: 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… ▽ More 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. △ Less

Submitted 15 May, 2015; v1 submitted 25 August, 2014; originally announced August 2014.

Comments: 12 pages, 10 figures, and 2 tables. Version accepted by MNRAS

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

Authors: Sarah R. Loebman, Zeljko Ivezic, Thomas R. Quinn, Jo Bovy, Charlotte R. Christensen, Mario Juric, Rok Roskar, Alyson M. Brooks, Fabio Governato

Abstract: 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 genera… ▽ More 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. △ Less

Submitted 24 November, 2014; v1 submitted 22 August, 2014; originally announced August 2014.

Comments: 15 pages, 21 figures, accepted to ApJ. If viewing using Preview, some figures may appear blurry. If possible, please view using Adobe PDF Viewer

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

Authors: Fabio Governato, Daniel Weisz, Andrew Pontzen, Sarah Loebman, Darren Reed, Alyson M. Brooks, Peter Behroozi, Charlotte Christensen, Piero Madau, Lucio Mayer, Sijing Shen, Matthew Walker, Thomas Quinn, Benjamin W. Keller, James Wadsley

Abstract: 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. In… ▽ More 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. △ Less

Submitted 23 November, 2014; v1 submitted 30 June, 2014; originally announced July 2014.

Comments: 13 pages including Appendix on Color Magnitude Diagrams. Accepted by MNRAS. Added one plot and details on ChaNGa implementation. Reduced number of citations after editorial request

Journal ref: MNRAS, 2015, 448, 792

Effects of Inclination on Measuring Velocity Dispersion and Implications for Black Holes

Authors: Jillian Bellovary, Kelly Holley-Bockelmann, Kayhan Gültekin, Charlotte Christensen, Fabio Governato, Alyson Brooks, Sarah Loebman, Ferah Munshi

Abstract: The relation of central black hole mass and stellar spheroid velocity dispersion (the M-$σ$ relation) is one of the best-known and tightest correlations linking black holes and their host galaxies. There has been much scrutiny concerning the difficulty of obtaining accurate black hole measurements, and rightly so; however, it has been taken for granted that measurements of velocity dispersion are… ▽ More The relation of central black hole mass and stellar spheroid velocity dispersion (the M-$σ$ relation) is one of the best-known and tightest correlations linking black holes and their host galaxies. There has been much scrutiny concerning the difficulty of obtaining accurate black hole measurements, and rightly so; however, it has been taken for granted that measurements of velocity dispersion are essentially straightforward. We examine five disk galaxies from cosmological SPH simulations and find that line-of-sight effects due to galaxy orientation can affect the measured $σ$ by 30%, and consequently black hole mass predictions by up to 1.0 dex. Face-on orientations correspond to systematically lower velocity dispersion measurements, while more edge-on orientations give higher velocity dispersions, due to contamination by disk stars when measuring line of sight quantities. We caution observers that the uncertainty of velocity dispersion measurements is at least 20 km/s, and can be much larger for moderate inclinations. This effect may account for some of the scatter in the locally measured M-$σ$ relation, particularly at the low-mass end. We provide a method for correcting observed $σ_{\rm los}$ values for inclination effects based on observable quantities. △ Less

Submitted 22 September, 2014; v1 submitted 1 May, 2014; originally announced May 2014.

Comments: 11 pages, 7 figures, replaced with accepted version

Simulating disc galaxy bulges that are consistent with observed scaling relations

Authors: C. R. Christensen, A. M. Brooks, D. B. Fisher, F. Governato, J. McCleary, T. R. Quinn, S. Shen, J. Wadsley

Abstract: We present a detailed comparison between the photometric properties of the bulges of two simulated galaxies and those of a uniform sample of observed galaxies. This analysis shows that the simulated galaxies have bulges with realistic surface brightnesses for their sizes and magnitude. These two field disc galaxies have rotational velocities ~ 100 km/s and were integrated to a redshift of zero in… ▽ More We present a detailed comparison between the photometric properties of the bulges of two simulated galaxies and those of a uniform sample of observed galaxies. This analysis shows that the simulated galaxies have bulges with realistic surface brightnesses for their sizes and magnitude. These two field disc galaxies have rotational velocities ~ 100 km/s and were integrated to a redshift of zero in a fully cosmological Lambda cold dark matter context as part of high-resolution smoothed particle hydrodynamic simulations. We performed bulge-disc decompositions of the galaxies using artificial observations, in order to conduct a fair comparison to observations. We also dynamically decomposed the galaxies and compared the star formation histories of the bulges to those of the entire galaxies. These star formation histories showed that the bulges were primarily formed before z = 1 and during periods of rapid star formation. Both galaxies have large amounts of early star formation, which is likely related to the relatively high bulge-to-disc ratios also measured for them. Unlike almost all previous cosmological simulations, the realistically concentrated bulges of these galaxies do not lead to unphysically high rotational velocities, causing them to naturally lie along the observed Tully-Fisher relation. △ Less

Submitted 11 February, 2014; originally announced February 2014.

Comments: 6 pages, 3 figures. Accepted for publication in MNRAS Letters

The Broadband XMM-Newton and NuSTAR X-ray Spectra of Two Ultraluminous X-ray Sources in the Galaxy IC 342

Authors: Vikram Rana, Fiona A. Harrison, Matteo Bachetti, Dominic J. Walton, Felix Furst, Didier Barret, Jon M. Miller, Andrew C. Fabian, Steven E. Boggs, Finn C. Christensen, William W. Craig, Brian W. Grefenstette, Charles J. Hailey, Kristin K. Madsen, Andrew F. Ptak, Daniel Stern, Natalie A. Webb, William W. Zhang

Abstract: We present results for two Ultraluminous X-ray Sources (ULXs), IC 342 X-1 and IC 342 X-2, using two epochs of XMM-Newton and NuSTAR observations separated by $\sim$7 days. We observe little spectral or flux variability above 1 keV between epochs, with unabsorbed 0.3--30 keV luminosities being $1.04^{+0.08}_{-0.06} \times 10^{40}$ erg s$^{-1}$ for IC 342 X-1 and $7.40\pm0.20 \times 10^{39}$ erg s… ▽ More We present results for two Ultraluminous X-ray Sources (ULXs), IC 342 X-1 and IC 342 X-2, using two epochs of XMM-Newton and NuSTAR observations separated by $\sim$7 days. We observe little spectral or flux variability above 1 keV between epochs, with unabsorbed 0.3--30 keV luminosities being $1.04^{+0.08}_{-0.06} \times 10^{40}$ erg s$^{-1}$ for IC 342 X-1 and $7.40\pm0.20 \times 10^{39}$ erg s$^{-1}$ for IC 342 X-2, so that both were observed in a similar, luminous state. Both sources have a high absorbing column in excess of the Galactic value. Neither source has a spectrum consistent with a black hole binary in low/hard state, and both ULXs exhibit strong curvature in their broadband X-ray spectra. This curvature rules out models that invoke a simple reflection-dominated spectrum with a broadened iron line and no cutoff in the illuminating power-law continuum. X-ray spectrum of IC 342 X-1 can be characterized by a soft disk-like black body component at low energies and a cool, optically thick Comptonization continuum at high energies, but unique physical interpretation of the spectral components remains challenging. The broadband spectrum of IC 342 X-2 can be fit by either a hot (3.8 keV) accretion disk, or a Comptonized continuum with no indication of a seed photon population. Although the seed photon component may be masked by soft excess emission unlikely to be associated with the binary system, combined with the high absorption column, it is more plausible that the broadband X-ray emission arises from a simple thin blackbody disk component. Secure identification of the origin of the spectral components in these sources will likely require broadband spectral variability studies. △ Less

Submitted 19 November, 2014; v1 submitted 19 January, 2014; originally announced January 2014.

Comments: 12 pages, 11 figures, 5 Tables, Accepted for publication in The Astrophysical Journal

The Pressure of the Star Forming ISM in Cosmological Simulations

Authors: Ferah Munshi, Charlotte Christensen, Thomas R. Quinn, Fabio Governato, James Wadsley, Sarah Loebman, Sijing Shen

Abstract: We examine the pressure of the star-forming interstellar medium (ISM) of Milky-Way sized disk galaxies using fully cosmological SPH+N-body, high resolution simulations. These simulations include explicit treatment of metal-line cooling in addition to dust and self-shielding, $\mathrm{H_{2}}$ based star formation. The 4 simulated halos have masses ranging from a few times $10^{10}$ to nearly… ▽ More We examine the pressure of the star-forming interstellar medium (ISM) of Milky-Way sized disk galaxies using fully cosmological SPH+N-body, high resolution simulations. These simulations include explicit treatment of metal-line cooling in addition to dust and self-shielding, $\mathrm{H_{2}}$ based star formation. The 4 simulated halos have masses ranging from a few times $10^{10}$ to nearly $10^{12}$ solar masses. Using a kinematic decomposition of these galaxies into present-day bulge and disk components, we find that the typical pressure of the star-forming ISM in the present-day bulge is higher than that in the present-day disk by an order of magnitude. We also find that pressure of the star-forming ISM at high redshift is on average, higher than ISM pressures at low redshift. This explains the why the bulge forms at higher pressures: the disk assembles at lower redshift, when the ISM is lower pressure and the bulge forms at high redshift, when the ISM is at higher pressure. If ISM pressure and IMF variation are tied together as suggested in studies like \cite{Conroy2012}, these results could indicate a time-dependent IMF in Milky-Way like systems, as well as a different IMF in the bulge and the disk. △ Less

Submitted 5 December, 2013; originally announced December 2013.

Comments: Accepted 12/3/2013 to ApJL

Galaxy Formation with local photoionisation feedback I. Methods

Authors: R. Kannan, G. S. Stinson, A. V. Macciò, J. F. Hennawi, R. Woods, J. Wadsley, S. Shen, T. Robitaille, S. Cantalupo, T. R. Quinn, C. Christensen

Abstract: We present a first study of the effect of local photoionising radiation on gas cooling in smoothed particle hydrodynamics simulations of galaxy formation. We explore the combined effect of ionising radiation from young and old stellar populations. The method computes the effect of multiple radiative sources using the same tree algorithm used for gravity, so it is computationally efficient and well… ▽ More We present a first study of the effect of local photoionising radiation on gas cooling in smoothed particle hydrodynamics simulations of galaxy formation. We explore the combined effect of ionising radiation from young and old stellar populations. The method computes the effect of multiple radiative sources using the same tree algorithm used for gravity, so it is computationally efficient and well resolved. The method foregoes calculating absorption and scattering in favour of a constant escape fraction for young stars to keep the calculation efficient enough to simulate the entire evolution of a galaxy in a cosmological context to the present day. This allows us to quantify the effect of the local photoionisation feedback through the whole history of a galaxy`s formation. The simulation of a Milky Way like galaxy using the local photoionisation model forms ~ 40 % less stars than a simulation that only includes a standard uniform background UV field. The local photoionisation model decreases star formation by increasing the cooling time of the gas in the halo and increasing the equilibrium temperature of dense gas in the disc. Coupling the local radiation field to gas cooling from the halo provides a preventive feedback mechanism which keeps the central disc light and produces slowly rising rotation curves without resorting to extreme feedback mechanisms. These preliminary results indicate that the effect of local photoionising sources is significant and should not be ignored in models of galaxy formation. △ Less

Submitted 28 October, 2013; v1 submitted 24 October, 2013; originally announced October 2013.

Comments: Accepted for Publication in MNRAS, 13 pages, 13 figures

The AGORA High-Resolution Galaxy Simulations Comparison Project

Authors: Ji-hoon Kim, Tom Abel, Oscar Agertz, Greg L. Bryan, Daniel Ceverino, Charlotte Christensen, Charlie Conroy, Avishai Dekel, Nickolay Y. Gnedin, Nathan J. Goldbaum, Javiera Guedes, Oliver Hahn, Alexander Hobbs, Philip F. Hopkins, Cameron B. Hummels, Francesca Iannuzzi, Dusan Keres, Anatoly Klypin, Andrey V. Kravtsov, Mark R. Krumholz, Michael Kuhlen, Samuel N. Leitner, Piero Madau, Lucio Mayer, Christopher E. Moody , et al. (21 additional authors not shown)

Abstract: We introduce the AGORA project, a comprehensive numerical study of well-resolved galaxies within the LCDM cosmology. Cosmological hydrodynamic simulations with force resolutions of ~100 proper pc or better will be run with a variety of code platforms to follow the hierarchical growth, star formation history, morphological transformation, and the cycle of baryons in and out of 8 galaxies with halo… ▽ More We introduce the AGORA project, a comprehensive numerical study of well-resolved galaxies within the LCDM cosmology. Cosmological hydrodynamic simulations with force resolutions of ~100 proper pc or better will be run with a variety of code platforms to follow the hierarchical growth, star formation history, morphological transformation, and the cycle of baryons in and out of 8 galaxies with halo masses M_vir ~= 1e10, 1e11, 1e12, and 1e13 Msun at z=0 and two different ("violent" and "quiescent") assembly histories. The numerical techniques and implementations used in this project include the smoothed particle hydrodynamics codes GADGET and GASOLINE, and the adaptive mesh refinement codes ART, ENZO, and RAMSES. The codes will share common initial conditions and common astrophysics packages including UV background, metal-dependent radiative cooling, metal and energy yields of supernovae, and stellar initial mass function. These are described in detail in the present paper. Subgrid star formation and feedback prescriptions will be tuned to provide a realistic interstellar and circumgalactic medium using a non-cosmological disk galaxy simulation. Cosmological runs will be systematically compared with each other using a common analysis toolkit, and validated against observations to verify that the solutions are robust - i.e., that the astrophysical assumptions are responsible for any success, rather than artifacts of particular implementations. The goals of the AGORA project are, broadly speaking, to raise the realism and predictive power of galaxy simulations and the understanding of the feedback processes that regulate galaxy "metabolism." The proof-of-concept dark matter-only test of the formation of a galactic halo with a z=0 mass of M_vir ~= 1.7e11 Msun by 9 different versions of the participating codes is also presented to validate the infrastructure of the project. △ Less

Submitted 24 December, 2013; v1 submitted 12 August, 2013; originally announced August 2013.

Comments: 21 pages, 6 figures, Accepted for publication in the Astrophysical Journal Supplement, Image resolution greatly reduced, High-resolution version of this article and more information about the AGORA Project including the science goals of the 13 Working Groups are available at http://www.AGORAsimulations.org/ and at http://sites.google.com/site/santacruzcomparisonproject/details/

Journal ref: Astrophysical Journal Supplement Series, 210, 14, 2014

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

Authors: Charlotte Christensen, Fabio Governato, Thomas Quinn, Alyson M. Brooks, David B. Fisher, Sijing Shen, Jacqueline McCleary, James Wadsley

Abstract: 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 o… ▽ More 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) △ Less

Submitted 24 January, 2014; v1 submitted 1 November, 2012; originally announced November 2012.

Comments: 19 pages, 14 figures, submitted to MNRAS

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

Authors: Sarah R. Loebman, Zeljko Ivezic, Thomas R. Quinn, Fabio Governato, Alyson M. Brooks, Charlotte R. Christensen, Mario Juric

Abstract: 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… ▽ More 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. △ Less

Submitted 12 September, 2012; originally announced September 2012.

Comments: accepted to ApJ Letters, six pages, five figures

Reproducing the Stellar Mass/Halo Mass Relation in Simulated LCDM Galaxies: Theory vs Observational Estimates

Authors: Ferah Munshi, F. Governato, A. M. Brooks, C. Christensen, S. Shen, S. Loebman, B. Moster, T. Quinn, J. Wadsley

Abstract: We examine the present-day total stellar-to-halo mass (SHM) ratio as a function of halo mass for a new sample of simulated field galaxies using fully cosmological, LCDM, high resolution SPH + N-Body simulations.These simulations include an explicit treatment of metal line cooling, dust and self-shielding, H2 based star formation and supernova driven gas outflows. The 18 simulated halos have masses… ▽ More We examine the present-day total stellar-to-halo mass (SHM) ratio as a function of halo mass for a new sample of simulated field galaxies using fully cosmological, LCDM, high resolution SPH + N-Body simulations.These simulations include an explicit treatment of metal line cooling, dust and self-shielding, H2 based star formation and supernova driven gas outflows. The 18 simulated halos have masses ranging from a few times 10^8 to nearly 10^12 solar masses. At z=0 our simulated galaxies have a baryon content and morphology typical of field galaxies. Over a stellar mass range of 2.2 x 10^3 to 4.5 x 10^10 solar masses, we find extremely good agreement between the SHM ratio in simulations and the present-day predictions from the statistical Abundance Matching Technique presented in Moster et al. (2012). This improvement over past simulations is due to a number systematic factors, each decreasing the SHM ratios: 1) gas outflows that reduce the overall SF efficiency but allow for the formation of a cold gas component 2) estimating the stellar masses of simulated galaxies using artificial observations and photometric techniques similar to those used in observations and 3) accounting for a systematic, up to 30 percent overestimate in total halo masses in DM-only simulations, due to the neglect of baryon loss over cosmic times. Our analysis suggests that stellar mass estimates based on photometric magnitudes can underestimate the contribution of old stellar populations to the total stellar mass, leading to stellar mass errors of up to 50 percent for individual galaxies. These results highlight the importance of using proper techniques to compare simulations with observations and reduce the perceived tension between the star formation efficiency in galaxy formation models and in real galaxies. △ Less

Submitted 6 September, 2012; originally announced September 2012.

Comments: Submitted to ApJ 9 pages, 5 figures

High Velocity Dispersion in A Rare Grand Design Spiral Galaxy at Redshift z=2.18

Authors: David R. Law, Alice E. Shapley, Charles C. Steidel, Naveen A. Reddy, Charlotte R. Christensen, Dawn K. Erb

Abstract: Although relatively common in the local Universe, only one grand-design spiral galaxy has been spectroscopically confirmed to lie at z>2 (HDFX 28; z=2.011), and may prove to be a major merger that simply resembles a spiral in projection. The rarity of spirals has been explained as a result of disks being dynamically 'hot' at z>2 which may instead favor the formation of commonly-observed clumpy str… ▽ More Although relatively common in the local Universe, only one grand-design spiral galaxy has been spectroscopically confirmed to lie at z>2 (HDFX 28; z=2.011), and may prove to be a major merger that simply resembles a spiral in projection. The rarity of spirals has been explained as a result of disks being dynamically 'hot' at z>2 which may instead favor the formation of commonly-observed clumpy structures. Alternatively, current instrumentation may simply not be sensitive enough to detect spiral structures comparable to those in the modern Universe. At redshifts <2, the velocity dispersion of disks decreases, and spiral galaxies are more numerous by z~1. Here we report observations of the grand design spiral galaxy Q2343-BX442 at z=2.18. Spectroscopy of ionized gas shows that the disk is dynamically hot, implying an uncertain origin for the spiral structure. The kinematics of the galaxy are consistent with a thick disk undergoing a minor merger, which can drive the formation of short-lived spiral structure. A duty cycle of < 100 Myr for such tidally-induced spiral structure in a hot massive disk is consistent with their rarity. △ Less

Submitted 17 July, 2012; originally announced July 2012.

Comments: Accepted for publication in Nature (July 19 2012). Includes 15-page supplement

Baryons Matter: Why Luminous Satellite Galaxies Have Reduced Central Masses

Authors: Adi Zolotov, Alyson M. Brooks, Beth Willman, Fabio Governato, Andrew Pontzen, Charlotte Christensen, Avishai Dekel, Tom Quinn, Sijing Shen, James Wadsley

Abstract: 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… ▽ More 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. △ Less

Submitted 4 December, 2012; v1 submitted 29 June, 2012; originally announced July 2012.

Comments: Version accepted for publication in ApJ. Now includes an appendix with a discussion of numerical convergence

Implementing Molecular Hydrogen in Hydrodynamic Simulations of Galaxy Formation

Authors: Charlotte Christensen, Thomas Quinn, Fabio Governato, Adrienne Stilp, Sijing Shen, James Wadsley

Abstract: Motivated by the observed connection between molecular hydrogen (H2) and star formation, we present a method for tracking the non-equilibrium abundance and cooling processes of H2 and H2-based star formation in Smoothed Particle Hydrodynamic simulations. The local abundances of H2 are calculated by integrating over the hydrogen chemical network. This calculation includes the gas-phase and dust gra… ▽ More Motivated by the observed connection between molecular hydrogen (H2) and star formation, we present a method for tracking the non-equilibrium abundance and cooling processes of H2 and H2-based star formation in Smoothed Particle Hydrodynamic simulations. The local abundances of H2 are calculated by integrating over the hydrogen chemical network. This calculation includes the gas-phase and dust grain formation of H2, shielding of HI and H2, and photodissociation of H2 by Lyman-Werner radiation from nearby stellar populations. Because this model does not assume equilibrium abundances, it is particularly well suited for simulations that model low-metallicity environments, such as dwarf galaxies and the early Universe. We further introduce an explicit link between star formation and local H2 abundance. This link limits star formation to "star-forming regions," represented by areas with abundant H2. With this implementation, we determine the effect of H2 on star formation in a cosmological simulation of a dwarf galaxy. This simulation is the first cosmological simulation with non-equilibrium H2 abundances to be integrated to a redshift of zero or to include efficient SN feedback. We find that our simulations are consistent with the observed Tully-Fisher, global Kennicutt-Schmidt, and resolved Kennicutt-Schmidt relations. We find that the inclusion of shielding of both the atomic and molecular hydrogen and, to a lesser extent, the additional cooling from H2 at temperatures between 200 and 5000 K increases the amount of cold gas in the galaxies. The changes to the ISM result in an increased amount of cold, dense gas in the disk of the galaxy and the formation of a clumpier interstellar media (ISM). The explicit link between star formation and H2 and the clumpier ISM results in a bluer galaxy with a greater spatial distribution of star formation at a redshift of zero. (abridged) △ Less

Submitted 24 May, 2012; originally announced May 2012.

Comments: 20 pages, 15 figures, resubmitted to MNRAS following referee report

Morphology of High-Multiplicity Events in Heavy Ion Collisions

Authors: P. Naselsky, C. H. Christensen, P. R. Christensen, P. H. Damgaard, A. Frejsel, J. J. Gaardhøje, A. Hansen, M. Hansen, J. Kim, O. Verkhodanov, U. A. Wiedemann

Abstract: We discuss opportunities that may arise from subjecting high-multiplicity events in relativistic heavy ion collisions to an analysis similar to the one used in cosmology for the study of fluctuations of the Cosmic Microwave Background (CMB). To this end, we discuss examples of how pertinent features of heavy ion collisions including global characteristics, signatures of collective flow and event-w… ▽ More We discuss opportunities that may arise from subjecting high-multiplicity events in relativistic heavy ion collisions to an analysis similar to the one used in cosmology for the study of fluctuations of the Cosmic Microwave Background (CMB). To this end, we discuss examples of how pertinent features of heavy ion collisions including global characteristics, signatures of collective flow and event-wise fluctuations are visually represented in a Mollweide projection commonly used in CMB analysis, and how they are statistically analyzed in an expansion over spherical harmonic functions. If applied to the characterization of purely azimuthal dependent phenomena such as collective flow, the expansion coefficients of spherical harmonics are seen to contain redundancies compared to the set of harmonic flow coefficients commonly used in heavy ion collisions. Our exploratory study indicates, however, that these redundancies may offer novel opportunities for a detailed characterization of those event-wise fluctuations that remain after subtraction of the dominant collective flow signatures. By construction, the proposed approach allows also for the characterization of more complex collective phenomena like higher-order flow and other sources of fluctuations, and it may be extended to the characterization of phenomena of non-collective origin such as jets. △ Less

Submitted 7 August, 2012; v1 submitted 2 April, 2012; originally announced April 2012.

Comments: Matches version accepted for publication in Physical Review C. 13 pages, 9 figures

Journal ref: Phys. Rev. C 86, 024916 (2012)

Cuspy No More: How Outflows Affect the Central Dark Matter and Baryon Distribution in Lambda CDM Galaxies

Authors: F. Governato, A. Zolotov, A. Pontzen, C. Christensen, S. H. Oh, A. M. Brooks, T. Quinn, S. Shen, J. Wadsley

Abstract: We examine the evolution of the inner dark matter (DM) and baryonic density profile of a new sample of simulated field galaxies using fully cosmological, Lambda CDM, high resolution SPH + N-Body simulations. These simulations include explicit H2 and metal cooling, star formation (SF) and supernovae (SNe) driven gas outflows. Starting at high redshift, rapid, repeated gas outflows following bursty… ▽ More We examine the evolution of the inner dark matter (DM) and baryonic density profile of a new sample of simulated field galaxies using fully cosmological, Lambda CDM, high resolution SPH + N-Body simulations. These simulations include explicit H2 and metal cooling, star formation (SF) and supernovae (SNe) driven gas outflows. Starting at high redshift, rapid, repeated gas outflows following bursty SF transfer energy to the DM component and significantly flatten the originally `cuspy' central DM mass profile of galaxies with present day stellar masses in the 10^4.5 -- 10^9.8 Msolar range. At z=0, the central slope of the DM density profile of our galaxies (measured between 0.3 and 0.7 kpc from their centre) is well fitted by rhoDM propto r^alpha with alpha \simeq -0.5 + 0.35 log_10(Mstar/10^8Msolar) where Mstar is the stellar mass of the galaxy and 4 < log_10 Mstar < 9.4. These values imply DM profiles flatter than those obtained in DM--only simulations and in close agreement with those inferred in galaxies from the THINGS and LITTLE THINGS survey. Only in very small halos, where by z=0 star formation has converted less than ~ 0.03% of the original baryon abundance into stars, outflows do not flatten the original cuspy DM profile out to radii resolved by our simulations. The mass (DM and baryonic) measured within the inner 500 pc of each simulated galaxy remains nearly constant over four orders of magnitudes in stellar mass for Mstar 10^9 Msolar. This finding is consistent with estimates for faint Local Group dwarfs and field galaxies. These results address one of the outstanding problems faced by the CDM model, namely the strong discrepancy between the original predictions of cuspy DM profiles and the shallower central DM distribution observed in galaxies. △ Less

Submitted 7 February, 2012; v1 submitted 2 February, 2012; originally announced February 2012.

Comments: MNRAS in press. Accepted version, a few references added. 12 pages. Animation at http://youtu.be/FbcgEovabDI?hd=1

Star Formation and Feedback in Smoothed Particle Hydrodynamic Simulations II: Resolution Effects

Authors: Charlotte Christensen, Thomas Quinn, Gregory Stinson, Jillian Bellovary, James Wadsley

Abstract: We examine the effect of mass and force resolution on a specific star formation (SF) recipe using a set of N-body/Smooth Particle Hydrodynamic simulations of isolated galaxies. Our simulations span halo masses from 10^9 to 10^13 solar masses, more than four orders of magnitude in mass resolution, and two orders of magnitude in the gravitational softening length, epsilon, representing the force res… ▽ More We examine the effect of mass and force resolution on a specific star formation (SF) recipe using a set of N-body/Smooth Particle Hydrodynamic simulations of isolated galaxies. Our simulations span halo masses from 10^9 to 10^13 solar masses, more than four orders of magnitude in mass resolution, and two orders of magnitude in the gravitational softening length, epsilon, representing the force resolution. We examine the total global star formation rate, the star formation history, and the quantity of stellar feedback and compare the disk structure of the galaxies. Based on our analysis, we recommend using at least 10^4 particles each for the dark matter and gas component and a force resolution of epsilon approximately equal to 10^-3 R_vir when studying global SF and feedback. When the spatial distribution of stars is important, the number of gas and dark matter particles must be increased to at least 10^5 of each. Low mass resolution simulations with fixed softening lengths show particularly weak stellar disks due to two-body heating. While decreasing spatial resolution in low mass resolution simulations limits two-body effects, density and potential gradients cannot be sustained. Regardless of the softening, low-mass resolution simulations contain fewer high density regions where SF may occur. Galaxies of approximately 10^10 solar masses display unique sensitivity to both mass and force resolution. This mass of galaxy has a shallow potential and is on the verge of forming a disk. The combination of these factors give this galaxy the potential for strong gas outflows driven by supernova feedback and make it particularly sensitive to any changes to the simulation parameters. △ Less

Submitted 10 May, 2010; originally announced May 2010.

Comments: 34 pages, 7 figures. Accepted for publication in ApJ.

The ACS Nearby Galaxy Survey Treasury

Authors: Julianne J. Dalcanton, Benjamin F. Williams, Anil C. Seth, Andrew Dolphin, Jon Holtzman, Keith Rosema, Evan D. Skillman, Andrew Cole, Leo Girardi, Stephanie M. Gogarten, Igor D. Karachentsev, Knut Olsen, Daniel Weisz, Charlotte Christensen, Ken Freeman, Karoline Gilbert, Carme Gallart, Jason Harris, Paul Hodge, Roelof S. de Jong, Valentina Karachentseva, Mario Mateo, Peter B. Stetson, Maritza Tavarez, Dennis Zaritsky , et al. (2 additional authors not shown)

Abstract: The ACS Nearby Galaxy Survey Treasury (ANGST) is a systematic survey to establish a legacy of uniform multi-color photometry of resolved stars for a volume-limited sample of nearby galaxies (D<4 Mpc). The survey volume encompasses 69 galaxies in diverse environments, including close pairs, small & large groups, filaments, and truly isolated regions. The galaxies include a nearly complete range o… ▽ More The ACS Nearby Galaxy Survey Treasury (ANGST) is a systematic survey to establish a legacy of uniform multi-color photometry of resolved stars for a volume-limited sample of nearby galaxies (D<4 Mpc). The survey volume encompasses 69 galaxies in diverse environments, including close pairs, small & large groups, filaments, and truly isolated regions. The galaxies include a nearly complete range of morphological types spanning a factor of ~10^4 in luminosity and star formation rate. The survey data consists of images taken with ACS on HST, supplemented with archival data and new WFPC2 imaging taken after the failure of ACS. Survey images include wide field tilings covering the full radial extent of each galaxy, and single deep pointings in uncrowded regions of the most massive galaxies in the volume. The new wide field imaging in ANGST reaches median 50% completenesses of m_F475W=28.0 mag, m_F606W=27.3 mag, and m_F814W=27.3 mag, several magnitudes below the tip of the red giant branch (TRGB). The deep fields reach magnitudes sufficient to fully resolve the structure in the red clump. The resulting photometric catalogs are publicly accessible and contain over 34 million photometric measurements of >14 million stars. In this paper we present the details of the sample selection, imaging, data reduction, and the resulting photometric catalogs, along with an analysis of the photometric uncertainties (systematic and random), for both the ACS and WFPC2 imaging. We also present uniformly derived relative distances measured from the apparent magnitude of the TRGB. △ Less

Submitted 25 May, 2009; v1 submitted 22 May, 2009; originally announced May 2009.

Comments: 54 pages, including 24 pages of figures and 16 pages of tables. Project website and data available at http://www.nearbygalaxies.org/ . Data is also available through MAST. Scheduled to appear in the Astrophysical Journal Supplements. (Replaced to fix several figures that were damaged during compression)

Journal ref: Astrophys.J.Supp.183:67-108,2009

Increasing the Number of Underrepresented Minorities in Astronomy Through K-12 Education and Public Outreach (Paper II)

Authors: Dara Norman, CSMA, NSHP, Marcel Agueros, Scott F. Anderson, Andrew Baker, Adam Burgasser, Kelle Cruz, David J. Ernst, Eric Gawiser, Anita Krishnamurthi, Hyun-chul Lee, Kenneth Mighell, Charles McGruder, Philip J. Sakimoto, Kartik Sheth, Dave Soderblom, Michael Strauss, Donald Walter, Andrew West, UW Pre-Map staff - Eric Agol, Jeremiah Murphy, Sarah Garner, Jill Bellovary, Sarah Schmidt , et al. (8 additional authors not shown)

Abstract: In order to attract, recruit and retain underrepresented minority students to pursue Astronomy and related fields, we must ensure that there continues to be a well qualified pool of graduate and undergraduate students from which to recruit. This required pool of people are today's elementary, middle and high school students. The Astronomy community must be proactive in demonstrating the importan… ▽ More In order to attract, recruit and retain underrepresented minority students to pursue Astronomy and related fields, we must ensure that there continues to be a well qualified pool of graduate and undergraduate students from which to recruit. This required pool of people are today's elementary, middle and high school students. The Astronomy community must be proactive in demonstrating the importance of pursing scientific study and careers to these students and their parents. Only by actively engaging these communities can U.S Astronomy hope to increase the numbers of minority PhDs and continue to be a leader in Astronomical discovery and knowledge. △ Less

Submitted 26 March, 2009; originally announced March 2009.

Comments: cover page + 8 pages, An Astro2010 Decadal Survey State of the Profession Position Paper

Increasing the Number of Underrepresented Minorities in Astronomy at the Undergraduate, Graduate, and Postdoctoral Levels (Paper I)

Authors: Dara Norman, CSMA, NSHP, Marcel Agueros, Scott F. Anderson, Andrew Baker, Adam Burgasser, Kelle Cruz, David J. Ernst, Eric Gawiser, Anita Krishnamurthi, Hyun-chul Lee, Kenneth Mighell, Charles McGruder, Philip J. Sakimoto, Kartik Sheth, Dave Soderblom, Michael Strauss, Donald Walter, Andrew West, UW Pre-Map staff - Eric Agol, Jeremiah Murphy, Sarah Garner, Jill Bellovary, Sarah Schmidt , et al. (8 additional authors not shown)

Abstract: If the ethnic makeup of the astronomy profession is to achieve parity with the general population within one generation (~30 years), the number of underrepresented minorities earning graduate degrees in astronomy and astrophysics must increase in the coming decade by a factor of 5 to 10. To accomplish this, the profession must develop and invest in mechanisms to more effectively move individuals… ▽ More If the ethnic makeup of the astronomy profession is to achieve parity with the general population within one generation (~30 years), the number of underrepresented minorities earning graduate degrees in astronomy and astrophysics must increase in the coming decade by a factor of 5 to 10. To accomplish this, the profession must develop and invest in mechanisms to more effectively move individuals across critical educational junctures to the PhD and beyond. Early and continuous research engagement starting in the undergraduate years is critical to this vision, in which the federally funded research internship programs (e.g. NSF REU, NASA GSRP) and national centers/observatories play a vital role. Regionally based partnerships with minority-serving institutions (MSIs) are crucial for tapping extant pools of minority talent, as are post-baccalaurate and/or masters degree "bridging" programs that provide critical stepping stones to the PhD. Because of the strong undergraduate physics, engineering, and computer science backgrounds of many students from MSIs, we suggest that instrument development and large scale computing/data-mining are particularly promising avenues for engagement in the coming decade. △ Less

Submitted 26 March, 2009; originally announced March 2009.

Comments: cover page + 15 pages, An Astro2010 Decadal Survey State of the Profession Position Paper

Increasing the Number of Underrepresented Minorities in Astronomy: Executive Summary

Authors: Dara Norman, CSMA, NSHP, Marcel Agueros, Scott F. Anderson, Andrew Baker, Adam Burgasser, Kelle Cruz, David J. Ernst, Eric Gawiser, Anita Krishnamurthi, Hyun-chul Lee, Kenneth Mighell, Charles McGruder, Philip J. Sakimoto, Kartik Sheth, Dave Soderblom, Michael Strauss, Donald Walter, Andrew West, UW Pre-Map staff - Eric Agol, Jeremiah Murphy, Sarah Garner, Jill Bellovary, Sarah Schmidt , et al. (8 additional authors not shown)

Abstract: Promoting racial and ethnic diversity is critically important to the future success and growth of the field of astronomy. The raw ability, drive and interest required to excel in the field is distributed without regard to race, gender, or socioeconomic background. By not actively promoting diversity in our field we risk losing talented people to other professions (or losing them entirely), which… ▽ More Promoting racial and ethnic diversity is critically important to the future success and growth of the field of astronomy. The raw ability, drive and interest required to excel in the field is distributed without regard to race, gender, or socioeconomic background. By not actively promoting diversity in our field we risk losing talented people to other professions (or losing them entirely), which means that there will be astronomical discoveries that simply won't get made. There is demonstrated evidence that STEM fields benefit from diverse perspectives on problems that require more complex thought processes. This is especially relevant to a field like astronomy where more and more work is being done collaboratively. The lack of notable growth in African American, Hispanic, and Native American representation in astronomy indicates that the 'pipeline' for these individuals is systemically leaky at critical junctures. Substantially more effort must be directed toward improving the educational and career development of minorities to insure that these potential colleagues are supported through the process. However, simply recognizing that the pipeline is faulty is woefully inadequate. There must be very specific, targeted solutions to help improve the situation. With this in mind, we offer two position papers addressing specific areas of improvement that we identify as (a) essential for any foreseeable progress in the field, and (b) attainable in the 2010-2020 decade. These position papers focus primarily on African Americans, Hispanics, and Native Americans. Although we do not directly address issues of Asian Americans, Pacific Islanders, and other groups, many of the recommendations made here can be adapted to address issues faced by these groups as well. △ Less

Submitted 26 March, 2009; originally announced March 2009.

Comments: 5 pages, An Astro2010 Decadal Survey State of the Profession Position Paper

Forming a Large Disc Galaxy from a z<1 Major Merger

Authors: F. Governato, C. B. Brook, A. M. Brooks, L. Mayer, B. Willman, P. Jonsson, A. M. Stilp, L. Pope, C. Christensen, J. Wadsley, T. Quinn

Abstract: Using high resolution SPH simulations in a fully cosmological Lambda CDM context we study the formation of a bright disk dominated galaxy that originates from a "wet" major merger at z=0.8. The progenitors of the disk galaxy are themselves disk galaxies that formed from early major mergers between galaxies with blue colors. A substantial thin stellar disk grows rapidly following the last major m… ▽ More Using high resolution SPH simulations in a fully cosmological Lambda CDM context we study the formation of a bright disk dominated galaxy that originates from a "wet" major merger at z=0.8. The progenitors of the disk galaxy are themselves disk galaxies that formed from early major mergers between galaxies with blue colors. A substantial thin stellar disk grows rapidly following the last major merger and the present day properties of the final remnant are typical of early type spiral galaxies, with an i band B/D ~0.65, a disk scale length of 7.2 kpc, g-r = 0.5 mag, an HI line width (W_{20}/2) of 238 km/sec and total magnitude i = -22.4. The key ingredients for the formation of a dominant stellar disk component after a major merger are: i) substantial and rapid accretion of gas through cold flows followed at late times by cooling of gas from the hot phase, ii) supernova feedback that is able to partially suppress star formation during mergers and iii) relative fading of the spheroidal component. The gas fraction of the progenitors' disks does not exceed 25% at z<3, emphasizing that the continuous supply of gas from the local environment plays a major role in the regrowth of disks and in keeping the galaxies blue. The results of this simulation alleviate the problem posed for the existence of disk galaxies by the high likelihood of interactions and mergers for galaxy sized halos at relatively low z. △ Less

Submitted 20 July, 2009; v1 submitted 1 December, 2008; originally announced December 2008.

Comments: MNRAS in press. Minor additions to text. Movie at: http://www.astro.washington.edu/fabio/movies/Merger.mpg

Journal ref: MNRAS, 398, pp. 312-320 (2009)