skip to main content


Search for: All records

Award ID contains: 1911111

Note: When clicking on a Digital Object Identifier (DOI) number, you will be taken to an external site maintained by the publisher. Some full text articles may not yet be available without a charge during the embargo (administrative interval).
What is a DOI Number?

Some links on this page may take you to non-federal websites. Their policies may differ from this site.

  1. ABSTRACT

    We use the GRUMPY galaxy formation model based on a suite of zoom-in, high-resolution, dissipationless Λ Cold Dark Matter (ΛCDM) simulations of the Milky Way (MW) sized haloes to examine total matter density within the half-mass radius of stellar distribution, ρtot(< r1/2), of satellite dwarf galaxies around the MW hosts and their mass assembly histories. We compare model results to ρtot(< r1/2) estimates for observed dwarf satellites of the Milky Way spanning their entire luminosity range. We show that observed MW dwarf satellites exhibit a trend of decreasing total matter density within a half-mass radius, ρtot(< r1/2), with increasing stellar mass. This trend is in general agreement with the trend predicted by the model. None of the observed satellites are overly dense compared to the results of our ΛCDM-based model. We also show that although the halo mass of many satellite galaxies is comparable to the halo mass of the MW progenitor at z ≳ 10, at these early epochs halos that survive as satellites to z = 0 are located many virial radii away from the MW progenitors and thus do not have a chance to merge with it. Our results show that neither the densities estimated in observed Milky Way satellites nor their mass assembly histories pose a challenge to the ΛCDM model. In fact, the broad agreement between density trends with the stellar mass of the observed and model galaxies can be considered as yet another success of the model.

     
    more » « less
  2. ABSTRACT

    Anomalously high nitrogen-to-oxygen abundance ratios [N/O] are observed in globular clusters (GCs), among the field stars of the Milky Way (MW), and even in the gas in a z ≈ 11 galaxy. Using data from the APOGEE Data Release 17 and the Gaia Data Release 3, we present several independent lines of evidence that most of the MW’s high-[N/O] stars were born in situ in massive bound clusters during the early, pre-disc evolution of the Galaxy. Specifically, we show that distributions of metallicity [Fe/H], energy, the angular momentum Lz, and distance of the low-metallicity high-[N/O] stars match the corresponding distributions of stars of the Aurora population and of the in situ GCs. We also show that the fraction of in situ field high-[N/O] stars, fN/O, increases rapidly with decreasing metallicity. During epochs when metallicity evolves from $\rm [Fe/H]=-1.5$ to $\rm [Fe/H]=-0.9$, the Galaxy spins up and transitions from a turbulent Aurora state to a coherently rotating disc. This transformation is accompanied by many qualitative changes. In particular, we show that high N/O abundances similar to those observed in GN-z11 were common before the spin-up ($\rm [Fe/H]\lesssim -1.5$) when up to $\approx 50~{{\ \rm per\ cent}}-70~{{\ \rm per\ cent}}$ of the in situ stars formed in massive bound clusters. The dramatic drop of fN/O at $\rm [Fe/H]\gtrsim -0.9$ indicates that after the disc emerges the fraction of stars forming in massive bound clusters decreases by two orders of magnitude.

     
    more » « less
  3. ABSTRACT

    We present a simple regulator-type framework designed specifically for modelling formation of dwarf galaxies. Despite its simplicity, when coupled with realistic mass accretion histories of haloes from simulations and reasonable choices for model parameter values, the framework can reproduce a remarkably broad range of observed properties of dwarf galaxies over seven orders of magnitude in stellar mass. In particular, we show that the model can simultaneously match observational constraints on the stellar mass–halo mass relation, as well as observed relations between stellar mass and gas phase and stellar metallicities, gas mass, size, and star formation rate, as well as general form and diversity of star formation histories of observed dwarf galaxies. The model can thus be used to predict photometric properties of dwarf galaxies hosted by dark matter haloes in N-body simulations, such as colours, surface brightnesses, and mass-to-light ratios and to forward model observations of dwarf galaxies. We present examples of such modelling and show that colours and surface brightness distributions of model galaxies are in good agreement with observed distributions for dwarfs in recent observational surveys. We also show that in contrast with the common assumption, the absolute magnitude–halo mass relation is generally predicted to have a non-power law form in the dwarf regime, and that the fraction of haloes that host detectable ultra-faint galaxies is sensitive to reionization redshift (zrei) and is predicted to be consistent with observations for zrei ≲ 9.

     
    more » « less
  4. ABSTRACT

    We use accurate estimates of aluminium abundance from the APOGEE Data Release 17 and Gaia Early Data Release 3 astrometry to select a highly pure sample of stars with metallicity −1.5 ≲ [Fe/H] ≲ 0.5 born in-situ in the Milky Way proper. The low-metallicity ([Fe/H]  ≲ −1.3) in-situ component we dub Aurora is kinematically hot with an approximately isotropic velocity ellipsoid and a modest net rotation. Aurora stars exhibit large scatter in metallicity and in many element abundance ratios. The median tangential velocity of the in-situ stars increases sharply with metallicity between [Fe/H] = −1.3 and −0.9, the transition that we call the spin-up. The observed and theoretically expected age–metallicity correlations imply that this increase reflects a rapid formation of the MW disc over ≈1–2 Gyr. The transformation of the stellar kinematics as a function of [Fe/H] is accompanied by a qualitative change in chemical abundances: the scatter drops sharply once the Galaxy builds up a disc during later epochs corresponding to [Fe/H] > −0.9. Results of galaxy formation models presented in this and other recent studies strongly indicate that the trends observed in the MW reflect generic processes during the early evolution of progenitors of MW-sized galaxies: a period of chaotic pre-disc evolution, when gas is accreted along cold narrow filaments and when stars are born in irregular configurations, and subsequent rapid disc formation. The latter signals formation of a stable hot gaseous halo around the MW progenitor, which changes the mode of gas accretion and allows development of coherently rotating disc.

     
    more » « less
  5. null (Ed.)
    We compare the performance of energy-based and entropy-conservative schemes for modeling nonthermal energy components, such as unresolved turbulence and cosmic rays, using idealized fluid dynamics tests and isolated galaxy simulations. While both methods are aimed to model advection and adiabatic compression or expansion of different energy components, the energy-based scheme numerically solves the non-conservative equation for the energy density evolution, while the entropy-conserving scheme uses a conservative equation for modified entropy. Using the standard shock tube and Zel'dovich pancake tests, we show that the energy-based scheme results in a spurious generation of nonthermal energy on shocks, while the entropy-conserving method evolves the energy adiabatically to machine precision. We also show that, in simulations of an isolated Lstar galaxy, switching between the schemes results in ~20-30% changes of the total star formation rate and a significant difference in morphology, particularly near the galaxy center. We also outline and test a simple method that can be used in conjunction with the entropy-conserving scheme to model the injection of nonthermal energies on shocks. Finally, we discuss how the entropy-conserving scheme can be used to capture the kinetic energy dissipated by numerical viscosity into the subgrid turbulent energy implicitly, without explicit source terms that require calibration and can be rather uncertain. Our results indicate that the entropy-conserving scheme is the preferred choice for modeling nonthermal energy components, a conclusion that is equally relevant for Eulerian and moving-mesh fluid dynamics codes. 
    more » « less
  6. null (Ed.)
    ABSTRACT A self-similar spherical collapse model predicts a dark matter (DM) splashback and accretion shock in the outskirts of galaxy clusters while missing a key ingredient of structure formation – processes associated with mergers. To fill this gap, we perform simulations of merging self-similar clusters and investigate their DM and gas evolution in an idealized cosmological context. Our simulations show that the cluster rapidly contracts during the major merger and the splashback radius rsp decreases, approaching the virial radius rvir. While in the self-similar model rsp depends on a smooth mass accretion rate parameter Γs, our simulations show that in the presence of mergers, rsp responds to the changes in the total mass accretion rate Γvir, which accounts for both mergers and smooth accretion. The scatter of the Γvir − rsp/rvir relation indicates a generally low Γs ∼ 1 in clusters in cosmological simulations. In contrast to the DM, the hot gaseous atmospheres significantly expand by the merger-accelerated (MA-) shocks formed when the runaway merger shocks overtake the outer accretion shock. After a major merger, the MA-shock radius is larger than rsp by a factor of up to ∼1.7 for Γs ≲ 1 and is ∼rsp for Γs ≳ 3. This implies that (1) mergers could easily generate the MA-shock-splashback offset measured in cosmological simulations, and (2) the smooth mass accretion rate is small in regions away from filaments where MA-shocks reside. We further discuss the shapes of the DM haloes, various shocks, and contact discontinuities formed at different epochs of the merger, and the ram-pressure stripping in cluster outskirts. 
    more » « less
  7. null (Ed.)
    ABSTRACT We examine the thermodynamic state and cooling of the low-z circumgalactic medium (CGM) in five FIRE-2 galaxy formation simulations of Milky Way-mass galaxies. We find that the CGM in these simulations is generally multiphase and dynamic, with a wide spectrum of largely non-linear density perturbations sourced by the accretion of gas from the intergalactic medium (IGM) and outflows from both the central and satellite galaxies. We investigate the origin of the multiphase structure of the CGM with a particle-tracking analysis and find that most of the low-entropy gas has cooled from the hot halo as a result of thermal instability triggered by these perturbations. The ratio of cooling to free-fall time-scales tcool/tff in the hot component of the CGM spans a wide range of ∼1−100 at a given radius but exhibits approximately constant median values of ∼5−20 at all radii 0.1Rvir < r < Rvir. These are similar to the ≈10−20 value typically adopted as the thermal instability threshold in ‘precipitation’ models of the ICM. Consequently, a one-dimensional model based on the assumption of a constant tcool/tff and hydrostatic equilibrium approximately reproduces the number density and entropy profiles of each simulation but only if it assumes the metallicity profile and temperature boundary condition taken directly from the simulation. We explicitly show that the tcool/tff value of a gas parcel in the hot component of the CGM does not predict its probability of subsequently accreting on to the central galaxy. This suggests that the value of tcool/tff is a poor predictor of thermal stability in gaseous haloes in which large-amplitude density perturbations are prevalent. 
    more » « less
  8. null (Ed.)
    The spatial decorrelation of dense molecular gas and young stars observed on ≲ 1 kiloparsec scales in nearby galaxies indicates rapid dispersal of star-forming regions by stellar feedback. We explore the sensitivity of this decorrelation to different processes controlling the structure of the interstellar medium, the abundance of molecular gas, star formation, and feedback in a suite of simulations of an isolated dwarf galaxy with structural properties similar to NGC300 that self-consistently model radiative transfer and molecular chemistry. Our fiducial simulation reproduces the magnitude of decorrelation and its scale dependence measured in NGC300, and we show that this agreement is due to different aspects of feedback, including H2 dissociation, gas heating by the locally variable UV field, early mechanical feedback, and supernovae. In particular, early radiative and mechanical feedback affect the correlation on ≲100 pc scales, while supernovae play a significant role on ≳100 pc scales. The correlation is also sensitive to the choice of the local star formation efficiency per freefall time, eps_ff, which provides a strong observational constraint on eps_ff when the global star formation rate is independent of its value. Finally, we explicitly show that the degree of correlation between the peaks of molecular gas and star formation density is directly related to the distribution of the lifetimes of star-forming regions. 
    more » « less