Category: Uncategorized

McKeown et al., available on arXiv

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

Patel et al., available on arXiv

Abstract: We investigate stellar elemental abundance patterns at z = 0 in 8 low-mass (M_* = 10^6 – 10^9 M_sun) galaxies in the Feedback in Realistic Environments (FIRE-2) cosmological simulations. Using magnesium (Mg) as a representative alpha-element, we explore stellar abundance patterns in [Mg/Fe] versus [Fe/H], which follow an overall monotonic trend that evolved slowly over time. Beyond this, we explore 3 notable secondary features in enrichment (found in three different case-study galaxies) that arise from a galaxy merger or bursty star formation. First, we observe a secondary track with a lower [Mg/Fe] than the main trend. At z = 0, stars from this track are predominantly found within 2-6 kpc of the center; they were accreted in a 1:3 total-mass-ratio merger ~ 0.4 Gyr ago. Second, we find a distinct elemental bi-modality that forms following a strong burst in star formation in a galaxy at t_lookback ~ 10 Gyr. This burst quenched star formation for ~ 0.66 Gyr, allowing Ia supernovae to enrich the system with iron before star formation resumed. Third, we examine stripes in enrichment that run roughly orthogonal to the dominant [Mg/Fe] versus [Fe/H] trend; these stripes correspond to short bursts of star formation during which core-collapse supernovae enrich the surrounding medium with Mg (and Fe) on short timescales. If observed, these features would substantiate the utility of elemental abundances in revealing the assembly and star formation histories of dwarf galaxies. We explore the observability of these features for upcoming spectroscopic studies. Our results show that precise measurements of elemental abundance patterns can reveal critical events in the formation histories of low-mass galaxies.

Chan et al., available on arXiv

Abstract: Cosmic rays (CRs) are an important component in the interstellar medium (ISM), but their effect on the dynamics of the disk-halo interface (< 10 kpc from the disk) is still unclear. We study the influence of CRs on the gas above the disk with high-resolution FIRE-2 cosmological simulations of late-type Lstar galaxies at redshift around zero. We compare runs with and without CR feedback (with constant anisotropic diffusion around 3e29 cm^2/s and streaming). Our simulations capture the relevant disk halo interactions, including outflows, inflows, and galactic fountains. Extra-planar gas in all of the runs satisfies dynamical balance, where total pressure balances the weight of the overlying gas. While the kinetic pressure from non-uniform motion (>1-kpc scale) dominates in the midplane, thermal and bulk pressures (or CR pressure if included) take over at large heights. We find that with CR feedback, (1) the warm (1e4K) gas is slowly accelerated by CRs; (2) the hot (> 5e5K) gas scale height is suppressed; (3) the warm-hot (2e4-5e5K) medium becomes the most volume-filling phase in the disk-halo interface. We develop a novel conceptual model of the near-disk gas dynamics in low-redshift Lstar galaxies: With CRs, the disk-halo interface is filled with CR-driven warm winds and hot super-bubbles that are propagating into the CGM with a small fraction falling back to the disk. Without CRs, most outflows from hot superbubbles are trapped by the existing hot halo and gravity, so typically they form galactic fountains.

Cunningham et al., available on arXiv

Abstract: In the era of large-scale spectroscopic surveys in the Local Group (LG), we can explore using chemical abundances of halo stars to study the star formation and chemical enrichment histories of the dwarf galaxy progenitors of the Milky Way (MW) and M31 stellar halos. In this paper, we investigate using the Chemical Abundance Ratio Distributions (CARDs) of seven stellar halos from the Latte suite of FIRE-2 simulations. We attempt to infer galaxies’ assembly histories by modelling the CARDs of the stellar halos of the Latte galaxies as a linear combination of template CARDs from disrupted dwarfs, with different stellar masses M* and quenching times t100. We present a method for constructing these templates using present-day dwarf galaxies. For four of the seven Latte halos studied in this work, we recover the mass spectrum of accreted dwarfs to a precision of <10%. For the fraction of mass accreted as a function of t100, we find residuals of 20−30% for five of the seven simulations. We discuss the failure modes of this method, which arise from the diversity of star formation and chemical enrichment histories dwarf galaxies can take. These failure cases can be robustly identified by the high model residuals. Though the CARDs modeling method does not successfully infer the assembly histories in these cases, the CARDs of these disrupted dwarfs contain signatures of their unusual formation histories. Our results are promising for using CARDs to learn more about the histories of the progenitors of the MW and M31 stellar halos.

Hopkins et al., available on arXiv

Abstract: We present the first simulations evolving resolved spectra of cosmic rays (CRs) from MeV-TeV energies (including electrons, positrons, (anti)protons, and heavier nuclei), in live kinetic-MHD galaxy simulations with star formation and feedback. We utilize new numerical methods including terms often neglected in historical models, comparing Milky Way analogues with phenomenological scattering coefficients $\nu$ to Solar-neighborhood (LISM) observations (spectra, B/C, e^+/e^-, p^-/p, 10^Be/9^Be, ionization). We show it is possible to reproduce observations with simple single-power-law injection and scattering coefficients (scaling with rigidity R), similar to previous (non-dynamical) calculations. We also find: (1) The circum-galactic medium in realistic galaxies necessarily imposes a ~10 kpc CR scattering halo, influencing the required nu(R). (2) Increasing the normalization of nu(R) re-normalizes CR secondary spectra but also changes primary spectral slopes, owing to source distribution and loss effects. (3) Diffusive/turbulent reacceleration is unimportant and generally sub-dominant to gyroresonant/streaming losses, which are sub-dominant to adiabatic/convective terms dominated by ~0.1-1 kpc turbulent/fountain motions. (4) CR spectra vary considerably across galaxies; certain features can arise from local structure rather than transport physics. (5) Systematic variation in CR ionization rates between LISM and molecular clouds (or Galactic position) arises naturally without invoking alternative sources. (6) Abundances of CNO nuclei require most CR acceleration occurs around when reverse shocks form in SNe, not in OB wind bubbles or later Sedov-Taylor stages of SNe remnants.

Kado-Fong et al., available on arXiv

Abstract: Extended, old, and round stellar halos appear to be ubiquitous around high-mass dwarf galaxies (10^8.5 < Mstar/Msun < 10^9.6) in the observed universe. However, it is unlikely that these dwarfs have undergone a sufficient number of minor mergers to form stellar halos that are composed of predominantly accreted stars. Here, we demonstrate that FIRE-2 (Feedback in Realistic Environments) cosmological zoom-in simulations are capable of producing dwarf galaxies with realistic structure, including both a thick disk and round stellar halo. Crucially, these stellar halos are formed in-situ, largely via the outward migration of disk stars. However, there also exists a large population of "non-disky" dwarfs in FIRE that lack a well-defined disk/halo and do not resemble the observed dwarf population. These non-disky dwarfs tend to be either more gas poor or to have burstier recent star formation histories than the disky dwarfs, suggesting that star formation feedback may be preventing disk formation. Both classes of dwarfs underscore the power of a galaxy's intrinsic shape -- which is a direct quantification of the distribution of the galaxy's stellar content -- to interrogate the feedback implementation in simulated galaxies.

Stern et al., available on arXiv

Abstract: Recent searches for the hosts of high-redshift (z~4) damped Ly-alpha absorbers (DLAs) have detected bright galaxies at distances of tens of kpc from the DLA. Using the FIRE-2 cosmological zoom simulations, we argue that these relatively large distances are due to a predominantly cool and neutral inner circumgalactic medium (CGM) surrounding high-redshift galaxies. The inner CGM is cool because of the short cooling time of hot gas in <~10^12 Msun halos, which implies that accretion and feedback energy are radiated quickly, while it is neutral due to the high volume densities and column densities at high redshift which shield cool gas from photoionization. Our analysis predicts large DLA covering factors (>50%) out to impact parameters ∼0.3((1+z)/5)^3/2 Rvir from the central galaxies at z>1, equivalent to a physical distance of ∼21 M12^(1/3)*((1+z)/5)^1/2 kpc (Rvir and M12 are the halo virial radius and mass in units of 10^12 Msun, respectively). This implies that DLA covering factors at z~4 may be comparable to unity out to a distance ~10 times larger than stellar half-mass radii. A predominantly neutral inner CGM in the early universe suggests that its mass and metallicity can be directly constrained by CGM absorption surveys, without resorting to large ionization corrections as required for ionized CGM.

Trapp et al., available on arXiv

Abstract: Observations indicate that a continuous supply of gas is needed to maintain observed star formation rates in large, disky galaxies. To fuel star formation, gas must reach the inner regions of such galaxies. Despite its crucial importance for galaxy evolution, how and where gas joins galaxies is poorly constrained observationally and is rarely explored in fully cosmological simulations. To investigate gas accretion in the vicinity of galaxies, we analyze the FIRE-2 cosmological zoom-in simulations for 4 Milky Way mass galaxies (M_halo ~ 10E12 solar masses), focusing on simulations with cosmic ray physics. We find that at z~0, gas approaches the disk with angular momentum similar to the gaseous disk edge and low radial velocities, piling-up near the edge and settling into full rotational support. Accreting gas moves predominantly parallel to the disk with small but nonzero vertical velocity components, and joins the disk largely in the outskirts as opposed to “raining” down onto the disk. Once in the disk, gas trajectories are complex, being dominated by spiral arm induced oscillations and feedback. However, time and azimuthal averages show clear but slow net radial infall with transport speeds of 1-3 km/s and net mass fluxes through the disk on the order of one solar mass per year, comparable to the star formation rates of the galaxies and decreasing towards galactic center as gas is sunk into star formation. These rates are slightly higher in simulations without cosmic rays (1-7 km/s, ~4-5 solar masses per year). We find overall consistency of our results with observational constraints and discuss prospects of future observations of gas flows in and around galaxies.

Panithanpaisal et al., available on arXiv

Abstract: Stellar streams record the accretion history of their host galaxy. We present a set of simulated streams from disrupted dwarf galaxies in 13 cosmological simulations of Milky Way (MW)-mass galaxies from the FIRE-2 suite at z=0, including 7 isolated Milky Way-mass systems and 6 hosts resembling the MW-M31 pair (full dataset at: this https URL). In total, we identify 106 simulated stellar streams, with no significant differences in the number of streams and masses of their progenitors between the isolated and paired environments. We resolve simulated streams with stellar masses ranging from ~5*10^5 up to ~10^9 Msun, similar to the mass range between the Orphan and Sagittarius streams in the MW. We confirm that present-day simulated satellite galaxies are good proxies for stellar stream progenitors, with similar properties including their stellar mass function, velocity dispersion, [Fe/H] and [alpha/H] evolution tracks, and orbital distribution with respect to the galactic disk plane. Each progenitor’s lifetime is marked by several important timescales: its infall, star-formation quenching, and stream-formation times. We show that the ordering of these timescales is different between progenitors with stellar masses higher and lower than ~2*10^6 Msun. Finally, we show that the main factor controlling the rate of phase-mixing, and therefore fading, of tidal streams from satellite galaxies in MW-mass hosts is non-adiabatic evolution of the host potential. Other factors commonly used to predict phase-mixing timescales, such as progenitor mass and orbital circularity, show virtually no correlation with the number of dynamical times required for a stream to become phase-mixed.

Nikakhtar et al., available on arXiv

Abstract: The standard picture of galaxy formation motivates the decomposition of the Milky Way into 3–4 stellar populations with distinct kinematic and elemental abundance distributions: the thin disk, thick disk, bulge, and stellar halo. To test this idea, we construct a Gaussian mixture model (GMM) for both simulated and observed stars in the Solar neighborhood, using measured velocities and iron abundances (i.e., an augmented Toomre diagram) as the distributions to be decomposed. We compare results for the Gaia-APOGEE DR16 crossmatch catalog of the Solar neighborhood with those from a suite of synthetic Gaia-APOGEE crossmatches constructed from FIRE-2 cosmological simulations of Milky Way-mass galaxies. We find that in both the synthetic and real data, the best-fit GMM uses five independent components, some of whose properties resemble the standard populations predicted by galaxy formation theory. Two components can be identified unambiguously as the thin disk and another as the halo. However, instead of a single counterpart to the thick disk, there are three intermediate components with different age and alpha abundance distributions (although these data are not used to construct the model). We use decompositions of the synthetic data to show that the classified components indeed correspond to stars with different origins. By analogy with the simulated data, we show that our mixture model of the real Gaia-APOGEE crossmatch distinguishes the following components: (1) a classic thin disk of young stars on circular orbits (46%), (2) thin disk stars heated by interactions with satellites (22%), (3, 4) two components representing the velocity asymmetry of the alpha-enhanced thick disk (27%), and (5) a stellar halo consistent with early, massive accretion (4%).