Search for: All records

ABSTRACT We use reduced-mass eigentensors to quantify the 3D ellipsoidal shape evolution of 13 Milky Way-mass galaxies simulated using zoom simulations with fire-2 physics; all but one form discs at $z=0$. We find that all of our Milky Way progenitors go through phases when they are elongated. They often oscillate between spheroidal and elongated shapes in the early Universe over billion-year time-scales, with $$\sim 25-45~{{\ \rm per\ cent}}$$ of the population having elongated luminosity-weighted shapes at any given time at $z = 0.5-8.5$. In contrast, all stellar populations in our $z=0$ Milky Way analogues are symmetric about their minor axes at $z=0$, even though the old and intermediate-age stellar populations were often arranged in the shape of elongated pickles or triaxial spheroids at the time they formed meaning these populations changed shape significantly over time. During their transient elongated phases, our galaxies have anisotropic velocity dispersion ellipsoids directed along their spatial major axis; however, their shapes do not correlate with their dark matter fraction nor with the shapes and orientations of their underlying dark matter haloes. We find that when treated as a population, the fraction of our galaxy progenitors that are elongated at $$z\gt 0.5$$ is roughly consistent with what is observed for systems of the same mass and redshift. Our results suggest that observed elongated galaxies seen in the early Universe with James Webb Space Telescope and Hubble Space Telescope are not stable structures, but rather transitory phases that are nevertheless statistically common. Some of these observed objects may evolve into Milky Way-like galaxies at $z=0$. 

Abstract While galaxy rotation curves (RCs) provide one of the most powerful methods for measuring dark matter profiles in the inner regions of rotation-supported galaxies, at the dwarf scale there are factors that can complicate this analysis. Given the expectation of a universal profile in dark-matter-only simulations, the diversity of the observed RCs has become an often-discussed issue in Lambda cold dark matter cosmology on galactic scales. We analyze a suite of Feedback in Realistic Environments simulations of 10¹⁰–10¹²M_⊙halos with standard cold dark matter and compare the true circular velocity to RC reconstructions. We find that, for galaxies with well-ordered gaseous disks, the measured RC may deviate from the true circular velocity by at most ∼10% within the radius of the disk. However, nonequilibrium behaviors, noncircular motions, and nonthermal and nonkinetic stresses may cause much larger discrepancies, of ∼50% or more. Most RC reconstructions underestimate the true circular velocity, while some reconstructions transiently overestimate it in the central few kiloparsecs, due to dynamical phenomena. We further demonstrate that the features that contribute to these failures are not always visibly obvious in HIobservations. If such dwarf galaxies are included in galaxy catalogs, they may give rise to the appearance of “artificial” RC diversity that does not reflect the true variation in underlying dark matter profiles. 

Abstract We study the intrinsic and observable shapes of approximately 700 star-forming galaxies with stellar masses of 10⁸–10¹¹M_⊙from the FIREbox simulation atz= 0. We calculate intrinsic axis ratios using inertia tensors weighted by three morphology types: “All Stars,” “Young Stars,” and “Luminosity-weighted Stars.” Young Stars shows mass-dependent 3D configurations, with spheroidal, elongated, and disky shapes dominant at stellar masses of 10^8.5, 10^9.5, and 10^10.5M_⊙, respectively. Using the radiative transfer codeSKIRT, we construct mock images for each galaxy and show that projected short-to-long axis ratios,q, inferred from 2D Sérsic fits are most closely related to Luminosity-weighted Stars tensor shapes and least resemble the All Stars’ shapes. This suggests observed 2D shape distributions should not be compared to predictions based on 3D stellar mass shapes. Next, we construct a sample of mock images projected in random orientations and compare them to observed axis ratio distributions from the GAMA survey. At stellar masses below 10¹⁰M_⊙, we produce too few galaxies with observedq < 0.4 and none withq < 0.2, suggesting that FIREbox does not produce enough low-mass disk galaxies. At higher masses, 10¹⁰–10¹¹M_⊙, we find that the predictedqdistribution is sensitive to the dust-to-metal ratio; using our fiducial model, the distribution ofqvalues is formally consistent with observations, but there is tension with our ability to produce enough very thin systems withq < 0.2. Future observational and theoretical programs aimed at understanding disk and thin-disk fractions will provide crucial tests of galaxy formation models. 

Abstract We explore how a realistic surface brightness (SB) detection limit ofμ_V≈ 32.5 mag arcsec⁻²for stars at the edges of ultrafaint galaxies affects our ability to infer their underlying properties. We use a sample of 19 galaxies with stellar masses ≈400–40,000M_⊙simulated with FIRE-2 physics and baryonic mass resolution of 30M_⊙. The SB cut leads to smaller sizes, lower stellar masses, and lower stellar velocity dispersions than the values inferred without the cut. However, by imposing this realistic limit, our inferred galaxy properties lie closer to observed populations in the mass-size plane, better match observed velocity dispersions as a function of stellar mass, and better reproduce derived circular velocities as a function of half-light radius. For the most massive galaxies in our sample, the SB cut leads to higher mean [Fe/H] values, but the increase is not enough to match the observed MZR. Finally, we demonstrate that the common J. Wolf et al. dynamical mass estimator is less accurate when the SB cut is applied. For our lowest-mass galaxies, in particular, excluding the low-surface brightness outskirts causes us to overestimate their central dark-matter densities and virial masses. This suggests that attempts to use mass estimates of ultrafaint galaxies to constrain dark-matter physics or to place constraints on the low-mass threshold of galaxy formation must take into account surface brightness limits or risk significant biases. 

ABSTRACT We study the morphology of hundreds of simulated central galaxies in the stellar mass range $$M_\star =$$ 107.5–1011  $$\rm M_\odot$$ from the FIREbox cosmological volume. We demonstrate that FIREbox is able to predict a wide variety of morphologies, spanning from disc-dominated objects to spheroidal galaxies supported by stellar velocity dispersion. However, the simulations predict a strong relation between morphology (degree of rotational support) and stellar mass: galaxies comparable to the Milky Way are often disc-dominated while the presence of stellar discs mostly vanishes for dwarfs with $$M_\star < 10^9 ~$$\rm M_\odot$$. This defines a ‘morphology transition’ regime for galaxies with $$10^9 < M_\star /\rm {M_\odot }< 10^{10}$$ in which discs become increasingly common, but below which discs are rare. We show that burstiness in the star formation history and the deepening of the gravitational potential strongly correlate in our simulations with this transition regime, with discs forming in objects with lower levels of burstiness in the last $$\sim 6$$ Gyr and haloes with mass $$\sim 10^{11} ~ \rm {{\rm M}_{\odot }}$$ and above. While observations support a transition towards thicker discs in the regime of dwarfs, our results are in partial disagreement with observations of at least some largely rotationally supported gas discs in dwarfs with $$M_\star < 10^9$$\rm M_\odot$$. This study highlights dwarf morphology as a fundamental benchmark for testing future galaxy formation models. 

Abstract We utilize the cosmological volume simulation FIREbox to investigate how a galaxy’s environment influences its size and dark matter content. Our study focuses on approximately 1200 galaxies (886 central and 332 satellite halos) in the low-mass regime, with stellar masses between 10⁶and 10⁹M_⊙. We analyze the size–mass relation (r₅₀–M_⋆), the inner dark matter mass–stellar mass ( $M_{\mathrm{DM}}^{50}$ –M_⋆) relation, and the halo mass–stellar mass (M_halo–M_⋆) relation. At fixed stellar mass, we find that galaxies experiencing stronger tidal influences, indicated by higher Perturbation Indices (PI > 1) are generally larger and have lower halo masses relative to their counterparts with lower Perturbation Indices (PI < 1). Applying a Random Forest regression model, we show that both the environment (PI) and halo mass (M_halo) are significant predictors of a galaxy’s relative size and dark matter content. Notably, becauseM_halois also strongly affected by the environment, our findings indicate that environmental conditions not only influence galactic sizes and relative inner dark matter content directly, but also indirectly, through their impact on halo mass. Our results highlight a critical interplay between environmental factors and halo mass in shaping galaxy properties, affirming the environment as a fundamental driver in galaxy formation and evolution. 

ABSTRACT The observationally inferred size versus stellar–mass relationship (SMR) for low-mass galaxies provides an important test for galaxy formation models. However, the relationship relies on assumptions that relate observed luminosity profiles to underlying stellar mass profiles. Here we use the Feedback in Realistic Environments simulations of low-mass galaxies to explore how the predicted SMR changes depending on whether one uses star-particle counts directly or mock observations. We reproduce the SMR found in The Exploration of Local Volume Satellites survey remarkably well only when we infer stellar masses and sizes using mock observations. However, when we use star particles to directly infer stellar masses and half-mass radii, we find that our galaxies are too large and obey an SMR with too little scatter compared to observations. This discrepancy between the ‘true’ galaxy size and mass and those derived in the mock observation approach is twofold. First, our simulated galaxies have higher and more varied mass-to-light ratios (MLR) at a fixed colour than those commonly adopted, which tends to underestimate their stellar masses compared to their true, simulated values. Second, our galaxies have radially increasing MLR gradients therefore using a single MLR tends to underpredict the mass in the outer regions. Similarly, the true half-mass radius is larger than the half-light radius because the light is more concentrated than the mass. If our simulations are accurate representations of the real Universe, then the relationship between galaxy size and stellar mass is even tighter for low-mass galaxies than is commonly inferred from observed relations. 

ABSTRACT The radial acceleration relation (RAR) connects the total gravitational acceleration of a galaxy at a given radius, atot(r), with that accounted for by baryons at the same radius, abar(r). The shape and tightness of the RAR for rotationally-supported galaxies have characteristics in line with MOdified Newtonian Dynamics (MOND) and can also arise within the cosmological constant + cold dark matter (ΛCDM) paradigm. We use zoom simulations of 20 galaxies with stellar masses of M⋆ ≃ 107–11 M⊙ to study the RAR in the FIRE-2 simulations. We highlight the existence of simulated galaxies with non-monotonic RAR tracks that ‘hook’ down from the average relation. These hooks are challenging to explain in Modified Inertia theories of MOND, but naturally arise in all of our ΛCDM-simulated galaxies that are dark-matter dominated at small radii and have feedback-induced cores in their dark matter haloes. We show, analytically and numerically, that downward hooks are expected in such cored haloes because they have non-monotonic acceleration profiles. We also extend the relation to accelerations below those traced by disc galaxy rotation curves. In this regime, our simulations exhibit ‘bends’ off of the MOND-inspired extrapolation of the RAR, which, at large radii, approach atot ≈ abar/fb, where fb is the cosmic baryon fraction. Future efforts to search for these hooks and bends in real galaxies will provide interesting tests for MOND and ΛCDM. 

Abstract We present ∼300 stellar metallicity measurements in two faint M31 dwarf galaxies, Andromeda XVI (M_V= −7.5) and Andromeda XXVIII (M_V= –8.8), derived using metallicity-sensitive calcium H and K narrowband Hubble Space Telescope imaging. These are the first individual stellar metallicities in And XVI (95 stars). Our And XXVIII sample (191 stars) is a factor of ∼15 increase over literature metallicities. For And XVI, we measure $〈\mathrm{[Fe/H]}〉=-{2.17}_{-0.05}^{+0.05}$ , ${\sigma }_{\mathrm{[Fe/H]}}={0.33}_{-0.07}^{+0.07}$ , and ∇_[Fe/H]= −0.23 ± 0.15 dex $R_e^{-1}$ . We find that And XVI is more metal-rich than Milky Way ultrafaint dwarf galaxies of similar luminosity, which may be a result of its unusually extended star formation history. For And XXVIII, we measure $〈\mathrm{[Fe/H]}〉=-{1.95}_{-0.04}^{+0.04}$ , ${\sigma }_{\mathrm{[Fe/H]}}={0.34}_{-0.05}^{+0.05}$ , and ∇_[Fe/H]= −0.46 ± 0.10 dex $R_e^{-1}$ , placing it on the dwarf galaxy mass–metallicity relation. Neither galaxy has a metallicity distribution function (MDF) with an abrupt metal-rich truncation, suggesting that star formation fell off gradually. The stellar metallicity gradient measurements are among the first for faint (L≲ 10⁶L_⊙) galaxies outside the Milky Way halo. Both galaxies’ gradients are consistent with predictions from the FIRE simulations, where an age–gradient strength relationship is the observational consequence of stellar feedback that produces dark matter cores. We include a catalog for community spectroscopic follow-up, including 19 extremely metal-poor ([Fe/H] < –3.0) star candidates, which make up 7% of And XVI’s MDF and 6% of And XXVIII’s. 

Abstract We measure the metallicities of 374 red giant branch (RGB) stars in the isolated, quenched dwarf galaxy Tucana using Hubble Space Telescope narrowband (F395N) calcium H and K imaging. Our sample is a factor of ∼7 larger than what is available from previous studies. Our main findings are as follows. (i) A global metallicity distribution function (MDF) with $〈\mathrm{[Fe/H]}〉=-{1.55}_{-0.04}^{+0.04}$ and ${\sigma }_{\mathrm{[Fe/H]}}={0.54}_{-0.03}^{+0.03}$ . (ii) A metallicity gradient of −0.54 ± 0.07 dex $R_e^{-1}$ (−2.1 ± 0.3 dex kpc⁻¹) over the extent of our imaging (∼2.5R_e), which is steeper than literature measurements. Our finding is consistent with predicted gradients from the publicly available FIRE-2 simulations, in which bursty star formation creates stellar population gradients and dark matter cores. (iii) Tucana’s bifurcated RGB has distinct metallicities: a blue RGB with $〈\mathrm{[Fe/H]}〉=-{1.78}_{-0.06}^{+0.06}$ and ${\sigma }_{\mathrm{[Fe/H]}}={0.44}_{-0.06}^{+0.07}$ and a red RGB with $〈\mathrm{[Fe/H]}〉=-{1.08}_{-0.07}^{+0.07}$ and ${\sigma }_{\mathrm{[Fe/H]}}={0.42}_{-0.06}^{+0.06}$ . (iv) At fixed stellar mass, Tucana is more metal-rich than Milky Way satellites by ∼0.4 dex, but its blue RGB is chemically comparable to the satellites. Tucana’s MDF appears consistent with star-forming isolated dwarfs, though MDFs of the latter are not as well populated. (v) About 2% of Tucana’s stars have [Fe/H] < −3% and 20% have [Fe/H] > −1. We provide a catalog for community spectroscopic follow-up.