skip to main content
US FlagAn official website of the United States government
dot gov icon
Official websites use .gov
A .gov website belongs to an official government organization in the United States.
https lock icon
Secure .gov websites use HTTPS
A lock ( lock ) or https:// means you've safely connected to the .gov website. Share sensitive information only on official, secure websites.


Title: Predictions for detecting a turndown in the baryonic Tully–Fisher relation
ABSTRACT The baryonic Tully–Fisher relation (bTFR) provides an empirical connection between baryonic mass and dynamical mass (measured by the maximum rotation velocity) for galaxies. Due to the impact of baryonic feedback in the shallower potential wells of dwarf galaxies, the bTFR is predicted to turn down at low masses from the extrapolated power-law relation at high masses. The low-mass end of the bTFR is poorly constrained due to small samples and difficulty in connecting the galaxy’s gas kinematics to its dark matter halo. Simulations can help us understand this connection and interpret observations. We measure the bTFR with 75 dwarf galaxies from the Marvel-ous and Marvelous Massive Dwarfs hydrodynamic simulations. Our sample has M$$_\star = 10^6-10^9$$ M$$_\odot$$, and is mostly gas dominated. We compare five velocity methods: V$$_\text{out,circ}$$ (spatially resolved mass-enclosed), V$$_\text{out,mid}$$ (spatially resolved mid-plane gravitational potential), and unresolved H i linewidths at different percentages of the peak flux (W$$_\text{10}$$, W$$_\text{20}$$, and W$$_\text{50}$$). We find an intrinsic turndown in the bTFR for maximum halo speeds $$\lesssim 50$$ km s$$^{-1}$$, or total baryonic mass M$$_\text{bary}\lesssim 10^{8.5}$$ M$$_\odot$$. We find that observing H i in lower-mass galaxies to the conventional surface density limit of 1 M$$_\odot$$ pc$$^{-2}$$ is not enough to detect a turndown in the bTFR; none of the H i velocity methods, spatially resolved or unresolved, recover the turndown, and we find bTFR slopes consistent with observations of higher-mass galaxies. However, we predict that the turndown can be recovered by resolved rotation curves if the H i limit is $$\lesssim 0.08$$ M$$_\odot$$ pc$$^{-2}$$, which is within the sensitivity of current H i surveys like FEASTS and MHONGOOSE.  more » « less
Award ID(s):
2008110
PAR ID:
10617002
Author(s) / Creator(s):
; ; ; ; ; ; ;
Publisher / Repository:
Oxford University Press
Date Published:
Journal Name:
Monthly Notices of the Royal Astronomical Society
Volume:
541
Issue:
3
ISSN:
0035-8711
Format(s):
Medium: X Size: p. 2180-2196
Size(s):
p. 2180-2196
Sponsoring Org:
National Science Foundation
More Like this
  1. Abstract The ratio of baryonic-to-dark matter in present-day galaxies constrains galaxy formation theories and can be determined empirically via the baryonic Tully–Fisher relation (BTFR), which compares a galaxy’s baryonic mass ( M bary ) to its maximum rotation velocity ( V max ). The BTFR is well determined at M bary > 10 8 M ⊙ , but poorly constrained at lower masses due to small samples and the challenges of measuring rotation velocities in this regime. For 25 galaxies with high-quality data and M bary ≲ 10 8 M ⊙ , we estimate M bary from infrared and H i observations and V max from the H i gas rotation. Many of the V max values are lower limits because the velocities are still rising at the edge of the detected H i disks ( R max ); consequently, most of our sample has lower velocities than expected from extrapolations of the BTFR at higher masses. To estimate V max , we map each galaxy to a dark matter halo assuming density profiles with and without cores. In contrast to noncored profiles, we find the cored profile rotation curves are still rising at R max values, similar to the data. When we compare the V max values derived from the cored density profiles to our M bary measurements, we find a turndown of the BTFR at low masses that is consistent with Λ cold dark matter predictions and implies baryon fractions of 1%–10% of the cosmic value. Although we are limited by the sample size and assumptions inherent in mapping measured rotational velocities to theoretical rotation curves, our results suggest that galaxy formation efficiency drops at masses below M bary ∼ 10 8 M ⊙ , corresponding to M 200 ∼ 10 10 M ⊙ . 
    more » « less
  2. Abstract Testing the standard cosmological model (ΛCDM) at small scales is challenging. Galaxies that inhabit low-mass dark matter halos provide an ideal test bed for dark matter models by linking observational properties of galaxies at small scales (low mass, low velocity) to low-mass dark matter halos. However, the observed kinematics of these galaxies do not align with the kinematics of the dark matter halos predicted to host them, obscuring our understanding of the low-mass end of the galaxy–halo connection. We use deep Hiobservations of low-mass galaxies at high spectral resolution in combination with cosmological simulations of dwarf galaxies to better understand the connection between dwarf galaxy kinematics and low-mass halos. Specifically, we use Hiline widths to directly compare to the maximum velocities in a dark matter halo and find that each deeper measurement approaches the expected one-to-one relationship between the observed kinematics and the predicted kinematics in ΛCDM. We also measure baryonic masses and place these on the baryonic Tully–Fisher relation (BTFR). Again, our deepest measurements approach the theoretical predictions for the low-mass end of this relation, a significant improvement on similar measurements based on line widths measured at 50% and 20% of the peak. Our data also hint at the rollover in the BTFR predicted by hydrodynamical simulations of ΛCDM for low-mass galaxies. 
    more » « less
  3. ABSTRACT We study a suite of extremely high-resolution cosmological Feedback in Realistic Environments simulations of dwarf galaxies ($$M_{\rm halo} \lesssim 10^{10}\rm \, M_{\odot }$$), run to z = 0 with $$30\, \mathrm{M}_{\odot }$$ resolution, sufficient (for the first time) to resolve the internal structure of individual supernovae remnants within the cooling radius. Every halo with $$M_{\rm halo} \gtrsim 10^{8.6}\, \mathrm{M}_{\odot }$$ is populated by a resolved stellar galaxy, suggesting very low-mass dwarfs may be ubiquitous in the field. Our ultra-faint dwarfs (UFDs; $$M_{\ast }\lt 10^{5}\, \mathrm{M}_{\odot }$$) have their star formation (SF) truncated early (z ≳ 2), likely by reionization, while classical dwarfs ($$M_{\ast }\gt 10^{5}\, \mathrm{M}_{\odot }$$) continue forming stars to z < 0.5. The systems have bursty star formation histories, forming most of their stars in periods of elevated SF strongly clustered in both space and time. This allows our dwarf with M*/Mhalo > 10−4 to form a dark matter core $${\gt}200\rm \, pc$$, while lower mass UFDs exhibit cusps down to $${\lesssim}100\rm \, pc$$, as expected from energetic arguments. Our dwarfs with $$M_{\ast }\gt 10^{4}\, \mathrm{M}_{\odot }$$ have half-mass radii (R1/2) in agreement with Local Group (LG) dwarfs (dynamical mass versus R1/2 and stellar rotation also resemble observations). The lowest mass UFDs are below surface brightness limits of current surveys but are potentially visible in next-generation surveys (e.g. LSST). The stellar metallicities are lower than in LG dwarfs; this may reflect pre-enrichment of the LG by the massive hosts or Pop-III stars. Consistency with lower resolution studies implies that our simulations are numerically robust (for a given physical model). 
    more » « less
  4. ABSTRACT Atomic hydrogen (H i) serves a crucial role in connecting galactic-scale properties such as star formation with the large-scale structure of the Universe. While recent numerical simulations have successfully matched the observed covering fraction of H i near Lyman Break Galaxies (LBGs) and in the foreground of luminous quasars at redshifts $$z \lesssim 3$$, the low-mass end remains as-of-yet unexplored in observational and computational surveys. We employ a cosmological, hydrodynamical simulation (FIREbox) supplemented with zoom-in simulations (MassiveFIRE) from the Feedback In Realistic Environments (FIRE) project to investigate the H i covering fraction of Lyman Limit Systems ($$N_{{\text{H}}\, \rm{{\small I}}} \gtrsim 10^{17.2}$$ cm$$^{-2}$$) across a wide range of redshifts ($z=0-6$) and halo masses ($$10^8-10^{13} \, \,\mathrm{ M}_{\odot }$$ at $z=0$, $$10^8-10^{11}\, \,\mathrm{ M}_{\odot }$$ at $z=6$) in the absence of feedback from active galactic nuclei. We find that the covering fraction inside haloes exhibits a strong increase with redshift, with only a weak dependence on halo mass for higher mass haloes. For massive haloes ($$M_{\mathrm{vir}} \sim 10^{11}-10^{12} \,\mathrm{ M}_{\odot }$$), the radial profiles showcase scale-invariance and remain independent of mass. The radial dependence is well captured by a fitting function. The covering fractions in our simulations are in good agreement with measurements of the covering fraction in LBGs. Our comprehensive analysis unveils a complex dependence with redshift and halo mass for haloes with $$M_{\mathrm{vir}} \lesssim 10^{10} \,\mathrm{ M}_{\odot }$$ that future observations aim to constrain, providing key insights into the physics of structure formation and gas assembly. 
    more » « less
  5. null (Ed.)
    ABSTRACT We present the first set of cosmological baryonic zoom-in simulations of galaxies including dissipative self-interacting dark matter (dSIDM). These simulations utilize the Feedback In Realistic Environments galaxy formation physics, but allow the dark matter to have dissipative self-interactions analogous to standard model forces, parametrized by the self-interaction cross-section per unit mass, (σ/m), and the dimensionless degree of dissipation, 0 < fdiss < 1. We survey this parameter space, including constant and velocity-dependent cross-sections, and focus on structural and kinematic properties of dwarf galaxies with $$M_{\rm halo} \sim 10^{10-11}{\, \rm M_\odot }$$ and $$M_{\ast } \sim 10^{5-8}{\, \rm M_\odot }$$. Central density profiles (parametrized as ρ ∝ rα) of simulated dwarfs become cuspy when $$(\sigma /m)_{\rm eff} \gtrsim 0.1\, {\rm cm^{2}\, g^{-1}}$$ (and fdiss = 0.5 as fiducial). The power-law slopes asymptote to α ≈ −1.5 in low-mass dwarfs independent of cross-section, which arises from a dark matter ‘cooling flow’. Through comparisons with dark matter only simulations, we find the profile in this regime is insensitive to the inclusion of baryons. However, when $$(\sigma /m)_{\rm eff} \ll 0.1\, {\rm cm^{2}\, g^{-1}}$$, baryonic effects can produce cored density profiles comparable to non-dissipative cold dark matter (CDM) runs but at smaller radii. Simulated galaxies with $$(\sigma /m) \gtrsim 10\, {\rm cm^{2}\, g^{-1}}$$ and the fiducial fdiss develop significant coherent rotation of dark matter, accompanied by halo deformation, but this is unlike the well-defined thin ‘dark discs’ often attributed to baryon-like dSIDM. The density profiles in this high cross-section model exhibit lower normalizations given the onset of halo deformation. For our surveyed dSIDM parameters, halo masses and galaxy stellar masses do not show appreciable difference from CDM, but dark matter kinematics and halo concentrations/shapes can differ. 
    more » « less