Showing 1–39 of 39 results for author: Planelles, S

Vortex-p: a Helmholtz-Hodge and Reynolds decomposition algorithm for particle-based simulations

Authors: David Vallés-Pérez, Susana Planelles, Vicent Quilis, Frederick Groth, Tirso Marin-Gilabert, Klaus Dolag

Abstract: Astrophysical turbulent flows display an intrinsically multi-scale nature, making their numerical simulation and the subsequent analyses of simulated data a complex problem. In particular, two fundamental steps in the study of turbulent velocity fields are the Helmholtz-Hodge decomposition (compressive+solenoidal; HHD) and the Reynolds decomposition (bulk+turbulent; RD). These problems are relativ… ▽ More Astrophysical turbulent flows display an intrinsically multi-scale nature, making their numerical simulation and the subsequent analyses of simulated data a complex problem. In particular, two fundamental steps in the study of turbulent velocity fields are the Helmholtz-Hodge decomposition (compressive+solenoidal; HHD) and the Reynolds decomposition (bulk+turbulent; RD). These problems are relatively simple to perform numerically for uniformly-sampled data, such as the one emerging from Eulerian, fix-grid simulations; but their computation is remarkably more complex in the case of non-uniformly sampled data, such as the one stemming from particle-based or meshless simulations. In this paper, we describe, implement and test vortex-p, a publicly available tool evolved from the vortex code, to perform both these decompositions upon the velocity fields of particle-based simulations, either from smoothed particle hydrodynamics (SPH), moving-mesh or meshless codes. The algorithm relies on the creation of an ad-hoc adaptive mesh refinement (AMR) set of grids, on which the input velocity field is represented. HHD is then addressed by means of elliptic solvers, while for the RD we adapt an iterative, multi-scale filter. We perform a series of idealised tests to assess the accuracy, convergence and scaling of the code. Finally, we present some applications of the code to various SPH and meshless finite-mass (MFM) simulations of galaxy clusters performed with OpenGadget3, with different resolutions and physics, to showcase the capabilities of the code. △ Less

Submitted 2 July, 2024; originally announced July 2024.

Comments: Accepted for publication in Computer Physics Communications. 20 pages, 13 figures + appendices. The code is available in https://github.com/dvallesp/vortex-p and the documentation can be accessed in https://vortex-particles.github.io

Journal ref: Computer Physics Communications, Volume 304, November 2024, 109305

Cosmic accretion shocks as a tool to measure the dark matter mass of galaxy clusters

Authors: David Vallés-Pérez, Vicent Quilis, Susana Planelles

Abstract: Cosmological accretion shocks created during the formation of galaxy clusters are a ubiquitous phenomenon all around the Universe. These shocks, and their features, are intimately related with the gravitational energy put into play during galaxy cluster formation. Studying a sample of simulated galaxy clusters and their associated accretion shocks, we show that objects in our sample sit in a plane… ▽ More Cosmological accretion shocks created during the formation of galaxy clusters are a ubiquitous phenomenon all around the Universe. These shocks, and their features, are intimately related with the gravitational energy put into play during galaxy cluster formation. Studying a sample of simulated galaxy clusters and their associated accretion shocks, we show that objects in our sample sit in a plane within the three dimensional-space of cluster total mass, shock radius, and Mach number (a measure of shock intensity). Using this relation, and considering that forthcoming new observations will be able to measure shock radii and intensities, we put forward the idea that the dark matter content of galaxy clusters could be indirectly measured with an error up to around 30 per cent at the $1σ$ confidence level. This procedure would be a new and independent method to measure the dark matter mass in cosmic structures, and a novel constraint to the accepted $Λ$CDM paradigm. △ Less

Submitted 1 July, 2024; originally announced July 2024.

Comments: Accepted for publication in Nature Astronomy. This version corresponds to the preprint, prior to any formal peer review. The published version, including peer-review improvements, is freely accessible through the SharedIt platform in https://rdcu.be/dMm5v. This version: 13 pages, 3 figures (main paper + methods), plus supplementary material (9 pages, 10 figures)

CAVITY, Calar Alto Void Integral-field Treasury surveY and project extension

Authors: I. Pérez, S. Verley, L. Sánchez-Menguiano, T. Ruiz-Lara, R. García-Benito, S. Duarte Puertas, A. Jiménez, J. Domínguez-Gómez, D. Espada, R. F. Peletier, J. Román, M. I. Rodríguez, P. Sánchez Alarcón, M. Argudo-Fernández, G. Torres-Ríos, B. Bidaran, M. Alcázar-Laynez, R. van de Weygaert, S. F. Sánchez, U. Lisenfeld, A. Zurita, E. Florido, J. M. van der Hulst, G. Blázquez-Calero, P. Villalba-González , et al. (36 additional authors not shown)

Abstract: We have learnt in the last decades that the majority of galaxies belong to high density regions interconnected in a sponge-like fashion. This large-scale structure is characterised by clusters, filaments, walls, where most galaxies concentrate, but also under-dense regions, called voids. The void regions and the galaxies within represent an ideal place for the study of galaxy formation and evoluti… ▽ More We have learnt in the last decades that the majority of galaxies belong to high density regions interconnected in a sponge-like fashion. This large-scale structure is characterised by clusters, filaments, walls, where most galaxies concentrate, but also under-dense regions, called voids. The void regions and the galaxies within represent an ideal place for the study of galaxy formation and evolution as they are largely unaffected by the complex physical processes that transform galaxies in high-density environments. These void galaxies can hold the key as well to answer current challenges to the $Λ$CDM paradigm. The Calar Alto Void Integral-field Treasury surveY (CAVITY) is a Legacy project approved by the Calar Alto Observatory to obtain spatially resolved spectroscopic information of $\sim300$ void galaxies in the Local Universe (0.005 < z < 0.050) covering from -17.0 to -21.5 in $\rm r$ band absolute magnitude. It officially started in January 2021 and has been awarded 110 useful dark observing nights at the 3.5 m telescope using the PMAS spectrograph. Complementary follow-up projects including deep optical imaging, integrated, as well as resolved CO data, and integrated HI spectra, have joint the PMAS observations and naturally complete the scientific aim of characterising galaxies in cosmic voids. The extension data has been denominated CAVITY+. The data will be available to the whole community in different data releases, the first of which is planned for July 2024, and it will provide the community with PMAS data cubes for around 100 void galaxies through a user friendly, and well documented, database platform. We present here the survey, sample selection, data reduction, quality control schemes, science goals, and some examples of the scientific power of the CAVITY and CAVITY+ data. △ Less

Submitted 24 May, 2024; v1 submitted 7 May, 2024; originally announced May 2024.

Comments: 19 pages, 11 figures. Accepted for publication in A&A - Replacement after A&A minor language edition

The imprints of galaxy cluster internal dynamics on the Sunyaev-Zeldovich effect

Authors: Óscar Monllor-Berbegal, David Vallés-Pérez, Susana Planelles, Vicent Quilis

Abstract: Forthcoming measurements of the Sunyaev-Zeldovich (SZ) effect in galaxy clusters will dramatically improve our understanding of the main intra-cluster medium (ICM) properties and how they depend on the particular thermal and dynamical state of the associated clusters. Using a sample of simulated galaxy clusters we assess the impact of the ICM internal dynamics on both the thermal and kinetic SZ ef… ▽ More Forthcoming measurements of the Sunyaev-Zeldovich (SZ) effect in galaxy clusters will dramatically improve our understanding of the main intra-cluster medium (ICM) properties and how they depend on the particular thermal and dynamical state of the associated clusters. Using a sample of simulated galaxy clusters we assess the impact of the ICM internal dynamics on both the thermal and kinetic SZ effects (tSZ and kSZ, respectively). We produce synthetic maps of the SZ effect, for the simulated clusters. For each galaxy cluster in the sample, its dynamical state is estimated by using a combination of well-established indicators. We use the correlations between SZ maps and cluster dynamical state, to look for the imprints of the evolutionary events, mainly mergers, on the SZ signals. The kinetic effect shows a remarkable correlation with the dynamical state: unrelaxed clusters present a higher radial profile and an overall stronger signal at all masses and radii. Furthermore, the kSZ signal is correlated with rotation for relaxed clusters, while for the disturbed systems the effect is dominated by other motions such as bulk flows, turbulence, etc. The kSZ effect shows a dipolar pattern when averaging over cluster dynamical classes, especially for the relaxed population. This feature can be exploited to stack multiple kSZ maps in order to recover a stronger dipole signal that would be correlated with the global rotation properties of the sample. The SZ effect can be used as a tool to estimate the dynamical state of galaxy clusters, especially to segregate those clusters with a quiescent evolution from those with a rich record of recent merger events. Our results suggest that the forthcoming observational data measuring the SZ signal in clusters could be used as a complementary strategy to classify the evolutionary history of galaxy clusters. △ Less

Submitted 25 March, 2024; originally announced March 2024.

Comments: 13 pages. 11 figures. Accepted for publication in A&A

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

Authors: David Vallés-Pérez, Susana Planelles, Óscar Monllor-Berbegal, Vicent Quilis

Abstract: 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 modera… ▽ More 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. △ Less

Submitted 5 January, 2023; originally announced January 2023.

Comments: 15 pages, 10 figures; accepted for publication in MNRAS

Sociology and hierarchy of voids: A study of seven nearby CAVITY galaxy voids and their dynamical CosmicFlows-3 environment

Authors: H. M. Courtois, R. van de Weygaert, M. Aubert, D. Pomarède, D. Guinet, J. Domínguez-Gómez, E. Florido, L. Galbany, R. García-Benito, J. M. van der Hulst, K. Kreckel, R. E. Miura, I. Pérez, S. Planelles, V. Quilis, J. Román, M. Sánchez-Portal

Abstract: Context. The present study addresses a key question related to our understanding of the relation between void galaxies and their environment: the relationship between luminous and dark matter in and around voids. Aims. To explore the extent to which local Universe voids are empty of matter, we study the full (dark+luminous) matter content of seven nearby cosmic voids that are fully contained withi… ▽ More Context. The present study addresses a key question related to our understanding of the relation between void galaxies and their environment: the relationship between luminous and dark matter in and around voids. Aims. To explore the extent to which local Universe voids are empty of matter, we study the full (dark+luminous) matter content of seven nearby cosmic voids that are fully contained within the CosmicFlows-3 volume. Methods. We obtained the matter-density profiles of seven cosmic voids using two independent methods. These were built from the galaxy redshift space two-point correlation function in conjunction with peculiar velocity gradients from the CosmicFlows-3 dataset. Results. The results are striking, because when the redshift survey is used, all voids show a radial positive gradient of galaxies, while based on the dynamical analysis, only three of these voids display a clear underdensity of matter in their center. Conclusions. This work constitutes the most detailed observational analysis of voids conducted so far, and shows that void emptiness should be derived from dynamical information. From this limited study, the Hercules void appears to be the best candidate for a local Universe pure "pristine volume", expanding in three directions with no dark matter located in that void. △ Less

Submitted 8 March, 2023; v1 submitted 29 November, 2022; originally announced November 2022.

Comments: Submitted A\&A Nov 29, 2022 - AA/2022/45578 / Accepted March 3rd, 2023

Journal ref: A&A 673, A38 (2023)

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

Authors: David Vallés-Pérez, Susana Planelles, Vicent Quilis

Abstract: 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 revisi… ▽ More 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. △ Less

Submitted 4 May, 2022; originally announced May 2022.

Comments: 21 pages, 15 figures. Accepted for publication in A&A. Find the code repository in https://github.com/dvallesp/ASOHF and the code documentation in https://asohf.github.io/

Journal ref: A&A 664, A42 (2022)

Void replenishment: how voids accrete matter over cosmic history

Authors: David Vallés-Pérez, Vicent Quilis, Susana Planelles

Abstract: Cosmic voids are underdense regions filling up most of the volume in the Universe. They are expected to emerge in regions comprising negative initial density fluctuations, and subsequently expand as the matter around them collapses and forms walls, filaments and clusters. We report results from the analysis of a cosmological simulation specially designed to accurately describe low-density regions,… ▽ More Cosmic voids are underdense regions filling up most of the volume in the Universe. They are expected to emerge in regions comprising negative initial density fluctuations, and subsequently expand as the matter around them collapses and forms walls, filaments and clusters. We report results from the analysis of a cosmological simulation specially designed to accurately describe low-density regions, such as cosmic voids. Contrary to the common expectation, we find that voids also experience significant mass inflows over cosmic history. On average, $10\%$ of the mass of voids in the sample at $z \sim 0$ is accreted from overdense regions, reaching values beyond $35\%$ for a significant fraction of voids. More than half of the mass entering the voids lingers on periods of time $\sim 10 \, \mathrm{Gyr}$ well inside them, reaching inner radii. This would imply that part of the gas lying inside voids at a given time proceeds from overdense regions (e.g., clusters or filaments), where it could have been pre-processed, thus challenging the scenario of galaxy formation in voids, and dissenting from the idea of being pristine environments. △ Less

Submitted 16 September, 2021; originally announced September 2021.

Comments: Accepted for publication in The Astrophysical Journal Letters (ApJL). 11 pages (main text: 6 pages), 3 figures and 3 supplementary figures included in the PDF; the movie on Supplementary Figure 3 can be accessed on https://www.uv.es/davape3/voids/figS3_animation.mp4

Journal ref: The Astrophysical Journal Letters, vol 920, L2 (2021)

Exploring the role of cosmological shock waves in the Dianoga simulations of galaxy clusters

Authors: Susana Planelles, Stefano Borgani, Vicent Quilis, Giuseppe Murante, Veronica Biffi, Elena Rasia, Klaus Dolag, Gian Luigi Granato, Cinthia Ragone-Figueroa

Abstract: Cosmological shock waves are ubiquitous to cosmic structure formation and evolution. As a consequence, they play a major role in the energy distribution and thermalization of the intergalactic medium (IGM). We analyze the Mach number distribution in the Dianoga simulations of galaxy clusters performed with the SPH code GADGET-3. The simulations include the effects of radiative cooling, star format… ▽ More Cosmological shock waves are ubiquitous to cosmic structure formation and evolution. As a consequence, they play a major role in the energy distribution and thermalization of the intergalactic medium (IGM). We analyze the Mach number distribution in the Dianoga simulations of galaxy clusters performed with the SPH code GADGET-3. The simulations include the effects of radiative cooling, star formation, metal enrichment, supernova and active galactic nuclei feedback. A grid-based shock-finding algorithm is applied in post-processing to the outputs of the simulations. This procedure allows us to explore in detail the distribution of shocked cells and their strengths as a function of cluster mass, redshift and baryonic physics. We also pay special attention to the connection between shock waves and the cool-core/non-cool core (CC/NCC) state and the global dynamical status of the simulated clusters. In terms of general shock statistics, we obtain a broad agreement with previous works, with weak (low-Mach number) shocks filling most of the volume and processing most of the total thermal energy flux. As a function of cluster mass, we find that massive clusters seem more efficient in thermalising the IGM and tend to show larger external accretion shocks than less massive systems. We do not find any relevant difference between CC and NCC clusters. However, we find a mild dependence of the radial distribution of the shock Mach number on the cluster dynamical state, with disturbed systems showing stronger shocks than regular ones throughout the cluster volume. △ Less

Submitted 22 August, 2021; originally announced August 2021.

Comments: 19 pages, 15 figures, accepted for publication in MNRAS

Troubled cosmic flows: turbulence, enstrophy and helicity from the assembly history of the intracluster medium

Authors: David Vallés-Pérez, Susana Planelles, Vicent Quilis

Abstract: Both simulations and observations have shown that turbulence is a pervasive phenomenon in cosmic scenarios, yet it is particularly difficult to model numerically due to its intrinsically multiscale character which demands high resolutions. Additionally, turbulence is tightly connected to the dynamical state and the formation history of galaxies and galaxy clusters, producing a diverse phenomenlogy… ▽ More Both simulations and observations have shown that turbulence is a pervasive phenomenon in cosmic scenarios, yet it is particularly difficult to model numerically due to its intrinsically multiscale character which demands high resolutions. Additionally, turbulence is tightly connected to the dynamical state and the formation history of galaxies and galaxy clusters, producing a diverse phenomenlogy which requires large samples of such structures to attain robust conclusions. In this work, we use an adaptive mesh refinement (AMR) cosmological simulation to explore the generation and dissipation of turbulence in galaxy clusters, in connection to its assembly history. We find that major mergers, and more generally accretion of gas, is the main process driving turbulence in the intracluster medium. We have especially focused on solenoidal turbulence, which can be quantified through enstrophy. Our results seem to confirm a scenario for its generation which involves baroclinicity and compression at the external (accretion) and internal (merger) shocks, followed by vortex stretching downstream of them. We have also looked at the infall of mass to the cluster beyond its virial boundary, finding that gas follows trajectories with some degree of helicity, as it has already developed some vorticity in the external shocks. △ Less

Submitted 7 April, 2021; v1 submitted 24 March, 2021; originally announced March 2021.

Comments: 19 pages, 12 figures; accepted for publication in MNRAS

Journal ref: MNRAS, 504-1 (2021), 510-527

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

Authors: David Vallés-Pérez, Susana Planelles, Vicent Quilis

Abstract: 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… ▽ More 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. △ Less

Submitted 11 February, 2021; originally announced February 2021.

Comments: 12 pages, 8 figures; accepted for publication in Computer Physics Communications

Journal ref: Computer Physics Communications, Volume 263, June 2021, 107892

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

Authors: David Vallés-Pérez, Susana Planelles, Vicent Quilis

Abstract: 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… ▽ More 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. △ Less

Submitted 29 September, 2020; originally announced September 2020.

Comments: accepted for publication in MNRAS

Journal ref: MNRAS, 499-2 (2020), 2303-2318

Cosmic magnetic fields with MASCLET: an application to galaxy clusters

Authors: Vicent Quilis, José Ma. Martí, Susana Planelles

Abstract: 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… ▽ More 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. △ Less

Submitted 27 March, 2020; originally announced March 2020.

Comments: 14 pages, 13 figures, accepted for publication in MNRAS

The large-scale environment from cosmological simulations II: The redshift evolution and distributions of baryons

Authors: Weiguang Cui, Alexander Knebe, Noam I. Libeskind, Susana Planelles, Xiaohu Yang, Wei Cui, Romeel Dave, Xi Kang, Robert Mostoghiu, Lister Staveley-Smith, Huiyuan Wang, Peng Wang, Gustavo Yepes

Abstract: Following Cui et al. 2018 (hereafter Paper I) on the classification of large-scale environments (LSE) at z = 0, we push our analysis to higher redshifts and study the evolution of LSE and the baryon distributions in them. Our aim is to investigate how baryons affect the LSE as a function of redshift. In agreement with Paper I, the baryon models have negligible effect on the LSE over all investigat… ▽ More Following Cui et al. 2018 (hereafter Paper I) on the classification of large-scale environments (LSE) at z = 0, we push our analysis to higher redshifts and study the evolution of LSE and the baryon distributions in them. Our aim is to investigate how baryons affect the LSE as a function of redshift. In agreement with Paper I, the baryon models have negligible effect on the LSE over all investigated redshifts. We further validate the conclusion obtained in Paper I that the gas web is an unbiased tracer of total matter -- even better at high redshifts. By separating the gas mainly by temperature, we find that about 40 per cent of gas is in the so-called warm-hot intergalactic medium (WHIM). This fraction of gas mass in the WHIM decreases with redshift, especially from z = 1 (29 per cent) to z = 2.1 (10 per cent). By separating the whole WHIM gas mass into the four large-scale environments (i.e. voids, sheets, filaments, and knots), we find that about half of the WHIM gas is located in filaments. Although the total gas mass in WHIM decreases with redshift, the WHIM mass fractions in the different LSE seem unchanged. △ Less

Submitted 12 March, 2019; v1 submitted 25 February, 2019; originally announced February 2019.

Comments: 13 pages, 8 figures, published version in MNRAS. Key results are shown in Figs. 3, 4, 8

Journal ref: Monthly Notices of the Royal Astronomical Society, Volume 485, Issue 2, 11 May 2019, Pages 2367-2379

BCG Mass Evolution in Cosmological Hydro-Simulations

Authors: C. Ragone-Figueroa, G. L. Granato, M. E. Ferraro, G. Murante, V. Biffi, S. Borgani, S. Planelles, E. Rasia

Abstract: We analyze the stellar growth of Brightest Cluster Galaxies (BCGs) produced by cosmological zoom-in hydrodynamical simulations of the formation of massive galaxy clusters. The evolution of the stellar mass content is studied considering different apertures, and tracking backwards either the main progenitor of the $z=0$ BCG or that of the cluster hosting the BCG at $z=0$. Both methods lead to simil… ▽ More We analyze the stellar growth of Brightest Cluster Galaxies (BCGs) produced by cosmological zoom-in hydrodynamical simulations of the formation of massive galaxy clusters. The evolution of the stellar mass content is studied considering different apertures, and tracking backwards either the main progenitor of the $z=0$ BCG or that of the cluster hosting the BCG at $z=0$. Both methods lead to similar results up to $z \simeq 1.5$. The simulated BCGs masses at $z=0$ are in agreement with recent observations. In the redshift interval from $z=1$ to $z=0$ we find growth factors 1.3, 1.6 and 3.6 for stellar masses within 30kpc, 50kpc and 10% of $R_{500}$ respectively. The first two factors, and in general the mass evolution in this redshift range, are in agreement with most recent observations. The last larger factor is similar to the growth factor obtained by a semi-analytical model (SAM). Half of the star particles that end up in the inner 50 kpc was typically formed by redshift $\sim$ 3.7, while the assembly of half of the BCGs stellar mass occurs on average at lower redshifts $\sim 1.5$. This assembly redshift correlates with the mass attained by the cluster at high $z \gtrsim 1.3$, due to the broader range of the progenitor clusters at high-$z$. The assembly redshift of BCGs decreases with increasing apertures. Our results are compatible with the {\it inside-out} scenario. Simulated BCGs could lack intense enough star formation (SF) at high redshift, while possibly exhibit an excess of residual SF at low redshift. △ Less

Submitted 20 June, 2018; v1 submitted 21 March, 2018; originally announced March 2018.

Comments: 13 pages, 10 figures, Accepted by MNRAS 2018 June 17

Mass-Metallicity Relation from Cosmological Hydrodynamical Simulations and X-ray Observations of Galaxy Groups and Clusters

Authors: Nhut Truong, Elena Rasia, Veronica Biffi, François Mernier, Norbert Werner, Massimo Gaspari, Stefano Borgani, Susana Planelles, Dunja Fabjan, Giuseppe Murante

Abstract: Recent X-ray observations of galaxy clusters show that the distribution of intra-cluster medium (ICM) metallicity is remarkably uniform in space and time. In this paper, we analyse a large sample of simulated objects, from poor groups to rich clusters, to study the dependence of the metallicity and related quantities on the mass of the systems. The simulations are performed with an improved versio… ▽ More Recent X-ray observations of galaxy clusters show that the distribution of intra-cluster medium (ICM) metallicity is remarkably uniform in space and time. In this paper, we analyse a large sample of simulated objects, from poor groups to rich clusters, to study the dependence of the metallicity and related quantities on the mass of the systems. The simulations are performed with an improved version of the Smoothed-Particle-Hydrodynamics \texttt{GADGET-3} code and consider various astrophysical processes including radiative cooling, metal enrichment and feedback from stars and active galactic nuclei (AGN). The scaling between the metallicity and the temperature obtained in the simulations agrees well in trend and evolution with the observational results obtained from two data samples characterised by a wide range of masses and a large redshift coverage. We find that the iron abundance in the cluster core ($r<0.1R_{500}$) does not correlate with the temperature nor presents a significant evolution. The scale invariance is confirmed when the metallicity is related directly to the total mass. The slope of the best-fitting relations is shallow ($β\sim-0.1$) in the innermost regions ($r<0.5R_{500}$) and consistent with zero outside. We investigate the impact of the AGN feedback and find that it plays a key role in producing a constant value of the outskirts metallicity from groups to clusters. This finding additionally supports the picture of early enrichment. △ Less

Submitted 15 January, 2019; v1 submitted 16 March, 2018; originally announced March 2018.

Comments: MNRAS, accepted

Multi-wavelength mock observations of the WHIM in a simulated galaxy cluster

Authors: Susana Planelles, Petar Mimica, Vicent Quilis, Carlos Cuesta-Martínez

Abstract: About half of the expected total baryon budget in the local Universe is `missing'. Hydrodynamical simulations suggest that most of the missing baryons are located in a mildly overdense, warm-hot intergalactic medium (WHIM), which is difficult to be detected at most wavelengths. In this paper we explore multi-wavelength synthetic observations of a massive galaxy cluster developed in a full Eulerian… ▽ More About half of the expected total baryon budget in the local Universe is `missing'. Hydrodynamical simulations suggest that most of the missing baryons are located in a mildly overdense, warm-hot intergalactic medium (WHIM), which is difficult to be detected at most wavelengths. In this paper we explore multi-wavelength synthetic observations of a massive galaxy cluster developed in a full Eulerian-AMR cosmological simulation. A novel numerical procedure is applied on the outputs of the simulation, which are post-processed with a full-radiative transfer code that allows to compute the change of the intensity at any frequency along the null-geodesic of photons. We compare the emission from the whole inter-galactic medium (IGM) and from the WHIM component (defined as the gas with a temperature in the range $10^5-10^7$ K) at three observational bands associated to thermal X-rays, thermal and kinematic Sunyaev-Zel'dovich effect, and radio emission. The synthetic maps produced by this procedure could be directly compared with existing observational maps and could be used as a guide for future observations with forthcoming instruments. The analysis of the different emissions associated to a high-resolution galaxy cluster is in broad agreement with previous simulated and observational estimates of both gas components. △ Less

Submitted 26 February, 2018; originally announced February 2018.

Comments: Accepted for publication in MNRAS

The origin of ICM enrichment in the outskirts of present-day galaxy clusters from cosmological hydrodynamical simulations

Authors: Veronica Biffi, Susana Planelles, Stefano Borgani, Elena Rasia, Giuseppe Murante, Dunja Fabjan, Massimo Gaspari

Abstract: The uniformity of the intra-cluster medium (ICM) enrichment level in the outskirts of nearby galaxy clusters suggests that chemical elements were deposited and widely spread into the intergalactic medium before the cluster formation. This observational evidence is supported by numerical findings from cosmological hydrodynamical simulations, as presented in Biffi et al. (2017), including the effect… ▽ More The uniformity of the intra-cluster medium (ICM) enrichment level in the outskirts of nearby galaxy clusters suggests that chemical elements were deposited and widely spread into the intergalactic medium before the cluster formation. This observational evidence is supported by numerical findings from cosmological hydrodynamical simulations, as presented in Biffi et al. (2017), including the effect of thermal feedback from active galactic nuclei. Here, we further investigate this picture, by tracing back in time the spatial origin and metallicity evolution of the gas residing at z=0 in the outskirts of simulated galaxy clusters. In these regions, we find a large distribution of iron abundances, including a component of highly-enriched gas, already present at z=2. At z>1, the gas in the present-day outskirts was distributed over tens of virial radii from the the main cluster and had been already enriched within high-redshift haloes. At z=2, about 40% of the most Fe-rich gas at z=0 was not residing in any halo more massive than 1e11 Msun/h in the region and yet its average iron abundance was already 0.4, w.r.t. the solar value by Anders & Grevesse (1989). This confirms that the in situ enrichment of the ICM in the outskirts of present-day clusters does not play a significant role, and its uniform metal abundance is rather the consequence of the accretion of both low-metallicity and pre-enriched (at z>2) gas, from the diffuse component and through merging substructures. These findings do not depend on the mass of the cluster nor on its core properties. △ Less

Submitted 12 February, 2018; v1 submitted 16 January, 2018; originally announced January 2018.

Comments: MNRAS, accepted

Shaken Snow Globes: Kinematic Tracers of the Multiphase Condensation Cascade in Massive Galaxies, Groups, and Clusters

Authors: M. Gaspari, M. McDonald, S. L. Hamer, F. Brighenti, P. Temi, M. Gendron-Marsolais, J. Hlavacek-Larrondo, A. C. Edge, N. Werner, P. Tozzi, M. Sun, J. M. Stone, G. R. Tremblay, M. T. Hogan, D. Eckert, S. Ettori, H. Yu, V. Biffi, S. Planelles

Abstract: We propose a novel method to constrain turbulence and bulk motions in massive galaxies, groups and clusters, exploring both simulations and observations. As emerged in the recent picture of the top-down multiphase condensation, the hot gaseous halos are tightly linked to all other phases in terms of cospatiality and thermodynamics. While hot halos (10^7 K) are perturbed by subsonic turbulence, war… ▽ More We propose a novel method to constrain turbulence and bulk motions in massive galaxies, groups and clusters, exploring both simulations and observations. As emerged in the recent picture of the top-down multiphase condensation, the hot gaseous halos are tightly linked to all other phases in terms of cospatiality and thermodynamics. While hot halos (10^7 K) are perturbed by subsonic turbulence, warm (10^4 K) ionized and neutral filaments condense out of the turbulent eddies. The peaks condense into cold molecular clouds (< 100 K) raining in the core via chaotic cold accretion (CCA). We show all phases are tightly linked via the ensemble (wide-aperture) velocity dispersion along the line of sight. The correlation arises in complementary long-term AGN feedback simulations and high-resolution CCA runs, and is corroborated by the combined Hitomi and new IFU measurements in Perseus cluster. The ensemble multiphase gas distributions are characterized by substantial spectral line broadening (100-200 km/s) with mild line shift. On the other hand, pencil-beam detections sample the small-scale clouds displaying smaller broadening and significant line shift up to several 100 km/s, with increased scatter due to the turbulence intermittency. We present new ensemble sigma_v of the warm Halpha+[NII] gas in 72 observed cluster/group cores: the constraints are consistent with the simulations and can be used as robust proxies for the turbulent velocities, in particular for the challenging hot plasma (otherwise requiring extremely long X-ray exposures). We show the physically motivated criterion C = t_cool/t_eddy ~ 1 best traces the condensation extent region and presence of multiphase gas in observed clusters/groups. The ensemble method can be applied to many available datasets and can substantially advance our understanding of multiphase halos in light of the next-generation multiwavelength missions. △ Less

Submitted 9 February, 2018; v1 submitted 19 September, 2017; originally announced September 2017.

Comments: 18 pages, 5 figures, 2 tables; ApJ in press - typos corrected, formatting improved (new emulateapj mod included here)

Is ram-pressure stripping an efficient mechanism to remove gas in galaxies?

Authors: Vicent Quilis, Susana Planelles, Elena Ricciardelli

Abstract: We study how the gas in a sample of galaxies (M* > 10e9 Msun) in clusters, obtained in a cosmological simulation, is affected by the interaction with the intra-cluster medium (ICM). The dynamical state of each elemental parcel of gas is studied using the total energy. At z ~ 2, the galaxies in the simulation are evenly distributed within clusters, moving later on towards more central locations. In… ▽ More We study how the gas in a sample of galaxies (M* > 10e9 Msun) in clusters, obtained in a cosmological simulation, is affected by the interaction with the intra-cluster medium (ICM). The dynamical state of each elemental parcel of gas is studied using the total energy. At z ~ 2, the galaxies in the simulation are evenly distributed within clusters, moving later on towards more central locations. In this process, gas from the ICM is accreted and mixed with the gas in the galactic halo. Simultaneously, the interaction with the environment removes part of the gas. A characteristic stellar mass around M* ~ 10e10 Msun appears as a threshold marking two differentiated behaviours. Below this mass, galaxies are located at the external part of clusters and have eccentric orbits. The effect of the interaction with the environment is marginal. Above, galaxies are mainly located at the inner part of clusters with mostly radial orbits with low velocities. In these massive systems, part of the gas, strongly correlated with the stellar mass of the galaxy, is removed. The amount of removed gas is sub-dominant compared with the quantity of retained gas which is continuously influenced by the hot gas coming from the ICM. The analysis of individual galaxies reveals the existence of a complex pattern of flows, turbulence and a constant fuelling of gas to the hot corona from the ICM that could make the global effect of the interaction of galaxies with their environment to be substantially less dramatic than previously expected. △ Less

Submitted 28 March, 2017; originally announced March 2017.

Comments: 17 pages, 12 figures, accepted for publication in MNRAS

The history of chemical enrichment in the intracluster medium from cosmological simulations

Authors: V. Biffi, S. Planelles, S. Borgani, D. Fabjan, E. Rasia, G. Murante, L. Tornatore, K. Dolag, G. L. Granato, M. Gaspari, A. M. Beck

Abstract: The distribution of metals in the intracluster medium (ICM) of galaxy clusters provides valuable information on their formation and evolution, on the connection with the cosmic star formation and on the effects of different gas processes. By analyzing a sample of simulated galaxy clusters, we study the chemical enrichment of the ICM, its evolution, and its relation with the physical processes incl… ▽ More The distribution of metals in the intracluster medium (ICM) of galaxy clusters provides valuable information on their formation and evolution, on the connection with the cosmic star formation and on the effects of different gas processes. By analyzing a sample of simulated galaxy clusters, we study the chemical enrichment of the ICM, its evolution, and its relation with the physical processes included in the simulation and with the thermal properties of the core. These simulations, consisting of re-simulations of 29 Lagrangian regions performed with an upgraded version of the SPH GADGET-3 code, have been run including two different sets of baryonic physics: one accounts for radiative cooling, star formation, metal enrichment and supernova (SN) feedback, and the other one further includes the effects of feedback from active galactic nuclei (AGN). In agreement with observations, we find an anti-correlation between entropy and metallicity in cluster cores, and similar radial distributions of heavy-element abundances and abundance ratios out to large cluster-centric distances (~R180). In the outskirts, namely outside of ~0.2R180, we find a remarkably homogeneous metallicity distribution, with almost flat profiles of the elements produced by either SNIa or SNII. We investigated the origin of this phenomenon and discovered that it is due to the widespread displacement of metal-rich gas by early (z>2-3) AGN powerful bursts, acting on small high-redshift haloes. Our results also indicate that the intrinsic metallicity of the hot gas for this sample is on average consistent with no evolution between z=2 and z=0, across the entire radial range. △ Less

Submitted 27 January, 2017; originally announced January 2017.

Comments: 18 pages; submitted to MNRAS

Pressure of the hot gas in simulations of galaxy clusters

Authors: S. Planelles, D. Fabjan, S. Borgani, G. Murante, E. Rasia, V. Biffi, N. Truong, C. Ragone-Figueroa, G. L. Granato, K. Dolag, E. Pierpaoli, A. M. Beck, Lisa K. Steinborn, M. Gaspari

Abstract: We analyze the radial pressure profiles, the ICM clumping factor and the Sunyaev-Zel'dovich (SZ) scaling relations of a sample of simulated galaxy clusters and groups identified in a set of hydrodynamical simulations based on an updated version of the TreePM-SPH GADGET-3 code. Three different sets of simulations are performed: the first assumes non-radiative physics, the others include, among othe… ▽ More We analyze the radial pressure profiles, the ICM clumping factor and the Sunyaev-Zel'dovich (SZ) scaling relations of a sample of simulated galaxy clusters and groups identified in a set of hydrodynamical simulations based on an updated version of the TreePM-SPH GADGET-3 code. Three different sets of simulations are performed: the first assumes non-radiative physics, the others include, among other processes, AGN and/or stellar feedback. Our results are analyzed as a function of redshift, ICM physics, cluster mass and cluster cool-coreness or dynamical state. In general, the mean pressure profiles obtained for our sample of groups and clusters show a good agreement with X-ray and SZ observations. Simulated cool-core (CC) and non-cool-core (NCC) clusters also show a good match with real data. We obtain in all cases a small (if any) redshift evolution of the pressure profiles of massive clusters, at least back to z=1. We find that the clumpiness of gas density and pressure increases with the distance from the cluster center and with the dynamical activity. The inclusion of AGN feedback in our simulations generates values for the gas clumping ($\sqrt C_ρ\sim 1.2$ at $R_{200}$) in good agreement with recent observational estimates. The simulated $Y_{SZ}-M$ scaling relations are in good accordance with several observed samples, especially for massive clusters. As for the scatter of these relations, we obtain a clear dependence on the cluster dynamical state, whereas this distinction is not so evident when looking at the subsamples of CC and NCC clusters. △ Less

Submitted 21 December, 2016; originally announced December 2016.

Comments: 22 pages, 14 figures, 2 tables. Submitted to MNRAS

Cosmological hydrodynamical simulations of galaxy clusters: X-ray scaling relations and their evolution

Authors: N. Truong, E. Rasia, P. Mazzotta, S. Planelles, V. Biffi, D. Fabjan, A. M. Beck, S. Borgani, K. Dolag, M. Gaspari, G. L. Granato, G. Murante, C. Ragone-Figueroa, L. K. Steinborn

Abstract: We analyse cosmological hydrodynamical simulations of galaxy clusters to study the X-ray scaling relations between total masses and observable quantities such as X-ray luminosity, gas mass, X-ray temperature, and $Y_{X}$. Three sets of simulations are performed with an improved version of the smoothed particle hydrodynamics GADGET-3 code. These consider the following: non-radiative gas, star forma… ▽ More We analyse cosmological hydrodynamical simulations of galaxy clusters to study the X-ray scaling relations between total masses and observable quantities such as X-ray luminosity, gas mass, X-ray temperature, and $Y_{X}$. Three sets of simulations are performed with an improved version of the smoothed particle hydrodynamics GADGET-3 code. These consider the following: non-radiative gas, star formation and stellar feedback, and the addition of feedback by active galactic nuclei (AGN). We select clusters with $M_{500} > 10^{14} M_{\odot} E(z)^{-1}$, mimicking the typical selection of Sunyaev-Zeldovich samples. This permits to have a mass range large enough to enable robust fitting of the relations even at $z \sim 2$. The results of the analysis show a general agreement with observations. The values of the slope of the mass-gas mass and mass-temperature relations at $z=2$ are 10 per cent lower with respect to $z=0$ due to the applied mass selection, in the former case, and to the effect of early merger in the latter. We investigate the impact of the slope variation on the study of the evolution of the normalization. We conclude that cosmological studies through scaling relations should be limited to the redshift range $z=0-1$, where we find that the slope, the scatter, and the covariance matrix of the relations are stable. The scaling between mass and $Y_X$ is confirmed to be the most robust relation, being almost independent of the gas physics. At higher redshifts, the scaling relations are sensitive to the inclusion of AGNs which influences low-mass systems. The detailed study of these objects will be crucial to evaluate the AGN effect on the ICM. △ Less

Submitted 25 January, 2018; v1 submitted 30 June, 2016; originally announced July 2016.

Comments: 24 pages, 11 figures, 5 tables, replaced to match accepted version

Journal ref: 2018, MNRAS, 474, 4089T

On the nature of hydrostatic equilibrium in galaxy clusters

Authors: V. Biffi, S. Borgani, G. Murante, E. Rasia, S. Planelles, G. L. Granato, C. Ragone-Figueroa, A. M. Beck, M. Gaspari, K. Dolag

Abstract: In this paper we investigate the level of hydrostatic equilibrium (HE) in the intra-cluster medium of simulated galaxy clusters, extracted from state-of-the-art cosmological hydrodynamical simulations performed with the Smoothed-Particle-Hydrodynamic code GADGET-3. These simulations include several physical processes, among which stellar and AGN feedback, and have been performed with an improved v… ▽ More In this paper we investigate the level of hydrostatic equilibrium (HE) in the intra-cluster medium of simulated galaxy clusters, extracted from state-of-the-art cosmological hydrodynamical simulations performed with the Smoothed-Particle-Hydrodynamic code GADGET-3. These simulations include several physical processes, among which stellar and AGN feedback, and have been performed with an improved version of the code that allows for a better description of hydrodynamical instabilities and gas mixing processes. Evaluating the radial balance between the gravitational and hydrodynamical forces, via the gas accelerations generated, we effectively examine the level of HE in every object of the sample, its dependence on the radial distance from the center and on the classification of the cluster in terms of either cool-coreness or dynamical state. We find an average deviation of 10-20% out to the virial radius, with no evident distinction between cool-core and non-cool-core clusters. Instead, we observe a clear separation between regular and disturbed systems, with a more significant deviation from HE for the disturbed objects. The investigation of the bias between the hydrostatic estimate and the total gravitating mass indicates that, on average, this traces very well the deviation from HE, even though individual cases show a more complex picture. Typically, in the radial ranges where mass bias and deviation from HE are substantially different, the gas is characterized by a significant amount of random motions (>~30 per cent), relative to thermal ones. As a general result, the HE-deviation and mass bias, at given interesting distance from the cluster center, are not very sensitive to the temperature inhomogeneities in the gas. △ Less

Submitted 7 June, 2016; originally announced June 2016.

Comments: 14 pages, 13 figures, 1 table; submitted

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

Authors: Francisco Villaescusa-Navarro, Susana Planelles, Stefano Borgani, Matteo Viel, Elena Rasia, Giuseppe Murante, Klaus Dolag, Lisa K. Steinborn, Veronica Biffi, Alexander M. Beck, Cinthia Ragone-Figueroa

Abstract: 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 a… ▽ More 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. △ Less

Submitted 8 December, 2015; v1 submitted 14 October, 2015; originally announced October 2015.

Comments: 19 pages, 12 figures. Accepted for publication in MNRAS

Journal ref: Monthly Notices of the Royal Astronomical Society 2016 456 (4): 3553-3570

Cool Core Clusters from Cosmological Simulations

Authors: E. Rasia, S. Borgani, G. Murante, S. Planelles, A. M. Beck, V. Biffi, C. Ragone-Figueroa, G. L. Granato, L. K. Steinborn, K. Dolag

Abstract: We present results obtained from a set of cosmological hydrodynamic simulations of galaxy clusters, aimed at comparing predictions with observational data on the diversity between cool-core (CC) and non-cool-core (NCC) clusters. Our simulations include the effects of stellar and AGN feedback and are based on an improved version of the smoothed particle hydrodynamics code GADGET-3, which ameliorate… ▽ More We present results obtained from a set of cosmological hydrodynamic simulations of galaxy clusters, aimed at comparing predictions with observational data on the diversity between cool-core (CC) and non-cool-core (NCC) clusters. Our simulations include the effects of stellar and AGN feedback and are based on an improved version of the smoothed particle hydrodynamics code GADGET-3, which ameliorates gas mixing and better captures gas-dynamical instabilities by including a suitable artificial thermal diffusion. In this Letter, we focus our analysis on the entropy profiles, the primary diagnostic we used to classify the degree of cool-coreness of clusters, and on the iron profiles. In keeping with observations, our simulated clusters display a variety of behaviors in entropy profiles: they range from steadily decreasing profiles at small radii, characteristic of cool-core systems, to nearly flat core isentropic profiles, characteristic of non-cool-core systems. Using observational criteria to distinguish between the two classes of objects, we find that they occur in similar proportions in both simulations and in observations. Furthermore, we also find that simulated cool-core clusters have profiles of iron abundance that are steeper than those of NCC clusters, which is also in agreement with observational results. We show that the capability of our simulations to generate a realistic cool-core structure in the cluster population is due to AGN feedback and artificial thermal diffusion: their combined action allows us to naturally distribute the energy extracted from super-massive black holes and to compensate for the radiative losses of low-entropy gas with short cooling time residing in the cluster core. △ Less

Submitted 31 October, 2015; v1 submitted 14 September, 2015; originally announced September 2015.

Comments: 6 pages, 4 figures, accepted in ApJL, v2 contains some modifications on the text (results unchanged)

Simulation-based marginal likelihood for cluster strong lensing cosmology

Authors: Madhura Killedar, Stefano Borgani, Dunja Fabjan, Klaus Dolag, Gian Luigi Granato, Massimo Meneghetti, Susana Planelles, Cinthia Ragone-Figueroa

Abstract: 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 th… ▽ More 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. △ Less

Submitted 1 September, 2017; v1 submitted 20 July, 2015; originally announced July 2015.

Comments: 15 pages, 6 figures, 1 table, accepted to MNRAS

The Spanish Square Kilometre Array White Book

Authors: J. A. Acosta-Pulido, I. Agudo, A. Alberdi, J. Alcolea, E. J. Alfaro, A. Alonso-Herrero, G. Anglada, P. Arnalte-Mur, Y. Ascasibar, B. Ascaso, R. Azulay, R. Bachiller, A. Baez-Rubio, E. Battaner, J. Blasco, C. B. Brook, V. Bujarrabal, G. Busquet, M. D. Caballero-Garcia, C. Carrasco-Gonzalez, J. Casares, A. J. Castro-Tirado, L. Colina, F. Colomer, I. de Gregorio-Monsalvo , et al. (94 additional authors not shown)

Abstract: The Square Kilometre Array (SKA) is called to revolutionise essentially all areas of Astrophysics. With a collecting area of about a square kilometre, the SKA will be a transformational instrument, and its scientific potential will go beyond the interests of astronomers. Its technological challenges and huge cost requires a multinational effort, and Europe has recognised this by putting the SKA on… ▽ More The Square Kilometre Array (SKA) is called to revolutionise essentially all areas of Astrophysics. With a collecting area of about a square kilometre, the SKA will be a transformational instrument, and its scientific potential will go beyond the interests of astronomers. Its technological challenges and huge cost requires a multinational effort, and Europe has recognised this by putting the SKA on the roadmap of the European Strategy Forum for Research Infrastructures (ESFRI). The Spanish SKA White Book is the result of the coordinated effort of 120 astronomers from 40 different research centers. The book shows the enormous scientific interest of the Spanish astronomical community in the SKA and warrants an optimum scientific exploitation of the SKA by Spanish researchers, if Spain enters the SKA project. △ Less

Submitted 17 June, 2015; v1 submitted 10 June, 2015; originally announced June 2015.

Comments: Editors: Miguel A. Perez-Torres (Editor-in-chief), L. Verdes-Montenegro, J. C. Guirado, A. Alberdi, J. Martin-Pintado, R. Bachiller, D. Herranz, J. M. Girart, J. Gorgas, C. Hernandez-Monteagudo, S. Migliari and J. M. Rodriguez Espinosa; 289 pages, published in June 2015 by Sociedad Española de Astronomía

An improved SPH scheme for cosmological simulations

Authors: A. M. Beck, G. Murante, A. Arth, R. -S. Remus, A. F. Teklu, J. M. F. Donnert, S. Planelles, M. C. Beck, P. Foerster, M. Imgrund, K. Dolag, S. Borgani

Abstract: We present an implementation of smoothed particle hydrodynamics (SPH) with improved accuracy for simulations of galaxies and the large-scale structure. In particular, we combine, implement, modify and test a vast majority of SPH improvement techniques in the latest instalment of the GADGET code. We use the Wendland kernel functions, a particle wake-up time-step limiting mechanism and a time-depend… ▽ More We present an implementation of smoothed particle hydrodynamics (SPH) with improved accuracy for simulations of galaxies and the large-scale structure. In particular, we combine, implement, modify and test a vast majority of SPH improvement techniques in the latest instalment of the GADGET code. We use the Wendland kernel functions, a particle wake-up time-step limiting mechanism and a time-dependent scheme for artificial viscosity, which includes a high-order gradient computation and shear flow limiter. Additionally, we include a novel prescription for time-dependent artificial conduction, which corrects for gravitationally induced pressure gradients and largely improves the SPH performance in capturing the development of gas-dynamical instabilities. We extensively test our new implementation in a wide range of hydrodynamical standard tests including weak and strong shocks as well as shear flows, turbulent spectra, gas mixing, hydrostatic equilibria and self-gravitating gas clouds. We jointly employ all modifications; however, when necessary we study the performance of individual code modules. We approximate hydrodynamical states more accurately and with significantly less noise than standard SPH. Furthermore, the new implementation promotes the mixing of entropy between different fluid phases, also within cosmological simulations. Finally, we study the performance of the hydrodynamical solver in the context of radiative galaxy formation and non-radiative galaxy cluster formation. We find galactic disks to be colder, thinner and more extended and our results on galaxy clusters show entropy cores instead of steadily declining entropy profiles. In summary, we demonstrate that our improved SPH implementation overcomes most of the undesirable limitations of standard SPH, thus becoming the core of an efficient code for large cosmological simulations. △ Less

Submitted 20 October, 2015; v1 submitted 25 February, 2015; originally announced February 2015.

Comments: 21 figures, 2 tables, accepted to MNRAS

Large-Scale Structure Formation: from the first non-linear objects to massive galaxy clusters

Authors: S. Planelles, D. R. G. Schleicher, A. M. Bykov

Abstract: The large-scale structure of the Universe formed from initially small perturbations in the cosmic density field, leading to galaxy clusters with up to 10^15 Msun at the present day. Here, we review the formation of structures in the Universe, considering the first primordial galaxies and the most massive galaxy clusters as extreme cases of structure formation where fundamental processes such as gr… ▽ More The large-scale structure of the Universe formed from initially small perturbations in the cosmic density field, leading to galaxy clusters with up to 10^15 Msun at the present day. Here, we review the formation of structures in the Universe, considering the first primordial galaxies and the most massive galaxy clusters as extreme cases of structure formation where fundamental processes such as gravity, turbulence, cooling and feedback are particularly relevant. The first non-linear objects in the Universe formed in dark matter halos with 10^5-10^8 Msun at redshifts 10-30, leading to the first stars and massive black holes. At later stages, larger scales became non-linear, leading to the formation of galaxy clusters, the most massive objects in the Universe. We describe here their formation via gravitational processes, including the self-similar scaling relations, as well as the observed deviations from such self-similarity and the related non-gravitational physics (cooling, stellar feedback, AGN). While on intermediate cluster scales the self-similar model is in good agreement with the observations, deviations from such self-similarity are apparent in the core regions, where numerical simulations do not reproduce the current observational results. The latter indicates that the interaction of different feedback processes may not be correctly accounted for in current simulations. Both in the most massive clusters of galaxies as well as during the formation of the first objects in the Universe, turbulent structures and shock waves appear to be common, suggesting them to be ubiquitous in the non-linear regime. △ Less

Submitted 15 April, 2014; originally announced April 2014.

Comments: Review article. Accepted for publication in Space Science Reviews. It will appear as a contribution to an ISSI book

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

Authors: Kai Hoffmann, Susana Planelles, Enrique Gaztanaga, Alexander Knebe, Frazer R. Pearce, Hanni Lux, Julian Onions, Stuart I. Muldrew, Pascal Elahi, Peter Behroozi, Yago Ascasibar, Jiaxin Han, Michal Maciejewski, Manuel E. Merchan, Mark Neyrinck, Andrés N. Ruiz, Mario A. Sgro

Abstract: 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 dist… ▽ More 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. △ Less

Submitted 10 September, 2014; v1 submitted 9 January, 2014; originally announced January 2014.

Comments: 15 pages, 11 figures, 2 tables. 1 figure removed to reduce paper extension, shorter discussion, accepted for publication in MNRAS

On the role of AGN feedback on the thermal and chemodynamical properties of the hot intra-cluster medium

Authors: S. Planelles, S. Borgani, D. Fabjan, M. Killedar, G. Murante, G. L. Granato, C. Ragone-Figueroa, K. Dolag

Abstract: We present an analysis of the properties of the ICM in an extended set of cosmological hydrodynamical simulations of galaxy clusters and groups performed with the TreePM+SPH GADGET-3 code. Besides a set of non-radiative simulations, we carried out two sets of simulations including radiative cooling, star formation, metal enrichment and feedback from supernovae, one of which also accounts for the e… ▽ More We present an analysis of the properties of the ICM in an extended set of cosmological hydrodynamical simulations of galaxy clusters and groups performed with the TreePM+SPH GADGET-3 code. Besides a set of non-radiative simulations, we carried out two sets of simulations including radiative cooling, star formation, metal enrichment and feedback from supernovae, one of which also accounts for the effect of feedback from AGN resulting from gas accretion onto super-massive black holes. These simulations are analysed with the aim of studying the relative role played by SN and AGN feedback on the general properties of the diffuse hot baryons in galaxy clusters and groups: scaling relations, temperature, entropy and pressure radial profiles, and ICM chemical enrichment. We find that simulations including AGN feedback produce scaling relations that are in good agreement with X-ray observations at all mass scales. However, our simulations are not able to account for the observed diversity between CC and NCC clusters: unlike for observations, we find that temperature and entropy profiles of relaxed and unrelaxed clusters are quite similar and resemble more the observed behaviour of NCC clusters. As for the pattern of metal enrichment, we find that an enhanced level of iron abundance is produced by AGN feedback with respect to the case of purely SN feedback. As a result, while simulations including AGN produce values of iron abundance in groups in agreement with observations, they over-enrich the ICM in massive clusters. The efficiency of AGN feedback in displacing enriched gas from halos into the inter-galactic medium at high redshift also creates a widespread enrichment in the outskirts of clusters and produces profiles of iron abundance whose slope is in better agreement with observations. △ Less

Submitted 4 November, 2013; originally announced November 2013.

Comments: 23 pages, 14 figures, 1 table, accepted for publication in MNRAS

The structure of cosmic voids in a LCDM Universe

Authors: E. Ricciardelli, V. Quilis, S. Planelles

Abstract: Eulerian cosmological codes are especially suited to properly describe the low density regions. This property makes this class of codes excellent tools to study the formation and evolution of cosmic voids. Following such ideas, we present the results of an Eulerian adaptive mesh refinement (AMR) hydrodynamical and N-body simulation, that contrary to the common practice, has been designed to refine… ▽ More Eulerian cosmological codes are especially suited to properly describe the low density regions. This property makes this class of codes excellent tools to study the formation and evolution of cosmic voids. Following such ideas, we present the results of an Eulerian adaptive mesh refinement (AMR) hydrodynamical and N-body simulation, that contrary to the common practice, has been designed to refine the computational grid in the underdense regions of the simulated volume. Thus, the void regions are better described due to the combined effect of the Eulerian character of the numerical technique and the use of high numerical resolution from the AMR approach. To analyse the outcome of this simulation, we have constructed a new void finder optimally suited to find the hierarchy of voids in AMR simulations. The algorithm identifies voids starting from the cells with least density and highest velocity divergence and then expanding the underdense volume until reaching the void walls, defined from the steepness of the density gradient. At redshift z=0, in a cosmological box of comoving side length 100 Mpc/h, we identify hundreds of voids with sizes up to 17 Mpc/h and typical density contrast of -0.8, which show a complex morphology and an intricate hierarchy of nested structures. The analysis of their mass density profile leads to the conclusion that a universal density profile can be applied to voids of any size, density, morphology and redshift. △ Less

Submitted 12 June, 2013; originally announced June 2013.

Comments: 14 pages, 13 figures. Accepted for publication in MNRAS

Structure Finding in Cosmological Simulations: The State of Affairs

Authors: Alexander Knebe, Frazer R. Pearce, Hanni Lux, Yago Ascasibar, Peter Behroozi, Javier Casado, Christine Corbett Moran, Juerg Diemand, Klaus Dolag, Rosa Dominguez-Tenreiro, Pascal Elahi, Bridget Falck, Stefan Gottloeber, Jiaxin Han, Anatoly Klypin, Zarija Lukic, Michal Maciejewski, Cameron K. McBride, Manuel E. Merchan, Stuart I. Muldrew, Mark Neyrinck, Julian Onions, Susana Planelles, Doug Potter, Vicent Quilis , et al. (10 additional authors not shown)

Abstract: 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 remain… ▽ More 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. △ Less

Submitted 26 July, 2013; v1 submitted 2 April, 2013; originally announced April 2013.

Comments: 28 pages containing 13 figures & 4 tables + 9 pages appendix containing another 4 tables + 4 pages of references, accepted for publication in MNRAS

Cosmological shock waves: clues to the formation history of haloes

Authors: Susana Planelles, Vicent Quilis

Abstract: 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… ▽ More 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) △ Less

Submitted 4 October, 2012; originally announced October 2012.

Comments: 14 pages, 9 figures, accepted for publication in MNRAS

Baryon Census in Hydrodynamical Simulations of Galaxy Clusters

Authors: Susana Planelles, Stefano Borgani, Klaus Dolag, Stefano Ettori, Dunja Fabjan, Giuseppe Murante, Luca Tornatore

Abstract: 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 coolin… ▽ More 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. △ Less

Submitted 11 February, 2013; v1 submitted 23 September, 2012; originally announced September 2012.

Comments: 17 pages, 8 figures, 4 tables. Additional figure, minor modifications and references added to match the accepted version

Haloes gone MAD: The Halo-Finder Comparison Project

Authors: Alexander Knebe, Steffen R. Knollmann, Stuart I. Muldrew, Frazer R. Pearce, Miguel Angel Aragon-Calvo, Yago Ascasibar, Peter S. Behroozi, Daniel Ceverino, Stephane Colombi, Juerg Diemand, Klaus Dolag, Bridget L. Falck, Patricia Fasel, Jeff Gardner, Stefan Gottloeber, Chung-Hsing Hsu, Francesca Iannuzzi, Anatoly Klypin, Zarija Lukic, Michal Maciejewski, Cameron McBride, Mark C. Neyrinck, Susana Planelles, Doug Potter, Vicent Quilis , et al. (12 additional authors not shown)

Abstract: [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 develope… ▽ More [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). △ Less

Submitted 5 April, 2011; originally announced April 2011.

Comments: 27 interesting pages, 20 beautiful figures, and 4 informative tables accepted for publication in MNRAS. The high-resolution version of the paper as well as all the test cases and analysis can be found at the web site http://popia.ft.uam.es/HaloesGoingMAD

ASOHF: a new adaptive spherical overdensity halo finder

Authors: Susana Planelles, Vicent Quilis

Abstract: 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 sta… ▽ More 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. △ Less

Submitted 16 June, 2010; originally announced June 2010.

Comments: Latex format, 16 pages, 11 postscript figures, accepted for publication in Astronomy and Astrophysics

Galaxy cluster mergers

Authors: Susana Planelles, Vicent Quilis

Abstract: 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.… ▽ More 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. △ Less

Submitted 22 June, 2009; originally announced June 2009.

Comments: 17 pages, 11 figures, accepted for publication in MNRAS