Source author record

Susana Planelles

Susana Planelles appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record

Catalog footprint

What is connected

15works
5topics
4close collaborators

Actions

Connect this record

Log in to claim

Research graph

See the researcher in context

Open full explorer

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

15 published item(s)

preprint2025arXiv

The signature of major mergers on the hydrostatic mass bias of galaxy clusters

While galaxy cluster masses are fundamental cosmological observables, estimates based on intra-cluster medium observations rely on hydrostatic equilibrium, introducing a systematic bias. We investigate how mergers drive the time evolution of this hydrostatic mass bias, identifying the dominant physical mechanisms and their dependence on dynamical state and merger history. Using a high-resolution AMR Eulerian+$N$-body cosmological simulation, we analyse a sample of cluster mergers within $1.5 \leq z \leq 0$, comparing true and hydrostatic masses derived from gas density and temperature profiles, and tracing their evolution. At $z=0$, the hydrostatic mass bias shows a mild correlation with dynamical state. During major mergers, the bias follows a characteristic trend: a sharp negative dip around the merger time, a transient positive peak, and a gradual return to pre-merger levels. This behaviour is primarily driven by morphological and dynamical reconfigurations of the gas density within the ICM, while thermodynamical processes play a secondary role. The pattern shows no strong dependence on secondary parameters, such as mass ratio or impact parameter, but it can be fitted to a simple time-dependent functional form. This trend is present at radii $r\le R_{\mathrm{vir}}$, although with reduced amplitude and shorter timescales as the radius decreases. Hydrostatic mass bias is closely linked, albeit in a non-trivial way, with the merging history of galaxy clusters. We find that the bias values are weakly correlated with the dynamical state of clusters. Nevertheless, our results give a robust estimation of the hydrostatic mass bias values in the pre-merger, merging, and post-merger phases. These findings highlight the importance of delving deeper into the observational assessment of cluster assembly state in order to improve mass estimations for cosmological analyses.

preprint2023arXiv

On the choice of the most suitable indicator for the assembly state of dark matter haloes through cosmic time

The dynamical state and morphological features of galaxies and galaxy clusters, and their high-redshift precursors, are tightly connected with their assembly history, encoding crucial information about the formation and evolution of such cosmic structures. As a first step towards finding an optimal indicator of the assembly state of observed structures, we use a cosmological simulation of a moderate volume to critically examine the best definition of an indicator that is able to discriminate dark matter haloes undergoing mergers and/or strong accretion from haloes experimenting a relaxed evolution. Using a combination of centre offset, virial ratio, mean radial velocity, sparsity and ellipticity of the dark matter halo, we study how the thresholds on these parameters, as well as their relative weights, should evolve with redshift to provide the best classification possible. This allows us to split a sample of haloes in a totally relaxed, a marginally relaxed and an unrelaxed subsamples. The resulting classification strongly correlates with the merging activity obtained from the analysis of complete merger trees extracted from whole simulation data. The results on how the different indicators depend on redshift and halo mass, and their optimal combination to better match the true assembly history of haloes, could constitute relevant hints to find a suitable set of indicators applicable to observational data.

preprint2022arXiv

The halo finding problem revisited: a deep revision of the ASOHF code

Context. New-generation cosmological simulations are providing huge amounts of data, whose analysis becomes itself a cutting-edge computational problem. In particular, the identification of gravitationally bound structures, known as halo finding, is one of the main analyses. A handful of codes developed to tackle this task have been presented during the last years. Aims. We present a deep revision of the already existing code ASOHF. The algorithm has been throughfully redesigned in order to improve its capabilities to find bound structures and substructures, both using dark matter particles and stars, its parallel performance, and its abilities to handle simulation outputs with vast amounts of particles. This upgraded version of ASOHF is conceived to be a publicly available tool. Methods. A battery of idealised and realistic tests are presented in order to assess the performance of the new version of the halo finder. Results. In the idealised tests, ASOHF produces excellent results, being able to find virtually all the structures and substructures placed within the computational domain. When applied to realistic data from simulations, the performance of our finder is fully consistent with the results from other commonly used halo finders, with remarkable performance in substructure detection. Besides, ASOHF turns out to be extremely efficient in terms of computational cost. Conclusions. We present a public, deeply revised version of the ASOHF halo finder. The new version of the code produces remarkable results finding haloes and subhaloes in cosmological simulations, with an excellent parallel performance and with a contained computational cost.

preprint2021arXiv

Unravelling cosmic velocity flows: a Helmholtz-Hodge decomposition algorithm for cosmological simulations

In the context of intra-cluster medium turbulence, it is essential to be able to split the turbulent velocity field in a compressive and a solenoidal component. We describe and implement a new method for this aim, i.e., performing a Helmholtz-Hodge decomposition, in multi-grid, multi-resolution descriptions, focusing on (but not being restricted to) the outputs of AMR cosmological simulations. The method is based on solving elliptic equations for a scalar and a vector potential, from which the compressive and the solenoidal velocity fields, respectively, are derived through differentiation. These equations are addressed using a combination of Fourier (for the base grid) and iterative (for the refinement grids) methods. We present several idealised tests for our implementation, reporting typical median errors in the order of $1\unicode{x2030}$-$1\%$, and with 95-percentile errors below a few percents. Additionally, we also apply the code to the outcomes of a cosmological simulation, achieving similar accuracy at all resolutions, even in the case of highly non-linear velocity fields. We finally take a closer look to the decomposition of the velocity field around a massive galaxy cluster.

preprint2020arXiv

Cosmic magnetic fields with MASCLET: an application to galaxy clusters

We describe and test a new version of the adaptive mesh refinement (AMR) cosmological code MASCLET. The new version of the code includes all the ingredients of its previous version plus a description of the evolution of the magnetic field under the approximation of the ideal magneto-hydrodynamics (MHD). To preserve the divergence-free condition of MHD, the original divergence cleaning algorithm of Dedner et al. (2002) is implemented. We present a set of well-known 1D and 2D tests, such as several shock-tube problems, the fast rotor and the Orszag-Tang vortex. The performance of the code in all the tests is excellent with estimated median relative errors of $\nabla \cdot {\bf B}$ in the 2D tests smaller than $5 \times 10^{-5}$ for the fast rotor test, and $5 \times 10^{-3}$ for the Orszag-Tang vortex. As an astrophysical application of the code, we present a simulation of a cosmological box of 40 comoving Mpc side length in which a primordial uniform comoving magnetic field of strength 0.1 nG is seeded. The simulation shows how the magnetic field is channelled along the filaments of gas and is concentrated and amplified within galaxy clusters. Comparison with the values expected from pure compression reveals an additional amplification of the magnetic field caused by turbulence in the central region of the cluster. Values of the order of $\sim 1 μ$G are obtained in clusters at $z\sim 0$ with median relative errors of $\nabla \cdot {\bf B}$ below 0.4\%. The implications of a proper description of the dynamics of the magnetic field and their possible observational counterparts in future facilities are discussed.

preprint2020arXiv

On the accretion history of galaxy clusters: temporal and spatial distribution

We analyse the results of an Eulerian AMR cosmological simulation in order to quantify the mass growth of galaxy clusters, exploring the differences between dark matter and baryons. We have determined the mass assembly histories (MAHs) of each of the mass components and computed several proxies for the instantaneous mass accretion rate (MAR). The mass growth of both components is clearly dominated by the contribution of major mergers, but high MARs can also occur during smooth accretion periods. We explored the correlations between MARs, merger events and clusters' environments, finding the mean densities in $1 \leq r/R_{200m} \leq 1.5$ to correlate strongly with $Γ_{200m}$ in massive clusters which undergo major mergers through their MAH. From the study of the dark matter velocity profiles, we find a strong anticorrelation between the MAR proxies $Γ_{200m}$ and $α_{200m}$. Last, we present a novel approach to study the angularly-resolved distribution of gas accretion flows in simulations, which allows to extract and interpret the main contributions to the accretion picture and to assess systematic differences between the thermodynamical properties of each of these contributions using multipolar analysis. We have preliminarily applied the method to the best numerically-resolved cluster in our simulation. Amongst the most remarkable results, we find that the gas infalling through the cosmic filaments has systematically lower entropy compared to the isotropic component, but we do not find a clear distinction in temperature.

preprint2017arXiv

Simulation-based marginal likelihood for cluster strong lensing cosmology

Comparisons between observed and predicted strong lensing properties of galaxy clusters have been routinely used to claim either tension or consistency with $Λ$CDM cosmology. However, standard approaches to such cosmological tests are unable to quantify the preference for one cosmology over another. We advocate approximating the relevant Bayes factor using a marginal likelihood that is based on the following summary statistic: the posterior probability distribution function for the parameters of the scaling relation between Einstein radii and cluster mass, $α$ and $β$. We demonstrate, for the first time, a method of estimating the marginal likelihood using the X-ray selected $z>0.5$ MACS clusters as a case in point and employing both N-body and hydrodynamic simulations of clusters. We investigate the uncertainty in this estimate and consequential ability to compare competing cosmologies, that arises from incomplete descriptions of baryonic processes, discrepancies in cluster selection criteria, redshift distribution, and dynamical state. The relation between triaxial cluster masses at various overdensities provide a promising alternative to the strong lensing test.

preprint2015arXiv

Neutral hydrogen in galaxy clusters: impact of AGN feedback and implications for intensity mapping

By means of zoom-in hydrodynamic simulations we quantify the amount of neutral hydrogen (HI) hosted by groups and clusters of galaxies. Our simulations, which are based on an improved formulation of smoothed particle hydrodynamics (SPH), include radiative cooling, star formation, metal enrichment and supernova feedback, and can be split in two different groups, depending on whether feedback from active galactic nuclei (AGN) is turned on or off. Simulations are analyzed to account for HI self-shielding and the presence of molecular hydrogen. We find that the mass in neutral hydrogen of dark matter halos monotonically increases with the halo mass and can be well described by a power-law of the form $M_{\rm HI}(M,z)\propto M^{3/4}$. Our results point out that AGN feedback reduces both the total halo mass and its HI mass, although it is more efficient in removing HI. We conclude that AGN feedback reduces the neutral hydrogen mass of a given halo by $\sim50\%$, with a weak dependence on halo mass and redshift. The spatial distribution of neutral hydrogen within halos is also affected by AGN feedback, whose effect is to decrease the fraction of HI that resides in the halo inner regions. By extrapolating our results to halos not resolved in our simulations we derive astrophysical implications from the measurements of $Ω_{\rm HI}(z)$: halos with circular velocities larger than $\sim25~{\rm km/s}$ are needed to host HI in order to reproduce observations. We find that only the model with AGN feedback is capable of reproducing the value of $Ω_{\rm HI}b_{\rm HI}$ derived from available 21cm intensity mapping observations.

preprint2014arXiv

Subhaloes gone Notts: Subhaloes as tracers of the dark matter halo shape

We study the shapes of subhalo distributions from four dark-matter-only simulations of Milky Way type haloes. Comparing the shapes derived from the subhalo distributions at high resolution to those of the underlying dark matter fields we find the former to be more triaxial if theanalysis is restricted to massive subhaloes. For three of the four analysed haloes the increased triaxiality of the distributions of massive subhaloes can be explained by a systematic effect caused by the low number of objects. Subhaloes of the fourth halo show indications for anisotropic accretion via their strong triaxial distribution and orbit alignment with respect to the dark matter field. These results are independent of the employed subhalo finder. Comparing the shape of the observed Milky Way satellite distribution to those of high-resolution subhalo samples from simulations, we find an agreement for samples of bright satellites, but significant deviations if faint satellites are included in the analysis. These deviations might result from observational incompleteness.

preprint2013arXiv

Baryon Census in Hydrodynamical Simulations of Galaxy Clusters

We carry out an analysis of a set of cosmological SPH hydrodynamical simulations of galaxy clusters and groups aimed at studying the total baryon budget in clusters, and how this budget is shared between the hot diffuse component and the stellar component. Using the TreePM+SPH GADGET-3 code, we carried out one set of non-radiative simulations, and two sets of simulations including radiative cooling, star formation and feedback from supernovae (SN), one of which also accounting for the effect of feedback from active galactic nuclei (AGN). The analysis is carried out with the twofold aim of studying the implication of stellar and hot gas content on the relative role played by SN and AGN feedback, and to calibrate the cluster baryon fraction and its evolution as a cosmological tool. We find that both radiative simulation sets predict a trend of stellar mass fraction with cluster mass that tends to be weaker than the observed one. However this tension depends on the particular set of observational data considered. Including the effect of AGN feedback alleviates this tension on the stellar mass and predicts values of the hot gas mass fraction and total baryon fraction to be in closer agreement with observational results. We further compute the ratio between the cluster baryon content and the cosmic baryon fraction, Y_b, as a function of cluster-centric radius and redshift. At R_500 we find for massive clusters with M_500>2\times10^{14} h^{-1} M_sun that Y_b is nearly independent of the physical processes included and characterized by a negligible redshift evolution: Y_{b,500}=0.85+/-0.03 with the error accounting for the intrinsic r.m.s. scatter within the set of simulated clusters. At smaller radii, R_2500, the typical value of Y_b slightly decreases, by an amount that depends on the physics included in the simulations, while its scatter increases by about a factor of two.

preprint2013arXiv

Structure Finding in Cosmological Simulations: The State of Affairs

The ever increasing size and complexity of data coming from simulations of cosmic structure formation demands equally sophisticated tools for their analysis. During the past decade, the art of object finding in these simulations has hence developed into an important discipline itself. A multitude of codes based upon a huge variety of methods and techniques have been spawned yet the question remained as to whether or not they will provide the same (physical) information about the structures of interest. Here we summarize and extent previous work of the "halo finder comparison project": we investigate in detail the (possible) origin of any deviations across finders. To this extent we decipher and discuss differences in halo finding methods, clearly separating them from the disparity in definitions of halo properties. We observe that different codes not only find different numbers of objects leading to a scatter of up to 20 per cent in the halo mass and Vmax function, but also that the particulars of those objects that are identified by all finders differ. The strength of the variation, however, depends on the property studied, e.g. the scatter in position, bulk velocity, mass, and the peak value of the rotation curve is practically below a few per cent, whereas derived quantities such as spin and shape show larger deviations. Our study indicates that the prime contribution to differences in halo properties across codes stems from the distinct particle collection methods and -- to a minor extent -- the particular aspects of how the procedure for removing unbound particles is implemented. We close with a discussion of the relevance and implications of the scatter across different codes for other fields such as semi-analytical galaxy formation models, gravitational lensing, and observables in general.

preprint2012arXiv

Cosmological shock waves: clues to the formation history of haloes

Shock waves developed during the formation and evolution of cosmic structures encode crucial information on the hierarchical formation of the Universe. We analyze an Eulerian AMR hydro + N-body simulation in a $Λ$CDM cosmology focused on the study of cosmological shock waves. The combination of a shock-capturing algorithm together with the use of a halo finder allows us to study the morphological structures of the shock patterns, the statistical properties of shocked cells, and the correlations between the cosmological shock waves appearing at different scales and the properties of the haloes harbouring them. The shocks in the simulation can be split into two broad classes: internal weak shocks related with evolutionary events within haloes, and external strong shocks associated with large-scale events. The shock distribution function contains information on the abundances and strength of the different shocks, and it can be fitted by a double power law with a break in the slope around a Mach number of 20. We introduce a generalised scaling relation that correlates the average Mach numbers within the virial radius of haloes and their virial masses. In this plane, haloes occupy different areas according to their early evolutionary histories: those with a quiet evolution have an almost constant Mach number independently of their masses, whereas haloes undergoing significant merger events very early in their evolution show a linear dependence with their masses. At high redshift, the halo distribution in this scaling relation forms a L-like pattern that changes due to the evolution of the haloes. The analysis of the propagation speed and size of the shock waves around haloes could give some hints on the formation time and main features of the haloes. (Abridged)

preprint2011arXiv

Haloes gone MAD: The Halo-Finder Comparison Project

[abridged] We present a detailed comparison of fundamental dark matter halo properties retrieved by a substantial number of different halo finders. These codes span a wide range of techniques including friends-of-friends (FOF), spherical-overdensity (SO) and phase-space based algorithms. We further introduce a robust (and publicly available) suite of test scenarios that allows halo finder developers to compare the performance of their codes against those presented here. This set includes mock haloes containing various levels and distributions of substructure at a range of resolutions as well as a cosmological simulation of the large-scale structure of the universe. All the halo finding codes tested could successfully recover the spatial location of our mock haloes. They further returned lists of particles (potentially) belonging to the object that led to coinciding values for the maximum of the circular velocity profile and the radius where it is reached. All the finders based in configuration space struggled to recover substructure that was located close to the centre of the host halo and the radial dependence of the mass recovered varies from finder to finder. Those finders based in phase space could resolve central substructure although they found difficulties in accurately recovering its properties. Via a resolution study we found that most of the finders could not reliably recover substructure containing fewer than 30-40 particles. However, also here the phase space finders excelled by resolving substructure down to 10-20 particles. By comparing the halo finders using a high resolution cosmological volume we found that they agree remarkably well on fundamental properties of astrophysical significance (e.g. mass, position, velocity, and peak of the rotation curve).

preprint2010arXiv

ASOHF: a new adaptive spherical overdensity halo finder

We present and test a new halo finder based on the spherical overdensity (SO) method. This new adaptive spherical overdensity halo finder (ASOHF) is able to identify dark matter haloes and their substructures (subhaloes) down to the scales allowed by the analysed simulations. The code has been especially designed for the adaptive mesh refinement cosmological codes, although it can be used as a stand-alone halo finder for N-body codes. It has been optimised for the purpose of building the merger tree of the haloes. In order to verify the viability of this new tool, we have developed a set of bed tests that allows us to estimate the performance of the finder. Finally, we apply the halo finder to a cosmological simulation and compare the results obtained to those given by other well known publicly available halo finders.

preprint2009arXiv

Galaxy cluster mergers

We present the results of an Eulerian adaptive mesh refinement (AMR) hydrodynamical and N-body simulation in a $Λ$CDM cosmology. The simulation incorporates common cooling and heating processes for primordial gas. A specific halo finder has been designed and applied in order to extract a sample of galaxy clusters directly obtained from the simulation without considering any resimulating scheme. We have studied the evolutionary history of the cluster halos, and classified them into three categories depending on the merger events they have undergone: major mergers, minor mergers, and relaxed clusters. The main properties of each one of these classes and the differences among them are discussed. The collisions among galaxy clusters are produced naturally by the non-linear evolution in the simulated cosmological volume, no controlled collisions have been considered. We pay special attention to discuss the role of merger events as a source of feedback and reheating, and their effects on the existence of cool cores in galaxy clusters, as well as in the scaling relations.