The MARVELous Dwarfs meet the : Constraints on Dwarf Galaxies using a heroically large simulated sample Ferah Munshi VIDA Fellow, Vanderbilt University (—> Assistant Prof, OU this fall)

In collaboration with: Alyson Brooks, Dan Weisz, Jillian Bellovary, Kelly Holley-Bockelmann , Elaad Applebaum, Charlotte Christensen+ UW N-body Shop*

Email: [email protected] The Abundance of Dwarf Galaxies

Ferah Munshi VIDA Fellow, Vanderbilt University (—> Assistant Prof, OU this fall)

In collaboration with: Alyson Brooks, Dan Weisz, Jillian Bellovary, Kelly Holley-Bockelmann , Elaad Applebaum, Charlotte Christensen + UW N-body Shop*

Email: [email protected] I use cosmological N- body + SPH galaxy simulations to figure out how galaxies form and evolve- i.e., how they come to look as they do today. What is a (N-body) Simulation?

Modeling a dynamical system of particles, usually under the influence of physical forces, in this case: gravity For me: stars + dark matter, acting under the influence of gravity, within a galaxy What is an N-body + SPH Simulation? • SPH= “smoothed particle hydrodynamics” • computational method used for simulating fluid flows- ie, gas • Gas is divided into a set of discrete elements, referred to as “particles” • “cosmological”= from early times all the way to present day Galaxies are made up of stars, gas and dark matter (the majority of a galaxy is in dark matter)

FEEDBACK can imprint its affects on all three components

In dwarf galaxies, feedback is key in understanding the dark matter profiles. (α)

FEEDBACK

Figure courtesy of F. Governato Galaxies are made up of stars, gas and dark matter (the majority of a galaxy is in dark matter)

FEEDBACK can imprint its affects on all three components

In dwarf galaxies, feedback is key in understanding the dark matter profiles.

How does feedback imprint itself on the other components? All feedback mechanisms have this in common:

They heat gas, drive outflows, and suppress star formation

In order to simulate a galaxy, you must be able to model feedback. Feedback is necessary to form realistic* galaxies. *realistic= look like observed galaxies in basic properties

Stellar e.g. winds from massive stars

Feedback Supernova

Black Hole e.g. AGN feedback

Depending on mass of galaxy, different sources have varying importance So how do you know you’re modeling feedback correctly? Compare to observations!

Big question #1: How do feedback and star formation affect the stellar to mass relationship (SMHM)? What is the abundance and scatter of low mass galaxies?

How do they populate their dark matter halos? MARVELous Dwarf Volumes Captain Marvel • run with ChaNGa a.k.a 40 Thieves • mgas=1.4e3 Msun, mstar=400Msun, mdark=6e3Msun 3 • Effective resolution (4096) COMPLETE TOTAL # RESOLVED DWARFS = 64

COMPLETE Rogue Storm Elektra

COMPLETE COMPLETE

13 MARVELous Dwarfs Cpt Marvel run encapsulating multiple subgrid models- Captain Marvel 1. High density threshold SF (MC run) [formerly 40 Thieves] 2. H2 based SF (H2 run) [formerly 40 Thieves] 3. H2 based SF + “SM”BHs (BH COMPLETE implementation from Tremmel+ 2015) TOTAL # RESOLVED DWARFS = 64 COMPLETE Rogue Storm Elektra

COMPLETE COMPLETE 14 Justice League Dwarfs • run with ChaNGa • 4 volumes containing a • mstar=3.9e3 Msun, Milky-Way sized halo- mgas=8.1e3Msun, Sandra, Ruth, Sonia,

mdark=1.3e5Msun Elena • Effective resolution (3072)3 TOTAL # RESOLVED DWARFS = 101

15 MARVELous Dwarf Volumes + Justice League Dwarfs = 165 High-resolution simulated dwarfs Charm Nbody GrAvity solver • Massively parallel SPH (smoothed particle hydrodynamics); fully cosmological • SNe feedback creating realistic outflows • SF linked to shielded gas • Optimized SF parameters • NEW SPH implementation • Previous gen code: Gasoline

Menon+ 2014, Governato+ 2014 z=0 DM density

z=0 Gas density Sandra: highest simulated redshift —> present day Mass-Metallicity matches observations Mass-Metallicity matches observations Star Formation consistent with local dwarfs Cumulative SFH

We cover a wide range of SFHs Cumulative SFH

We cover a wide range of SFHs Simulations can predict fraction of halos that remain dark till present day

Only ~20% of 107 solar mass halos actually host a galaxy Munshi, Brooks + submitted25 ~Atomic Cooling Limit

As lower masses are probed, fewer and fewer halos are occupied!

The occupied halos have differently shaped cumulative SFHs!26 Low mass end of SMHM is poorly constrained…

Scatter ~ 0.1 dex

Scatter ~ 0.9 dex

Satellites Centrals

Munshi+ submitted 40 Thieves: vMC Low mass end of SMHM is poorly constrained…

MARVEL + Justice League Dwarfs

H2 based SF SMBH physics Higher SN feedback …regardless of subgrid physics Single halo- single luminosity assumption breaks down at low mass end

29 Dark halos and extremely low mass halos contribute significantly to scatter scatter in stellar mass for a given halo mass scatter in halo mass for a given stellar mass

Garrison-Kimmel et al. (2017) With one star halos ()

With dark halos () Munshi+ 2017 The scatter has observational ramifications for the stellar mass function

In essence, Munshi+ 2017 is a *tool* to populate low mass halos stochastically, in order to make predictions Munshi+ 2017 Two runs, same initial conditions, different SF/ Feedback prescription -predict similar satellite mass functions

However…

Munshi+, in prep They predict vastly differing numbers of luminous satellites

Munshi+, in prep What’s the difference? Conditions of gas where stars are forming (density!)

Munshi+, in prep, Christensen+ 2012 Dif. SF physics predicts different frequencies of dwarf satellites in LG

Newest FIRE results consistent with H2 run- ask Coral

Munshi+ in prep Dif. SF physics predicts different faint end slope of mass function

Low mass end is sensitive to subgrid physics

Need observations to constrain models!

Munshi+ in prep Summary

• Abundance of dwarf galaxies is largely unconstrained- abundance matching breaks down here! • Simulations, like MARVEL + Justice League, can begin to constrain fraction of populated galaxies and the scatter • BUT! the physics of your simulation changes your predictions- we have to be very careful

Munshi+ 2017 Feedback is necessary to form realistic* galaxies. *realistic= look like observed galaxies in basic properties

Stellar e.g. winds from massive stars

Feedback Supernova

Black Hole e.g. AGN feedback

Depending on mass of galaxy, different sources have varying importance Stars are dwarfs which host a MBH

Bellovary, Cleary, Munshi + in prep • MBH seeds form via direct collapse • Probabilistic approach- similar conditions to SF • MBH formation prescription as in Tremmel+ 2015,2016

• MBHs form at high -z • Truncation due to propagation of metals (need pristine gas) • Formation halo mass as expected from models of direct collapse (e.g. Lodato & Natarajan 2006) (virialized halo gas reaches 104K)

Bellovary, Cleary, Munshi + in prep These MBHs are not AGNs; are not distinguishable from other X-ray sources MBHs preferentially form in denser environments

They are not necessarily in the center of their halos!

Bellovary, Cleary, Munshi + in prep Most common MBH merger ratio is dwarf’s MBH + SMBH of the MW progenitor nearby

MBH mergers happen at all z's. LISA will be sensitive to those at high-z. LISA is a way to probe structure formation at high-z through these MBH mergers! Bellovary, Cleary, Munshi + in prep Part 2 Summary

• MBHs in dwarfs preferentially form in higher density environments • MBH mergers at high-z observed by LISA will be a tracer of early structure formation • MBHs are not active (now or in the past)- hard to observe in light Talk Summary • Combination of Justice League dwarfs + MARVELous dwarfs = 165 high-res dwarfs • Can study broad range of Captain Marvel dwarf properties including SMHM, SFHs, MBHs, radial gradients, resolved SFHs Justice League Elektra Rogue

Storm 45