IllustrisTNG Simulations Research Articles

Warped Disk Galaxies. I. Linking U-type Warps in Groups/Clusters to Jellyfish Galaxies

Warped disk galaxies are classified into two morphologies: S and U types. Conventional theories routinely attribute both types to galactic tidal interaction and/or gas accretion, but reproducing U types in simulations is extremely challenging. Here we investigate whether both types are governed by the same mechanisms using the most extensive sample of ∼8000 nearby (0.02 < z < 0.06) massive (M */M ⊙ > 109) edge-on disks from the Sloan Digital Sky Survey. We find that U types show on average bluer optical colors and a higher specific star formation rate (sSFR) than S types, with more strongly warped U types having a higher sSFR. We also find that while the S-type warp properties correlate with the tidal force by the nearest neighbor regardless of the environment, there is no such correlation for U types in groups/clusters, suggesting a nontidal mechanism could be at play for U types, such as ram pressure stripping (RPS). Indeed, U types are more common in groups/clusters than in fields and they have stellar mass, gas fraction, sSFR enhancement, and phase-space distribution closely analogous to RPS-induced jellyfish galaxies in clusters. We furthermore show that the stellar disks of most RPS galaxies in the IllustrisTNG simulation are warped in a U shape and bent in the opposite direction of stripped gas tails, satisfying theoretical expectations for stellar warps embedded in jellyfishes. We therefore suggest that despite the majority of U types that live in fields being still less explained, RPS can be an alternative origin for those in groups/clusters.

Continuous Wavelet Analysis of Matter Clustering Using the Gaussian-derived Wavelet

Continuous wavelet analysis has been increasingly employed in various fields of science and engineering due to its remarkable ability to maintain optimal resolution in both space and scale. Here, we introduce wavelet-based statistics, including the wavelet power spectrum, wavelet cross correlation, and wavelet bicoherence, to analyze the large-scale clustering of matter. For this purpose, we perform wavelet transforms on the density distribution obtained from the one-dimensional Zel’dovich approximation and then measure the wavelet power spectra and wavelet bicoherences of this density distribution. Our results suggest that the wavelet power spectrum and wavelet bicoherence can identify the effects of local environments on the clustering at different scales. Moreover, we apply the statistics based on the three-dimensional isotropic wavelet to the IllustrisTNG simulation at z = 0, and investigate the environmental dependence of the matter clustering. We find that the clustering strength of the total matter increases with increasing local density except on the largest scales. Besides, we notice that the gas traces dark matter better than stars on large scales in all environments. On small scales, the cross correlation between the dark matter and gas first decreases and then increases with increasing density. This is related to the impacts of the active galactic nucleus feedback on the matter distribution, which also varies with the density environment in a similar trend to the cross correlation between dark matter and gas. Our findings are qualitatively consistent with previous studies on matter clustering.

Exploring the effect of baryons on the radial distribution of satellite galaxies with GAMA and IllustrisTNG

ABSTRACT We explore the radial distribution of satellite galaxies in groups in the Galaxy and Mass Assembly (GAMA) survey and the IllustrisTNG simulations. Considering groups with masses $12.0 \le \log _{10} (\mathcal {M}_h / h^{-1} \, \mathrm{M}_{\odot }) \lt 14.8$ at z &amp;lt; 0.267, we find a good agreement between GAMA and a sample of TNG300 groups and galaxies designed to match the GAMA selection. Both display a flat profile in the centre of groups, followed by a decline that becomes steeper towards the group edge, and normalized profiles show no dependence on group mass. Using matched satellites from TNG and dark matter-only TNG-Dark runs we investigate the effect of baryons on satellite radial location. At z = 0, we find that the matched subhaloes from the TNG-Dark runs display a much flatter radial profile: namely, satellites selected above a minimum stellar mass exhibit both smaller halocentric distances and longer survival times in the full-physics simulations compared to their dark-matter only analogues. We then divide the TNG satellites into those which possess TNG-Dark counterparts and those which do not, and develop models for the radial positions of each. We find the satellites with TNG-Dark counterparts are displaced towards the halo centre in the full-physics simulations, and this difference has a power-law behaviour with radius. For the ‘orphan’ galaxies without TNG-Dark counterparts, we consider the shape of their radial distribution and provide a model for their motion over time, which can be used to improve the treatment of satellite galaxies in semi-analytic and semi-empirical models of galaxy formation.

IMaNGA: mock MaNGA galaxies based on IllustrisTNG and MaStar SSPs – I. Construction and analysis of the mock data cubes

ABSTRACT Galaxy formation and evolution simulations are essential tools to probe poorly known astrophysics processes, but particular care is needed to compare simulations with galaxy observations, as observed data need to be modelled as well. We present a method to generate mock galaxies from the hydro-dynamical IllustrisTNG simulations which are suited to compare with integral field spectroscopic observation of galaxies from the SDSS-IV/Mapping Nearby Galaxies at Apache Point Observatory (MaNGA) survey. First, we include the same instrumental effects and procedures as adopted in the acquisition and analysis of real data. Furthermore, we generate the galaxy spectra from the simulations using new stellar population models based on the MaNGA stellar library (MaStar). In this way, our mock data cubes have the same spatial sampling, cover the same wavelength range (3600–10 300 Å), and share the same spectral resolution (R ≈ 1800) and flux calibration of real MaNGA galaxy spectra. In this first paper, we demonstrate the method over an early- and a late-type simulated galaxy from TNG50. We analyse the correspondent mock MaNGA-like data cubes with the same full spectral fitting code, firefly, which was used for the observed spectra. We find that the intrinsic and recovered age and metallicity gradients are consistent within 1σ, with residuals over all tassels consistent with 0 at the 68 ${{\ \rm per\ cent}}$ confidence level. We also perform the challenging test at comparing intrinsic and recovered star formation histories, finding a close resemblance between input and output. In follow-up papers, we will present a full simulated MaNGA-like catalogue (≈10 000 galaxies) with a comprehensive comparison of TNG50 simulations to MaNGA observational results.

Star formation characteristics of CNN-identified post-mergers in the Ultraviolet Near Infrared Optical Northern Survey (UNIONS)

ABSTRACT The importance of the post-merger epoch in galaxy evolution has been well documented, but post-mergers are notoriously difficult to identify. While the features induced by mergers can sometimes be distinctive, they are frequently missed by visual inspection. In addition, visual classification efforts are highly inefficient because of the inherent rarity of post-mergers (~1 per cent in the low-redshift Universe), and non-parametric statistical merger selection methods do not account for the diversity of post-mergers or the environments in which they appear. To address these issues, we deploy a convolutional neural network (CNN) that has been trained and evaluated on realistic mock observations of simulated galaxies from the IllustrisTNG simulations, to galaxy images from the Canada France Imaging Survey, which is part of the Ultraviolet Near Infrared Optical Northern Survey. We present the characteristics of the galaxies with the highest CNN-predicted post-merger certainties, as well as a visually confirmed subset of 699 post-mergers. We find that post-mergers with high CNN merger probabilities [p(x) &amp;gt; 0.8] have an average star formation rate that is 0.1 dex higher than a mass- and redshift-matched control sample. The SFR enhancement is even greater in the visually confirmed post-merger sample, a factor of 2 higher than the control sample.

Lopsided galactic disks in IllustrisTNG

A significant fraction of nearby late-type galaxies are lopsided. We study the asymmetry of the stellar component in a sample of well-resolved disky galaxies selected from the last snapshot of the Illustris TNG100 simulation based on their flatness and rotational support. Among 1912 disks, we identify 161 objects with significant asymmetry in terms of them = 1 Fourier mode of the stellar component within (1−2) stellar half-mass radii and describe their properties using three representative examples. The profiles of them = 1 mode typically increase with radius, and the corresponding phase is constant in the asymmetric region, signifying a global distortion. Following the evolution of the lopsided disks over time, we find that their history is rather uneventful and the occurrence of the asymmetry is fairly recent. Only about 1/3 of the lopsided disks experienced any strong interaction recently that could have led to the distortion of their shape: 24% were affected by a more massive object and 9% underwent a gas-rich merger. Still, a majority of lopsided disks show a significant increase in their recent star formation rate. The most frequent mechanism for the formation of lopsided disks thus seems to be asymmetric star formation probably related to gas accretion, although the distortions in the gas and stars are not strongly correlated. This picture is supported by the finding that the lopsided population on average contains more gas, has higher star formation rate, lower metallicity and bluer color than the remaining disks. These correlations are similar to those seen in real galaxies, even though the fraction of simulated lopsided disks (8%) is much lower than in observations (30%). The observed correlation between the presence of the asymmetry and a bar is not reproduced either. These discrepancies may be due to overquenching or insufficient resolution of IllustrisTNG simulations.

Disc instability and bar formation: view from the IllustrisTNG simulations

ABSTRACT We make use of $z\, {=}\, 0$ samples of strongly barred and unbarred disc galaxies from the TNG100 and TNG50 cosmological hydrodynamical simulations to assess the performance of the simple disc instability criterion proposed by Efstathiou, Lake &amp; Negroponte (1982) (ELN-criterion). We find that strongly barred galaxies generally assemble earlier, are more star dominated in their central regions, and have more massive and more compact discs than unbarred galaxies. The ELN-criterion successfully identifies ${\sim }\, 75$ and ${\sim }\, 80{{ per\ cent}}$ of the strongly barred and the unbarred galaxies, respectively. Strongly barred galaxies that the criterion fails to identify tend to have more extended discs, higher spin values and bars that assembled later than is typical for the bulk of the barred population. The bars in many of these cases appear to be produced by an interaction with a close neighbour (i.e. to be externally triggered) rather than to result from secular growth in the disc. On the other hand, we find that unbarred galaxies misclassified as barred by the ELN-criterion typically have stellar discs similar to those of barred galaxies, although more extended in the vertical direction and less star-dominated in their central regions, possibly reflecting later formation times. In addition, the bulge component of these galaxies is significantly more prominent at early times than in the strongly barred sample. Thus, the ELN-criterion robustly identifies secular bar instabilities in most simulated disc galaxies, but additional environmental criteria are needed to account for interaction-induced bar formation.

3D intrinsic shapes of quiescent galaxies in observations and simulations

ABSTRACT We study the intrinsic 3D shapes of quiescent galaxies over the last half of cosmic history based on their axial ratio distribution. To this end, we construct a sample of unprecedented size, exploiting multiwavelength u-to-Ks photometry from the deep wide-area surveys KiDS+VIKING paired with high-quality i-band imaging from HSC-SSP. The dependences of the shapes on mass, redshift, photometric bulge prominence and environment are considered. For comparison, the intrinsic shapes of quenched galaxies in the IllustrisTNG simulations are analysed and contrasted with their formation history. We find that over the full 0 &amp;lt; z &amp;lt; 0.9 range, and in both simulations and observations, spheroidal 3D shapes become more abundant at $M_* \gt 10^{11}\, \mathrm{M}_{\odot }$, with the effect being most pronounced at lower redshifts. In TNG, the most massive galaxies feature the highest ex situ stellar mass fractions, pointing to violent relaxation via mergers as the mechanism responsible for their 3D shape transformation. Larger differences between observed and simulated shapes are found at low to intermediate masses. At any mass, the most spheroidal quiescent galaxies in TNG feature the highest bulge mass fractions, and, conversely, observed quiescent galaxies with the highest bulge-to-total ratios are found to be intrinsically the roundest. Finally, we detect an environmental influence on galaxy shape, at least at the highest masses, such that at fixed mass and redshift, quiescent galaxies tend to be rounder in denser environments.

Multi-epoch machine learning 1: Unravelling nature versus nurture for galaxy formation

Abstract We present a novel machine learning method for predicting the baryonic properties of dark matter only subhalos from N-body simulations. Our model is built using the extremely randomized tree (ERT) algorithm and takes subhalo properties over a wide range of redshifts as its input features. We train our model using the IllustrisTNG simulations to predict blackhole mass, gas mass, magnitudes, star formation rate, stellar mass, and metallicity. We compare the results of our method with a baseline model from previous works, and against a model that only considers the mass history of the subhalo. We find that our new model significantly outperforms both of the other models. We then investigate the predictive power of each input by looking at feature importance scores from the ERT algorithm. We produce feature importance plots for each baryonic property, and find that they differ significantly. We identify low redshifts as being most important for predicting star formation rate and gas mass, with high redshifts being most important for predicting stellar mass and metallicity, and consider what this implies for nature versus nurture. We find that the physical properties of galaxies investigated in this study are all driven by nurture and not nature. The only property showing a somewhat stronger impact of nature is the present-day star formation rate of galaxies. Finally we verify that the feature importance plots are discovering physical patterns, and that the trends shown are not an artefact of the ERT algorithm.

Realistic galaxy images and improved robustness in machine learning tasks from generative modelling

ABSTRACTWe examine the capability of generative models to produce realistic galaxy images. We show that mixing generated data with the original data improves the robustness in downstream machine learning tasks. We focus on three different data sets: analytical Sérsic profiles, real galaxies from the COSMOS survey, and galaxy images produced with the SKIRT code, from the IllustrisTNG simulation. We quantify the performance of each generative model, using the Wasserstein distance between the distributions of morphological properties (e.g. the Gini-coefficient, the asymmetry, and ellipticity), the surface brightness distribution on various scales (as encoded by the power spectrum), the bulge statistic, and the colour for the generated and source data sets. With an average Wasserstein distance (Fréchet Inception Distance) of 7.19 × 10−2 (0.55), 5.98 × 10−2 (1.45), and 5.08 × 10−2 (7.76) for the Sérsic, COSMOS and SKIRT data set, respectively, our best models convincingly reproduce even the most complicated galaxy properties and create images that are visually indistinguishable from the source data. We demonstrate that by supplementing the training data set with generated data, it is possible to significantly improve the robustness against domain-shifts and out-of-distribution data. In particular, we train a convolutional neural network to denoise a data set of mock observations. By mixing generated images into the original training data, we obtain an improvement of 11 and 45 per cent in the model performance regarding domain-shifts in the physical pixel size and background noise level, respectively.

The Redshift Evolution of the Binary Black Hole Merger Rate: A Weighty Matter

Gravitational-wave detectors are starting to reveal the redshift evolution of the binary black hole (BBH) merger rate, R BBH(z). We make predictions for R BBH(z) as a function of black hole mass for systems originating from isolated binaries. To this end, we investigate correlations between the delay time and black hole mass by means of the suite of binary population synthesis simulations, COMPAS. We distinguish two channels: the common envelope (CE), and the stable Roche-lobe overflow (RLOF) channel, characterized by whether the system has experienced a common envelope or not. We find that the CE channel preferentially produces BHs with masses below about 30 M ⊙ and short delay times (t delay ≲ 1 Gyr), while the stable RLOF channel primarily forms systems with BH masses above 30 M ⊙ and long delay times (t delay ≳ 1 Gyr). We provide a new fit for the metallicity-dependent specific star formation rate density based on the Illustris TNG simulations, and use this to convert the delay time distributions into a prediction of R BBH(z). This leads to a distinct redshift evolution of R BBH(z) for high and low primary BH masses. We furthermore find that, at high redshift, R BBH(z) is dominated by the CE channel, while at low redshift, it contains a large contribution (∼40%) from the stable RLOF channel. Our results predict that, for increasing redshifts, BBHs with component masses above 30 M ⊙ will become increasingly scarce relative to less massive BBH systems. Evidence of this distinct evolution of R BBH(z) for different BH masses can be tested with future detectors.

Relative distribution of dark matter, gas, and stars around cosmic filaments in the IllustrisTNG simulation

We present a comprehensive study of the distribution of matter around different populations of large-scale cosmic filaments, using the IllustrisTNG simulation atz = 0. We computed the dark matter (DM), gas, and stellar radial density profiles of filaments, and we characterise the distribution of the baryon fraction in these structures. We find that baryons exactly follow the underlying DM distribution only down tor ∼ 7 Mpc to the filament spines. At shorter distances (r &lt; 7 Mpc), the baryon fraction profile of filaments departs from the cosmic value Ωb/Ωm. While in ther ∼ 0.7−7 Mpc radial domain this departure is due to the radial accretion of the warm-hot intergalactic medium (WHIM) towards the filament cores (creating an excess of baryons with respect to the cosmic fraction), the cores of filaments (r &lt; 0.7 Mpc) show a clear baryon depletion instead. The analysis of the efficiency of active galactic nuclei (AGN) feedback events in filaments reveals that they are potentially powerful enough to eject gas outside of the gravitational potential wells of filaments. We show that the large-scale environment (i.e. denser versus less dense, hotter versus colder regions) has a non-negligible effect on the absolute values of the DM, gas, and stellar densities around filaments. Nevertheless, the relative distribution of baryons with respect to the underlying DM density field is found to be independent of the filament population. Finally, we provide scaling relations between the gas density, temperature, and pressure for the different populations of cosmic filaments. We compare these relations to those pertaining to clusters of galaxies, and find that these cosmic structures occupy separate regions of the density-temperature and density-pressure planes.

Coevolution of Brightest Cluster Galaxies and Their Host Clusters in IllustrisTNG

We use IllustrisTNG simulations to explore the dynamic scaling relation between massive clusters and their—central—brightest cluster galaxies (BCGs). The IllustrisTNG-300 simulation we use includes 280 massive clusters from the z = 0 snapshot with M 200 > 1014 M ⊙, enabling a robust statistical analysis. We derive the line-of-sight velocity dispersion of the stellar particles of the BCGs (σ *,BCG), analogous to the observed BCG stellar velocity dispersion. We also compute the subhalo velocity dispersion to measure the cluster velocity dispersion (σ cl). Both σ *,BCG and σ cl are proportional to the cluster halo mass, but the slopes differ slightly. Thus, like the observed relation, σ *,BCG/σ cl declines as a function of σ cl, but the scatter is large. We explore the redshift evolution of the σ *,BCG − σ cl scaling relation for z ≲ 1 in a way that can be compared directly with observations. The scaling relation has a similar slope at high redshift, but the scatter increases because of the large scatter in σ *,BCG. The simulations imply that high-redshift BCGs are dynamically more complex than their low-redshift counterparts.

Mock HUBS observations of hot gas with IllustrisTNG

The lack of adequate X-ray observing capability is seriously impeding the progress in understanding the hot phase of circumgalactic medium (CGM), which is predicted to extend to the virial radius of a galaxy or beyond, and thus in acquiring key boundary conditions for studying galaxy evolution. To this end, the Hot Universe Baryon Surveyor (HUBS) is proposed. HUBS is designed to probe hot CGM by detecting its emission or absorption lines with a non-dispersive X-ray spectrometer of high resolution and high throughput. The spectrometer consists of a 60 × 60 array of microcalorimeters, with each detector providing an energy resolution of 2 eV, and is placed in the focal plane of an X-ray telescope of 1∘ field-of-view. With such a design, the spectrometer is highly optimized for detecting X-ray-emitting hot gas in the CGM of local galaxies, as well as in intra-group medium (IGrM), intra-cluster medium (ICM), or intergalactic medium (IGM). To assess the scientific potential of HUBS, in this work, we created mock observations of galaxies, groups, and clusters at different redshifts with the IllustrisTNG simulation. Focusing exclusively on emission studies, we took into account the effects of light cone, Galactic foreground emission, and background AGN contribution in the mock observations. From the observations, we made mock X-ray images and spectra, analyzed them to derive the properties of the emitting gas in each case, and compared the results with the input parameters from the simulation. The results show that HUBS is well suited for studying hot CGM at low redshifts. The redshift range is significantly extended for measuring IGrM and ICM. The sensitivity limits are also presented for detecting extended emission of low surface brightness.

A systematic search for galaxy protocluster cores at the transition epoch of their star formation activity

ABSTRACT The redshift of z ∼ 1.5 is the transition epoch of protoclusters (PCs) from the star-forming phase into the quenching phase, and hence an appropriate era to investigate the build up of the quenched population. We define a ‘core’ as the most massive halo in a given PC, where environmental effects are likely to work most effectively, and search for cores at 1 &amp;lt; z &amp;lt; 1.5. We use a photometric redshift catalogue of a wide (effective area of ${\sim}22.2\, \mathrm{deg}^{2}$) and deep ($i\sim 26.8\, \mathrm{mag}$) optical survey with Subaru Hyper-Suprime Cam. Regarding galaxies with log (M*/M⊙) &amp;gt; 11.3 as the central galaxies of PC cores, we estimate their average halo mass by clustering analysis and find it to be log (Mh/M⊙) ∼ 13.7. An expected mass growth by the IllustrisTNG simulation and the observed overdensities around them suggest that the PC cores we find are progenitors of present-day clusters. Classifying our galaxy sample into red and blue galaxies, we calculate the stellar mass function (SMF) and the red galaxy fraction. The SMFs in the PC cores are more-top heavy than field, implying early high-mass galaxy formation and disruption of low-mass galaxies. We also find that the red fraction increases with stellar mass, consistent with stellar mass dependent environmental quenching recently found at z &amp;gt; 1. Interestingly, although the cores with red and blue centrals have similar halo masses, only those with red centrals show a significant red fraction excess compared to the field, suggesting a conformity effect. Some observational features of PC cores may imply that the conformity is caused by assembly bias.

Dust spectral energy distributions in Milky Way-like galaxies in the IllustrisTNG simulations based on the evolution of grain size distribution

ABSTRACT To understand how the evolution of grain size distribution in galaxies affects observed dust properties, we apply a post-processing dust evolution model to galaxy merger trees from the IllustrisTNG cosmological hydrodynamical simulation. Our dust model includes stellar dust production, sputtering in hot gas, dust growth by accretion and coagulation in the dense interstellar medium (ISM), and shattering in the diffuse ISM. We decompose the grain size distribution into different dust species depending on the elemental abundances and the dense ISM fraction given by the simulation. In our previous work, we focused on Milky Way (MW) analogues and reproduced the observed MW extinction curve. In this study, we compute dust spectral energy distributions (SEDs) for the MW analogues. Our simulated SEDs broadly reproduce the observed MW SED within their dispersion and so does the observational data of nearby galaxies, although they tend to underpredict the MW SED at short wavelengths where emission is dominated by polycyclic aromatic hydrocarbons. We find that metallicity and dense gas fraction are the most critical factors for the SED shape, through their influence on coagulation and shattering. The overall success of our models in reproducing the MW SED further justifies the dust evolution processes included in the model and predicts the dispersion in the SEDs caused by the variety in the assembly history. We also show that the most significant increase in the dust SED occurs between redshifts z ∼ 3 and 2 in the progenitors of the simulated MW-like galaxies.

Quenching of Massive Disk Galaxies in the IllustrisTNG Simulation

Abstract A rare population of massive disk galaxies have been found to invade the red sequence dominated by early-type galaxies. The formation and origins of these red/quenched massive disk galaxies have recently gained great interest. The quenching mechanisms that are usually proposed, such as bar quenching and environment quenching, do not seem to be suitable for those bulgeless quenched disks in a low-density environment. In this paper, we use the TNG300 simulation to investigate the formation of massive quenched central disk galaxies. It is found that these galaxies contain less gas than their star-forming counterparts and harbor giant super massive black holes (SMBHs; above 108 M ⊙). By tracing their formation history, we found that quenched disk galaxies formed early and preserved disk morphology for cosmological timescales. They have experienced less than one major merger on average, and mini-mergers (mass ratio &lt;1/10) have mainly contributed to the growth of their SMBHs. In the IllustrisTNG simulation, the black hole feedback mode switches from thermal to kinetic feedback when the black hole mass is more massive than ∼108 M ⊙, which is more efficient to eject gas outside of the galaxy and to suppress further cooling of the hot gaseous halo. We conclude that the dominant quenching mechanism in massive red/quenched disk galaxies is kinetic active galactic nuclei feedback.

Scatter in the satellite galaxy SHMR: fitting functions, scaling relations, and physical processes from the IllustrisTNG simulation

ABSTRACT The connection between galaxies and their dark matter haloes is often described with the stellar-to-halo mass relation (SHMR). Satellite galaxies in clusters follow an SHMR distinct from central galaxies because of the environmental processes that they are subject to, and the variety of accretion histories leads to an important scatter in this relation. In this work, we use the suite of magnetohydrodynamical simulations IllustrisTNG to study the scatter in the satellite galaxy SHMR, and extract the parameters that can best allow to understand it. Active galaxies, that represent a very small fraction of cluster galaxies, follow a very different relation than their passive counterparts, mainly because they were accreted much more recently. For this latter population, we find that the distance to the cluster centre is a good predictor of variations in the SHMR, but some information on the galaxy orbital history, such as the distance of closest approach to the host centre, is an even better one, although it is in practice more difficult to measure. In addition, we found that galaxy compactness is also correlated with the SHMR, while the host cluster properties (mass and concentration, formation redshift, mass and size of BCG) do not play a significant role. We provide accurate fitting functions and scaling relations to the scientific community, useful to predict the subhalo mass given a set of observable parameters. Finally, we connect the scatter in the SHMR to the physical processes affecting galaxies in clusters, and how they impact the different satellite subpopulations.

Intrinsic alignments in IllustrisTNG and their implications for weak lensing: Tidal shearing and tidal torquing mechanisms put to the test

ABSTRACT Accurate measurements of the cosmic shear signal require a separation of the true weak gravitational lensing signal from intrinsic shape correlations of galaxies. These ‘intrinsic alignments’ of galaxies originate from galaxy formation processes and are expected to be correlated with the gravitational field through tidal processes affecting the galaxies, such as tidal shearing for elliptical galaxies and tidal torquing for spiral galaxies. In this study, we use morphologically selected samples of elliptical and spiral galaxies from the illustrisTNG simulation at z = 0 and z = 1 to test the commonly employed linear (tidal shearing) and quadratic (tidal torquing) models for intrinsic alignments. We obtain local measurements of the linear and quadratic alignment parameters, including corrections for large-scale anisotropies of the cosmologically small simulation volume, and study their dependence on galaxy and environmental properties. We find a significant alignment signal for elliptical galaxies (linear model), that increases with mass and redshift. Spiral galaxies (quadratic model), on the other hand, exhibit a significant signal only for the most massive objects at z = 1. We show the quadratic model for spiral galaxies to break down at its fundamental assumptions, and simultaneously obtain a significant signal of spiral galaxies to align according to the linear model. We use the derived alignment parameters to compute intrinsic alignment spectra and estimate the expected contamination in the weak lensing signal obtained by Euclid.

Emulating Sunyaev–Zeldovich images of galaxy clusters using autoencoders

ABSTRACT We develop a machine-learning (ML) algorithm that generates high-resolution thermal Sunyaev–Zeldovich (SZ) maps of novel galaxy clusters given only halo mass and mass accretion rate (MAR). The algorithm uses a conditional variational autoencoder (CVAE) in the form of a convolutional neural network and is trained with SZ maps generated from the IllustrisTNG simulation. Our method can reproduce many of the details of galaxy clusters that analytical models usually lack, such as internal structure and aspherical distribution of gas created by mergers, while achieving the same computational feasibility, allowing us to generate mock SZ maps for over 105 clusters in 30 s on a laptop. We show that the model is capable of generating novel clusters (i.e. not found in the training set) and that the model accurately reproduces the effects of mass and MAR on the SZ images, such as scatter, asymmetry, and concentration, in addition to modelling merging sub-clusters. This work demonstrates the viability of ML-based methods for producing the number of realistic, high-resolution maps of galaxy clusters necessary to achieve statistical constraints from future SZ surveys.