Cosmological Hydrodynamical Simulations Research Articles

Simulating the LOcal Web (SLOW)

Context. This is the second paper in a series presenting the results from a 500 h−1Mpc large constrained simulation of the local Universe (SLOW). The initial conditions for this cosmological hydro-dynamical simulation are based on peculiar velocities derived from the CosmicFlows-2 catalog. The simulation follows cooling, star formation, and the evolution of super-massive black holes. This allows one to directly predict observable properties of the intracluster medium (ICM) within galaxy clusters, including X-ray luminosity, temperatures, and the Compton-y signal. Aims. Comparing the properties of observed galaxy clusters within the local Universe with the properties of their simulated counterparts enables us to assess the effectiveness of the initial condition constraints in accurately replicating the mildly nonlinear properties of the largest, collapsed objects within the simulation. Methods. Based on the combination of several, publicly available surveys we compiled a sample of galaxy clusters within the local Universe, of which we were able to cross-identify 46 of them with an associated counterpart within the SLOW simulation. We then derived the probability of the cross identification based on mass, X-ray luminosity, temperature, and Compton-y by comparing it to a random selection. Results. Our set of 46 cross-identified local Universe clusters contains the 13 most massive clusters from the Planck SZ catalog as well as 70% of clusters with M500 larger than 2 × 1014 M⊙. Compared to previous constrained simulations of the local volume, we found in SLOW a much larger amount of replicated galaxy clusters, where their simulation-based mass prediction falls within the uncertainties of the observational mass estimates. Comparing the median observed and simulated masses of our cross-identified sample allows us to independently deduce a hydrostatic mass bias of (1 − b)≈0.87. Conclusions. The SLOW constrained simulation of the local Universe faithfully reproduces numerous fundamental characteristics of a sizable number of galaxy clusters within our local neighborhood, opening a new avenue for studying the formation and evolution of a large set of individual galaxy clusters as well as testing our understanding of physical processes governing the ICM.

Helium reionization from empirical quasar luminosity functions before and after JWST

ABSTRACT Recently, models of the quasar luminosity function (QLF) rooted on large observational compilations have been produced that, unlike their predecessors, feature a smooth evolution with time. This bypasses the need to assume an ionizing emissivity evolution when simulating helium reionization with observations-based QLF, thus yielding more robust constraints. We combine one such QLF with a cosmological hydrodynamical simulation and 3D multifrequency radiative transfer. The simulated reionization history is consistently delayed in comparison to most other models in the literature. The predicted intergalactic medium temperature is larger than the observed one at $z \lesssim 3$. Through forward modelling of the He ii Lyman $\alpha$ forest, we show that our model produces an extended helium reionization and successfully matches the bulk of the observed effective optical depth distribution, although it overionizes the Universe at $z\lesssim 2.8$ as the effect of small-scale Lyman Limit Systems not being resolved. We thoroughly characterize transmission regions and dark gaps in He ii Lyman $\alpha$ forest sightlines. We quantify their sensitivity to the helium reionization, opening a new avenue for further observational studies of this epoch. Finally, we explore the implications for helium reionization of the large number of active galactic nuclei revealed at $z\gtrsim 5$ by JWST. We find that such modifications do not affect any observable at $z\le 4$, except in our most extreme model, indicating that the observed abundance of high-z AGNs does not bear consequences for helium reionization.

Tracing the evolutionary pathways of dust and cold gas in high-z quiescent galaxies with SIMBA

Recent discoveries of copious amounts of dust in quiescent galaxies (QGs) at high redshifts ($z 1-2$) challenge the conventional view that these objects have a negligible interstellar medium (ISM) in proportion to their stellar mass. We made use of the SIMBA hydrodynamic cosmological simulation to explore how dust and cold gas evolve in QGs and are linked to the quenching processes affecting them. We applied a novel method for tracking the changes in the ISM dust abundance across the evolutionary history of QGs identified at $0 < z in both cluster and field environments. The QGs transition from a diversity of quenching pathways, both rapidly and slowly, and they exhibit a wide range of times that elapsed between the quenching event and cold gas removal (from $ Myr to $ Gyr). Contrary to some claims, we find that quenching modes attributed to the feedback from active galactic nuclei (AGNs) do not affect dust and cold gas within the same timescales. Remarkably, QGs may replenish their dust content in the quenched phase primarily due to internal processes and marginally by external factors such as minor mergers. Prolonged grain growth on gas-phase metals appears to be the key mechanism for dust re-formation, which is effective within $ Myr after the quenching event and rapidly increases the dust-to-gas mass ratio in QGs above the standard values ($ DGR Consequently, despite heavily depleted cold gas reservoirs, roughly half of QGs maintain little evolution of their ISM dust with stellar age within the first 2 Gyr following the quenching. Overall, we predict that relatively dusty QGs dust /M_ star $) arise from both fast and slow quenchers, and they are prevalent in quenched systems of intermediate and low stellar masses ($9< star /M_ odot )<10.5$). This strong prediction poses an immediate quest for observational synergy between, for example, the James Webb Space Telescope (JWST) and the Atacama Large Millimetre Array (ALMA).

Applying machine learning to Galactic Archaeology: how well can we recover the origin of stars in Milky Way-like galaxies?

ABSTRACT We present several machine learning (ML) models developed to efficiently separate stars formed in situ in Milky Way-type galaxies from those that were formed externally and later accreted. These models, which include examples from artificial neural networks, decision trees, and dimensionality reduction techniques, are trained on a sample of disc-like, Milky Way-mass galaxies drawn from the artemis cosmological hydrodynamical zoom-in simulations. We find that the input parameters which provide an optimal performance for these models consist of a combination of stellar positions, kinematics, chemical abundances ([Fe/H] and [α/Fe]), and photometric properties. Models from all categories perform similarly well, with area under the precision–recall curve (PR-AUC) scores of ≃ 0.6. Beyond a galactocentric radius of 5 kpc, models retrieve $\gt 90~{{\ \rm per\ cent}}$ of accreted stars, with a sample purity close to 60 per cent, however the purity can be increased by adjusting the classification threshold. For one model, we also include host galaxy-specific properties in the training, to account for the variability of accretion histories of the hosts, however this does not lead to an improvement in performance. The ML models can identify accreted stars even in regions heavily dominated by the in-situ component (e.g. in the disc), and perform well on an unseen suite of simulations (the auriga simulations). The general applicability bodes well for application of such methods on observational data to identify accreted substructures in the Milky Way without the need to resort to selection cuts for minimizing the contamination from in-situ stars.

The Environmental Dependence of the Stellar Mass - Gas Metallicity Relation in Horizon Run 5

Abstract Metallicity offers a unique window into the baryonic history of the cosmos, being instrumental in probing evolutionary processes in galaxies between different cosmic environments. We aim to quantify the contribution of these environments to the scatter in the mass-metallicity relation (MZR) of galaxies. By analysing the galaxy distribution within the cosmic skeleton of the Horizon Run 5 cosmological hydrodynamical simulation at redshift z = 0.625, computed using a careful calibration of the T-ReX filament finder, we identify galaxies within three main environments: nodes, filaments and voids. We also classify galaxies based on the dynamical state of the clusters and the length of the filaments in which they reside. We find that the cosmic environment significantly contributes to the scatter in the MZR; in particular, both the gas metallicity and its average relative standard deviation increase when considering denser large-scale environments. The difference in the average metallicity between galaxies within relaxed and unrelaxed clusters is ≈0.1dex, with both populations displaying positive residuals, δZg, from the averaged MZR. Moreover, the difference in metallicity between node and void galaxies accounts for ≈0.14 dex in the scatter of the MZR at stellar mass M⋆ ≈ 109.35 M⊙. Finally, both the average [O/Fe] in the gas and the galaxy gas fraction decrease when moving to higher large-scale densities in the simulation, suggesting that the cores of cosmic environments host – on average – older and more massive galaxies, whose enrichment is affected by a larger number of Type Ia Supernova events.

On the Origin of the Variety of Velocity Dispersion Profiles of Galaxies

Observed and simulated galaxies exhibit a significant variation in their velocity dispersion profiles. We examine the inner and outer slopes of stellar velocity dispersion profiles using integral field spectroscopy data from two surveys, SAMI (for z < 0.115) and CALIFA (for z < 0.03), comparing them with results from two cosmological hydrodynamic simulations: Horizon-AGN (for z = 0.017) and NewHorizon (for z ≲ 1). The simulated galaxies closely reproduce the variety of velocity dispersion slopes and stellar mass dependence of both inner and outer radii (0.5 r 50 and 3 r 50) as observed, where r 50 stands for half-light radius. The inner slopes are mainly influenced by the relative radial distribution of the young and old stars formed in situ: a younger center shows a flatter inner profile. The presence of accreted (ex situ) stars has two effects on the velocity dispersion profiles. First, because they are more dispersed in spatial and velocity distributions compared to in situ formed stars, it increases the outer slope of the velocity dispersion profile. It also causes the velocity anisotropy to be more radial. More massive galaxies have a higher fraction of stars formed ex situ and hence show a higher slope in outer velocity dispersion profile and a higher degree of radial anisotropy. The diversity in the outer velocity dispersion profiles reflects the diverse assembly histories among galaxies.

Zooming by in the CARPoolGP Lane: New CAMELS-TNG Simulations of Zoomed-in Massive Halos

Galaxy formation models within cosmological hydrodynamical simulations contain numerous parameters with nontrivial influences over the resulting properties of simulated cosmic structures and galaxy populations. It is computationally challenging to sample these high dimensional parameter spaces with simulations, in particular for halos in the high-mass end of the mass function. In this work, we develop a novel sampling and reduced variance regression method, CARPoolGP, which leverages built-in correlations between samples in different locations of high dimensional parameter spaces to provide an efficient way to explore parameter space and generate low-variance emulations of summary statistics. We use this method to extend the Cosmology and Astrophysics with machinE Learning Simulations to include a set of 768 zoom-in simulations of halos in the mass range of 1013–1014.5 M ⊙ h −1 that span a 28-dimensional parameter space in the IllustrisTNG model. With these simulations and the CARPoolGP emulation method, we explore parameter trends in the Compton Y–M, black hole mass–halo mass, and metallicity–mass relations, as well as thermodynamic profiles and quenched fractions of satellite galaxies. We use these emulations to provide a physical picture of the complex interplay between supernova and active galactic nuclei feedback. We then use emulations of the Y–M relation of massive halos to perform Fisher forecasts on astrophysical parameters for future Sunyaev–Zeldovich observations and find a significant improvement in forecasted constraints. We publicly release both the simulation suite and CARPoolGP software package.

A stream come true: Connecting tidal tails, shells, streams, and planes with galaxy kinematics and formation history

Context. The rapidly improving quality and resolution of both low surface brightness observations and cosmological simulations of galaxies enable us to address the important question of how the formation history is imprinted in the outer unrelaxed regions of galaxies, and to inspect the correlations of these imprints with another tracer of galaxy formation, the internal kinematics. Aims. Using the hydrodynamical cosmological simulation called Magneticum Pathfinder, we identified tidal tails, shells, streams, and satellite planes, and connected them to the amount of rotational support and the formation histories of the host galaxies. This presents the first combined statistical census considering all these four types of features in hydrodynamical cosmological simulations. Methods. Tidal features were visually classified from a three-dimensional rendering of the simulated galaxies by several scientists independently. Only features that were identified by at least half of the participating individuals were considered to be existing features. The data on satellite planes and kinematic properties of the simulated galaxies were taken from previous work. The results were compared to observations, especially from the MATLAS survey. Results. Generally, prominent features are much more common around elliptical than around disk galaxies. Shells are preferentially found around kinematically slowly rotating galaxies in both simulations and observations, while streams can be found around all types of galaxies, with a slightly higher probability to be present around less rotationally supported galaxies. Tails and satellite planes, however, appear independently of the internal kinematics of the central galaxy, indicating that they are formed through processes that have not (yet) affected the internal kinematics. Prolate rotators have the overall highest probability to exhibit tidal features, but the highest likelihood for a specific type of feature is found for galaxies with kinematically distinct cores (KDCs), nearly 20% of which exhibit streams. Conclusions. As shells are formed through radial merger events while streams are remnants of circular merger infall, this suggests that the orbital angular momentum of the merger event plays a more crucial role in transforming the host galaxy than previously anticipated. The existence of a shell around a given slow rotator furthermore is a sign of a radial merger formation for this particular slow rotator because one-third of the galaxies with a shell were transformed into slow rotators by the merger event that also caused the shells. The appearance of a stream around a KDC is a direct indicator for the multiple merger formation pathway of that KDC as opposed to the major merger pathway.

Escape Velocity Mass of A1063

We measure the radius–velocity phase-space edge profile of A1063 using galaxy redshifts from Karman et al. and Mercurio et al. Combined with a cosmological model and after accounting for interlopers and sampling effects, we infer the escape velocity profile. Using the Poisson equation, we then directly constrain the gravitational potential profile and find excellent agreement between three different density models. For the Navarro–Frenk–White profile, we find log10(M 200,crit) = 15.40−0.12+0.06 M ⊙, consistent to within 1σ of six recently published lensing masses. We argue that this consistency is due to the fact that the escape technique shares no common systematics with lensing other than radial binning. These masses are 2–4σ lower than estimates made using X-ray data, in addition to the Gómezet al. velocity dispersion estimate. We measure the 1D velocity dispersion within r 200 to be σv=1477−99+87 km s−1, which combined with our escape velocity mass, brings the dispersion for A1063 in-line with hydrodynamic cosmological simulations for the first time.

Effectiveness of halo and galaxy properties in reducing the scatter in the stellar-to-halo mass relation

ABSTRACT The stellar-to-halo mass relation (SHMR) is a fundamental relationship between galaxies and their host dark matter haloes. In this study, we examine the scatter in this relation for primary galaxies in the semi-analytic l-galaxies model and two cosmological hydrodynamical simulations, EAGLE (Evolution and Assembly of Galaxies and their Environments) and TNG100-1. We find that in low-mass haloes, more massive galaxies tend to reside in haloes with higher concentration, earlier formation time, greater environmental density, earlier major mergers, and, to have older stellar populations, which is consistent with findings in various studies. Quantitative analysis reveals the varying significance of halo and galaxy properties in determining SHMR scatter across simulations and models. In EAGLE and TNG100-1, halo concentration and formation time primarily influence SHMR scatter for haloes with $M_{\rm h}\lt 10^{12}{\, \rm M_\odot }$, but the influence diminishes at high mass. Baryonic processes play a more significant role in LGal. For haloes with $M_{\rm h} \lt 10^{11}{\, \rm M_\odot }$ and $10^{12}{\, \rm M_\odot }\lt M_{\rm h}\lt 10^{13}{\, \rm M_\odot }$, the main drivers of scatter are galaxy star formation rate and age. In the $10^{11.5}{\, \rm M_\odot }\lt M_{\rm h} \lt 10^{12}{\, \rm M_\odot }$ range, halo concentration and formation time are the primary factors. And for haloes with $M_{\rm h} \gt 10^{13}{\, \rm M_\odot }$, supermassive black hole mass becomes more important. Interestingly, it is found that active galactic nucleus feedback may increase the amplitude of the scatter and decrease the dependence on halo properties at high masses.

A self-similar model of galaxy formation and dark halo relaxation

We develop a spherical self-similar model for the formation of a galaxy through gas collapsing in an isolated self-gravitating dark matter halo. As is well known, the self-similarity assumption makes the problem eminently tractable by reducing it to a system of ordinary differential equations. We improve upon the existing literature on self-similar collapse in two ways. First, we include the effects of radiative cooling and the formation of a pseudo-disk at the center of collapse, in a parametrised manner. More importantly, we solve for the evolution of gas and dark matter simultaneously and self-consistently using a novel iterative approach. As a result, our model produces shell trajectories of both gas and dark matter that qualitatively agree with the results of full hydrodynamical simulations of self-gravitating systems. We discuss the impact of various ingredients such as the accretion rate, gas equation of state, disk radius and cooling rate amplitude on the evolution of the gas shells, although we leave the inclusion of stellar and black hole activity to future work. The self-consistent evolution of gas and dark matter allows us to study the response (or `quasi-adiabatic relaxation') of the dark matter trajectories to the presence of collapsing gas, an effect that has gained increasing importance recently in the context of precision estimates of small-scale statistics like the matter power spectrum. Our default configuration produces a relaxation relation in qualitative agreement with that seen in cosmological hydrodynamical simulations, and further allows us to easily study the impact of the model ingredients mentioned above. As an initial application, we vary one ingredient at a time and find that the accretion rate and gas equation of state have the largest impact on the relaxation relation, while the cooling amplitude plays only a minor role. Our model thus provides a convenient framework to rapidly explore the coupled nonlinear impact of multiple astrophysical processes on the mass and velocity profiles of dark matter in galactic halos, and consequently on observables such as rotation curves and gravitational lensing signals.

Dissipative Dark Substructure: The Consequences of Atomic Dark Matter on Milky Way Analog Subhalos

Abstract Using cosmological hydrodynamical zoom-in simulations, we explore the properties of subhalos in Milky Way analogs that contain a subcomponent of atomic dark matter (ADM). ADM differs from cold dark matter (CDM) due to the presence of self-interactions that lead to energy dissipation, analogous to standard model baryons. This model can arise in dark sectors that are natural and theoretically motivated extensions to the standard model. The simulations used in this work were carried out using GIZMO and utilize the FIRE-2 galaxy formation physics in the standard model baryonic sector. For the parameter points we consider, the ADM gas cools efficiently, allowing it to collapse to the center of subhalos. This increases a subhalo’s central density and affects its orbit, with more subhalos surviving small pericentric passages. The subset of subhalos that host satellite galaxies have cuspier density profiles and smaller stellar half-mass radii relative to CDM. The entire population of dwarf galaxies produced in the ADM simulations is more compact than those seen in CDM simulations, unable to reproduce the entire diversity of observed dwarf galaxy structures. Additionally, we also identify a population of highly compact subhalos that consist nearly entirely of ADM and form in the central region of the host, where they can leave distinctive imprints in the baryonic disk. This work presents the first detailed exploration of subhalo properties in a strongly dissipative dark matter scenario, providing intuition for how other regions of ADM parameter space, as well as other dark sector models, would impact galactic-scale observables.

Measurement of the small-scale 3D Lyman-α forest power spectrum

Small-scale correlations measured in the Lyman-α (Lyα) forest encode information about the intergalactic medium and the primordial matter power spectrum. In this article, we present and implement a simple method to measure the 3-dimensional power spectrum, P 3D, of the Lyα forest at wavenumbers k corresponding to small, ∼ Mpc scales. In order to estimate P 3D from sparsely and unevenly distributed data samples, we rely on averaging 1-dimensional Fourier Transforms, as previously carried out to estimate the 1-dimensional power spectrum of the Lyα forest, P 1D. This methodology exhibits a very low computational cost. We confirm the validity of this approach through its application to Nyx cosmological hydrodynamical simulations. Subsequently, we apply our method to the eBOSS DR16 Lyα forest sample, providing as a proof of principle, a first P 3D measurement averaged over two redshift bins z = 2.2 and z = 2.4. This work highlights the potential for forthcoming P 3D measurements, from upcoming large spectroscopic surveys, to untangle degeneracies in the cosmological interpretation of P 1D.

Euclid preparation

The Euclid photometric survey of galaxy clusters stands as a powerful cosmological tool, with the capacity to significantly propel our understanding of the Universe. Despite being subdominant to dark matter and dark energy, the baryonic component of our Universe holds substantial influence over the structure and mass of galaxy clusters. This paper presents a novel model that can be used to precisely quantify the impact of baryons on the virial halo masses of galaxy clusters using the baryon fraction within a cluster as a proxy for their effect. Constructed on the premise of quasi-adiabaticity, the model includes two parameters, which are calibrated using non-radiative cosmological hydrodynamical simulations, and a single large-scale simulation from the Magneticum set, which includes the physical processes driving galaxy formation. As a main result of our analysis, we demonstrate that this model delivers a remarkable 1% relative accuracy in determining the virial dark matter-only equivalent mass of galaxy clusters starting from the corresponding total cluster mass and baryon fraction measured in hydrodynamical simulations. Furthermore, we demonstrate that this result is robust against changes in cosmological parameters and against variation of the numerical implementation of the subresolution physical processes included in the simulations. Our work substantiates previous claims regarding the impact of baryons on cluster cosmology studies. In particular, we show how neglecting these effects would lead to biased cosmological constraints for a Euclid-like cluster abundance analysis. Importantly, we demonstrate that uncertainties associated with our model arising from baryonic corrections to cluster masses are subdominant when compared to the precision with which mass–observable (i.e. richness) relations will be calibrated using Euclid and to our current understanding of the baryon fraction within galaxy clusters.

Active galactic nucleus X-ray luminosity function and absorption function in the Early Universe (3 ≤ z ≤ 6)

The X-ray luminosity function (XLF) of active galactic nuclei (AGN) offers a robust tool to study the evolution and the growth of the supermassive black-hole population over cosmic time. Owing to the limited area probed by X-ray surveys, optical surveys are routinely used to probe the accretion in the high-redshift Universe z ≥ 3. However, optical surveys may be incomplete because they are strongly affected by dust redenning. In this work we derive the XLF and its evolution at high redshifts (z ≥ 3) using a large sample of AGN selected in different fields with various areas and depths covering a wide range of luminosities. Additionally, we put the tightest yet constraints on the absorption function in this redshift regime. In particular, we used more than 600 soft X-ray selected (0.5 − 2 keV) high-z sources in the Chandra deep fields, the Chandra COSMOS Legacy survey, and the XMM-XXL northern field. We derived the X-ray spectral properties for all sources via spectral fitting, using a consistent technique and model. To model the parametric form of the XLF and the absorption function, we used a Bayesian methodology, allowing us to correctly propagate the uncertainties for the observed X-ray properties of our sources and also the absorption effects. The evolution of XLF is in agreement with a pure density evolution model similar to what is witnessed at optical wavelengths, although a luminosity-dependent density evolution model cannot be securely ruled out. A large fraction (∼60%) of our sources are absorbed by column densities of NH ≥ 1023 cm−2, while ∼17% of the sources are Compton-Thick. Our results favour a scenario where both the interstellar medium of the host and the AGN torus contribute to the obscuration. The derived black hole accretion rate density is roughly in agreement with the large-scale cosmological hydrodynamical simulations, if one takes into account the results that the X-ray AGN are hosted by massive galaxies, while it differs from that derived using JWST data. The latter could be due to the differences in the AGN and host-galaxy properties.

The origin of cold gas in the circumgalactic medium

Context. The presence of cold gas (T ≲ 104 K) in the circumgalactic medium (CGM) of galaxies has been confirmed in observations and in high-resolution simulations, but its origin is still a puzzle. Possible mechanisms are cold accretion from the intergalactic medium (IGM), clumps embedded in outflows and transported from the disk, and gas detaching from the hot CGM phase via thermal instabilities. Aims. In this work we characterize the history of cold CGM gas in order to identify the dominant origin channels at different evolutionary stages of the main galaxy. Methods. To this end, we tracked gas particles in different snapshots of the smoothed particle hydrodynamics (SPH) cosmological zoom-in simulation Eris2k. We performed a backward tracking of cold gas, starting from different redshifts until we could identify one of the followings origins for the particle: cold inflow, ejection from the disk, cooling down in situ, or stripping from a satellite. We also performed a forward tracking of gas in different components of the galaxy (such as the disk and outflows). Results. We find a clear transition between two epochs. For z &gt; 2, most cold gas (up to 80%) in the CGM comes from cold accretion streams as the galaxy is accreting in the cold mode from the IGM. At lower z, gas either cools down in situ after several recycles (with 10–20% of the gas cooling in outflow), or it is ejected directly from the disk (up to 30%). Outflows have a major contribution to the cold CGM gas budget at z &lt; 1, with almost 50% of the hot gas cooling in outflow. Finally, we discuss possible mechanisms for CGM cooling, showing that the thermally unstable gas with tcool/tff &lt; 1 (precipitation-regulated feedback) is abundant up to r ∼ 100 kpc and cooling times are shorter than 50 Myr for densities n &gt; 10−2 cm−3.

Predicting Stellar Masses of the First Galaxies Using Graph Neural Networks

Theoretical models of galaxy formation and evolution are primarily investigated through cosmological simulations and semi-analytical models. The former method consumes O(106) core-hours explicitly modeling the dynamics of the galaxies, whereas the latter method only requires O(103) core-hours foregoing directly simulating internal structure for computational efficiency. In this work, we present a proof-of-concept machine learning regression model, using a graph neural network architecture, to predict the stellar mass of high-redshift galaxies solely from their dark matter merger trees, trained from a radiation hydrodynamics cosmological simulation of the first galaxies.

Characterizing tidal features around galaxies in cosmological simulations

ABSTRACT Tidal features provide signatures of recent mergers and offer a unique insight into the assembly history of galaxies. The Vera C. Rubin Observatory’s Legacy Survey of Space and Time (LSST) will enable an unprecedentedly large survey of tidal features around millions of galaxies. To decipher the contributions of mergers to galaxy evolution it will be necessary to compare the observed tidal features with theoretical predictions. Therefore, we use cosmological hydrodynamical simulations NewHorizon, eagle, IllustrisTNG, and Magneticum to produce LSST-like mock images of z ∼ 0 galaxies (z ∼ 0.2 for NewHorizon) with $M_{\scriptstyle \star ,\text{ 30 pkpc}}\ge 10^{9.5}$ M$_{\scriptstyle \odot }$. We perform a visual classification to identify tidal features and classify their morphology. We find broadly good agreement between the simulations regarding their overall tidal feature fractions: $f_{{\small NewHorizon}}=0.40\pm 0.06$, $f_{{\small EAGLE}}=0.37\pm 0.01$, $f_{{\small TNG}}=0.32\pm 0.01$, and $f_{{\small Magneticum}}=0.32\pm 0.01$, and their specific tidal feature fractions. Furthermore, we find excellent agreement regarding the trends of tidal feature fraction with stellar and halo mass. All simulations agree in predicting that the majority of central galaxies of groups and clusters exhibit at least one tidal feature, while the satellite members rarely show such features. This agreement suggests that gravity is the primary driver of the occurrence of visually identifiable tidal features in cosmological simulations, rather than subgrid physics or hydrodynamics. All predictions can be verified directly with LSST observations.

Dense stellar clump formation driven by strong quasar winds in the FIRE cosmological hydrodynamic simulations

ABSTRACT We investigate the formation of dense stellar clumps in a suite of high-resolution cosmological zoom-in simulations of a massive, star-forming galaxy at z ∼ 2 under the presence of strong quasar winds. Our simulations include multiphase ISM physics from the Feedback In Realistic Environments (FIRE) project and a novel implementation of hyper-refined accretion disc winds. We show that powerful quasar winds can have a global negative impact on galaxy growth while in the strongest cases triggering the formation of an off-centre clump with stellar mass ${\rm M}_{\star }\sim 10^{7}\, {\rm M}_{\odot }$, effective radius ${\rm R}_{\rm 1/2\, \rm Clump}\sim 20\, {\rm pc}$, and surface density $\Sigma _{\star } \sim 10^{4}\, {\rm M}_{\odot }\, {\rm pc}^{-2}$. The clump progenitor gas cloud is originally not star-forming, but strong ram pressure gradients driven by the quasar winds (orders of magnitude stronger than experienced in the absence of winds) lead to rapid compression and subsequent conversion of gas into stars at densities much higher than the average density of star-forming gas. The AGN-triggered star-forming clump reaches ${\rm SFR} \sim 50\, {\rm M}_{\odot }\, {\rm yr}^{-1}$ and $\Sigma _{\rm SFR} \sim 10^{4}\, {\rm M}_{\odot }\, {\rm yr}^{-1}\, {\rm kpc}^{-2}$, converting most of the progenitor gas cloud into stars in ∼2 Myr, significantly faster than its initial free-fall time and with stellar feedback unable to stop star formation. In contrast, the same gas cloud in the absence of quasar winds forms stars over a much longer period of time (∼35 Myr), at lower densities, and losing spatial coherency. The presence of young, ultra-dense, gravitationally bound stellar clumps in recently quenched galaxies could thus indicate local positive feedback acting alongside the strong negative impact of powerful quasar winds, providing a plausible formation scenario for globular clusters.

The history and mass content of cluster galaxies in the EAGLE simulation

We explore the mass content of galaxies residing in galaxy clusters at $z=0$ in the cosmological hydrodynamical simulation. We also explore the galaxies' mass build-up through cosmic time. We used a galaxy catalogue generated with the algorithm, which identifies subhaloes consistently over time by tracking their dynamical evolution throughout the simulation. The satellite subhalo-to-stellar mass relation (SHSMR) is well described by a double power law, becoming increasingly steeper with stellar mass. At stellar masses $9&lt; satellites have 20-25 the subhalo mass of central galaxies at fixed stellar mass. At high stellar masses, $ the satellite SHSMR is consistent with that of centrals. The satellite SHSMR decreases steeply for satellites closer to the cluster centre, even in projection, broadly consistent with recent weak lensing measurements. The scatter in the satellite SHSMR is larger than that of central galaxies at all cluster masses and cluster-centric distances $R&lt;R_ 200m $. The SHSMR scatter decreases with stellar mass by about 12 over an order of magnitude, but this dependence can be explained by the mixing of infall times when binning by stellar mass. By splitting satellites into direct and indirect infallers (those that fell into their current host as a central galaxy or as the satellite of an infalling group, respectively) we clearly show the impact of pre-processing separately on satellite galaxies' dark and stellar mass. There is significant dark matter pre-processing; the most recent infallers into massive clusters ($ had already lost up to 50 of their dark matter by the time of infall, particularly if they fell in indirectly as satellites of another host. On the contrary, on average, satellite galaxies are still gaining stellar mass at the time of infall and they do so for another 2 Gyr afterwards, although we see evidence of a slowing growth for indirect infallers. How much and for how long they continue to gain stellar mass depends primarily on the gas mass fraction available at infall. Overall, pre- and post-processing have similar impacts on the satellite SHSMR. Finally, we provide a simple prescription to infer the mean mass loss experienced by satellites as a function of cluster-centric distance based on a comparison to central galaxies, convenient for observational weak lensing measurements.