Search

How mergers affect galaxy formation depends on both feedback processes, and on the geometry and strength of the mergers themselves. We introduce the PARADIGM project, where we study the response of a simulated Milky-Way-mass galaxy ($M_{\rm 200c} \sim 10^{12}M_{\odot}$ at $z=0$) forming in a cosmological setting to differing merger histories, using genetically modified initial conditions, each simulated with the VINTERGATAN and IllustrisTNG codes. While VINTERGATAN has been developed with an emphasis on resolving the cold interstellar medium, IllustrisTNG uses a subgrid two-phase model and consequently scales to large volume simulations, making them ideal to examine complementary views on how merger histories and feedback interact. Our genetic modifications alter the mass ratio of an important $z \approx 2$ merger while maintaining the halo's $z=0$ mass. Whether simulated with VINTERGATAN or IllustrisTNG, smaller mass ratios for this early merger result in larger galaxies at $z=0$, due to a greater build up of a kinematically cold disc. We conclude that such broad trends are robustly reproducible; however, the normalization of the resulting stellar sizes is substantially different in the two codes (ranging between $0.5-1.7\ \rm{kpc}$ for VINTERGATAN but $1.3-7.0\ \rm{kpc}$ for IllustrisTNG). The VINTERGATAN galaxies systematically form stars earlier, leading to a larger bulge component. Despite the difference in size normalization, both simulation suites lie on the observed size-mass relation for their respective morphological types. In light of these results, we discuss the interplay between internal processes and large scale gravitational interactions and gas accretion, and how the two galaxy models converge on similar emergent trends but along different evolutionary pathways., Comment: Accepted to MNRAS

We investigate the impact of a galaxy's merger history on its system of satellites using the new \textsc{vintergatan-gm} suite of zoom-in hydrodynamical simulations of Milky Way-mass systems. The suite simulates five realizations of the same halo with targeted `genetic modifications' (GMs) of a $z \approx 2$ merger, but resulting in the same halo mass at $z=0$. We find that differences in the satellite stellar mass functions last for $2.25-4.25$ Gyr after the $z \approx 2$ merger; specifically, the haloes that have undergone smaller mergers host up to 60\% more satellites than those of the larger merger scenarios. However, by $z=0$ these differences in the satellite stellar mass functions have been erased. The differences in satellite numbers seen soon after the mergers are driven by several factors, including the timings of significant mergers (with $M_{\rm 200c}$ mass ratios $>1:30$ and bringing in $M_{\rm *} \geq 10^{8}{\rm M}_{\odot}$ at infall), the masses and satellite populations of the central and merging systems, and the subsequent extended history of smaller mergers. The results persist when measured at fixed central stellar mass rather than fixed time, implying that a host's recent merger history can be a significant source of scatter when reconstructing its dynamical properties from its satellite population., Comment: Accepted to MNRAS

We present a new suite of cosmological zoom-in hydrodynamical ($\approx 20\, \mathrm{pc}$ spatial resolution) simulations of Milky-Way mass galaxies to study how a varying mass ratio for a Gaia-Sausage-Enceladus (GSE) progenitor impacts the $z=0$ chemodynamics of halo stars. Using the genetic modification approach, we create five cosmological histories for a Milky-Way-mass dark matter halo ($M_{200} \approx 10^{12} \, M_\mathrm{\odot}$), incrementally increasing the stellar mass ratio of a $z\approx2$ merger from 1:25 to 1:2, while fixing the galaxy's final dynamical, stellar mass and large-scale environment. We find markedly different morphologies at $z=0$ following this change in early history, with a growing merger resulting in increasingly compact and bulge-dominated galaxies. Despite this structural diversity, all galaxies show a radially-biased population of inner halo stars like the Milky-Way's GSE which, surprisingly, has a similar magnitude, age, $\rm [Fe/H]$ and $\rm [\alpha/Fe]$ distribution whether the $z\approx2$ merger is more minor or major. This arises because a smaller ex-situ population at $z\approx2$ is compensated by a larger population formed in an earlier merger-driven starburst whose contribution to the GES can grow dynamically over time, with both populations strongly overlapping in the $\rm [Fe/H]-\rm [\alpha/Fe]$ plane. Our study demonstrates that multiple high-redshift histories can lead to similar $z=0$ chemodynamical features in the halo, highlighting the need for additional constraints to distinguish them, and the importance of considering the full spectrum of progenitors when interpreting $z=0$ data to reconstruct our Galaxy's past., Comment: Matching accepted version after minor changes

We analyse the quenched fractions, gas content, and star formation histories of ~1200 satellite galaxies with $M_* \geq 5 \times 10^6~{\rm M}_\odot$ around 198 Milky Way- (MW) and Andromeda-like (M31) hosts in TNG50, the highest-resolution simulation of IllustrisTNG. Satellite quenched fractions are larger for smaller masses, for smaller distances to their host galaxy, and in the more massive M31-like compared to MW-like hosts. As satellites cross their host's virial radius, their gas content drops: most satellites within 300 kpc lack detectable gas reservoirs at $z=0$, unless they are massive like the Magellanic Clouds and M32. Nevertheless, their stellar assembly exhibits a large degree of diversity. On average, the cumulative star formation histories are more extended for brighter, more massive satellites with a later infall, and for those in less massive hosts. Based on these relationships, we can even infer infall periods for observedMWand M31 dwarfs: e.g. 0-4 Gyr ago for the Magellanic Clouds and Leo I, 4-8 and 0-2 Gyr ago for M32 and IC 10, respectively. Ram pressure stripping (in combination with tidal stripping) deprives TNG50 satellites of their gas reservoirs and ultimately quenches their star formation, even though only a few per cent of the present-day satellites around the 198 TNG50 MW/M31-like hosts appear as jellyfish. The typical time since quenching for currently quenched TNG50 satellites is $6.9^{+2.5}_{-3.3}~{\rm Gyr}$ ago. The TNG50 results are consistent with the quenched fractions and stellar assembly of observed MW and M31 satellites, however, satellites of the SAGA survey with $M_* \sim 10^{8-9}~{\rm M}_\odot$ exhibit lower quenched fractions than TNG50 and other, observed analogues., Comment: Accepted for publication in MNRAS

We present the cumulative star-formation histories (SFHs) of >15000 dwarf galaxies ($M_{*}=10^{7-10}M_{\odot}$) from the TNG50 run of the IllustrisTNG suite across a vast range of environments. The key factors determining the dwarfs' SFHs are their status as central or satellite and their stellar mass, with centrals and more massive dwarfs assembling their stellar mass at later times on average compared to satellites and lower mass dwarfs. The satellites (in hosts of total mass $M_{200c,\,host}=10^{12-14.3}M_{\odot}$) assembled 90% of their z=0 stellar mass ~$7.0_{-5.5}^{+3.3}$ Gyr ago, while the centrals did so only ~$1.0_{-0.5}^{+4.0}$ Gyr ago. TNG50 predicts a large diversity in SFHs for both centrals and satellites, so that the stacked cumulative SFHs are representative of the TNG50 dwarf populations only in an average sense and individual dwarfs can have significantly different cumulative SFHs. Satellite dwarfs with the highest stellar mass to host mass ratios have the latest stellar mass assembly. Satellites at fixed stellar and host halo mass, found closer to the cluster centre, or accreted at earlier times, show significantly earlier stellar mass assembly. These trends, as well as the shapes of the SFHs themselves, are a manifestation of the varying proportions within a given subsample of quenched vs. star-forming galaxies, which exhibit markedly distinct SFH shapes. We also find a subtle effect whereby satellite dwarfs in the most massive hosts at z=0 have higher SFRs at early times, well before final infall into their z=0 host, compared to a control sample of centrals mass-matched at the time of accretion. This suggests that the large-scale environment can have a mild effect even on future satellites by providing the conditions for enhanced SF at early epochs. Our results are useful theoretical predictions for comparison to future resolved-stellar-population observations., Comment: Submitted to MNRAS. Comments are welcome. The TNG50 data is now fully publicly available at https://www.tng-project.org

We use the IllustrisTNG simulations to show how the fractions of quenched galaxies vary across different environments and cosmic time, and to quantify the role AGN feedback and preprocessing play in quenching group and cluster satellites. At $z=0$, we select galaxies with $M_* = 10^{9-12} M_{\odot}$ residing within ($\leq R_{200c}$) groups and clusters of total host mass $M_{200c}=10^{13-15.2} M_{\odot}$. TNG predicts a quenched fraction of $\sim70-90\%$ (on average) for centrals and satellites $\gtrsim 10^{10.5} M_{\odot}$, regardless of host mass, cosmic time ($0\leq z\leq0.5$), clustercentric distance and time since infall in the $z=0$ host. Low-mass centrals ($\lesssim 10^{10} M_{\odot}$), instead, are rarely quenched unless they become members of groups ($10^{13-14} M_{\odot}$) or clusters ($\geq10^{14} M_{\odot}$), where the quenched fraction rises to $\sim80\%$. The fraction of low-mass passive galaxies is higher closer to the host center and for more massive hosts. The population of low-mass satellites accreted $\gtrsim$4-6 Gyr ago in massive hosts is almost entirely passive, thus suggesting an upper limit for the time needed for environmental quenching to occur. In fact, $\sim30\%$ of group and cluster satellites that are quenched at $z=0$ were already quenched before falling into their current host, and the bulk of them quenched as early as 4 to 10 billion years ago. For low-mass galaxies ($\lesssim10^{10-10.5}M_{\odot}$), this is due to preprocessing, whereby current satellites may have been members of other hosts, and hence have undergone environmental processes, before falling into their final host, this mechanism being more common and more effective for the purposes of quenching for satellites found today in more massive hosts. On the other hand, massive galaxies quench on their own and because of AGN feedback, regardless of whether they are centrals or satellites., Comment: 21 Pages, 11 figures, 2 tables, 1 appendix, published in MNRAS

We study the stellar morphological evolution of disc galaxies within clusters in the TNG50 and TNG100 runs from the IllustrisTNG simulation suite. We select satellites of masses $10^{9.7} \leq M_{*,z=0}/\text{M}_{\odot} \leq 10^{11.6}$ residing in clusters of masses $10^{14} \lesssim M_{\text{200c,z=0}}/\text{M}_{\odot} \leq 10^{14.6}$ at $z=0$ and that were discs at accretion according to a kinematic morphology indicator (the circularity fraction). These are traced from the time of accretion to $z=0$ and compared to a control sample of central galaxies mass-matched at accretion. Most cluster discs become non-discy by $z=0$, in stark contrast with the control discs, of which a significant fraction remains discy over the same timescales. Cluster discs become non-discy accompanied by gas removal and star formation quenching, loss of dark matter and little growth or a loss of stellar mass. In contrast, control discs transform while also losing gas mass and quenching, but growing significantly in dark matter and stellar mass. Most cluster satellites change morphologies on similar timescales regardless of stellar mass, in $\sim0.5-4$ Gyr after accretion. Cluster discs that experienced more numerous and closer pericentric passages show the largest change in morphology. Morphological change in all cases requires the presence of a gravitational perturbation to drive stellar orbits to non-discy configurations, along with gas removal/heating to prevent replenishment of the disc through continued star-formation. For cluster discs, the perturbation is impulsive tidal shocking at pericentres and not tidal stripping of outer disc stellar material, whereas for control discs, a combination of mergers and AGN feedback appears to be the key driving force behind morphological transformations., Comment: Accepted to MNRAS

We study the stellar-to-halo mass relation (SHMR) for central and satellite galaxies with total dynamical masses above 10^10.5 Msun using the suite of cosmological magneto-hydrodynamical simulations IllustrisTNG. In particular, we quantify environmental effects on satellite populations from TNG50, TNG100, and TNG300 located within the virial radius of group- and cluster-like hosts with total masses of 10^12-15.2 Msun. At fixed stellar mass, the satellite SHMR exhibits a distinct shift towards lower dynamical mass compared to the SHMR of centrals. Conversely, at fixed dynamical mass, satellite galaxies appear to have larger stellar-to-total mass fractions than centrals by up to a factor of a few. The systematic deviation from the central SHMR is larger for satellites in more massive hosts, at smaller cluster-centric distances, with earlier infall times, and that inhabit higher local density environments; moreover, it is in place already at early times (z < 2). Systematic environmental effects might contribute to the perceived galaxy-to-galaxy variation in the measured SHMR when galaxies cannot be separated into satellites and centrals. The SHMR of satellites exhibits a larger scatter than centrals, over the whole range of dynamical mass (by up to 0.8 dex). The shift of the satellite SHMR results mostly from tidal stripping of their dark matter, which affects satellites in an outside-in fashion: the departure of the satellite SHMR from the centrals' relation diminishes for measurements of dynamical mass in progressively smaller apertures. Finally, we provide a family of fitting functions for the SHMR predicted by IllustrisTNG., Comment: Published in MNRAS. Key figures: 1 & 8. Corrected Equation 1 for the SHMR fitting formula

We compile a sample of SDSS galaxy clusters with high-quality Chandra X-ray data to directly study the influence of the dense intra-cluster medium (ICM) on the quenching of satellite galaxies. We study the quenched fractions of satellite galaxies as a function of ICM density for low- ($10^9 \lesssim M_\star \lesssim 10^{10}\,\mathrm{M_\odot}$), intermediate- ($10^{10} \lesssim M_\star \lesssim 10^{10.5}\,\mathrm{M_\odot}$), and high-mass ($M_\star \gtrsim 10^{10.5}\,\mathrm{M_\odot}$) satellite galaxies with $>\!3000$ satellite galaxies across 24 low-redshift ($z < 0.1$) clusters. For low-mass galaxies we find evidence for a broken powerlaw trend between satellite quenched fraction and local ICM density. The quenched fraction increases modestly at ICM densities below a threshold before increasing sharply beyond this threshold toward the cluster center. We show that this increase in quenched fraction at high ICM density is well matched by a simple, analytic model of ram pressure stripping. These results are consistent with a picture where low-mass cluster galaxies experience an initial, slow-quenching mode driven by steady gas depletion, followed by rapid quenching associated with ram pressure of cold-gas stripping near (one quarter of the virial radius, on average) the cluster center., Comment: 17 pages plus appendix, 10 figures, accepted to ApJ

We present a study of environmental effects and preprocessing in a large galaxy group using a high-resolution, zoom-in simulation run with the GASOLINE2 hydrodynamics code. We categorize galaxies that were always in distinct haloes as unaccreted, galaxies that were distinct before accretion onto the main group as single, and galaxies that were in external sub-groups before accretion onto the main group as grouped. The unaccreted galaxy population experiences steady growth in dark matter, gas and stellar mass. Both single- and group-accreted galaxies begin to lose dark matter and gas after first accretion onto any host but continue to grow in stellar mass. Individual trajectories show that galaxies cease mass growth within roughly three virial radii of the main group. Single galaxies continue to form stars until the group virial radius is crossed, when they begin to lose both dark matter and gas. Grouped galaxies peak in mass when joining their external sub-group, indicating that they experience preprocessing. Most accreted galaxies retain their accumulated stellar mass. The total mass loss is dominated by tidal stripping, with evidence for additional gas stripping via ram pressure. Most accreted galaxies are quenched $\sim$(0.5-2.5) Gyr after accretion onto any group. These differing histories place unaccreted, single and grouped galaxies in distinct regions of the stellar mass-to-halo mass (SMHM) relation. This suggests that preprocessed galaxies are a key source of scatter in the SMHM relation for mixed galaxy populations., Comment: Accepted for publication in MNRAS, 16 pages, 15 figures

We investigate the mass loss of galaxies in groups and clusters with high-resolution DM simulations. We detect weak mass segregation in the inner regions of group/cluster haloes, consistent with observational findings. This applies to samples of galaxy analogues selected using either their present-day mass or past maximum (peak) mass. We find a strong radial trend in the fractional mass lost by the galaxies since peak, independent of their mass. This suggests that segregation is due to massive galaxies having formed closer to the halo centres and not the preferential destruction of smaller galaxies near halo centres. We divide our sample into galaxies that were accreted as a group vs. as a single, distinct halo. We find strong evidence for preprocessing -- the grouped galaxies lose $\sim 35-45\%$ of their peak mass before being accreted onto their final host haloes, compared to single galaxies which lose $\sim12\%$. After accretion, however, the single galaxies lose more mass compared to the grouped ones. These results are consistent with a scenario in which grouped galaxies are preprocessed in smaller haloes while single galaxies `catch up' in terms of total mass loss once they are accreted onto the final host halo. The fractional mass loss is mostly independent of the galaxy mass and host mass, and increases with amount of time spent in a dense environment., Comment: 11 pages, 7 figures, accepted for publication in MNRAS

We investigate mass segregation in group and cluster environments by identifying galaxy analogues in high-resolution dark matter simulations. Subhalos identified by the AHF and ROCKSTAR halo finders have similar mass functions, independent of resolution, but different radial distributions due to significantly different subhalo hierarchies. We propose a simple way to classify subhalos as galaxy analogues. The radial distributions of galaxy analogues agree well at large halo-centric radii for both AHF and ROCKSTAR but disagree near parent halo centres where the phase-space information used by ROCKSTAR is essential. We see clear mass segregation at small radii (within $0.5\,r_{vir}$) with average galaxy analogue mass decreasing with radius. Beyond the virial radius, we find a mild trend where the average galaxy analogue mass increases with radius. These mass segregation trends are strongest in small groups and dominated by the segregation of low mass analogues. The lack of mass segregation in massive galaxy analogues suggests that the observed trends are driven by the complex accretion histories of the parent halos rather than dynamical friction., Comment: 19 pages, 13 figures, accepted for publication in MNRAS

It has been shown that galaxy properties depend strongly on their host environment. In order to understand the relevant physical processes driving galaxy evolution it is important to study the observed properties of galaxies in different environments. Mass segregation in bound galaxy structures is an important indicator of evolutionary history and dynamical friction timescales. Using group catalogues derived from the Sloan Digital Sky Survey Data Release 7 (SDSS DR7) we investigate mass segregation trends in galaxy groups at low redshift. We investigate average galaxy stellar mass as a function of group-centric radius and find evidence for weak mass segregation in SDSS groups. The magnitude of the mass segregation depends on both galaxy stellar mass limits and group halo mass. We show that the inclusion of low mass galaxies tends to strengthen mass segregation trends, and that the strength of mass segregation tends to decrease with increasing group halo mass. We find the same trends if we use the fraction of massive galaxies as a function of group-centric radius as an alternative probe of mass segregation. The magnitude of mass segregation that we measure, particularly in high-mass haloes, indicates that dynamical friction is not acting efficiently., Comment: 6 pages, 2 figures, accepted for publication in MNRAS Letters