-
The JWST Resolved Stellar Populations Early Release Science Program V. DOLPHOT Stellar Photometry for NIRCam and NIRISS
Authors:
Daniel R. Weisz,
Andrew E. Dolphin,
Alessandro Savino,
Kristen B. W. McQuinn,
Max J. B. Newman,
Benjamin F. Williams,
Nitya Kallivayalil,
Jay Anderson,
Martha L. Boyer,
Matteo Correnti,
Marla C. Geha,
Karin M. Sandstrom,
Andrew A. Cole,
Jack T. Warfield,
Evan D. Skillman,
Roger E. Cohen,
Rachael Beaton,
Alessandro Bressan,
Alberto Bolatto,
Michael Boylan-Kolchin,
Alyson M. Brooks,
James S. Bullock,
Charlie Conroy,
Michael C. Cooper,
Julianne J. Dalcanton
, et al. (16 additional authors not shown)
Abstract:
We present NIRCam and NIRISS modules for DOLPHOT, a widely-used crowded field stellar photometry package. We describe details of the modules including pixel masking, astrometric alignment, star finding, photometry, catalog creation, and artificial star tests (ASTs). We tested these modules using NIRCam and NIRISS images of M92 (a Milky Way globular cluster), Draco II (an ultra-faint dwarf galaxy),…
▽ More
We present NIRCam and NIRISS modules for DOLPHOT, a widely-used crowded field stellar photometry package. We describe details of the modules including pixel masking, astrometric alignment, star finding, photometry, catalog creation, and artificial star tests (ASTs). We tested these modules using NIRCam and NIRISS images of M92 (a Milky Way globular cluster), Draco II (an ultra-faint dwarf galaxy), and WLM (a star-forming dwarf galaxy). DOLPHOT's photometry is highly precise and the color-magnitude diagrams are deeper and have better definition than anticipated during original program design in 2017. The primary systematic uncertainties in DOLPHOT's photometry arise from mismatches in the model and observed point spread functions (PSFs) and aperture corrections, each contributing $\lesssim0.01$ mag to the photometric error budget. Version 1.2 of WebbPSF models, which include charge diffusion and interpixel capacitance effects, significantly reduced PSF-related uncertainties. We also observed minor ($\lesssim0.05$ mag) chip-to-chip variations in NIRCam's zero points, which will be addressed by the JWST flux calibration program. Globular cluster observations are crucial for photometric calibration. Temporal variations in the photometry are generally $\lesssim0.01$ mag, although rare large misalignment events can introduce errors up to 0.08 mag. We provide recommended DOLPHOT parameters, guidelines for photometric reduction, and advice for improved observing strategies. Our ERS DOLPHOT data products are available on MAST, complemented by comprehensive online documentation and tutorials for using DOLPHOT with JWST imaging data.
△ Less
Submitted 5 February, 2024;
originally announced February 2024.
-
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.
-
A Census from JWST of Extreme Emission Line Galaxies Spanning the Epoch of Reionization in CEERS
Authors:
Kelcey Davis,
Jonathan R. Trump,
Raymond C. Simons,
Elizabeth J. Mcgrath,
Stephen M. Wilkins,
Pablo Arrabal Haro,
Micaela B. Bagley,
Mark Dickinson,
Vital FernÁndez,
Ricardo O. AmorÍn,
Bren E. Backhaus,
Nikko J. Cleri,
Mario Llerena,
Samantha W. Brunker,
Guillermo Barro,
Laura Bisigello,
Madisyn Brooks,
Luca Costantin,
Alexander De La Vega,
Avishai Dekel,
Steven L. Finkelstein,
Nimish P. Hathi,
Michaela Hirschmann,
Jeyhan S. Kartaltepe,
Anton M. Koekemoer
, et al. (7 additional authors not shown)
Abstract:
We present a sample of 1165 extreme emission-line galaxies (EELGs) at 4<z<9 selected using James Webb Space Telescope (JWST) NIRCam photometry in the Cosmic Evolution Early Release Science (CEERS) program. We use a simple method to photometrically identify EELGs with Hb + [OIII] (combined) or Ha emission of observed-frame equivalent width EW >5000 AA. JWST/NIRSpec spectroscopic observations of a s…
▽ More
We present a sample of 1165 extreme emission-line galaxies (EELGs) at 4<z<9 selected using James Webb Space Telescope (JWST) NIRCam photometry in the Cosmic Evolution Early Release Science (CEERS) program. We use a simple method to photometrically identify EELGs with Hb + [OIII] (combined) or Ha emission of observed-frame equivalent width EW >5000 AA. JWST/NIRSpec spectroscopic observations of a subset (34) of the photometrically selected EELGs validate our selection method: all spectroscopically observed EELGs confirm our photometric identification of extreme emission, including some cases where the SED-derived photometric redshifts are incorrect. We find that the medium-band F410M filter in CEERS is particularly efficient at identifying EELGs, both in terms of including emission lines in the filter and in correctly identifying the continuum between Hb + [OIII] and Ha in the neighboring broad-band filters. We present examples of EELGs that could be incorrectly classified at ultra-high redshift (z>12) as a result of extreme Hb + [OIII] emission blended across the reddest photometric filters. We compare the EELGs to the broader (sub-extreme) galaxy population in the same redshift range and find that they are consistent with being the bluer, high equivalent width tail of a broader population of emission-line galaxies. The highest-EW EELGs tend to have more compact emission-line sizes than continuum sizes, suggesting that active galactic nuclei are responsible for at least some of the most extreme EELGs. Photometrically inferred emission-line ratios are consistent with ISM conditions with high ionization and moderately low metallicity, consistent with previous spectroscopic studies.
△ Less
Submitted 12 December, 2023;
originally announced December 2023.
-
Closing the Gap between Observed Low-Mass Galaxy HI Kinematics and CDM Predictions
Authors:
Amy Sardone,
Annika H. G. Peter,
Alyson M. Brooks,
Jane Kaczmarek
Abstract:
Testing the standard cosmological model ($Λ$CDM) at small scales is challenging. Galaxies that inhabit low-mass dark matter halos provide an ideal test bed for dark matter models by linking observational properties of galaxies at small scales (low mass, low velocity) to low-mass dark matter halos. However, the observed kinematics of these galaxies do not align with the kinematics of the dark matte…
▽ More
Testing the standard cosmological model ($Λ$CDM) at small scales is challenging. Galaxies that inhabit low-mass dark matter halos provide an ideal test bed for dark matter models by linking observational properties of galaxies at small scales (low mass, low velocity) to low-mass dark matter halos. However, the observed kinematics of these galaxies do not align with the kinematics of the dark matter halos predicted to host them, obscuring our understanding of the low-mass end of the galaxy-halo connection. We use deep HI observations of low-mass galaxies at high spectral resolution in combination with cosmological simulations of dwarf galaxies to better understand the connection between dwarf galaxy kinematics and low-mass halos. Specifically, we use HI line widths to directly compare to the maximum velocities in a dark matter halo, and find that each deeper measurement approaches the expected one-to-one relationship between the observed kinematics and the predicted kinematics in $Λ$CDM. We also measure baryonic masses and place these on the Baryonic Tully-Fisher relation (BTFR). Again, our deepest measurements approach the theoretical predictions for the low-mass end of this relation, a significant improvement on similar measurements based on line widths measured at 50\% and 20\% of the peak. Our data also hints at the rollover in the BTFR predicted by hydrodynamical simulations of $Λ$CDM for low-mass galaxies.
△ Less
Submitted 12 June, 2023;
originally announced June 2023.
-
Examining Uranus' zeta ring in Voyager 2 Wide-Angle-Camera Observations: Quantifying the Ring's Structure in 1986 and its Modifications prior to the Year 2007
Authors:
M. M. Hedman,
I. Regan,
T. Becker,
S. M. Brooks,
I. de Pater,
M. Showalter
Abstract:
The zeta ring is the innermost component of the Uranian ring system. It is of scientific interest because its morphology changed significantly between the Voyager 2 encounter in 1986 and subsequent Earth-based observations around 2007. It is also of practical interest because some Uranus mission concepts have the spacecraft pass through the inner flank of this ring. Recent re-examinations of the V…
▽ More
The zeta ring is the innermost component of the Uranian ring system. It is of scientific interest because its morphology changed significantly between the Voyager 2 encounter in 1986 and subsequent Earth-based observations around 2007. It is also of practical interest because some Uranus mission concepts have the spacecraft pass through the inner flank of this ring. Recent re-examinations of the Voyager 2 images have revealed additional information about this ring that provide a more complete picture of the ring's radial brightness profile and phase function. These data reveal that this ring's brightness varies with phase angle in a manner similar to other tenuous rings, consistent with it being composed primarily of sub-millimeter-sized particles. The total cross section of particles within this ring can also be estimated from these data, but translating that number into the actual risk to a spacecraft flying through this region depends on a number of model-dependent parameters. Fortunately, comparisons with Saturn's G and D rings allows the zeta-ring's particle number density to be compared with regions previously encountered by the Voyager and Cassini spacecraft. Finally, these data indicate that the observed changes in the zeta-ring's structure between 1986 and 2007 are primarily due to a substantial increase in the amount of dust at distances between 38,000 km and 40,000 km from Uranus' center.
△ Less
Submitted 31 May, 2023; v1 submitted 11 May, 2023;
originally announced May 2023.
-
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.
-
Catalog of Ultraviolet Bright Stars (CUBS): Strategies for UV occultation measurements, planetary illumination modeling, and sky map analyses using hybrid IUE-Kurucz spectra
Authors:
M. A. Velez,
K. D. Retherford,
V. Hue,
J. A. Kammer,
T. M. Becker,
G. R. Gladstone,
M. W. Davis,
T. K. Greathouse,
P. M. Molyneux,
S. M. Brooks,
U. Raut,
M. H. Versteeg
Abstract:
Ultraviolet spectroscopy is a powerful method to study planetary surface composition through reflectance measurements and atmospheric composition through stellar/solar occultations, transits of other planetary bodies, and direct imaging of airglow and auroral emissions. The next generation of ultraviolet spectrographs (UVS) on board ESA's JUICE (Jupiter Icy Moons Explorer) and NASA's Europa Clippe…
▽ More
Ultraviolet spectroscopy is a powerful method to study planetary surface composition through reflectance measurements and atmospheric composition through stellar/solar occultations, transits of other planetary bodies, and direct imaging of airglow and auroral emissions. The next generation of ultraviolet spectrographs (UVS) on board ESA's JUICE (Jupiter Icy Moons Explorer) and NASA's Europa Clipper missions will perform such measurements of Jupiter and its moons in the early 2030's. This work presents a compilation of a detailed UV stellar catalog, named CUBS, of targets with high intensity in the 50-210 nm wavelength range with applications relevant to planetary spectroscopy. These applications include: 1) Planning and simulating occultations, including calibration measurements; 2) Modeling starlight illumination of dark, nightside planetary surfaces primarily lit by the sky; and 3) Studying the origin of diffuse galactic UV light as mapped by existing datasets from Juno-UVS and others. CUBS includes information drawn from resources such as the International Ultraviolet Explorer (IUE) catalog and SIMBAD. We have constructed model spectra at 0.1 nm resolution for almost 90,000 targets using Kurucz models and, when available, IUE spectra. CUBS also includes robust checks for agreement between the Kurucz models and the IUE data. We also present a tool for which our catalog can be used to identify the best candidates for stellar occultation observations, with applications for any UV instrument. We report on our methods for producing CUBS and discuss plans for its implementation during ongoing and upcoming planetary missions.
△ Less
Submitted 7 March, 2023;
originally announced March 2023.
-
The JWST Resolved Stellar Populations Early Release Science Program II. Survey Overview
Authors:
Daniel R. Weisz,
Kristen B. W. McQuinn,
Alessandro Savino,
Nitya Kallivayalil,
Jay Anderson,
Martha L. Boyer,
Matteo Correnti,
Marla C. Geha,
Andrew E. Dolphin,
Karin M. Sandstrom,
Andrew A. Cole,
Benjamin F. Williams,
Evan D. Skillman,
Roger E. Cohen,
Max J. B. Newman,
Rachael Beaton,
Alessandro Bressan,
Alberto Bolatto,
Michael Boylan-Kolchin,
Alyson M. Brooks,
James S. Bullock,
Charlie Conroy,
M. C. Cooper,
Julianne J. Dalcanton,
Aaron L. Dotter
, et al. (17 additional authors not shown)
Abstract:
We present the JWST Resolved Stellar Populations Early Release Science (ERS) science program. We obtained 27.5 hours of NIRCam and NIRISS imaging of three targets in the Local Group (Milky Way globular cluster M92, ultra-faint dwarf galaxy Draco II, star-forming dwarf galaxy WLM), which span factors of $\sim10^5$ in luminosity, $\sim10^4$ in distance, and $\sim10^5$ in surface brightness. We descr…
▽ More
We present the JWST Resolved Stellar Populations Early Release Science (ERS) science program. We obtained 27.5 hours of NIRCam and NIRISS imaging of three targets in the Local Group (Milky Way globular cluster M92, ultra-faint dwarf galaxy Draco II, star-forming dwarf galaxy WLM), which span factors of $\sim10^5$ in luminosity, $\sim10^4$ in distance, and $\sim10^5$ in surface brightness. We describe the survey strategy, scientific and technical goals, implementation details, present select NIRCam color-magnitude diagrams (CMDs), and validate the NIRCam exposure time calculator (ETC). Our CMDs are among the deepest in existence for each class of target. They touch the theoretical hydrogen burning limit in M92 ($<0.08$ $M_{\odot}$; SNR $\sim5$ at $m_{F090W}\sim28.2$; $M_{F090W}\sim+13.6$), include the lowest-mass stars observed outside the Milky Way in Draco II (0.09 $M_{\odot}$; SNR $=10$ at $m_{F090W}\sim29$; $M_{F090W}\sim+12.1$), and reach $\sim1.5$ magnitudes below the oldest main sequence turnoff in WLM (SNR $=10$ at $m_{F090W}\sim29.5$; $M_{F090W}\sim+4.6$). The PARSEC stellar models provide a good qualitative match to the NIRCam CMDs, though are $\sim0.05$ mag too blue compared to M92 F090W$-$F150W data. The NIRCam ETC (v2.0) matches the SNRs based on photon noise from DOLPHOT stellar photometry in uncrowded fields, but the ETC may not be accurate in more crowded fields, similar to what is known for HST. We release beta versions of DOLPHOT NIRCam and NIRISS modules to the community. Results from this ERS program will establish JWST as the premier instrument for resolved stellar populations studies for decades to come.
△ Less
Submitted 11 January, 2023;
originally announced January 2023.
-
AGN quenching in simulated dwarf galaxies
Authors:
Ray S. Sharma,
Alyson M. Brooks,
Michael Tremmel,
Jillian Bellovary,
Thomas R. Quinn
Abstract:
We examine the quenching characteristics of $328$ isolated dwarf galaxies $\left(10^{8} < M_{\rm star}/M_\odot < 10^{10} \right)$ within the \Rom{} cosmological hydrodynamic simulation. Using mock observation methods, we identify isolated dwarf galaxies with quenched star formation and make direct comparisons to the quenched fraction in the NASA Sloan Atlas (NSA). Similar to other cosmological sim…
▽ More
We examine the quenching characteristics of $328$ isolated dwarf galaxies $\left(10^{8} < M_{\rm star}/M_\odot < 10^{10} \right)$ within the \Rom{} cosmological hydrodynamic simulation. Using mock observation methods, we identify isolated dwarf galaxies with quenched star formation and make direct comparisons to the quenched fraction in the NASA Sloan Atlas (NSA). Similar to other cosmological simulations, we find a population of quenched, isolated dwarf galaxies below $M_{\rm star} < 10^{9} M_\odot$ not detected within the NSA. We find that the presence of massive black holes (MBHs) in \Rom{} is largely responsible for the quenched, isolated dwarfs, while isolated dwarfs without an MBH are consistent with quiescent fractions observed in the field. Quenching occurs between $z=0.5-1$, during which the available supply of star-forming gas is heated or evacuated by MBH feedback. Mergers or interactions seem to play little to no role in triggering the MBH feedback. At low stellar masses, $M_{\rm star} \lesssim 10^{9.3} M_\odot$, quenching proceeds across several Gyr as the MBH slowly heats up gas in the central regions. At higher stellar masses, $M_{\rm star} \gtrsim 10^{9.3} M_\odot$, quenching occurs rapidly within $1$ Gyr, with the MBH evacuating gas from the central few kpc of the galaxy and driving it to the outskirts of the halo. Our results indicate the possibility of substantial star formation suppression via MBH feedback within dwarf galaxies in the field. On the other hand, the apparent over-quenching of dwarf galaxies due to MBH suggests higher resolution and/or better modeling is required for MBHs in dwarfs, and quenched fractions offer the opportunity to constrain current models.
△ Less
Submitted 9 November, 2022;
originally announced November 2022.
-
A hidden population of massive black holes in simulated dwarf galaxies
Authors:
Ray S. Sharma,
Alyson M. Brooks,
Michael Tremmel,
Jillian Bellovary,
Angelo Ricarte,
Thomas R. Quinn
Abstract:
We explore the characteristics of actively accreting MBHs within dwarf galaxies in the \textsc{Romulus25} cosmological hydrodynamic simulation. We examine the MBH occupation fraction, x-ray active fractions, and AGN scaling relations within dwarf galaxies of stellar mass $10^{8} < M_{\rm star} < 10^{10} M_\odot$ out to redshift $z=2$. In the local universe, the MBH occupation fraction is consisten…
▽ More
We explore the characteristics of actively accreting MBHs within dwarf galaxies in the \textsc{Romulus25} cosmological hydrodynamic simulation. We examine the MBH occupation fraction, x-ray active fractions, and AGN scaling relations within dwarf galaxies of stellar mass $10^{8} < M_{\rm star} < 10^{10} M_\odot$ out to redshift $z=2$. In the local universe, the MBH occupation fraction is consistent with observed constraints, dropping below unity at $M_{\rm star} < 3\times10^{10} M_{\odot}$, $M_{\rm 200} < 3\times10^{11} M_\odot$. Local dwarf AGN in \textsc{Romulus25} follow observed scaling relations between AGN x-ray luminosity, stellar mass, and star formation rate, though they exhibit slightly higher active fractions and number densities than comparable x-ray observations. Since $z=2$, the MBH occupation fraction has decreased, the population of dwarf AGN has become overall less luminous, and as a result, the overall number density of dwarf AGN has diminished. We predict the existence of a large population of MBHs in the local universe with low x-ray luminosities and high contamination from x-ray binaries and the hot interstellar medium that are undetectable by current x-ray surveys. These hidden MBHs make up $76\%$ of all MBHs in local dwarf galaxies, and include many MBHs that are undermassive relative to their host galaxy's stellar mass. Their detection relies not only on greater instrument sensitivity but on better modeling of x-ray contaminants or multi-wavelength surveys. Our results indicate dwarf AGN were substantially more active in the past despite being low-luminosity today, and indicate future deep x-ray surveys may uncover many hidden MBHs in dwarf galaxies out to at least $z=2$.
△ Less
Submitted 1 August, 2022; v1 submitted 10 March, 2022;
originally announced March 2022.
-
What's in a name? Quantifying the interplay between the Definition, Orientation, and Shape of Ultra-diffuse Galaxies Using the Romulus Simulations
Authors:
Jordan D. Van Nest,
F. Munshi,
A. C. Wright,
M. Tremmel,
A. M. Brooks,
D. Nagai,
T. Quinn
Abstract:
We explore populations of ultra-diffuse galaxies (UDGs) in isolated, satellite, and cluster environments using the Romulus25 and RomulusC simulations, including how the populations vary with UDG definition and viewing orientation. Using a fiducial definition of UDGs, we find that isolated UDGs have notably larger semi-major (b/a) and smaller semi-minor (c/a) axis ratios than their non-UDG counterp…
▽ More
We explore populations of ultra-diffuse galaxies (UDGs) in isolated, satellite, and cluster environments using the Romulus25 and RomulusC simulations, including how the populations vary with UDG definition and viewing orientation. Using a fiducial definition of UDGs, we find that isolated UDGs have notably larger semi-major (b/a) and smaller semi-minor (c/a) axis ratios than their non-UDG counterparts, i.e., they are more oblate, or diskier. This is in line with previous results that adopted the same UDG definition and showed that isolated UDGs form via early, high-spin mergers. However, the choice of UDG definition can drastically affect what subsets of a dwarf population are classified as UDGs, changing the number of UDGs by up to approximately 45% of the dwarf population. We also find that a galaxy's classification as a UDG is dependent on its viewing orientation, and this dependence decreases as environmental density increases. Overall, we conclude that some definitions for UDGs used in the literature manage to isolate a specific formation mechanism for isolated dwarfs, while less restrictive definitions erase a link to the formation mechanism. Thus, how we define UDG populations must be considered if we want to understand the formation and evolution of UDGs.
△ Less
Submitted 25 September, 2023; v1 submitted 30 August, 2021;
originally announced August 2021.
-
Joint Gas and Stellar Dynamical Models of WLM: An isolated dwarf galaxy within a cored, prolate DM halo
Authors:
Gigi Y. C. Leung,
Ryan Leaman,
Giuseppina Battaglia,
Glenn van de Ven,
Alyson M. Brooks,
Jorge Peñarrubia,
Kim A. Venn
Abstract:
We present multi-tracer dynamical models of the low mass ($M_{*} \sim 10^{7}$), isolated dwarf irregular galaxy WLM in order to simultaneously constrain the inner slope of the dark matter (DM) halo density profile ($γ$) and flattening ($q_\mathrm{DM}$), and the stellar orbital anisotropy ($β_{z}, β_{r}$). For the first time, we show how jointly constraining the mass distribution from the HI gas ro…
▽ More
We present multi-tracer dynamical models of the low mass ($M_{*} \sim 10^{7}$), isolated dwarf irregular galaxy WLM in order to simultaneously constrain the inner slope of the dark matter (DM) halo density profile ($γ$) and flattening ($q_\mathrm{DM}$), and the stellar orbital anisotropy ($β_{z}, β_{r}$). For the first time, we show how jointly constraining the mass distribution from the HI gas rotation curve and solving the Jeans' equations with discrete stellar kinematics leads to a factor of $\sim2$ reduction in the uncertainties on $γ$. The mass-anisotropy degeneracy is also partially broken, leading to reductions on uncertainty by $\sim 30\%$ on $M_\mathrm{vir}$ (and $\sim 70\%$ at the half-light radius) and $\sim 25\%$ on anisotropy. Our inferred value of $γ= 0.3 \pm 0.1$ is robust to the halo geometry, and in excellent agreement with predictions of stellar feedback driven DM core creation. The derived prolate geometry of the DM halo with $q_\mathrm{DM} = 2 \pm 1$ is consistent with $Λ$CDM simulations of dwarf galaxy halos. While self-interacting DM (SIDM) models with $σ/m_{X} \sim 0.6$ can reproduce this cored DM profile, the interaction events may sphericalise the halo. The simultaneously cored and prolate DM halo may therefore present a challenge for SIDM. Finally we find that the radial profile of stellar anisotropy in WLM ($β_{r}$) follows a nearly identical trend of increasing tangential anisotropy to the classical dSphs, Fornax and Sculptor. Given WLM's orbital history, this result may call into question whether such anisotropy is a consequence of tidal stripping in only one pericentric passage or if it instead is a feature of the largely self-similar formation and evolutionary pathways for some dwarf galaxies.
△ Less
Submitted 9 April, 2021;
originally announced April 2021.
-
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.
-
Frontiers in Planetary Rings Science
Authors:
Shawn M. Brooks,
Tracy M. Becker,
K. Baillié,
H. N. Becker,
E. T. Bradley,
J. E. Colwell,
J. N. Cuzzi,
I. de Pater,
S. Eckert,
M. El Moutamid,
S. G. Edgington,
P. R. Estrada,
M. W. Evans,
A. Flandes,
R. G. French,
Á. García,
M. K. Gordon,
M. M. Hedman,
H. -W. Hsu,
R. G. Jerousek,
E. A. Marouf,
B. K. Meinke,
P. D. Nicholson,
S. H. Pilorz,
M. R. Showalter
, et al. (3 additional authors not shown)
Abstract:
We now know that the outer solar system is host to at least six diverse planetary ring systems, each of which is a scientifically compelling target with the potential to inform us about the evolution, history and even the internal structure of the body it adorns. These diverse ring systems represent a set of distinct local laboratories for understanding the physics and dynamics of planetary disks,…
▽ More
We now know that the outer solar system is host to at least six diverse planetary ring systems, each of which is a scientifically compelling target with the potential to inform us about the evolution, history and even the internal structure of the body it adorns. These diverse ring systems represent a set of distinct local laboratories for understanding the physics and dynamics of planetary disks, with applications reaching beyond our Solar System. We highlight the current status of planetary rings science and the open questions before the community to promote continued Earth-based and spacecraft-based investigations into planetary rings. As future spacecraft missions are launched and more powerful telescopes come online in the decades to come, we urge NASA for continued support of investigations that advance our understanding of planetary rings, through research and analysis of data from existing facilities, more laboratory work and specific attention to strong rings science goals during future mission selections.
△ Less
Submitted 27 August, 2020;
originally announced August 2020.
-
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.
-
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.
-
The Saturn Ring Skimmer Mission Concept: The next step to explore Saturn's rings, atmosphere, interior, and inner magnetosphere
Authors:
Matthew S. Tiscareno,
Mar Vaquero,
Matthew M. Hedman,
Hao Cao,
Paul R. Estrada,
Andrew P. Ingersoll,
Kelly E. Miller,
Marzia Parisi,
David. H. Atkinson,
Shawn M. Brooks,
Jeffrey N. Cuzzi,
James Fuller,
Amanda R. Hendrix,
Robert E. Johnson,
Tommi Koskinen,
William S. Kurth,
Jonathan I. Lunine,
Philip D. Nicholson,
Carol S. Paty,
Rebecca Schindhelm,
Mark R. Showalter,
Linda J. Spilker,
Nathan J. Strange,
Wendy Tseng
Abstract:
The innovative Saturn Ring Skimmer mission concept enables a wide range of investigations that address fundamental questions about Saturn and its rings, as well as giant planets and astrophysical disk systems in general. This mission would provide new insights into the dynamical processes that operate in astrophysical disk systems by observing individual particles in Saturn's rings for the first t…
▽ More
The innovative Saturn Ring Skimmer mission concept enables a wide range of investigations that address fundamental questions about Saturn and its rings, as well as giant planets and astrophysical disk systems in general. This mission would provide new insights into the dynamical processes that operate in astrophysical disk systems by observing individual particles in Saturn's rings for the first time. The Ring Skimmer would also constrain the origin, history, and fate of Saturn's rings by determining their compositional evolution and material transport rates. In addition, the Ring Skimmer would reveal how the rings, magnetosphere, and planet operate as an inter-connected system by making direct measurements of the ring's atmosphere, Saturn's inner magnetosphere and the material owing from the rings into the planet. At the same time, this mission would clarify the dynamical processes operating in the planet's visible atmosphere and deep interior by making extensive high-resolution observations of cloud features and repeated measurements of the planet's extremely dynamic gravitational field. Given the scientific potential of this basic mission concept, we advocate that it be studied in depth as a potential option for the New Frontiers program.
△ Less
Submitted 16 September, 2020; v1 submitted 30 July, 2020;
originally announced July 2020.
-
The Diversity and Variability of Star Formation Histories in Models of Galaxy Evolution
Authors:
Kartheik G. Iyer,
Sandro Tacchella,
Shy Genel,
Christopher C. Hayward,
Lars Hernquist,
Alyson M. Brooks,
Neven Caplar,
Romeel Davé,
Benedikt Diemer,
John C. Forbes,
Eric Gawiser,
Rachel S. Somerville,
Tjitske K. Starkenburg
Abstract:
Understanding the variability of galaxy star formation histories (SFHs) across a range of timescales provides insight into the underlying physical processes that regulate star formation within galaxies. We compile the SFHs of galaxies at $z=0$ from an extensive set of models, ranging from cosmological hydrodynamical simulations (Illustris, IllustrisTNG, Mufasa, Simba, EAGLE), zoom simulations (FIR…
▽ More
Understanding the variability of galaxy star formation histories (SFHs) across a range of timescales provides insight into the underlying physical processes that regulate star formation within galaxies. We compile the SFHs of galaxies at $z=0$ from an extensive set of models, ranging from cosmological hydrodynamical simulations (Illustris, IllustrisTNG, Mufasa, Simba, EAGLE), zoom simulations (FIRE-2, g14, and Marvel/Justice League), semi-analytic models (Santa Cruz SAM) and empirical models (UniverseMachine), and quantify the variability of these SFHs on different timescales using the power spectral density (PSD) formalism. We find that the PSDs are well described by broken power-laws, and variability on long timescales ($\gtrsim1$ Gyr) accounts for most of the power in galaxy SFHs. Most hydrodynamical models show increased variability on shorter timescales ($\lesssim300$ Myr) with decreasing stellar mass. Quenching can induce $\sim0.4-1$ dex of additional power on timescales $>1$ Gyr. The dark matter accretion histories of galaxies have remarkably self-similar PSDs and are coherent with the in-situ star formation on timescales $>3$ Gyr. There is considerable diversity among the different models in their (i) power due to SFR variability at a given timescale, (ii) amount of correlation with adjacent timescales (PSD slope), (iii) evolution of median PSDs with stellar mass, and (iv) presence and locations of breaks in the PSDs. The PSD framework is a useful space to study the SFHs of galaxies since model predictions vary widely. Observational constraints in this space will help constrain the relative strengths of the physical processes responsible for this variability.
△ Less
Submitted 15 July, 2020;
originally announced July 2020.
-
The Formation of Isolated Ultra-Diffuse Galaxies in Romulus25
Authors:
Anna C. Wright,
Michael Tremmel,
Alyson M. Brooks,
Ferah Munshi,
Daisuke Nagai,
Ray S. Sharma,
Thomas R. Quinn
Abstract:
We use the \textsc{Romulus25} cosmological simulation volume to identify the largest-ever simulated sample of {\it field} ultra-diffuse galaxies (UDGs). At $z=0$, we find that isolated UDGs have average star formation rates, colors, and virial masses for their stellar masses and environment. UDGs have moderately elevated HI masses, being 70\% (300\%) more HI-rich than typical isolated dwarf galaxi…
▽ More
We use the \textsc{Romulus25} cosmological simulation volume to identify the largest-ever simulated sample of {\it field} ultra-diffuse galaxies (UDGs). At $z=0$, we find that isolated UDGs have average star formation rates, colors, and virial masses for their stellar masses and environment. UDGs have moderately elevated HI masses, being 70\% (300\%) more HI-rich than typical isolated dwarf galaxies at luminosities brighter (fainter) than M$_\mathrm{B}$=-14. However, UDGs are consistent with the general isolated dwarf galaxy population and make up $\sim$20\% of all field galaxies with 10$^7$<M$_\star$/M$_\odot$<10$^{9}$. The HI masses, effective radii, and overall appearances of our UDGs are consistent with existing observations of field UDGs, but we predict that many isolated UDGs have been missed by current surveys. Despite their isolation at $z=0$, the UDGs in our sample are the products of major mergers. Mergers are no more common in UDG than non-UDG progenitors, but mergers that create UDGs tend to happen earlier - almost never occurring after $z=1$, produce a temporary boost in spin, and cause star formation to be redistributed to the outskirts of galaxies, resulting in lower central star formation rates. The centers of the galaxies fade as their central stellar populations age, but their global star formation rates are maintained through bursts of star formation at larger radii, producing steeper negative g-r color gradients. This formation channel is unique relative to other proposals for UDG formation in isolated galaxies, demonstrating that UDGs can potentially be formed through multiple mechanisms.
△ Less
Submitted 7 January, 2021; v1 submitted 15 May, 2020;
originally announced May 2020.
-
Self-Interacting Dark Matter and the Delay of Super-Massive Black Hole Growth
Authors:
Akaxia Cruz,
Andrew Pontzen,
Marta Volonteri,
Thomas R. Quinn,
Michael Tremmel,
Alyson M. Brooks,
Nicole N. Sanchez,
Ferah Munshi,
Arianna Di Cintio
Abstract:
Using cosmological hydrodynamic simulations with physically motivated models of super-massive black hole (SMBH) formation and growth, we compare the assembly of Milky Way-mass ( $M_{\mathrm{vir}} \approx 7 \times 10^{11}$ $M_{\odot}$ at $z = 0$) galaxies in cold dark matter (CDM) and self-interacting dark matter (SIDM) models. Our SIDM model adopts a constant cross-section of 1 cm$^2$/g. We find t…
▽ More
Using cosmological hydrodynamic simulations with physically motivated models of super-massive black hole (SMBH) formation and growth, we compare the assembly of Milky Way-mass ( $M_{\mathrm{vir}} \approx 7 \times 10^{11}$ $M_{\odot}$ at $z = 0$) galaxies in cold dark matter (CDM) and self-interacting dark matter (SIDM) models. Our SIDM model adopts a constant cross-section of 1 cm$^2$/g. We find that SMBH formation is suppressed in the early universe due to SIDM coring. SMBH-SMBH mergers are also suppressed in SIDM as a consequence of the lower number of SMBHs formed. Lack of initial merger-driven SMBH growth in turn delays SMBH growth by billions of years in SIDM compared to CDM. Further, we find that this delayed growth suppresses SMBH accretion in the largest progenitors of the main SIDM galaxies during the first 5 Gyrs of their evolution. Nonetheless, by $z = 0.8$ the CDM and SIDM SMBH masses differ only by around 0.2 dex, so that both remain compatible with the $M_{BH}-M_{*}$ relation. We show that the reduced accretion causes the SIDM SMBHs to less aggressively regulate star formation in their host galaxies than their CDM counterparts, resulting in a factor of 3 or more stars being produced over the lifetime of the SIDM galaxies compared to the CDM galaxies. Our results highlight a new way in which SIDM can affect the growth and merger history of SMBHs and ultimately give rise to very different galaxy evolution compared to the classic CDM model.
△ Less
Submitted 17 April, 2020;
originally announced April 2020.
-
The Formation of Ultra-Diffuse Galaxies in the RomulusC Galaxy Cluster Simulation
Authors:
Michael Tremmel,
Anna C. Wright,
Alyson M. Brooks,
Ferah Munshi,
Daisuke Nagai,
Thomas R. Quinn
Abstract:
We study the origins of 122 ultra-diffuse galaxies (UDGs) in the {\sc RomulusC} zoom-in cosmological simulation of a galaxy cluster (M$_{200} = 1.15\times10^{14}$ M$_{\odot}$), one of the only such simulations capable of resolving the evolution and structure of dwarf galaxies (M$_{\star} < 10^9$ M$_{\odot}$). We find broad agreement with observed cluster UDGs and predict that they are not separate…
▽ More
We study the origins of 122 ultra-diffuse galaxies (UDGs) in the {\sc RomulusC} zoom-in cosmological simulation of a galaxy cluster (M$_{200} = 1.15\times10^{14}$ M$_{\odot}$), one of the only such simulations capable of resolving the evolution and structure of dwarf galaxies (M$_{\star} < 10^9$ M$_{\odot}$). We find broad agreement with observed cluster UDGs and predict that they are not separate from the overall cluster dwarf population. UDGs in cluster environments form primarily from dwarf galaxies that experienced early cluster in-fall and subsequent quenching due to ram pressure. The ensuing dimming of these dwarf galaxies due to passive stellar evolution results in a population of very low surface brightness galaxies that are otherwise typical dwarfs. UDGs and non-UDGs alike are affected by tidal interactions with the cluster potential. Tidal stripping of dark matter, as well as mass loss from stellar evolution, results in the adiabatic expansion of stars, particularly in the lowest mass dwarfs. High mass dwarf galaxies show signatures of tidal heating while low mass dwarfs that survive until $z=0$ typically have not experienced such impulsive interactions. There is little difference between UDGs and non-UDGs in terms of their dark matter halos, stellar morphology, colors, and location within the cluster. In most respects cluster UDG and non-UDGs alike are similar to isolated dwarf galaxies, except for the fact that they are typically quenched.
△ Less
Submitted 13 July, 2020; v1 submitted 15 August, 2019;
originally announced August 2019.
-
Dark Matter Science in the Era of LSST
Authors:
Keith Bechtol,
Alex Drlica-Wagner,
Kevork N. Abazajian,
Muntazir Abidi,
Susmita Adhikari,
Yacine Ali-Haïmoud,
James Annis,
Behzad Ansarinejad,
Robert Armstrong,
Jacobo Asorey,
Carlo Baccigalupi,
Arka Banerjee,
Nilanjan Banik,
Charles Bennett,
Florian Beutler,
Simeon Bird,
Simon Birrer,
Rahul Biswas,
Andrea Biviano,
Jonathan Blazek,
Kimberly K. Boddy,
Ana Bonaca,
Julian Borrill,
Sownak Bose,
Jo Bovy
, et al. (155 additional authors not shown)
Abstract:
Astrophysical observations currently provide the only robust, empirical measurements of dark matter. In the coming decade, astrophysical observations will guide other experimental efforts, while simultaneously probing unique regions of dark matter parameter space. This white paper summarizes astrophysical observations that can constrain the fundamental physics of dark matter in the era of LSST. We…
▽ More
Astrophysical observations currently provide the only robust, empirical measurements of dark matter. In the coming decade, astrophysical observations will guide other experimental efforts, while simultaneously probing unique regions of dark matter parameter space. This white paper summarizes astrophysical observations that can constrain the fundamental physics of dark matter in the era of LSST. We describe how astrophysical observations will inform our understanding of the fundamental properties of dark matter, such as particle mass, self-interaction strength, non-gravitational interactions with the Standard Model, and compact object abundances. Additionally, we highlight theoretical work and experimental/observational facilities that will complement LSST to strengthen our understanding of the fundamental characteristics of dark matter.
△ Less
Submitted 11 March, 2019;
originally announced March 2019.
-
Probing the Fundamental Nature of Dark Matter with the Large Synoptic Survey Telescope
Authors:
Alex Drlica-Wagner,
Yao-Yuan Mao,
Susmita Adhikari,
Robert Armstrong,
Arka Banerjee,
Nilanjan Banik,
Keith Bechtol,
Simeon Bird,
Kimberly K. Boddy,
Ana Bonaca,
Jo Bovy,
Matthew R. Buckley,
Esra Bulbul,
Chihway Chang,
George Chapline,
Johann Cohen-Tanugi,
Alessandro Cuoco,
Francis-Yan Cyr-Racine,
William A. Dawson,
Ana Díaz Rivero,
Cora Dvorkin,
Denis Erkal,
Christopher D. Fassnacht,
Juan García-Bellido,
Maurizio Giannotti
, et al. (75 additional authors not shown)
Abstract:
Astrophysical and cosmological observations currently provide the only robust, empirical measurements of dark matter. Future observations with Large Synoptic Survey Telescope (LSST) will provide necessary guidance for the experimental dark matter program. This white paper represents a community effort to summarize the science case for studying the fundamental physics of dark matter with LSST. We d…
▽ More
Astrophysical and cosmological observations currently provide the only robust, empirical measurements of dark matter. Future observations with Large Synoptic Survey Telescope (LSST) will provide necessary guidance for the experimental dark matter program. This white paper represents a community effort to summarize the science case for studying the fundamental physics of dark matter with LSST. We discuss how LSST will inform our understanding of the fundamental properties of dark matter, such as particle mass, self-interaction strength, non-gravitational couplings to the Standard Model, and compact object abundances. Additionally, we discuss the ways that LSST will complement other experiments to strengthen our understanding of the fundamental characteristics of dark matter. More information on the LSST dark matter effort can be found at https://lsstdarkmatter.github.io/ .
△ Less
Submitted 24 April, 2019; v1 submitted 4 February, 2019;
originally announced February 2019.
-
Understanding Dwarf Galaxies in order to Understand Dark Matter
Authors:
Alyson M. Brooks
Abstract:
Much progress has been made in recent years by the galaxy simulation community in making realistic galaxies, mostly by more accurately capturing the effects of baryons on the structural evolution of dark matter halos at high resolutions. This progress has altered theoretical expectations for galaxy evolution within a Cold Dark Matter (CDM) model, reconciling many earlier discrepancies between theo…
▽ More
Much progress has been made in recent years by the galaxy simulation community in making realistic galaxies, mostly by more accurately capturing the effects of baryons on the structural evolution of dark matter halos at high resolutions. This progress has altered theoretical expectations for galaxy evolution within a Cold Dark Matter (CDM) model, reconciling many earlier discrepancies between theory and observations. Despite this reconciliation, CDM may not be an accurate model for our Universe. Much more work must be done to understand the predictions for galaxy formation within alternative dark matter models.
△ Less
Submitted 30 November, 2018;
originally announced December 2018.
-
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.
-
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.
-
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.
-
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.
-
Going, going, gone dark: Quantifying the scatter in the faintest dwarf galaxies
Authors:
Ferah Munshi,
Alyson M. Brooks,
Elaad Applebaum,
Daniel R. Weisz,
Fabio Governato,
Thomas R. Quinn
Abstract:
We predict the stellar mass-halo mass (SMHM) relationship for dwarf galaxies and their satellites residing in halos down to M$_{halo} =$ 10$^7$ M$_{\odot}$ with 10$^4$ M$_{\odot} <$ M$_{star}$($z=0$) $< 10^8$ M$_{\odot}$, and quantify the predicted scatter in the relation at the low mass end, using cosmological simulations. The galaxies were drawn from a cosmological simulation of dwarf galaxies,…
▽ More
We predict the stellar mass-halo mass (SMHM) relationship for dwarf galaxies and their satellites residing in halos down to M$_{halo} =$ 10$^7$ M$_{\odot}$ with 10$^4$ M$_{\odot} <$ M$_{star}$($z=0$) $< 10^8$ M$_{\odot}$, and quantify the predicted scatter in the relation at the low mass end, using cosmological simulations. The galaxies were drawn from a cosmological simulation of dwarf galaxies, run with the N-body + SPH code, ChaNGA, at a high resolution of 60 pc. For M$_{halo} > 10^9$ M$_{\odot}$, the simulated SMHM relationship agrees with literature determinations, including exhibiting a small scatter. However, the scatter in the SMHM relation increases dramatically for lower-mass halos. We find that some of this scatter is due to {\em dark dwarfs}, halos devoid of stars. However, even when only considering well-resolved halos that contain a stellar population, the scatter in stellar mass reaches nearly 1 dex for M$_{halo}$($z=0$) 10$^7$ M$_{\odot}$. Much of this scatter is due to including satellites of the dwarf galaxies that have had their halo masses reduced through tidal stripping. The fraction of dark dwarfs (those that contain no stars) increases substantially with decreasing halo mass. When these dark halos are considered, the true scatter in the SMHM at low masses is even larger. At the faintest end of the SMHM relation probed by our simulations, a galaxy cannot be assigned a unique halo mass based solely on its luminosity. We provide a formula to stochastically populate low-mass halos following our results. Our predicted large scatter at low halo masses increases the slope of the resulting stellar mass function on the ultra-faint dwarf galaxy scales currently being probed by such surveys as the Dark Energy Survey or the Hyper-Suprime Cam Subaru Strategic Program, and in the future by the Large Synoptic Survey Telescope.
△ Less
Submitted 17 May, 2017;
originally announced May 2017.
-
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.
-
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.
-
Search for Gamma-ray Emission from Dark Matter Annihilation in the Small Magellanic Cloud with the Fermi Large Area Telescope
Authors:
Regina Caputo,
Matthew R. Buckley,
Pierrick Martin,
Eric Charles,
Alyson M. Brooks,
Alex Drlica-Wagner,
Jennifer M. Gaskins,
Matthew Wood
Abstract:
The Small Magellanic Cloud (SMC) is the second-largest satellite galaxy of the Milky Way and is only 60 kpc away. As a nearby, massive, and dense object with relatively low astrophysical backgrounds, it is a natural target for dark matter indirect detection searches. In this work, we use six years of Pass 8 data from the Fermi Large Area Telescope to search for gamma-ray signals of dark matter ann…
▽ More
The Small Magellanic Cloud (SMC) is the second-largest satellite galaxy of the Milky Way and is only 60 kpc away. As a nearby, massive, and dense object with relatively low astrophysical backgrounds, it is a natural target for dark matter indirect detection searches. In this work, we use six years of Pass 8 data from the Fermi Large Area Telescope to search for gamma-ray signals of dark matter annihilation in the SMC. Using data-driven fits to the gamma-ray backgrounds, and a combination of N-body simulations and direct measurements of rotation curves to estimate the SMC DM density profile, we found that the SMC was well described by standard astrophysical sources, and no signal from dark matter annihilation was detected. We set conservative upper limits on the dark matter annihilation cross section. These constraints are in agreement with stronger constraints set by searches in the Large Magellanic Cloud and approach the canonical thermal relic cross section at dark matter masses lower than 10 GeV in the $b\bar{b}$ and $τ^+τ^-$ channels.
△ Less
Submitted 29 March, 2016; v1 submitted 2 March, 2016;
originally announced March 2016.
-
Assessing Astrophysical Uncertainties in Direct Detection with Galaxy Simulations
Authors:
Jonathan D. Sloane,
Matthew R. Buckley,
Alyson M. Brooks,
Fabio Governato
Abstract:
We study the local dark matter velocity distribution in simulated Milky Way-mass galaxies, generated at high resolution with both dark matter and baryons. We find that the dark matter in the Solar neighborhood is influenced appreciably by the inclusion of baryons, increasing the speed of dark matter particles compared to dark matter-only simulations. The gravitational potential due to the presence…
▽ More
We study the local dark matter velocity distribution in simulated Milky Way-mass galaxies, generated at high resolution with both dark matter and baryons. We find that the dark matter in the Solar neighborhood is influenced appreciably by the inclusion of baryons, increasing the speed of dark matter particles compared to dark matter-only simulations. The gravitational potential due to the presence of a baryonic disk increases the amount of high velocity dark matter, resulting in velocity distributions which are more similar to the Maxwellian Standard Halo Model than predicted from dark matter-only simulations. Further, the velocity structures present in baryonic simulations possess a greater diversity than expected from dark matter-only simulation. We show the impact on the direct detection experiments LUX, DAMA/Libra, and CoGeNT using our simulated velocity distributions, and explore how resolution and halo mass within the Milky Way's estimated mass range impact the results. A Maxwellian fit to the velocity distribution tends to overpredict the amount of dark matter in the high velocity tail, even with baryons, and thus leads to overly optimistic direct detection bounds on models which are dependent on this region of phase space for an experimental signal. Our work further demonstrates that it is critical to transform simulated velocity distributions to the lab frame of reference, due to the fact that velocity structure in the Solar neighborhood appears when baryons are included. There is more velocity structure present when baryons are included than in dark matter-only simulations. Even when baryons are included, the importance of the velocity structure is not as apparent in the Galactic frame of reference as in the Earth frame.
△ Less
Submitted 29 July, 2016; v1 submitted 20 January, 2016;
originally announced January 2016.
-
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.
-
Particle tagging and its implications for stellar population dynamics
Authors:
Theo Le Bret,
Andrew Pontzen,
Andrew P. Cooper,
Carlos Frenk,
Adi Zolotov,
Alyson M. Brooks,
Fabio Governato,
Owen H. Parry
Abstract:
We establish a controlled comparison between the properties of galactic stellar halos obtained with hydrodynamical simulations and with `particle tagging'. Tagging is a fast way to obtain stellar population dynamics: instead of tracking gas and star formation, it `paints' stars directly onto a suitably defined subset of dark matter particles in a collisionless, dark-matter-only simulation.Our stud…
▽ More
We establish a controlled comparison between the properties of galactic stellar halos obtained with hydrodynamical simulations and with `particle tagging'. Tagging is a fast way to obtain stellar population dynamics: instead of tracking gas and star formation, it `paints' stars directly onto a suitably defined subset of dark matter particles in a collisionless, dark-matter-only simulation.Our study shows that there are conditions under which particle tagging generates good fits to the hydrodynamical stellar density profiles of a central Milky-Way-like galaxy and its most prominent substructure. Phase-space diffusion processes are crucial to reshaping the distribution of stars in infalling spheroidal systems and hence the final stellar halo. We conclude that the success of any particular tagging scheme hinges on this diffusion being taken into account, at a minimum by making use of `live' tagging schemes, in which particles are regularly tagged throughout the evolution of a galaxy.
△ Less
Submitted 23 February, 2015;
originally announced February 2015.
-
Search for Gamma-ray Emission from Dark Matter Annihilation in the Large Magellanic Cloud with the Fermi Large Area Telescope
Authors:
Matthew R. Buckley,
Eric Charles,
Jennifer M. Gaskins,
Alyson M. Brooks,
Alex Drlica-Wagner,
Pierrick Martin,
Geng Zhao
Abstract:
At a distance of 50 kpc and with a dark matter mass of $\sim10^{10}$ M$_{\odot}$, the Large Magellanic Cloud (LMC) is a natural target for indirect dark matter searches. We use five years of data from the Fermi Large Area Telescope (LAT) and updated models of the gamma-ray emission from standard astrophysical components to search for a dark matter annihilation signal from the LMC. We perform a rot…
▽ More
At a distance of 50 kpc and with a dark matter mass of $\sim10^{10}$ M$_{\odot}$, the Large Magellanic Cloud (LMC) is a natural target for indirect dark matter searches. We use five years of data from the Fermi Large Area Telescope (LAT) and updated models of the gamma-ray emission from standard astrophysical components to search for a dark matter annihilation signal from the LMC. We perform a rotation curve analysis to determine the dark matter distribution, setting a robust minimum on the amount of dark matter in the LMC, which we use to set conservative bounds on the annihilation cross section. The LMC emission is generally very well described by the standard astrophysical sources, with at most a $1-2σ$ excess identified near the kinematic center of the LMC once systematic uncertainties are taken into account. We place competitive bounds on the dark matter annihilation cross section as a function of dark matter particle mass and annihilation channel.
△ Less
Submitted 15 April, 2015; v1 submitted 3 February, 2015;
originally announced February 2015.
-
All about baryons: revisiting SIDM predictions at small halo masses
Authors:
A. Bastidas Fry,
F. Governato,
A. Pontzen,
T. Quinn,
M. Tremmel,
L. Anderson,
H. Menon,
A. M. Brooks,
J. Wadsley
Abstract:
We use cosmological hydrodynamic simulations to consistently compare the assembly of dwarf galaxies in both $Λ$ dominated, Cold (CDM) and Self--Interacting (SIDM) dark matter models. The SIDM model adopts a constant cross section of 2 $cm^{2}/g$, a relatively large value to maximize its effects. These are the first SIDM simulations that are combined with a description of stellar feedback that natu…
▽ More
We use cosmological hydrodynamic simulations to consistently compare the assembly of dwarf galaxies in both $Λ$ dominated, Cold (CDM) and Self--Interacting (SIDM) dark matter models. The SIDM model adopts a constant cross section of 2 $cm^{2}/g$, a relatively large value to maximize its effects. These are the first SIDM simulations that are combined with a description of stellar feedback that naturally drives potential fluctuations able to create dark matter cores. Remarkably, SIDM fails to significantly lower the central dark matter density at halo peak velocities V$_{max}$ $<$ 30 Km/s. This is due to the fact that the central regions of very low--mass field halos have relatively low central velocity dispersion and densities, leading to time scales for SIDM collisions greater than a Hubble time. CDM halos with V$_{max}$ $<$ 30 km/s have inefficient star formation, and hence weak supernova feedback. At a fixed 2 cm2/g SIDM cross section, the DM content of very low mass CDM and SIDM halos differs by no more than a factor of two within 100-200pc. At larger halo masses ($\sim$ 10$^{10}$ solar masses), the introduction of baryonic processes creates field dwarf galaxies with dark matter cores and central DM$+$baryon distributions that are effectively indistinguishable between CDM and SIDM. Both models are in broad agreement with observed Local Group field galaxies across the range of masses explored. To significantly differentiate SIDM from CDM at the scale of faint dwarf galaxies, a velocity dependent cross section that rapidly increases to values larger than 2 $cm^{2}/g$ for halos with V$_{max}$ < 25-30 Km/s needs to be introduced.
△ Less
Submitted 16 June, 2015; v1 submitted 2 January, 2015;
originally announced January 2015.
-
Delayed Star Formation in Isolated Dwarf Galaxies: HST Star Formation History of the Aquarius Dwarf Irregular
Authors:
Andrew A. Cole,
Daniel R. Weisz,
Andrew E. Dolphin,
Evan D. Skillman,
Alan W. McConnachie,
Alyson M. Brooks,
Ryan Leaman
Abstract:
We have obtained deep images of the highly isolated (d = 1 Mpc) Aquarius dwarf irregular galaxy (DDO 210) with the Hubble Space Telescope (HST) Advanced Camera for Surveys (ACS). The resulting color-magnitude diagram (CMD) reaches more than a magnitude below the oldest main-sequence turnoff, allowing us to derive the star formation history (SFH) over the entire lifetime of the galaxy with a timing…
▽ More
We have obtained deep images of the highly isolated (d = 1 Mpc) Aquarius dwarf irregular galaxy (DDO 210) with the Hubble Space Telescope (HST) Advanced Camera for Surveys (ACS). The resulting color-magnitude diagram (CMD) reaches more than a magnitude below the oldest main-sequence turnoff, allowing us to derive the star formation history (SFH) over the entire lifetime of the galaxy with a timing precision of ~10% of the lookback time. Using a maximum likelihood fit to the CMD we find that only ~10% of all star formation in Aquarius took place more than 10 Gyr ago (lookback time equivalent to redshift z ~2). The star formation rate increased dramatically ~6-8 Gyr ago (z ~ 0.7-1.1) and then declined until the present time. The only known galaxy with a more extreme confirmed delay in star formation is Leo A, a galaxy of similar M(HI)/M(stellar), dynamical mass, mean metallicity, and degree of isolation. The delayed stellar mass growth in these galaxies does not track the mean dark matter accretion rate from CDM simulations. The similarities between Leo A and Aquarius suggest that if gas is not removed from dwarf galaxies by interactions or feedback, it can linger for several gigayears without cooling in sufficient quantity to form stars efficiently. We discuss possible causes for the delay in star formation including suppression by reionization and late-time mergers. We find reasonable agreement between our measured SFHs and select cosmological simulations of isolated dwarfs. Because star formation and merger processes are both stochastic in nature, delayed star formation in various degees is predicted to be a characteristic (but not a universal) feature of isolated small galaxies.
△ Less
Submitted 4 September, 2014;
originally announced September 2014.
-
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.
-
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.
-
The SAMI Galaxy Survey: instrument specification and target selection
Authors:
J. J. Bryant,
M. S. Owers,
A. S. G. Robotham,
S. M. Croom,
S. P. Driver,
M. J. Drinkwater,
N. P. F. Lorente,
L. Cortese,
N. Scott,
M. Colless,
A. Schaefer,
E. N. Taylor,
I. S. Konstantopoulos,
J. T. Allen,
I. Baldry,
L. Barnes,
A. E. Bauer,
J. Bland-Hawthorn,
J. V. Bloom,
A. M. Brooks,
S. Brough,
G. Cecil,
W. Couch,
D. Croton,
R. Davies
, et al. (32 additional authors not shown)
Abstract:
The SAMI Galaxy Survey will observe 3400 galaxies with the Sydney-AAO Multi-object Integral-field spectrograph (SAMI) on the Anglo-Australian Telescope (AAT) in a 3-year survey which began in 2013. We present the throughput of the SAMI system, the science basis and specifications for the target selection, the survey observation plan and the combined properties of the selected galaxies. The survey…
▽ More
The SAMI Galaxy Survey will observe 3400 galaxies with the Sydney-AAO Multi-object Integral-field spectrograph (SAMI) on the Anglo-Australian Telescope (AAT) in a 3-year survey which began in 2013. We present the throughput of the SAMI system, the science basis and specifications for the target selection, the survey observation plan and the combined properties of the selected galaxies. The survey includes four volume limited galaxy samples based on cuts in a proxy for stellar mass, along with low-stellar mass dwarf galaxies all selected from the Galaxy And Mass Assembly (GAMA) survey. The GAMA regions were selected because of the vast array of ancillary data available, including ultraviolet through to radio bands. These fields are on the celestial equator at 9, 12, and 14.5 hours, and cover a total of 144 square degrees (in GAMA-I). Higher density environments are also included with the addition of eight clusters. The clusters have spectroscopy from 2dFGRS and SDSS and photometry in regions covered by the Sloan Digital Sky Survey (SDSS) and/or VLT Survey Telescope/ATLAS. The aim is to cover a broad range in stellar mass and environment, and therefore the primary survey targets cover redshifts 0.004 < z < 0.095, magnitudes r$_{pet}$ < 19.4, stellar masses $10^{7} - 10^{12}$ M$_{sol}$, and environments from isolated field galaxies through groups to clusters of $10^{15}$ M$_{sol}$.
△ Less
Submitted 17 February, 2015; v1 submitted 28 July, 2014;
originally announced July 2014.
-
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.
-
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.
-
The Comparative Chemical Evolution of an Isolated Dwarf Galaxy: A VLT and Keck Spectroscopic Survey of WLM
Authors:
Ryan Leaman,
Kim A. Venn,
Alyson M. Brooks,
Giuseppina Battaglia,
Andrew A. Cole,
Rodrigo A. Ibata,
Mike J. Irwin,
Alan W. McConnachie,
J. Trevor Mendel,
Else Starkenburg,
Eline Tolstoy
Abstract:
Building on our previous spectroscopic and photometric analysis of the isolated Local Group dwarf irregular (dIrr) galaxy WLM, we present a comparison of the metallicities of its RGB stars with respect to the well studied Local Group dwarf spheroidal galaxies (dSphs) and Magellanic Clouds. We calculate a mean metallicity of [Fe/H]$ = -1.28 \pm 0.02$, and intrinsic spread in metallicity of…
▽ More
Building on our previous spectroscopic and photometric analysis of the isolated Local Group dwarf irregular (dIrr) galaxy WLM, we present a comparison of the metallicities of its RGB stars with respect to the well studied Local Group dwarf spheroidal galaxies (dSphs) and Magellanic Clouds. We calculate a mean metallicity of [Fe/H]$ = -1.28 \pm 0.02$, and intrinsic spread in metallicity of $σ= 0.38 \pm 0.04$ dex, similar to the mean and spread observed in the massive dSph Fornax and the Small Magellanic Cloud. Thus, despite its isolated environment the global metallicity still follows expectations for WLM's mass and its global chemical evolution is similar to other nearby luminous dwarf galaxies (gas-rich or gas-poor). The data also show a radial gradient in [Fe/H] of $\rm{d[Fe/H]/dr_{c}} = -0.04 \pm 0.04$ dex $\rm{r_{c}^{-1}}$, which is flatter than that seen in the unbiased and spatially extended surveys of dSphs. Comparison of the spatial distribution of [Fe/H] in WLM, the Magellanic Clouds, and a sample of Local Group dSphs, shows an apparent dichotomy in the sense that the dIrrs have statistically flatter radial [Fe/H] gradients than the low angular momentum dSphs. The correlation between angular momentum and radial metallicity gradient is further supported when considering the Local Group dEs. This chemodynamic relationship offers a new and useful constraint for environment driven dwarf galaxy evolution models in the Local Group.
△ Less
Submitted 7 February, 2013;
originally announced February 2013.
-
Angular Momentum Acquisition in Galaxy Halos
Authors:
Kyle R. Stewart,
Alyson M. Brooks,
James S. Bullock,
Ariyeh H. Maller,
Juerg Diemand,
James Wadsley,
Leonidas A. Moustakas
Abstract:
We use high-resolution cosmological hydrodynamic simulations to study the angular momentum acquisition of gaseous halos around Milky Way sized galaxies. We find that cold mode accreted gas enters a galaxy halo with ~70% more specific angular momentum than dark matter averaged over cosmic time (though with a very large dispersion). In fact, we find that all matter has a higher spin parameter when m…
▽ More
We use high-resolution cosmological hydrodynamic simulations to study the angular momentum acquisition of gaseous halos around Milky Way sized galaxies. We find that cold mode accreted gas enters a galaxy halo with ~70% more specific angular momentum than dark matter averaged over cosmic time (though with a very large dispersion). In fact, we find that all matter has a higher spin parameter when measured at accretion than when averaged over the entire halo lifetime, and is well characterized by λ~0.1, at accretion. Combined with the fact that cold flow gas spends a relatively short time (1-2 dynamical times) in the halo before sinking to the center, this naturally explains why cold flow halo gas has a specific angular momentum much higher than that of the halo and often forms "cold flow disks". We demonstrate that the higher angular momentum of cold flow gas is related to the fact that it tends to be accreted along filaments.
△ Less
Submitted 25 July, 2013; v1 submitted 14 January, 2013;
originally announced January 2013.
-
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.
-
A Baryonic Solution to the Missing Satellites Problem
Authors:
Alyson M. Brooks,
Michael Kuhlen,
Adi Zolotov,
Dan Hooper
Abstract:
It has been demonstrated that the inclusion of baryonic physics can alter the dark matter densities in the centers of low-mass galaxies, making the central dark matter slope more shallow than predicted in pure cold dark matter simulations. This flattening of the dark matter profile can occur in the most luminous subhalos around Milky Way-mass galaxies. Zolotov et al. (2012) have suggested a correc…
▽ More
It has been demonstrated that the inclusion of baryonic physics can alter the dark matter densities in the centers of low-mass galaxies, making the central dark matter slope more shallow than predicted in pure cold dark matter simulations. This flattening of the dark matter profile can occur in the most luminous subhalos around Milky Way-mass galaxies. Zolotov et al. (2012) have suggested a correction to be applied to the central masses of dark matter-only satellites in order to mimic the affect of (1) the flattening of the dark matter cusp due to supernova feedback in luminous satellites, and (2) enhanced tidal stripping due to the presence of a baryonic disk. In this paper, we apply this correction to the z=0 subhalo masses from the high resolution, dark matter-only Via Lactea II (VL2) simulation, and find that the number of massive subhalos is dramatically reduced. After adopting a stellar mass to halo mass relationship for the VL2 halos, and identifying subhalos that are (1) likely to be destroyed by stripping and (2) likely to have star formation suppressed by photo-heating, we find that the number of massive, luminous satellites around a Milky Way-mass galaxy is in agreement with the number of observed satellites around the Milky Way or M31. We conclude that baryonic processes have the potential to solve the missing satellites problem.
△ Less
Submitted 29 April, 2013; v1 submitted 24 September, 2012;
originally announced September 2012.
-
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.
-
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.
-
Why Baryons Matter: The Kinematics of Dwarf Spheroidal Satellites
Authors:
Alyson M. Brooks,
Adi Zolotov
Abstract:
We use high resolution cosmological simulations of Milky Way-mass galaxies that include both baryons and dark matter to show that baryonic physics (energetic feedback from supernovae and subsequent tidal stripping) significantly reduces the dark matter mass in the central regions of luminous satellite galaxies. The reduced central masses of the simulated satellites reproduce the observed internal…
▽ More
We use high resolution cosmological simulations of Milky Way-mass galaxies that include both baryons and dark matter to show that baryonic physics (energetic feedback from supernovae and subsequent tidal stripping) significantly reduces the dark matter mass in the central regions of luminous satellite galaxies. The reduced central masses of the simulated satellites reproduce the observed internal dynamics of Milky Way and M31 satellites as a function of luminosity. We use these realistic satellites to update predictions for the observed velocity and luminosity functions of satellites around Milky Way-mass galaxies when baryonic effects are accounted for. We also predict that field dwarf galaxies in the same luminosity range as the Milky Way classical satellites should not exhibit velocities as low as the satellites, since the field dwarfs do not experience tidal stripping. Additionally, the early formation times of the satellites compared to field galaxies at the same luminosity may be apparent in the star formation histories of the two populations. Including baryonic physics in Cold Dark Matter models naturally explains the observed low dark matter densities in the Milky Way's dwarf spheroidal population. Our simulations therefore resolve the tension between kinematics predicted in Cold Dark Matter theory and observations of satellites, without invoking alternative forms of dark matter.
△ Less
Submitted 28 April, 2014; v1 submitted 10 July, 2012;
originally announced July 2012.