Search for: All records

Quantifying the connection between galaxies and their host dark matter halos has been key for testing cosmological models on various scales. BelowM_⋆∼ 10⁹M_⊙, such studies have primarily relied on the satellite galaxy population orbiting the Milky Way (MW). Here we present new constraints on the connection between satellite galaxies and their host dark matter subhalos using the largest sample of satellite galaxies in the Local Volume (D≲ 12 Mpc) to date. We use 250 confirmed and 71 candidate dwarf satellites around 27 MW-like hosts from the Exploration of Local VolumE Satellites (ELVES) Survey and use the semianalyticalSatGenmodel for predicting the population of dark matter subhalos expected in the same volume. Through a Bayesian model comparison of the observed and the forward-modeled satellite stellar mass functions (SSMFs), we infer the satellite stellar-to-halo mass relation. We find that the observed SSMF is best reproduced when subhalos at the low-mass end are populated by a relation of the form$M_{\star }\propto M_{\mathrm{peak}}^{\alpha }$, with a moderate slope of${\alpha }_{\mathrm{const}}=2.10\pm 0.01$and a low scatter, constant as a function of the peak halo mass, of${\sigma }_{\mathrm{const}}={0.06}_{-0.05}^{+0.07}$. A model with a steeper slope (α_grow= 2.39 ± 0.06) and a scatter that grows with decreasingM_peakis also consistent with the observed SSMF but is not required. Our new model for the satellite–subhalo connection, based on hundreds of Local Volume satellite galaxies, is in line with what was previously derived using only MW satellites.

 

Modelling self-gravity of collisionless fluids (e.g. ensembles of dark matter, stars, black holes, dust, and planetary bodies) in simulations is challenging and requires some force softening. It is often desirable to allow softenings to evolve adaptively, in any high-dynamic range simulation, but this poses unique challenges of consistency, conservation, and accuracy, especially in multiphysics simulations where species with different ‘softening laws’ may interact. We therefore derive a generalized form of the energy-and-momentum conserving gravitational equations of motion, applicable to arbitrary rules used to determine the force softening, together with consistent associated time-step criteria, interaction terms between species with different softening laws, and arbitrary maximum/minimum softenings. We also derive new methods to maintain better accuracy and conservation when symmetrizing forces between particles. We review and extend previously discussed adaptive softening schemes based on the local neighbour particle density, and present several new schemes for scaling the softening with properties of the gravitational field, i.e. the potential or acceleration or tidal tensor. We show that the ‘tidal softening’ scheme not only represents a physically motivated, translation and Galilean invariant and equivalence-principle respecting (and therefore conservative) method but also imposes negligible time-step or other computational penalties, ensuring that pairwise two-body scattering is small compared to smooth background forces and can resolve outstanding challenges in properly capturing tidal disruption of substructures (minimizing artificial destruction) while also avoiding excessive N-body heating. We make all of this public in the GIZMO code.

 

Isolated dwarf galaxies usually exhibit robust star formation but satellite dwarf galaxies are often devoid of young stars, even in Milky Way–mass groups. Dwarf galaxies thus offer an important laboratory of the environmental processes that cease star formation. We explore the balance of quiescent and star-forming galaxies (quenched fractions) for a sample of ∼400 satellite galaxies around 30 Local Volume hosts from the Exploration of Local VolumE Satellites (ELVES) Survey. We present quenched fractions as a function of satellite stellar mass, projected radius, and host halo mass, to conclude that overall, the quenched fractions are similar to the Milky Way, dropping below 50% at satelliteM_*≈ 10⁸M_⊙. We may see hints that quenching is less efficient at larger radii. Through comparison with the semianalytic modeling codeSatGen, we are also able to infer average quenching times as a function of satellite mass in host halo-mass bins. There is a gradual increase in quenching time with satellite stellar mass rather than the abrupt change from rapid to slow quenching that has been inferred for the Milky Way. We also generally infer longer average quenching times than recent hydrodynamical simulations. Our results are consistent with models that suggest a wide range of quenching times are possible via ram pressure stripping, depending on the clumpiness of the circumgalactic medium, the orbits of the satellites, and the degree of earlier preprocessing.

 

ABSTRACT We combine the isothermal Jeans model and the model of adiabatic halo contraction into a semi-analytic procedure for computing the density profile of self-interacting dark-matter (SIDM) haloes with the gravitational influence from the inhabitant galaxies. The model agrees well with cosmological SIDM simulations over the entire core-forming stage up to the onset of gravothermal core-collapse. Using this model, we show that the halo response to baryons is more diverse in SIDM than in CDM and depends sensitively on galaxy size, a desirable feature in the context of the structural diversity of bright dwarfs. The fast speed of the method facilitates analyses that would be challenging for numerical simulations – notably, we quantify the SIDM halo response as functions of the baryonic properties, on a fine mesh grid spanned by the baryon-to-total-mass ratio, Mb/Mvir, and galaxy compactness, r1/2/Rvir; we show with high statistical precision that for typical Milky-Way-like systems, the SIDM profiles are similar to their CDM counterparts; and we delineate the regime of core-collapse in the Mb/Mvir − r1/2/Rvir space, for a given cross section and concentration. Finally, we compare the isothermal Jeans model with the more sophisticated gravothermal fluid model, and show that the former yields faster core formation and agrees better with cosmological simulations. We attribute the difference to whether the target CDM halo is used as a boundary condition or as the initial condition for the gravothermal evolution, and thus comment on possible improvements of the fluid model. We have made our model publicly available at https://github.com/JiangFangzhou/SIDM. 

We utilize high-resolution cosmological simulations to reveal that high-redshift galaxies tend to undergo a robust ‘wet compaction’ event when near a ‘golden’ stellar mass of $\sim \!\!10^{10}\, \rm M_\odot$ . This is a gaseous shrinkage to a compact star-forming phase, a ‘blue nugget’ (BN), followed by central quenching of star formation to a compact passive stellar bulge, a ‘red nugget’ (RN), and a buildup of an extended gaseous disc and ring. Such nuggets are observed at cosmic noon and seed today’s early-type galaxies. The compaction is triggered by a drastic loss of angular momentum due to, e.g. wet mergers, counter-rotating cold streams, or violent disc instability. The BN phase marks drastic transitions in the galaxy structural, compositional, and kinematic properties. The transitions are from star forming to quenched inside-out, from diffuse to compact with an extended disc or ring and a stellar envelope, from dark matter to baryon central dominance, from prolate to oblate stellar shape, from pressure to rotation support, from low to high metallicity, and from supernova to AGN feedback. The central black hole growth, first suppressed by supernova feedback when below the golden mass, is boosted by the compaction, and the black hole keeps growing once the halo is massive enough to lock in the supernova ejecta.

 

ABSTRACT Self-interacting dark matter (SIDM) cosmologies admit an enormous diversity of dark matter (DM) halo density profiles, from low-density cores to high-density core-collapsed cusps. The possibility of the growth of high central density in low-mass haloes, accelerated if haloes are subhaloes of larger systems, has intriguing consequences for small-halo searches with substructure lensing. However, following the evolution of ${\lesssim}10^8 \, \mathrm{M}_\odot$ subhaloes in lens-mass systems (${\sim}10^{13}\, \mathrm{M}_\odot$) is computationally expensive with traditional N-body simulations. In this work, we develop a new hybrid semi-analytical + N-body method to study the evolution of SIDM subhaloes with high fidelity, from core formation to core-collapse, in staged simulations. Our method works best for small subhaloes (≲1/1000 host mass), for which the error caused by dynamical friction is minimal. We are able to capture the evaporation of subhalo particles by interactions with host halo particles, an effect that has not yet been fully explored in the context of subhalo core-collapse. We find three main processes drive subhalo evolution: subhalo internal heat outflow, host-subhalo evaporation, and tidal effects. The subhalo central density grows only when the heat outflow outweighs the energy gain from evaporation and tidal heating. Thus, evaporation delays or even disrupts subhalo core-collapse. We map out the parameter space for subhaloes to core-collapse, finding that it is nearly impossible to drive core-collapse in subhaloes in SIDM models with constant cross-sections. Any discovery of ultracompact dark substructures with future substructure lensing observations favours additional degrees of freedom, such as velocity-dependence, in the cross-section. 

As they grow, galaxies can transition from irregular/spheroidal with ‘bursty’ star formation histories (SFHs), to discy with smooth SFHs. But even in simulations, the direct physical cause of such transitions remains unclear. We therefore explore this in a large suite of numerical experiments re-running portions of cosmological simulations with widely varied physics, further validated with existing FIRE simulations. We show that gas supply, cooling/thermodynamics, star formation model, Toomre scale, galaxy dynamical times, and feedback properties do not have a direct causal effect on these transitions. Rather, both the formation of discs and cessation of bursty star formation are driven by the gravitational potential, but in different ways. Disc formation is promoted when the mass profile becomes sufficiently centrally concentrated in shape (relative to circularization radii): we show that this provides a well-defined dynamical centre, ceases to support the global ‘breathing modes’ that can persist indefinitely in less-concentrated profiles and efficiently destroy discs, promotes orbit mixing to form a coherent angular momentum, and stabilizes the disc. Smooth SF is promoted by the potential or escape velocity Vesc (not circular velocity Vc) becoming sufficiently large at the radii of star formation that cool, mass-loaded (momentum-conserving) outflows are trapped/confined near the galaxy, as opposed to escaping after bursts. We discuss the detailed physics, how these conditions arise in cosmological contexts, their relation to other correlated phenomena (e.g. inner halo virialization, vertical disc ‘settling’), and observations.

 

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.