Search

The shape and orientation of dark matter (DM) halos are sensitive to the micro-physics of the DM particle, yet in many mass models, the symmetry axes of the Milky Way's DM halo are often assumed to be aligned with the symmetry axes of the stellar disk. This is well-motivated for the inner DM halo but not for the outer halo. We use zoomed cosmological-baryonic simulations from the Latte suite of FIRE-2 Milky Way-mass galaxies to explore the evolution of the DM halo's orientation with radius and time, with or without a major merger with a Large Magellanic Cloud (LMC) analog, and when varying the DM model. In three of the four CDM halos we examine, the orientation of the halo minor axis diverges from the stellar disk vector by more than 20 degrees beyond about 30 galactocentric kpc, reaching a maximum of 30--90 degrees depending on the individual halo's formation history. In identical simulations using a model of self-interacting DM with $\sigma = 1 \, \mathrm{cm}^2 \, \mathrm{g}^{-1}$, the halo remains aligned with the stellar disk out to $\sim$200--400 kpc. Interactions with massive satellites ($M \gtrsim 4 \times 10^{10} \, \rm{M_\odot}$ at pericenter; $M \gtrsim 3.3 \times 10^{10} \, \rm{M_\odot}$ at infall) affect the orientation of the halo significantly, aligning the halo's major axis with the satellite galaxy from the disk to the virial radius. The relative orientation of the halo and disk beyond 30 kpc is a potential diagnostic of SIDM if the effects of massive satellites can be accounted for., Comment: Submitted to ApJ, 19 pages, 12 figures, 3 tables

We study how supersonic streaming velocities of baryons relative to dark matter -- a large-scale effect imprinted at recombination and coherent over $\sim 3$ Mpc scales -- affects the formation of dwarf galaxies at $z \gtrsim 5$. We perform cosmological hydrodynamic simulations, including and excluding streaming velocities, in regions centered on halos with $M_{\rm vir}(z=0) \approx 10^{10}$ M$_{\odot}$; the simulations are part of the Feedback In Realistic Environments (FIRE) project and run with FIRE-3 physics. Our simulations comprise many thousands of systems with halo masses between $M_{\rm vir} = 2\times10^{5}$ M$_{\odot}$ and $2\times10^9$ M$_{\odot}$ in the redshift range $z=20-5$. A few hundred of these galaxies form stars and have stellar masses ranging from 100 to $10^7$ M$_{\odot}$. While star formation is globally delayed by approximately 50 Myr in the streaming relative to non-streaming simulations and the number of luminous galaxies is correspondingly suppressed at high redshift in the streaming runs, these effects decay with time. By $z=5$, the properties of the simulated galaxies are nearly identical in the streaming versus non-streaming runs, indicating that any effects of streaming velocities on the properties of galaxies at the mass scale of classical dwarfs and larger do not persist to $z=0$., Comment: 12 pages, 8 figures, submitted to ApJ

We perform cosmological hydrodynamical simulations to study the formation of proto-globular cluster candidates in progenitors of present-day dwarf galaxies $(M_{\rm vir} \approx 10^{10}\, {\rm M}_\odot$ at $z=0$) as part of the "Feedback in Realistic Environment" (FIRE) project. Compact ($r_{1/2}<30$ pc), relatively massive ($0.5 \times 10^5 \lesssim M_{\star}/{\rm M}_\odot \lesssim 5\times10^5$), self-bound stellar clusters form at $11\gtrsim z \gtrsim 5$ in progenitors with $M_{\rm vir} \approx 10^9\,{\rm M}_\odot$. Cluster formation is triggered when at least $10^7\,{\rm M}_\odot$ of dense, turbulent gas reaches $\Sigma_{\rm gas} \approx 10^4\, {\rm M}_\odot\, {\rm pc}^{-2}$ as a result of the compressive effects of supernova feedback or from cloud-cloud collisions. The clusters can survive for $2-3\,{\rm Gyr}$; absent numerical effects, they would likely survive substantially longer, perhaps to $z=0$. The longest-lived clusters are those that form at significant distance -- several hundreds of pc -- from their host galaxy. We therefore predict that globular clusters forming in progenitors of present-day dwarf galaxies will be offset from any pre-existing stars within their host dark matter halos as opposed to deeply embedded within a well-defined galaxy. Properties of the nascent clusters are consistent with observations of some of the faintest and most compact high-redshift sources in \textit{Hubble Space Telescope} lensing fields and are at the edge of what will be detectable as point sources in deep imaging of non-lensed fields with the \textit{James Webb Space Telescope}. By contrast, the star clusters' host galaxies will remain undetectable., Comment: 14 pages, 5 figures, submitted to MNRAS

Self-interacting dark matter (SIDM) models have received great attention over the past decade as solutions to the small-scale puzzles of astrophysics. Though there are different implementations of dark matter (DM) self-interactions in N-body codes of structure formation, there has not been a systematic study to compare the predictions of these different implementations. We investigate the implementation of dark matter self-interactions in two simulation codes: Gizmo and Arepo. We begin with identical initial conditions for an isolated $10^{10}$ M$_\odot$ dark matter halo and investigate the evolution of the density and velocity dispersion profiles in Gizmo and Arepo for SIDM cross-section over mass of 1, 5, and 50 $\rm cm^2 g^{-1}$. Our tests are restricted to the core expansion phase where the core density decreases and core radius increases with time. We find better than 30% agreement between the codes for the density profile in this phase of evolution, with the agreement improving at higher resolution. We find that varying code-specific SIDM parameters changes the central halo density by less than 10% outside of the convergence radius. We argue that SIDM core formation is robust across the two different schemes and conclude that these codes can reliably differentiate between cross-sections of 1, 5, and 50 $\rm cm^2 g^{-1}$ but finer distinctions would require further investigation.

Increasingly, uncertainties in predictions from galaxy formation simulations (at sub-Milky Way masses) are dominated by uncertainties in stellar evolution inputs. In this paper, we present the full set of updates from the FIRE-2 version of the Feedback In Realistic Environments (FIRE) project code, to the next version, FIRE-3. While the transition from FIRE-1 to FIRE-2 focused on improving numerical methods, here we update the stellar evolution tracks used to determine stellar feedback inputs, e.g. stellar mass-loss (O/B and AGB), spectra (luminosities and ionization rates), and supernova rates (core-collapse and Ia), as well as detailed mass-dependent yields. We also update the low-temperature cooling and chemistry, to enable improved accuracy at $T \lesssim 10^{4}\,$K and densities $n\gg 1\,{\rm cm^{-3}}$, and the meta-galactic ionizing background. All of these synthesize newer empirical constraints on these quantities and updated stellar evolution and yield models from a number of groups, addressing different aspects of stellar evolution. To make the updated models as accessible as possible, we provide fitting functions for all of the relevant updated tracks, yields, etc, in a form specifically designed so they can be directly 'plugged in' to existing galaxy formation simulations. We also summarize the default FIRE-3 implementations of 'optional' physics, including spectrally-resolved cosmic rays and supermassive black hole growth and feedback., Comment: 26 pages, 12 figures. Submitted to MNRAS. Comments welcome

We describe a public data release of the FIRE-2 cosmological zoom-in simulations of galaxy formation, available at http://flathub.flatironinstitute.org/fire, from the Feedback In Realistic Environments (FIRE) project. FIRE-2 simulations achieve parsec-scale resolution to explicitly model the multi-phase interstellar medium while implementing direct models for stellar evolution and feedback, including stellar winds, core-collapse and Ia supernovae, radiation pressure, photoionization, and photoelectric heating. We release complete snapshots from 3 suites of simulations. The first comprises 20 simulations that zoom in on 14 Milky Way-mass galaxies, 5 SMC/LMC-mass galaxies, and 4 lower-mass galaxies including 1 ultra-faint; we release 39 snapshots across z = 0 - 10. The second comprises 4 massive galaxies, with 19 snapshots across z = 1 - 10. Finally, a high-redshift suite comprises 22 simulations, with 11 snapshots across z = 5 - 10. Each simulation also includes dozens of resolved lower-mass (satellite) galaxies in its zoom-in region. Snapshots include all stored properties for all dark matter, gas, and star particles, including 11 elemental abundances for stars and gas, and formation times (ages) of star particles. We also release accompanying (sub)halo catalogs, which include galaxy properties and member star particles. For the simulations to z = 0, including all Milky Way-mass galaxies, we release the formation coordinates and an "ex-situ" flag for all star particles, pointers to track particles across snapshots, catalogs of stellar streams, and multipole basis expansions for the halo mass distributions. We describe publicly available python packages for reading and analyzing these simulations., Comment: 22 pages. Accepted for publication in ApJS, matches the published version, with updated references. Now includes up to 39 snapshots per simulation and formation coordinates for all star particles. Data available at http://flathub.flatironinstitute.org/fire