-
Building models of the Universe with hydrodynamic simulations
Authors:
Edoardo Altamura
Abstract:
Hydrodynamic simulations have become irreplaceable in modern cosmology for exploring complex systems and making predictions to steer future observations. In Chapter 1, we begin with a philosophical discussion on the role of simulations in science. We argue that simulations can bridge the gap between empirical and fundamental knowledge. The validation of simulations stresses the importance of achie…
▽ More
Hydrodynamic simulations have become irreplaceable in modern cosmology for exploring complex systems and making predictions to steer future observations. In Chapter 1, we begin with a philosophical discussion on the role of simulations in science. We argue that simulations can bridge the gap between empirical and fundamental knowledge. The validation of simulations stresses the importance of achieving a balance between trustworthiness and scepticism. Next, Chapter 2 introduces the formation of structures and comparisons between synthetic and observational data. Chapter 3 describes the production pipeline of zoom-in simulations used to model individual objects and novel methods to mitigate known shortcomings. Then, we assessed the weak scaling of the SWIFT code and found it to be one of the hydrodynamic codes with the highest parallel efficiency. In Chapter 4, we study the rotational kinetic Sunyaev-Zeldovich (rkSZ) effect for high-mass galaxy clusters from the MACSIS simulations. We find a maximum signal greater than 100 $μ$K, 30 times stronger than early predictions from self-similar models, opening prospects for future detection. In Chapter 5, we address a tension between the distribution of entropy measured from observations and predicted by simulations of groups and clusters of galaxies. We find that most recent hydrodynamic simulations systematically over-predict the entropy profiles by up to one order of magnitude, leading to profiles that are shallower and higher than the power-law-like entropy profiles that have been observed. We discuss the dependence on different hydrodynamic and sub-grid parameters using variations of the EAGLE model. Chapter 6 explores the evolution of the profiles as a function of cosmic time. We report power-law-like entropy profiles at high redshift for both objects. However, at late times, an entropy plateau develops and alters the shape of the profile.
△ Less
Submitted 18 December, 2023;
originally announced December 2023.
-
Inferring the dark matter splashback radius from cluster gas and observable profiles in the FLAMINGO simulations
Authors:
Imogen Towler,
Scott T. Kay,
Joop Schaye,
Roi Kugel,
Matthieu Schaller,
Joey Braspenning,
Willem Elbers,
Carlos S. Frenk,
Juliana Kwan,
Jaime Salcido,
Marcel P. van Daalen,
Bert Vandenbroucke,
Edoardo Altamura
Abstract:
The splashback radius, coinciding with the minimum in the dark matter radial density gradient, is thought to be a universal definition of the edge of a dark matter halo. Observational methods to detect it have traced the dark matter using weak gravitational lensing or galaxy number counts. Recent attempts have also claimed the detection of a similar feature in Sunyaev-Zel'dovich (SZ) observations…
▽ More
The splashback radius, coinciding with the minimum in the dark matter radial density gradient, is thought to be a universal definition of the edge of a dark matter halo. Observational methods to detect it have traced the dark matter using weak gravitational lensing or galaxy number counts. Recent attempts have also claimed the detection of a similar feature in Sunyaev-Zel'dovich (SZ) observations of the hot intracluster gas. Here, we use the FLAMINGO simulations to investigate whether an extremum gradient in a similar position to the splashback radius is predicted to occur in the cluster gas profiles. We find that the minimum in the gradient of the stacked 3D gas density and pressure profiles, and the maximum in the gradient of the entropy profile, broadly align with the splashback feature though there are significant differences. While the dark matter splashback radius varies with specific mass accretion rate, in agreement with previous work, the radial position of the deepest minimum in the log-slope of the gas density is more sensitive to halo mass. In addition, we show that a similar minimum is also present in projected 2D pseudo-observable profiles: emission measure (X-ray); Compton-$y$ (SZ) and surface mass density (weak lensing). We find that the latter traces the dark matter results reasonably well albeit the minimum occurs at a slightly smaller radius. While results for the gas profiles are largely insensitive to accretion rate and various observable proxies for dynamical state, they do depend on the strength of the feedback processes.
△ Less
Submitted 8 March, 2024; v1 submitted 8 December, 2023;
originally announced December 2023.
-
Galaxy cluster rotation revealed in the MACSIS simulations with the kinetic Sunyaev-Zeldovich effect
Authors:
Edoardo Altamura,
Scott T. Kay,
Jens Chluba,
Imogen Towler
Abstract:
The kinetic Sunyaev-Zeldovich (kSZ) effect has now become a clear target for ongoing and future studies of the cosmic microwave background (CMB) and cosmology. Aside from the bulk cluster motion, internal motions also lead to a kSZ signal. In this work, we study the rotational kSZ effect caused by coherent large-scale motions of the cluster medium using cluster hydrodynamic cosmological simulation…
▽ More
The kinetic Sunyaev-Zeldovich (kSZ) effect has now become a clear target for ongoing and future studies of the cosmic microwave background (CMB) and cosmology. Aside from the bulk cluster motion, internal motions also lead to a kSZ signal. In this work, we study the rotational kSZ effect caused by coherent large-scale motions of the cluster medium using cluster hydrodynamic cosmological simulations. To utilise the rotational kSZ as a cosmological probe, simulations offer some of the most comprehensive data sets that can inform the modelling of this signal. In this work, we use the MACSIS data set to investigate the rotational kSZ effect in massive clusters specifically. Based on these models, we test stacking approaches and estimate the amplitude of the combined signal with varying mass, dynamical state, redshift and map-alignment geometry. We find that the dark matter, galaxy and gas spins are generally misaligned, an effect that can cause a sub-optimal estimation of the rotational kSZ effect when based on galaxy motions. Furthermore, we provide halo-spin-mass scaling relations that can be used to build a statistical model of the rotational kSZ. The rotational kSZ contribution, which is largest in massive unrelaxed clusters ($\gtrsim$100 $μ$K), could be relevant to studies of higher-order CMB temperature signals, such as the moving lens effect. The limited mass range of the MACSIS sample strongly motivates an extended investigation of the rotational kSZ effect in large-volume simulations to refine the modelling, particularly towards lower mass and higher redshift, and provide forecasts for upcoming cosmological CMB experiments (e.g. Simons Observatory, SKA-2) and X-ray observations (e.g. Athena/X-IFU).
△ Less
Submitted 16 June, 2023; v1 submitted 15 February, 2023;
originally announced February 2023.
-
The impact of stochastic modeling on the predictive power of galaxy formation simulations
Authors:
Josh Borrow,
Matthieu Schaller,
Yannick M. Bahe,
Joop Schaye,
Aaron D. Ludlow,
Sylvia Ploeckinger,
Folkert S. J. Nobels,
Edoardo Altamura
Abstract:
All modern galaxy formation models employ stochastic elements in their sub-grid prescriptions to discretise continuous equations across the time domain. In this paper, we investigate how the stochastic nature of these models, notably star formation, black hole accretion, and their associated feedback, that act on small ($<$ kpc) scales, can back-react on macroscopic galaxy properties (e.g. stellar…
▽ More
All modern galaxy formation models employ stochastic elements in their sub-grid prescriptions to discretise continuous equations across the time domain. In this paper, we investigate how the stochastic nature of these models, notably star formation, black hole accretion, and their associated feedback, that act on small ($<$ kpc) scales, can back-react on macroscopic galaxy properties (e.g. stellar mass and size) across long ($>$ Gyr) timescales. We find that the scatter in scaling relations predicted by the EAGLE model implemented in the SWIFT code can be significantly impacted by random variability between re-simulations of the same object, even when galaxies are resolved by tens of thousands of particles. We then illustrate how re-simulations of the same object can be used to better understand the underlying model, by showing how correlations between galaxy stellar mass and black hole mass disappear at the highest black hole masses ($M_{\rm BH} > 10^8$ M$_\odot$), indicating that the feedback cycle may be interrupted by external processes. We find that although properties that are collected cumulatively over many objects are relatively robust against random variability (e.g. the median of a scaling relation), the properties of individual galaxies (such as galaxy stellar mass) can vary by up to 25\%, even far into the well-resolved regime, driven by bursty physics (black hole feedback) and mergers between galaxies. We suggest that studies of individual objects within cosmological simulations be treated with caution, and that any studies aiming to closely investigate such objects must account for random variability within their results.
△ Less
Submitted 22 September, 2023; v1 submitted 15 November, 2022;
originally announced November 2022.
-
Gas clumping and its effect on hydrostatic bias in the MACSIS simulations
Authors:
Imogen Towler,
Scott Kay,
Edoardo Altamura
Abstract:
We use the MACSIS hydrodynamical simulations to estimate the extent of gas clumping in the intracluster medium of massive galaxy clusters and how it affects the hydrostatic mass bias. By comparing the clumping to the azimuthal scatter in the emission measure, an observational proxy, we find that they both increase with radius and are larger in higher-mass and dynamically perturbed systems. Similar…
▽ More
We use the MACSIS hydrodynamical simulations to estimate the extent of gas clumping in the intracluster medium of massive galaxy clusters and how it affects the hydrostatic mass bias. By comparing the clumping to the azimuthal scatter in the emission measure, an observational proxy, we find that they both increase with radius and are larger in higher-mass and dynamically perturbed systems. Similar trends are also seen for the azimuthal temperature scatter and non-thermal pressure fraction, both of which correlate with density fluctuations, with these values also increasing with redshift. However, in agreement with recent work, we find only a weak correlation between the clumping, or its proxies, and the hydrostatic mass bias. To reduce the effect of clumping in the projected profiles, we compute the azimuthal median following recent observational studies, and find this reduces the scatter in the bias. We also attempt to correct the cluster masses by using a non-thermal pressure term and find over-corrected mass estimates ($1-b=0.86$ to $1-b=1.15$) from 3D gas profiles but improved mass estimates ($1-b=0.75$ to $1-b=0.85$) from projected gas profiles, with the caveat of systematically increased scatter. We conclude that the cluster-averaged mass bias is minimised from applying a non-thermal pressure correction ($1-b=0.85$) with more modest reductions from selecting clusters that have low clumping ($1-b=0.79$) or are dynamically relaxed ($1-b=0.80$). However, the latter selection is most effective at minimising the scatter for individual objects. Such results can be tested with next generation X-ray missions equipped with high-resolution spectrometers such as Athena.
△ Less
Submitted 1 March, 2023; v1 submitted 2 November, 2022;
originally announced November 2022.
-
EAGLE-like simulation models do not solve the entropy core problem in groups and clusters of galaxies
Authors:
Edoardo Altamura,
Scott T. Kay,
Richard G. Bower,
Matthieu Schaller,
Yannick M. Bahé,
Joop Schaye,
Josh Borrow,
Imogen Towler
Abstract:
Recent high-resolution cosmological hydrodynamic simulations run with a variety of codes systematically predict large amounts of entropy in the intra-cluster medium at low redshift, leading to flat entropy profiles and a suppressed cool-core population. This prediction is at odds with X-ray observations of groups and clusters. We use a new implementation of the EAGLE galaxy formation model to inve…
▽ More
Recent high-resolution cosmological hydrodynamic simulations run with a variety of codes systematically predict large amounts of entropy in the intra-cluster medium at low redshift, leading to flat entropy profiles and a suppressed cool-core population. This prediction is at odds with X-ray observations of groups and clusters. We use a new implementation of the EAGLE galaxy formation model to investigate the sensitivity of the central entropy and the shape of the profiles to changes in the sub-grid model applied to a suite of zoom-in cosmological simulations of a group of mass $M_{500} = 8.8 \times 10^{12}~{\rm M}_\odot$ and a cluster of mass $2.9 \times 10^{14}~{\rm M}_\odot$. Using our reference model, calibrated to match the stellar mass function of field galaxies, we confirm that our simulated groups and clusters contain hot gas with too high entropy in their cores. Additional simulations run without artificial conduction, metal cooling or AGN feedback produce lower entropy levels but still fail to reproduce observed profiles. Conversely, the two objects run without supernova feedback show a significant entropy increase which can be attributed to excessive cooling and star formation. Varying the AGN heating temperature does not greatly affect the profile shape, but only the overall normalisation. Finally, we compared runs with four AGN heating schemes and obtained similar profiles, with the exception of bipolar AGN heating, which produces a higher and more uniform entropy distribution. Our study leaves open the question of whether the entropy core problem in simulations, and particularly the lack of power-law cool-core profiles, arise from incorrect physical assumptions, missing physical processes, or insufficient numerical resolution.
△ Less
Submitted 9 February, 2023; v1 submitted 18 October, 2022;
originally announced October 2022.
-
Serendipitous discovery of a physical binary quasar at z=1.76
Authors:
E. Altamura,
S. Brennan,
A. Leśniewska,
V. Pintér,
S. N. dos Reis,
T. Pursimo,
J. P. U. Fynbo,
S. Geier,
K. E. Heintz,
P. Møller
Abstract:
Binary quasars are extremely rare objects, used to investigate clustering on very small scales at different redshifts. The cases where the two quasar components are gravitationally bound, known as physical binary quasars, can also exhibit enhanced astrophysical activity and therefore are of particular scientific interest. Here we present the serendipitous discovery of a physical pair of quasars wi…
▽ More
Binary quasars are extremely rare objects, used to investigate clustering on very small scales at different redshifts. The cases where the two quasar components are gravitationally bound, known as physical binary quasars, can also exhibit enhanced astrophysical activity and therefore are of particular scientific interest. Here we present the serendipitous discovery of a physical pair of quasars with an angular separation of $Δθ= (8.76 \pm 0.11)$ arcsec. The redshifts of the two quasars are consistent within the errors and measured as $z = (1.76\pm 0.01)$. Under the motivated assumption that the pair does not arise from a single gravitationally lensed quasar, the resulting projected physical separation was estimated as $(76 \pm 1)$ kpc. For both targets we detected Si VI, C VI, C III], and Mg II emission lines. However, the two quasars show significantly different optical colours, one being among the most reddened quasars at $z > 1.5$ and the other with colours consistent with typical quasar colours at the same redshift. Therefore it is ruled out that the sources are a lensed system. This is our second serendipitous discovery of a pair of two quasars with different colours, having a separation $\lesssim 10$ arcsec, which extends the very limited catalogue of known quasar pairs. We ultimately argue that the number of binary quasars may have been significantly underestimated in previous photometric surveys, due to the bias arising from paired quasars with very different colours.
△ Less
Submitted 10 October, 2019;
originally announced October 2019.
