Search for: All records

Context.Disc-dominated galaxies can be difficult to accommodate in a hierarchical formation scenario such as Λ cold dark matter (ΛCDM), where mergers are an important growth mechanism. However, observational evidence indicates that these galaxies are common in the Universe. Aims.We seek to characterise the conditions that lead to the formation of disc-dominated galaxies within ΛCDM. Methods.We used dynamical decomposition of the stellar particles in all galaxies with stellar massM_∗= [10¹⁰− 10¹¹] M_⊙within the cosmological hydrodynamical simulation Illustris TNG100. We selected a sample of 43 mostly-disc galaxies that have less than ∼10% of their mass in a bulge component. For comparison, we also studied two additional stellar-mass matched samples: 43 intermediate galaxies having ∼30% of their stellar mass in the bulge and 43 with a purely spheroidal-like morphology. Results.We find that the selection purely based on stellar dynamics is able to reproduce the expected stellar population trends of different morphological types, with higher star-formation rates and younger stars in disc-dominated galaxies. Halo spin seems to play no role in the morphology of the galaxies, in agreement with previous works. At a fixedM_*, our mostly-disc and intermediate samples form in dark matter haloes that are two to ten times less massive than the spheroidal sample, highlighting a higher efficiency in disc galaxies to retain and condensate their baryons. On average, mergers are less prevalent in the buildup of discs than in spheroidal galaxies, but there is a large scatter, including the existence of mostly-disc galaxies, with 15%–30% of their stars coming from accreted origin. Discs start to form early on, settling their low vertical velocity dispersion as early as 9–10 Gyr ago, although the dominance of the disc over the spheroid was established more recently (3–4 Gyr lookback time). The most rotationally supported discs form in haloes with the lowest virial mass in the sample and the best aligned distribution of angular momentum in the gas. 

Abstract We study the stellar properties of a sample of simulated ultradiffuse galaxies (UDGs) with stellar massM_⋆= 10^7.5–10⁹M_⊙, selected from the TNG50 simulation, where UDGs form mainly in high-spin dwarf-mass halos. We divide our sample into star-forming and quenched UDGs, finding good agreement with the stellar assembly history measured in observations. Star-forming UDGs and quenched UDGs withM_⋆≥ 10⁸M_⊙in our sample are particularly inefficient at forming stars, having 2–10 times less stellar mass than non-UDGs for the same virial mass halo. These results are consistent with recent mass inferences in UDG samples and suggest that the most inefficient UDGs arise from a late assembly of the dark matter mass followed by a stellar growth that is comparatively slower (for star-forming UDGs) or that was interrupted due to environmental removal of the gas (for quenched UDGs). Regardless of efficiency, UDGs are 60% poorer in [Fe/H] than the population of non-UDGs at a fixed stellar mass, with the most extreme objects having metal content consistent with the simulated mass–metallicity relation atz∼ 2. Quenched UDGs stop their star formation in shorter timescales than non-UDGs of similar mass and are, as a consequence, alpha enhanced with respect to non-UDGs. We identify metallicity profiles in UDGs as a potential avenue to distinguish between different formation paths for these galaxies, where gentle formation as a result of high-spin halos would present well-defined declining metallicity radial profiles while powerful-outflows or tidal stripping formation models would lead to flatter or constant metallicity as a function of radius due to the inherent mixing of stellar orbits. 

Abstract Simulations of galaxy formation are mostly unable to resolve the energy-conserving phase of individual supernova events, having to resort to subgrid models to distribute the energy and momentum resulting from stellar feedback. However, the properties of these simulated galaxies, including the morphology, stellar mass formed, and the burstiness of the star formation history, are highly sensitive to the numerical choices adopted in these subgrid models. Using the SMUGGLE stellar feedback model, we carry out idealized simulations of anM_vir∼ 10¹⁰M_⊙dwarf galaxy, a regime where most simulation codes predict significant burstiness in star formation, resulting in strong gas flows that lead to the formation of dark matter cores. We find that by varying only the directional distribution of momentum imparted from supernovae to the surrounding gas, while holding the total momentum per supernova constant, bursty star formation may be amplified or completely suppressed, and the total stellar mass formed can vary by as much as a factor of ∼3. In particular, when momentum is primarily directed perpendicular to the gas disk, less bursty and lower overall star formation rates result, yielding less gas turbulence, more disky morphologies, and a retention of cuspy dark matter density profiles. An improved understanding of the nonlinear coupling of stellar feedback into inhomogeneous gaseous media is thus needed to make robust predictions for stellar morphologies and dark matter core formation in dwarfs independent of uncertain numerical choices in the baryonic treatment. 

ABSTRACT We use the TNG50 from the IllustrisTNG suite of cosmological hydrodynamical simulation, complemented by a catalogue of tagged globular clusters, to investigate the properties and build up of two extended luminous components: the intra-cluster light (ICL) and the intra-cluster globular clusters (ICGCs). We select the 39 most massive groups and clusters in the box, spanning the range of virial masses $$5 \times 10^{12} \lt \rm M_{200}/\rm {\rm M}_{\odot } \lt 2 \times 10^{14}$$. We find good agreement between predictions from the simulations and current observational estimates of the fraction of mass in the ICL and its radial extension. The stellar mass of the ICL is only $$\sim 10~{{\ \rm per\ cent}}$$–20 per cent of the stellar mass in the central galaxy but encodes useful information on the assembly history of the group or cluster. About half the ICL in all our systems is brought in by galaxies in a narrow stellar mass range, M* = 1010–1011 M⊙. However, the contribution of low-mass galaxies (M* < 1010 M⊙) to the build up of the ICL varies broadly from system to system, $$\sim 5~{{\ \rm per\ cent}}-45~{{\ \rm per\ cent}}$$, a feature that might be recovered from the observable properties of the ICL at z = 0. At fixed virial mass, systems where the accretion of dwarf galaxies plays an important role have shallower metallicity profiles, less metal content, and a lower stellar mass in the ICL than systems where the main contributors are more massive galaxies. We show that intra-cluster GCs are also good tracers of this history, representing a valuable alternative when diffuse light is not detectable. 

We study the fraction of the intra-cluster light (ICL) formed in-situ in the three most massive clusters of the TNG50 simulation, with virial masses $\sim {10}^{14}$. We find that a significant fraction of ICL stars ( $8\%$- $28\%$) are born in-situ. This amounts to a total stellar mass comparable to the central galaxy itself. Contrary to simple expectations, only a sub-dominant fraction of these in-situ ICL stars are born in the central regions and later re-distributed to more energetic orbits during mergers. Instead, many in-situ ICL stars form directly hundreds of kiloparsecs away from the central galaxy, in clouds condensing out of the circum-cluster medium. The simulations predict a present-date diffuse star formation rate of $$1 ${\mathrm{M}}_{\odot }$/yr, with higher rates at higher redshifts. The diffuse star forming component of the ICL is filamentary in nature, extends for hundreds of kiloparsecs and traces the distribution of neutral gas in the cluster host halo. We discuss briefly how numerical details of the baryonic treatment in the simulation, in particular the density threshold for star formation and the equation of state, may play a role in this result. We conclude that a sensitivity of $1.6\times {10}^{-19}-2.6\times {10}^{-18}$erg s ${}^{-1}$cm ${}^{-2}$arcsec ${}^{-2}$in H ${}_{\alpha }$flux (beyond current observational capabilities) would be necessary to detect this diffuse star-forming component in galaxy clusters. 

ABSTRACT Stellar feedback plays a crucial role in regulating baryon cycles of a galactic ecosystem, and may manifest itself in the formation of superbubbles in the interstellar medium. In this work, we used a set of high-resolution simulations to systematically study the properties and evolution of superbubbles in galactic environments. The simulations were based on the SMUGGLE galaxy formation framework using the hydrodynamical moving-mesh code arepo, reaching a spatial resolution of $$\sim 4 \, \rm pc$$ and mass resolution of $$\sim 10^3 \, \rm M_{\odot }$$. We identified superbubbles and tracked their time evolution using the parent stellar associations within the bubbles. The X-ray luminosity-size distribution of superbubbles in the fiducial run is largely consistent with the observations of nearby galaxies. The size of superbubbles shows a double-peaked distribution, with the peaks attributed to early feedback (radiative and stellar wind feedback) and supernova feedback. The early feedback tends to suppress the subsequent supernova feedback, and it is strongly influenced by star formation efficiency, which regulates the environmental density. Our results show that the volume filling factor of hot gas (T > 105.5 K) is about $$12~{{\ \rm per\ cent}}$$ averaged over a region of 4 kpc in height and 20 kpc in radius centred on the disc of the galaxy. Overall, the properties of superbubbles are sensitive to the choice of subgrid galaxy formation models and can, therefore, be used to constrain these models. 

ABSTRACT In this study, we modify the semi-analytic model galacticus in order to accurately reproduce the observed properties of dwarf galaxies in the Milky Way. We find that reproducing observational determinations of the halo occupation fraction and mass–metallicity relation for dwarf galaxies requires us to include H2 cooling, an updated ultraviolet background radiation model, and to introduce a model for the metal content of the intergalactic medium. By fine-tuning various model parameters and incorporating empirical constraints, we have tailored the model to match the statistical properties of Milky Way dwarf galaxies, such as their luminosity function and size–mass relation. We have validated our modified semi-analytic framework by undertaking a comparative analysis of the resulting galaxy–halo connection. We predict a total of $$300 ^{+75} _{-99}$$ satellites with an absolute V-band magnitude (MV) less than 0 within 300 kpc from our Milky Way analogues. The fraction of subhaloes that host a galaxy at least this bright drops to 50 per cent by a halo peak mass of ∼8.9 × 107 M⊙, consistent with the occupation fraction inferred from the latest observations of Milky Way satellite population. 

ABSTRACT The velocity dispersion of globular clusters (GCs) around ultra-diffuse galaxies (UDGs) in the Virgo cluster spans a wide range, including cases where GC kinematics suggest haloes as massive as (or even more massive than) that of the Milky Way around these faint dwarfs. We analyse the catalogues of GCs derived in post-processing from the TNG50 cosmological simulation to study the GC system kinematics and abundance of simulated UDGs in galaxy groups and clusters. UDGs in this simulation reside exclusively in dwarf-mass haloes with M200 ≲ 1011.2 M⊙. When considering only GCs gravitationally bound to simulated UDGs, we find GCs properties that overlap well with several observational measurements for UDGs. In particular, no bias towards overly massive haloes is inferred from the study of bound GCs, confirming that GCs are good tracers of UDG halo mass. However, we find that contamination by intracluster GCs may, in some cases, substantially increase velocity dispersion estimates when performing projected mock observations of our sample. We caution that targets with less than 10 GC tracers are particularly prone to severe uncertainties. Measuring the stellar kinematics of the host galaxy should help confirm the unusually massive haloes suggested by GC kinematics around some UDGs. 

The study of dynamically cold stellar streams reveals information about the gravitational potential where they reside and provides important constraints on the properties of dark matter. However, the intrinsic faintness of these streams makes their detection beyond Local environments highly challenging. Here, we report the detection of an extremely faint stellar stream (μ_g, max= 29.5 mag arcsec⁻²) with an extraordinarily coherent and thin morphology in the Coma Galaxy Cluster. This Giant Coma Stream spans ∼510 kpc in length and appears as a free-floating structure located at a projected distance of 0.8 Mpc from the center of Coma. We do not identify any potential galaxy remnant or core, and the stream structure appears featureless in our data. We interpret the Giant Coma Stream as being a recently accreted, tidally disrupting passive dwarf. Using the Illustris-TNG50 simulation, we identify a case with similar characteristics, showing that, although rare, these types of streams are predicted to exist in Λ-CDM. Our work unveils the presence of free-floating, extremely faint and thin stellar streams in galaxy clusters, widening the environmental context in which these objects are found ahead of their promising future application in the study of the properties of dark matter. 

We present a new methodology for simulating mid-infrared emission from polycyclic aromatic hydrocarbons (PAHs) in galaxy evolution simulations. To do this, we combine theoretical models of PAH emission features as they respond to varying interstellar radiation fields, grain-size distributions, and ionization states with a new model for dust evolution in galaxy simulations. We apply these models to three idealized arepo galaxy evolution simulations within the smuggle physics framework. We use these simulations to develop numerical experiments investigating the buildup of PAH masses and luminosities in galaxies in idealized analogs of the Milky Way, a dwarf galaxy, and a starburst disk. Our main results are as follows. Galaxies with high specific star formation rates have increased feedback energy per unit mass, and are able to shatter grains efficiently, driving up the fraction of ultrasmall grains. At the same time, in our model large radiation fields per unit gas density convert aliphatic grains into aromatics. The fraction of dust grains in the form of PAHs (q_PAH) can be understood as a consequence of these processes, and in our model PAHs form primarily from interstellar processing (shattering) of larger grains rather than from the growth of smaller grains. We find that the hardness of the radiation field plays a larger role than variations in the grain-size distribution in setting the total integrated PAH luminosities, though cosmological simulations are necessary to investigate fully the complex interplay of processes that drive PAH band luminosities in galaxies.