Search for: All records

We study the formation of ultradiffuse galaxies (UDGs) using the cosmological hydrodynamical simulation TNG50 of the Illustris-TNG suite. We define UDGs as dwarf galaxies in the stellar mass range $\rm {7.5 \le log (M_{\star } / {\rm M}_{\odot }) \le 9 }$ that are in the 5 per cent most extended tail of the simulated mass–size relation. This results in a sample of UDGs with half-mass radii $\rm {r_{h \star } \gtrsim 2 \ kpc}$ and surface brightness between $\rm {24.5}$ and $\rm {28 \ mag \ arcsec^{-2}}$, similar to definitions of UDGs in observations. The large cosmological volume in TNG50 allows for a comparison of UDGs properties in different environments, from the field to galaxy clusters with virial mass $\rm {M_{200} \sim 2 \times 10^{14} ~ {\rm M}_{\odot }}$. All UDGs in our sample have dwarf-mass haloes ($\rm {M_{200}\sim 10^{11} ~ {\rm M}_{\odot } }$) and show the same environmental trends as normal dwarfs: field UDGs are star-forming and blue while satellite UDGs are typically quiescent and red. The TNG50 simulation predicts UDGs that populate preferentially higher spin haloes and more massive haloes at fixed $\rm {M_{\star }}$ compared to non-UDG dwarfs. This applies also to most satellite UDGs, which are actually ‘born’ UDGs in the field and infall into groups and clusters without significant changes to their size. We find, however, a small subset of satellite UDGs ($\lesssim 10~{{\ \rm per\ cent}}$) with present-day stellar size a factor ≥1.5 larger than at infall, confirming that tidal effects, particularly in the lower mass dwarfs, are also a viable formation mechanism for some of these dwarfs, although sub-dominant in this simulation.

The tension between the diverging density profiles in Lambda cold dark matter simulations and the constant-density inner regions of observed galaxies is a long-standing challenge known as the ‘core–cusp’ problem. We demonstrate that the SMUGGLE galaxy formation model implemented in the arepo moving mesh code forms constant-density cores in idealized dwarf galaxies of M⋆ ≈ 8 × 107 Msun with initially cuspy dark matter (DM) haloes of M200 ≈ 1010 Msun. Identical initial conditions run with an effective equation of state interstellar medium model preserve cuspiness. Literature on the subject has pointed to the low density threshold for star formation, ρth, in such effective models as an obstacle to baryon-induced core formation. Using a SMUGGLE run with equal ρth, we demonstrate that core formation can proceed at low density thresholds, indicating that ρth is insufficient on its own to determine whether a galaxy develops a core. We reaffirm that the ability to resolve a multiphase interstellar medium at sufficiently high densities is a more reliable indicator of core formation than any individual model parameter. In SMUGGLE, core formation is accompanied by large degrees of non-circular motion, with gas rotational velocity profiles that consistently fall below the circular velocity $v_\text{circ} = \sqrt{GM/R}$ out to ∼2 kpc. Asymmetric drift corrections help recover the average underlying DM potential for some of our less efficient feedback runs, but time-variations in the instantaneous azimuthal gas velocity component are substantial, highlighting the need for careful modelling in the inner regions of dwarfs to infer the true distribution of DM.

We study gas inflows on to supermassive black holes using hydrodynamics simulations of isolated galaxies and idealized galaxy mergers with an explicit, multiphase interstellar medium (ISM). Our simulations use the recently developed ISM and stellar evolution model called Stars and MUltiphase Gas in GaLaxiEs (SMUGGLE). We implement a novel super-Lagrangian refinement scheme that increases the gas mass resolution in the immediate neighbourhood of the black holes (BHs) to accurately resolve gas accretion. We do not include black hole feedback in our simulations. We find that the complex and turbulent nature of the SMUGGLE ISM leads to highly variable BH accretion. BH growth in SMUGGLE converges at gas mass resolutions ≲3 × 103 M⊙. We show that the low resolution simulations combined with the super-Lagrangian refinement scheme are able to produce central gas dynamics and BH accretion rates very similar to that of the uniform high resolution simulations. We further explore BH fueling by simulating galaxy mergers. The interaction between the galaxies causes an inflow of gas towards the galactic centres and results in elevated and bursty star formation. The peak gas densities near the BHs increase by orders of magnitude resulting in enhanced accretion. Our results support the idea that galaxy mergers can trigger AGN activity, although the instantaneous accretion rate depends strongly on the local ISM. We also show that the level of merger-induced enhancement of BH fueling predicted by the SMUGGLE model is much smaller compared to the predictions by simulations using an effective equation of state model of the ISM.

We present a post-processing catalogue of globular clusters (GCs) for the 39 most massive groups and clusters in the TNG50 simulation of the IlllustrisTNG project (virial masses $M_{200} =[5\times 10^{12} \rm {\!-\!} 2 \times 10^{14}$] M⊙). We tag GC particles to all galaxies with stellar mass M* ≥ 5 × 106 M⊙, and we calibrate their masses to reproduce the observed power-law relation between GC mass and halo mass for galaxies with M200 ≥ 1011 M⊙ (corresponding to M* ∼ 109 M⊙). Here, we explore whether an extrapolation of this MGC–M200 relation to lower mass dwarfs is consistent with current observations. We find a good agreement between our predicted number and specific frequency of GCs in dwarfs with $\rm {\it M}_*=[5 \times 10^6 \rm {\!-\!} 10^9]$ M⊙ and observations. Moreover, we predict a steep decline in the GC occupation fraction for dwarfs with M* < 109 M⊙ that agrees well with current observational constraints. This declining occupation fraction is due to a combination of tidal stripping in all dwarfs plus a stochastic sampling of the GC mass function for dwarfs with M* < 107.5 M⊙. Our simulations also reproduce available constraints on the abundance of intracluster GCs in Virgo and Centaurus A. These successes provide support to the hypothesis that the MGC–M200 relation holds, albeit with more scatter, all the way down to the regime of classical dwarf spheroidals in these environments. Our GC catalogues are publicly available as part of the IllustrisTNG data release.

The abundance of the faintest galaxies provides insight into the nature of dark matter and the process of dwarf galaxy formation. In the LCDM scenario, low-mass haloes are so numerous that the efficiency of dwarf formation must decline sharply with decreasing halo mass in order to accommodate the relative scarcity of observed dwarfs and satellites in the Local Group. The nature of this decline contains important clues to the mechanisms regulating the onset of galaxy formation in the faintest systems. We explore here two possible models for the stellar mass (M*)–halo mass (M200) relation at the faint end, motivated by some of the latest LCDM cosmological hydrodynamical simulations. One model includes a sharp mass threshold below which no luminous galaxies form, as expected if galaxy formation proceeds only in systems above the hydrogen-cooling limit. In the second model, M* scales as a steep power law of M200 with no explicit cut-off, as suggested by recent semi-analytical work. Although both models predict satellite numbers around Milky Way-like galaxies consistent with current observations, they predict vastly different numbers of ultrafaint dwarfs and of satellites around isolated dwarf galaxies. Our results illustrate how the satellite mass function around dwarfs may be used to probe the M*–M200 relation at the faint end and to elucidate the mechanisms that determine which low-mass haloes ‘light up’ or remain dark in the LCDM scenario.

While it is generally believed that supermassive black holes (SMBHs) lie in most galaxies with bulges, few SMBHs have been confirmed in bulgeless galaxies. Identifying such a population could provide important insights to the BH seed population and secular BH growth. To this end, we obtained near-infrared (NIR) spectroscopic observations of a sample of low-redshift bulgeless galaxies with mid-infrared colors suggestive of active galactic nuclei (AGNs). We find additional evidence of AGN activity (such as coronal lines and broad permitted lines) in 69% (9/13) of the sample, demonstrating that mid-infrared selection is a powerful tool to detect AGNs. More than half of the galaxies with confirmed AGN activity show fast outflows in [Oiii] in the optical and/or [Sivi] in the NIR, with the latter generally having much faster velocities that are also correlated to their spatial extent. We are also able to obtain virial BH masses for some targets and find they fall within the scatter of other late-type galaxies in theM_BH–M_stellarrelation. The fact that they lack a significant bulge component indicates that secular processes, likely independent of major mergers, grew these BHs to supermassive sizes. Finally, we analyze the rotational gas kinematics and find two notable exceptions: two AGN hosts with outflows that appear to be rotating faster than expected. There is an indication that these two galaxies have stellar masses significantly lower than expected from their dark matter halo masses. This, combined with the observed AGN activity and strong gas outflows, may be evidence of the effects of AGN feedback.

Hydrogen emission lines can provide extensive information about star-forming galaxies in both the local and high-redshift Universe. We present a detailed Lyman continuum (LyC), Lyman-α (Lyα), and Balmer line (Hα and Hβ) radiative transfer study of a high-resolution isolated Milky Way simulation using the state-of-the-art Arepo-RT radiation hydrodynamics code with the SMUGGLE galaxy formation model. The realistic framework includes stellar feedback, non-equilibrium thermochemistry accounting for molecular hydrogen, and dust grain evolution in the interstellar medium (ISM). We extend our publicly available Cosmic Lyα Transfer (COLT) code with photoionization equilibrium Monte Carlo radiative transfer and various methodology improvements for self-consistent end-to-end (non-)resonant line predictions. Accurate LyC reprocessing to recombination emission requires modelling pre-absorption by dust ($f_\text{abs} \approx 27.5\,\rm{per\,\,cent}$), helium ionization ($f_\text{He} \approx 8.7\,\rm{per\,\,cent}$), and anisotropic escape fractions ($f_\text{esc} \approx 7.9\,\rm{per\,\,cent}$), as these reduce the available budget for hydrogen line emission ($f_\text{H} \approx 55.9\,\rm{per\,\,cent}$). We investigate the role of the multiphase dusty ISM, disc geometry, gas kinematics, and star formation activity in governing the physics of emission and escape, focusing on the time variability, gas-phase structure, and spatial spectral, and viewing angle dependence of the emergent photons. Isolated disc simulations are well-suited for comprehensive observational comparisons with local Hα surveys, but would require a proper cosmological circumgalactic medium (CGM) environment as well as less dust absorption and rotational broadening to serve as analogs for high-redshift Lyα emitting galaxies. Future applications of our framework to next-generation cosmological simulations of galaxy formation including radiation-hydrodynamics that resolve ≲10 pc multiphase ISM and ≲1 kpc CGM structures will provide crucial insights and predictions for current and upcoming Lyα observations.

As the Milky Way and its satellite system become more entrenched in near field cosmology efforts, the need for an accurate mass estimate of the Milky Way’s dark matter halo is increasingly critical. With the second and early third data releases of stellar proper motions from Gaia, several groups calculated full 6D phase-space information for the population of Milky Way satellite galaxies. Utilizing these data in comparison to subhalo properties drawn from the Phat ELVIS simulations, we constrain the Milky Way dark matter halo mass to be ∼1–1.2 × 1012 M⊙. We find that the kinematics of subhaloes drawn from more- or less-massive hosts (i.e. >1.2 × 1012 M⊙ or <1012 M⊙) are inconsistent, at the 3σ confidence level, with the observed velocities of the Milky Way satellites. The preferred host halo mass for the Milky Way is largely insensitive to the exclusion of systems associated with the Large Magellanic Cloud, changes in galaxy formation thresholds, and variations in observational completeness. As more Milky Way satellites are discovered, their velocities (radial, tangential, and total) plus Galactocentric distances will provide further insight into the mass of the Milky Way dark matter halo.

We present a suite of 16 high-resolution hydrodynamic simulations of an isolated dwarf galaxy (gaseous and stellar disc plus a stellar bulge) within an initially cuspy dark matter (DM) halo, including self-interactions between the DM particles; as well as stochastic star formation and subsequent supernova feedback (SNF), implemented using the stellar feedback model SMUGGLE. The simulations start from identical initial conditions, and we regulate the strength of DM self-interactions and SNF by systematically varying the self-interacting DM (SIDM) momentum transfer cross-section and the gas density threshold for star formation. The DM halo forms a constant density core of similar size and shape for several combinations of those two parameters. Haloes with cores that are formed due to SIDM (adiabatic cusp-core transformation) have velocity dispersion profiles that are closer to isothermal than those of haloes with cores that are formed due to SNF in simulations with bursty star formation (impulsive cusp-core transformation). Impulsive SNF can generate positive stellar age gradients and increase random motion in the gas at the centre of the galaxy. Simulated galaxies in haloes with cores that were formed adiabatically are spatially more extended, with stellar metallicity gradients that are shallower (at late times) than those of galaxies in other simulations. Such observable properties of the gas and the stars, which indicate either an adiabatic or an impulsive evolution of the gravitational potential, may be used to determine whether observed cores in DM haloes are formed through DM self-interactions or in response to impulsive SNF.

Galaxy mergers are known to host abundant young massive cluster (YMC) populations, whose formation mechanism is still not well-understood. Here, we present a high-resolution galaxy merger simulation with explicit star formation and stellar feedback prescriptions to investigate how mergers affect the properties of the interstellar medium and YMCs. Compared with a controlled simulation of an isolated galaxy, the mass fraction of dense and high-pressure gas is much higher in mergers. Consequently, the mass function of both molecular clouds and YMCs becomes shallower and extends to higher masses. Moreover, cluster formation efficiency is significantly enhanced and correlates positively with the star formation rate surface density and gas pressure. We track the orbits of YMCs and investigate the time evolution of tidal fields during the course of the merger. At an early stage of the merger, the tidal field strength correlates positively with YMC mass, λtid ∝ M0.71, which systematically affects the shape of the mass function and age distribution of the YMCs. At later times, most YMCs closely follow the orbits of their host galaxies, gradually sinking into the centre of the merger remnant due to dynamical friction, and are quickly dissolved via efficient tidal disruption. Interestingly, YMCs formed during the first passage, mostly in tidal tails and bridges, are distributed over a wide range of galactocentric radii, greatly increasing their survivability because of the much weaker tidal field in the outskirts of the merger system. These YMCs are promising candidates for globular clusters that survive to the present day.