ABSTRACT We study the spatially resolved (sub-kpc) gas velocity dispersion (σ)–star formation rate (SFR) relation in the FIRE-2 (Feedback in Realistic Environments) cosmological simulations. We specifically focus on Milky Way-mass disc galaxies at late times (z ≈ 0). In agreement with observations, we find a relatively flat relationship, with σ ≈ 15–30 km s−1 in neutral gas across 3 dex in SFRs. We show that higher dense gas fractions (ratios of dense gas to neutral gas) and SFRs are correlated at constant σ. Similarly, lower gas fractions (ratios of gas to stellar mass) are correlated with higher σ at constant SFR. The limits of the σ–ΣSFR relation correspond to the onset of strong outflows. We see evidence of ‘on-off’ cycles of star formation in the simulations, corresponding to feedback injection time-scales of 10–100 Myr, where SFRs oscillate about equilibrium SFR predictions. Finally, SFRs and velocity dispersions in the simulations agree well with feedback-regulated and marginally stable gas disc (Toomre’s Q = 1) model predictions, and the simulation data effectively rule out models assuming that gas turns into stars at (low) constant efficiency (i.e. 1 per cent per free-fall time). And although the simulation data do not entirely exclude gas accretion/gravitationally powered turbulence as a driver of σ, it appears to be subdominant to stellar feedback in the simulated galaxy discs at z ≈ 0.
more »
« less
Dynamical data mining captures disc–halo couplings that structure galaxies
Studying coupling between different galactic components is a challenging problem in galactic dynamics. Using basis function expansions (BFEs) and multichannel singular spectrum analysis (mSSA) as a means of dynamical data mining, we discover evidence for two multicomponent disc–halo dipole modes in a Milky-Way-like simulated galaxy. One of the modes grows throughout the simulation, while the other decays throughout the simulation. The multicomponent disc–halo modes are driven primarily by the halo, and have implications for the structural evolution of galaxies, including observations of lopsidedness and other non-axisymmetric structure. In our simulation, the modes create surface density features up to 10 per cent relative to the equilibrium model stellar disc. While the simulated galaxy was constructed to be in equilibrium, BFE + mSSA also uncovered evidence of persistent periodic signals incited by aphysical initial conditions disequilibrium, including rings and weak two-armed spirals, both at the 1 per cent level. The method is sensitive to distinct evolutionary features at and even below the 1 per cent level of surface density variation. The use of mSSA produced clean signals for both modes and disequilibrium, efficiently removing variance owing to estimator noise from the input BFE time series. The discovery of multicomponent halo–disc modes is strong motivation for application of BFE + mSSA to the rich zoo of dynamics of multicomponent interacting galaxies.
more » « less- NSF-PAR ID:
- 10401697
- Publisher / Repository:
- Oxford University Press
- Date Published:
- Journal Name:
- Monthly Notices of the Royal Astronomical Society
- Volume:
- 521
- Issue:
- 2
- ISSN:
- 0035-8711
- Page Range / eLocation ID:
- p. 1757-1774
- Format(s):
- Medium: X
- Sponsoring Org:
- National Science Foundation
More Like this
-
-
ABSTRACT We analyse the cold dark matter density profiles of 54 galaxy haloes simulated with Feedback In Realistic Environments (FIRE)-2 galaxy formation physics, each resolved within $0.5{{\ \rm per\ cent}}$ of the halo virial radius. These haloes contain galaxies with masses that range from ultrafaint dwarfs ($M_\star \simeq 10^{4.5}\, \mathrm{M}_{\odot }$) to the largest spirals ($M_\star \simeq 10^{11}\, \mathrm{M}_{\odot }$) and have density profiles that are both cored and cuspy. We characterize our results using a new, analytic density profile that extends the standard two-parameter Einasto form to allow for a pronounced constant density core in the resolved innermost radius. With one additional core-radius parameter, rc, this three-parameter core-Einasto profile is able to characterize our feedback-impacted dark matter haloes more accurately than other three-parameter profiles proposed in the literature. To enable comparisons with observations, we provide fitting functions for rc and other profile parameters as a function of both M⋆ and M⋆/Mhalo. In agreement with past studies, we find that dark matter core formation is most efficient at the characteristic stellar-to-halo mass ratio M⋆/Mhalo ≃ 5 × 10−3, or $M_{\star } \sim 10^9 \, \mathrm{M}_{\odot }$, with cores that are roughly the size of the galaxy half-light radius, rc ≃ 1−5 kpc. Furthermore, we find no evidence for core formation at radii $\gtrsim 100\ \rm pc$ in galaxies with M⋆/Mhalo < 5 × 10−4 or $M_\star \lesssim 10^6 \, \mathrm{M}_{\odot }$. For Milky Way-size galaxies, baryonic contraction often makes haloes significantly more concentrated and dense at the stellar half-light radius than DMO runs. However, even at the Milky Way scale, FIRE-2 galaxy formation still produces small dark matter cores of ≃ 0.5−2 kpc in size. Recent evidence for a ∼2 kpc core in the Milky Way’s dark matter halo is consistent with this expectation.more » « less
-
null (Ed.)ABSTRACT We present deep far-infrared observations of the nearby edge-on galaxy NGC 891 obtained with the Herschel Space Observatory and the Spitzer Space Telescope. The maps confirm the detection of thermal emission from the inner circumgalactic medium (halo) and spatially resolve a dusty superbubble and a dust spur (filament). The dust temperature of the halo component is lower than that of the disc but increases across a region of diameter ≈8.0 kpc extending at least 7.7 kpc vertically from one side of the disc, a region we call a superbubble because of its association with thermal X-ray emission and a minimum in the synchrotron scale height. This outflow is breaking through the thick disc and developing into a galactic wind, which is of particular interest because NGC 891 is not considered a starburst galaxy; the star formation rate surface density, 0.03 M⊙ yr−1 kpc−2, and gas fraction, just $10{{\ \rm per\ cent}}$ in the inner disc, indicate the threshold for wind formation is lower than previous work has suggested. We conclude that the star formation surface density is sufficient for superbubble blowout into the halo, but the cosmic ray electrons may play a critical role in determining whether this outflow develops into a fountain or escapes from the gravitational potential. The high dust-to-gas ratio in the dust spur suggests the material was pulled out of NGC 891 through the collision of a minihalo with the disc of NGC 891. We conclude that NGC 891 offers an example of both feedback and satellite interactions transporting dust into the halo of a typical galaxy.more » « less
-
more » « less
ABSTRACT We revisit the classical KZ problem – determination of the vertical force and implied total mass density distribution of the Milky Way disc – for a wide range of Galactocentric radius and vertical height using chemically selected thin and thick disc samples based on Apache Point Observatory Galactic Evolution Experiment spectroscopy combined with the Gaia astrometry. We derived the velocity dispersion profiles in Galactic cylindrical coordinates, and solved the Jeans equation for the two samples separately. The result is surprising that the total surface mass density as a function of vertical height as derived for these two chemically distinguished populations is different. The discrepancies are larger in the inner compared to the outer Galaxy, with the density calculated from thick disc being larger, independent of the Galactic radius. Furthermore, while there is an overall good agreement between the total mass density derived for the thick disc population and the standard halo model for vertical heights larger than 1 kpc, close to the mid-plane the mass density observed using the thick disc population is larger than that predicted from the standard halo model. We explore various implications of these discrepancies, and speculate their sources, including problems associated with the assumed density laws, velocity dispersion profiles, and the Galactic rotation curve, potential non-equilibrium of the Galactic disc, or a failure of the Navarro-Frenk-White (NFW) dark matter halo profile for the Milky Way. We conclude that the growing detail in hand on the chemodynamical distributions of Milky Way stars challenges traditional analytical treatments of the KZ problem.
-
more » « less
ABSTRACT The combination of galaxy–galaxy lensing (GGL) and galaxy clustering is a powerful probe of low-redshift matter clustering, especially if it is extended to the non-linear regime. To this end, we use an N-body and halo occupation distribution (HOD) emulator method to model the redMaGiC sample of colour-selected passive galaxies in the Dark Energy Survey (DES), adding parameters that describe central galaxy incompleteness, galaxy assembly bias, and a scale-independent multiplicative lensing bias Alens. We use this emulator to forecast cosmological constraints attainable from the GGL surface density profile ΔΣ(rp) and the projected galaxy correlation function wp, gg(rp) in the final (Year 6) DES data set over scales $r_p=0.3\!-\!30.0\, h^{-1} \, \mathrm{Mpc}$. For a $3{{\ \rm per\ cent}}$ prior on Alens we forecast precisions of $1.9{{\ \rm per\ cent}}$, $2.0{{\ \rm per\ cent}}$, and $1.9{{\ \rm per\ cent}}$ on Ωm, σ8, and $S_8 \equiv \sigma _8\Omega _m^{0.5}$, marginalized over all halo occupation distribution (HOD) parameters as well as Alens. Adding scales $r_p=0.3\!-\!3.0\, h^{-1} \, \mathrm{Mpc}$ improves the S8 precision by a factor of ∼1.6 relative to a large scale ($3.0\!-\!30.0\, h^{-1} \, \mathrm{Mpc}$) analysis, equivalent to increasing the survey area by a factor of ∼2.6. Sharpening the Alens prior to $1{{\ \rm per\ cent}}$ further improves the S8 precision to $1.1{{\ \rm per\ cent}}$, and it amplifies the gain from including non-linear scales. Our emulator achieves per cent-level accuracy similar to the projected DES statistical uncertainties, demonstrating the feasibility of a fully non-linear analysis. Obtaining precise parameter constraints from multiple galaxy types and from measurements that span linear and non-linear clustering offers many opportunities for internal cross-checks, which can diagnose systematics and demonstrate the robustness of cosmological results.
