Showing 1–38 of 38 results for author: Ragone-Figueroa, C

Intertwined Formation of $\rm{H_2}$, Dust, and Stars in Cosmological Simulations

Authors: Cinthia Ragone-Figueroa, Gian Luigi Granato, Massimiliano Parente, Giuseppe Murante, Milena Valentini, Stefano Borgani, Umberto Maio

Abstract: Context: Molecular hydrogen ($\rm{H_2}$) is crucial in galaxy formation and evolution, serving as the main fuel for star formation (SF). In metal-enriched environments, $\rm{H_2}$ primarily forms on interstellar dust grain surfaces. However, due to the complexities of modelling this process, SF in cosmological simulations often relies on empirical or theoretical frameworks validated only in the Lo… ▽ More Context: Molecular hydrogen ($\rm{H_2}$) is crucial in galaxy formation and evolution, serving as the main fuel for star formation (SF). In metal-enriched environments, $\rm{H_2}$ primarily forms on interstellar dust grain surfaces. However, due to the complexities of modelling this process, SF in cosmological simulations often relies on empirical or theoretical frameworks validated only in the Local Universe to estimate the abundance of $\rm{H_2}$. Aims: This study aims to model the connection between star, dust, and $\rm{H_2}$ formation processes in cosmological simulations. Methods: We include $\rm{H_2}$ formation on dust grain surfaces and account for molecule destruction and radiation shielding into the SF and feedback model MUPPI. Results: The model reproduces key properties of observed galaxies for stellar, dust, and $\rm{H_2}$ components. The cosmic density of $\rm{H_2}$ ($ρ_{\rm{H2}}$) peaks around $z=1.5$, then decreases by half towards $z=0$, showing milder evolution than observed. The $\rm{H_2}$ mass function since $z=2$ also shows gentler evolution. Our model successfully recovers the integrated molecular Kennicutt-Schmidt (mKS) law between surface star formation rate ($Σ_{\rm SFR}$) and surface $\rm{H_2}$ density ($Σ_{\rm H2}$) at $z=0$, already evident at $z=2$ with a higher normalization. We find hints of a broken power law with a steeper slope at higher $Σ_{\rm H2}$, aligning with some observational findings. Additionally, the $\rm{H_2}$-to-dust mass ratio in galaxies shows a decreasing trend with gas metallicity and stellar mass. The $\rm{H_2}$-to-dust mass fraction for the global galaxy population is higher at higher redshifts. The analysis of the atomic-to-molecular transition on a particle-by-particle basis suggests that gas metallicity cannot reliably substitute the dust-to-gas ratio in models simulating dust-promoted $\rm{H_2}$. △ Less

Submitted 8 July, 2024; originally announced July 2024.

Comments: submitted to A&A. Comments are welcome

Dynamical friction and evolution of black holes in cosmological simulations: a new implementation in OpenGadget3

Authors: Alice Damiano, Milena Valentini, Stefano Borgani, Luca Tornatore, Giuseppe Murante, Antonio Ragagnin, Cinthia Ragone-Figueroa, Klaus Dolag

Abstract: We implement a sub-resolution prescription for the unresolved dynamical friction onto black holes (BHs) in the OpenGadget3 code. We carry out cosmological simulations of a volume of 16 cMpc3 and zoom-ins of a galaxy group and of a galaxy cluster. The advantages of our new technique are assessed in comparison to commonly adopted methods to hamper spurious BH displacements, i.e. repositioning onto a… ▽ More We implement a sub-resolution prescription for the unresolved dynamical friction onto black holes (BHs) in the OpenGadget3 code. We carry out cosmological simulations of a volume of 16 cMpc3 and zoom-ins of a galaxy group and of a galaxy cluster. The advantages of our new technique are assessed in comparison to commonly adopted methods to hamper spurious BH displacements, i.e. repositioning onto a local minimum of the gravitational potential and ad-hoc boosting of the BH particle dynamical mass. The newly-introduced dynamical friction correction provides centering of BHs on host halos which is at least comparable with the other techniques. It predicts half as many merger events with respect to the repositioning prescription, with the advantage of being less prone to leave sub-structures without any central BH. Simulations featuring our dynamical friction prescription produce a smaller (by up to 50% with respect to repositioning) population of wandering BHs and final BH masses in good agreement with observations. As for individual BH-BH interactions, our dynamical friction model captures the gradual inspiraling of orbits before the merger occurs. By contrast, the repositioning scheme, in its most classical renditions considered, describes extremely fast mergers, while the dynamical mass misrepresents the BHs' dynamics, introducing numerical scattering between the orbiting BHs. Given its performances in describing the centering of BHs within host galaxies and the orbiting of BH pair before their merging, our dynamical friction correction opens interesting applications for an accurate description of the evolution of BH demography within cosmological simulations of galaxy formation at different cosmic epochs and within different environments. △ Less

Submitted 19 March, 2024; originally announced March 2024.

Comments: 19 pages, 10 figures, submitted to A&A

Star Formation and Dust in the Cosmic Web

Authors: Massimiliano Parente, Cinthia Ragone-Figueroa, Pablo López, Héctor J. Martínez, Andrés N. Ruiz, Laura Ceccarelli, Valeria Coenda, Facundo Rodriguez, Gian Luigi Granato, Andrea Lapi, Rien van de Weygaert

Abstract: The large-scale environment of the cosmic web is believed to impact galaxy evolution, but there is still no consensus regarding the mechanisms. We use a semi-analytic model (SAM) galaxy catalog to study the star formation and dust content of local galaxies in different cosmic environments of the cosmic web, namely voids, filaments, walls, and nodes. We find a strong impact of the environment only… ▽ More The large-scale environment of the cosmic web is believed to impact galaxy evolution, but there is still no consensus regarding the mechanisms. We use a semi-analytic model (SAM) galaxy catalog to study the star formation and dust content of local galaxies in different cosmic environments of the cosmic web, namely voids, filaments, walls, and nodes. We find a strong impact of the environment only for galaxies with $M_{\rm stars}\lesssim10^{10.8}\, M_\odot$: the less dense the environment, the larger the star formation rate and dust content at fixed stellar mass. This is attributed to the fact that galaxies in less dense environments typically feature younger stellar populations, a slower evolution of their stellar mass and a delayed star formation compared to galaxies in denser environments. As for galaxies with $M_{\rm stars}\gtrsim 10^{10.8}\, M_\odot$ differences among environments are milder due to the disc instability (DI) driven supermassive black hole (SMBH) growth implemented in the SAM, which makes SMBH growth, and thus galaxy quenching, environment insensitive. We qualitatively test our predictions against observations by identifying environments in the SDSS-DR16 using dust masses derived from the GAMA survey. The agreement is encouraging, particularly at ${\rm log} \, M_{\rm stars}/M_\odot\gtrsim 10.5-11$, where sSFRs and dust masses appear quite environment-insensitive. This result confirms the importance of in situ growth channels of SMBHs. △ Less

Submitted 11 March, 2024; v1 submitted 15 December, 2023; originally announced December 2023.

Comments: 21 pages, 12 figures, 3 tables. Accepted for publication in ApJ

The $z \lesssim 1$ drop of cosmic dust abundance in a semi-analytic framework

Authors: Massimiliano Parente, Cinthia Ragone-Figueroa, Gian Luigi Granato, Andrea Lapi

Abstract: Observations suggest that the amount of galactic dust in the Universe decreased by a factor $\sim 2-3$ during the last $\sim 8$ Gyr. However, cosmological models of galaxy evolution usually struggle to explain this decrease. Here we use the semi-analytic model (SAM) L-Galaxies2020 to show that this drop may be reproduced assuming standard prescriptions for dust production and evolution. We extend… ▽ More Observations suggest that the amount of galactic dust in the Universe decreased by a factor $\sim 2-3$ during the last $\sim 8$ Gyr. However, cosmological models of galaxy evolution usually struggle to explain this decrease. Here we use the semi-analytic model (SAM) L-Galaxies2020 to show that this drop may be reproduced assuming standard prescriptions for dust production and evolution. We extend the SAM with i) a state-of-the-art dust model which adopts the two-size approximation and ii) a new disc instability criterion which triggers bulge and central black hole growth. The model reproduces some fundamental properties of the local galaxy population, such as the fraction of spheroid-dominated galaxies and some scaling relations involving dust. Moreover, the model predicts a galactic dust drop from $z \sim 1 \rightarrow 0$, which becomes closer to the observed one when adopting the new treatment of disc instabilities. This result is related to the newly implemented super-massive black hole growth during disc instabilities, which enhances the quenching of massive galaxies. Consequently, these objects feature a lower gas and dust content. We provide a census of the contribution of all the processes affecting the galactic dust content. Accretion is the dominant dust mass growth process. Destruction by supernovae, astration and ejection by winds have all a non-negligible role in decreasing the overall dust content in galaxies below $z \sim 1$. We also discuss predictions concerning extra-galactic dust, confirming that a sputtering efficiency lower than the canonical one is required to match the few available observations. △ Less

Submitted 23 March, 2023; v1 submitted 6 February, 2023; originally announced February 2023.

Comments: 20 pages, 20 figures, 1 table. Accepted for publication in MNRAS

Impact of H$_{\rm 2}$-driven star formation and stellar feedback from low-enrichment environments on the formation of spiral galaxies

Authors: Milena Valentini, Klaus Dolag, Stefano Borgani, Giuseppe Murante, Umberto Maio, Luca Tornatore, Gian Luigi Granato, Cinthia Ragone-Figueroa, Andreas Burkert, Antonio Ragagnin, Elena Rasia

Abstract: The reservoir of molecular gas (H$_{\rm 2}$) represents the fuel for the star formation (SF) of a galaxy. Connecting the star formation rate (SFR) to the available H$_{\rm 2}$ is key to accurately model SF in cosmological simulations of galaxy formation. We investigate how modifying the underlying modelling of H$_{\rm 2}$ and the description of stellar feedback in low-metallicity environments (LMF… ▽ More The reservoir of molecular gas (H$_{\rm 2}$) represents the fuel for the star formation (SF) of a galaxy. Connecting the star formation rate (SFR) to the available H$_{\rm 2}$ is key to accurately model SF in cosmological simulations of galaxy formation. We investigate how modifying the underlying modelling of H$_{\rm 2}$ and the description of stellar feedback in low-metallicity environments (LMF, i.e. low-metallicity stellar feedback) in cosmological, zoomed-in simulations of a Milky Way-size halo influences the formation history of the forming, spiral galaxy and its final properties. We exploit two different models to compute the molecular fraction of cold gas (f$_{\rm H_{\rm 2}}$): $i)$ the theoretical model by Krumholz et al. (2009b) and $ii)$ the phenomenological prescription by Blitz & Rosolowsky (2006). We find that the model adopted to estimate f$_{\rm H_{\rm 2}}$ plays a key role in determining final properties and in shaping the morphology of the galaxy. The clumpier interstellar medium (ISM) and the more complex H$_{\rm 2}$ distribution that the Krumholz et al. (2009b) model predicts result in better agreement with observations of nearby disc galaxies. This shows how crucial it is to link the SFR to the physical properties of the star-forming, molecular ISM. The additional source of energy that LMF supplies in a metal-poor ISM is key in controlling SF at high redshift and in regulating the reservoir of SF across cosmic time. Not only is LMF able to regulate cooling properties of the ISM, but it also reduces the stellar mass of the galaxy bulge. These findings can foster the improvement of the numerical modelling of SF in cosmological simulations. △ Less

Submitted 28 October, 2022; v1 submitted 27 July, 2022; originally announced July 2022.

Comments: accepted for publication in MNRAS

Dust evolution with MUPPI in cosmological volumes

Authors: Massimiliano Parente, Cinthia Ragone-Figueroa, Gian Luigi Granato, Stefano Borgani, Giuseppe Murante, Milena Valentini, Alessandro Bressan, Andrea Lapi

Abstract: We study the evolution of dust in a cosmological volume using a hydrodynamical simulation in which the dust production is coupled with the MUPPI (MUlti Phase Particle Integrator) sub-resolution model of star formation and feedback. As for the latter, we keep as reference the model setup calibrated previously to match the general properties of Milky Way like galaxies in zoom-in simulations. However… ▽ More We study the evolution of dust in a cosmological volume using a hydrodynamical simulation in which the dust production is coupled with the MUPPI (MUlti Phase Particle Integrator) sub-resolution model of star formation and feedback. As for the latter, we keep as reference the model setup calibrated previously to match the general properties of Milky Way like galaxies in zoom-in simulations. However, we suggest that an increase of the star formation efficiency with the local dust to gas ratio would better reproduce the observed evolution of the cosmic star formation density. Moreover, the paucity of quenched galaxies at low redshift demands a stronger role of AGN feedback. We tune the parameters ruling direct dust production from evolved stars and accretion in the inter stellar medium to get scaling relations involving dust, stellar mass and metallicity in good agreement with observations. In low mass galaxies the accretion process is inefficient. As a consequence, they remain poorer in silicate and small grains than higher mass ones. We reproduce reasonably well the few available data on the radial distribution of dust outside the galactic region, supporting the assumption that the dust and gas dynamics are well coupled at galactic scales. △ Less

Submitted 6 July, 2022; v1 submitted 25 April, 2022; originally announced April 2022.

Comments: MNRAS accepted

Galaxies in the central regions of simulated galaxy clusters

Authors: Antonio Ragagnin, Massimo Meneghetti, Luigi Bassini, Cinthia Ragone-Figueroa, Gian Luigi Granato, Giulia Despali, Carlo Giocoli, Giovanni Granata, Lauro Moscardini, Pietro Bergamini, Elena Rasia, Milena Valentini, Stefano Borgani, Francesco Calura, Klaus Dolag, Claudio Grillo, Amata Mercurio, Giuseppe Murante, Priyamvada Natarajan, Piero Rosati, Giuliano Taffoni, Luca Tornatore, Luca Tortorelli

Abstract: In this paper, we assess the impact of numerical resolution and of the implementation of energy input from AGN feedback models on the inner structure of cluster sub-haloes in hydrodynamic simulations. We compare several zoom-in re-simulations of a sub-sample of the cluster-sized haloes studied in Meneghetti et al. (2020), obtained by varying mass resolution, softening length and AGN energy feedbac… ▽ More In this paper, we assess the impact of numerical resolution and of the implementation of energy input from AGN feedback models on the inner structure of cluster sub-haloes in hydrodynamic simulations. We compare several zoom-in re-simulations of a sub-sample of the cluster-sized haloes studied in Meneghetti et al. (2020), obtained by varying mass resolution, softening length and AGN energy feedback scheme. We study the impact of these different setups on the subhalo abundances, their radial distribution, their density and mass profiles and the relation between the maximum circular velocity, which is a proxy for subhalo compactness. Regardless of the adopted numerical resolution and feedback model, subhaloes with masses Msub < 1e11Msun/h, the most relevant mass-range for galaxy-galaxy strong lensing, have maximum circular velocities ~30% smaller than those measured from strong lensing observations of Bergamini et al. (2019). We also find that simulations with less effective AGN energy feedback produce massive subhaloes (Msub> 1e11 Msun/h ) with higher maximum circular velocity and that their Vmax - Msub relation approaches the observed one. However the stellar-mass number count of these objects exceeds the one found in observations and we find that the compactness of these simulated subhaloes is the result of an extremely over-efficient star formation in their cores, also leading to larger-than-observed subhalo stellar mass. We conclude that simulations are unable to simultaneously reproduce the observed stellar masses and compactness (or maximum circular velocities) of cluster galaxies. Thus, the discrepancy between theory and observations that emerged from the analysis of Meneghetti et al. (2020) persists. It remains an open question as to whether such a discrepancy reflects limitations of the current implementation of galaxy formation models or the LCDM paradigm. △ Less

Submitted 19 April, 2022; originally announced April 2022.

Comments: 11 pages, 10 figures, abstract is redacted to fit arXiv character count limit

Journal ref: A&A 665, A16 (2022)

The probability of galaxy-galaxy strong lensing events in hydrodynamical simulations of galaxy clusters

Authors: Massimo Meneghetti, Antonio Ragagnin, Stefano Borgani, Francesco Calura, Giulia Despali, Carlo Giocoli, Gian Luigi Granato, Claudio Grillo, Lauro Moscardini, Elena Rasia, Piero Rosati, Giuseppe Angora, Luigi Bassini, Pietro Bergamini, Gabriel B. Caminha, Giovanni Granata, Amata Mercurio, Robert Benton Metcalf, Priyamvada Natarajan, Mario Nonino, Giada Venusta Pignataro, Cinthia Ragone-Figueroa, Eros Vanzella, Ana Acebron, Klaus Dolag , et al. (5 additional authors not shown)

Abstract: Meneghetti et al. (2020) recently reported an excess of galaxy-galaxy strong lensing (GGSL) in galaxy clusters compared to expectations from the LCDM cosmological model. Theoretical estimates of the GGSL probability are based on the analysis of numerical hydrodynamical simulations in the LCDM cosmology. We quantify the impact of the numerical resolution and AGN feedback scheme adopted in cosmologi… ▽ More Meneghetti et al. (2020) recently reported an excess of galaxy-galaxy strong lensing (GGSL) in galaxy clusters compared to expectations from the LCDM cosmological model. Theoretical estimates of the GGSL probability are based on the analysis of numerical hydrodynamical simulations in the LCDM cosmology. We quantify the impact of the numerical resolution and AGN feedback scheme adopted in cosmological simulations on the predicted GGSL probability and determine if varying these simulation properties can alleviate the gap with observations. We repeat the analysis of Meneghetti et al. (2020) on cluster-size halos simulated with different mass and force resolutions and implementing several independent AGN feedback schemes. We find that improving the mass resolution by a factor of ten and twenty-five, while using the same galaxy formation model that includes AGN feedback, does not affect the GGSL probability. We find similar results regarding the choice of gravitational softening. On the contrary, adopting an AGN feedback scheme that is less efficient at suppressing gas cooling and star formation leads to an increase in the GGSL probability by a factor between three and six. However, we notice that such simulations form overly massive subhalos whose contribution to the lensing cross-section would be significant while their Einstein radii are too large to be consistent with the observations. The primary contributors to the observed GGSL cross-sections are subhalos with smaller masses, that are compact enough to become critical for lensing. The population with these required characteristics appears to be absent in simulations. △ Less

Submitted 19 April, 2022; originally announced April 2022.

Comments: 13 pages, 11 figures. Submitted for publication on Astronomy and Astrophysics

Journal ref: A&A 668, A188 (2022)

Machine Learning to identify ICL and BCG in simulated galaxy clusters

Authors: I. Marini, S. Borgani, A. Saro, G. Murante, G. L. Granato, C. Ragone-Figueroa, G. Taffoni

Abstract: Nowadays, Machine Learning techniques offer fast and efficient solutions for classification problems that would require intensive computational resources via traditional methods. We examine the use of a supervised Random Forest to classify stars in simulated galaxy clusters after subtracting the member galaxies. These dynamically different components are interpreted as the individual properties of… ▽ More Nowadays, Machine Learning techniques offer fast and efficient solutions for classification problems that would require intensive computational resources via traditional methods. We examine the use of a supervised Random Forest to classify stars in simulated galaxy clusters after subtracting the member galaxies. These dynamically different components are interpreted as the individual properties of the stars in the Brightest Cluster Galaxy (BCG) and IntraCluster Light (ICL). We employ matched stellar catalogues (built from the different dynamical properties of BCG and ICL) of 29 simulated clusters from the DIANOGA set to train and test the classifier. The input features are cluster mass, normalized particle cluster-centric distance, and rest-frame velocity. The model is found to correctly identify most of the stars, while the larger errors are exhibited at the BCG outskirt, where the differences between the physical properties of the two components are less obvious. We investigate the robustness of the classifier to numerical resolution, redshift dependence (up to $z=1$), and included astrophysical models. We claim that our classifier provides consistent results in simulations for $z<1$, at different resolution levels and with significantly different subgrid models. The phase-space structure is examined to assess whether the general properties of the stellar components are recovered: (i) the transition radius between BCG-dominated and ICL-dominated region is identified at $0.04$ \r200; (ii) the BCG outskirt ($> 0.1$ \r200) is significantly affected by uncertainties in the classification process. In conclusion, this work suggests the importance of employing Machine Learning to speed up a computationally expensive classification in simulations. △ Less

Submitted 3 June, 2022; v1 submitted 7 March, 2022; originally announced March 2022.

Comments: 19 pages, accepted to MNRAS

Velocity dispersion of the brightest cluster galaxies in cosmological simulations

Authors: I. Marini, S. Borgani, A. Saro, G. L. Granato, C. Ragone-Figueroa, B. Sartoris, K. Dolag, G. Murante, A. Ragagnin, Y. Wang

Abstract: Using the DIANOGA hydrodynamical zoom-in simulation set of galaxy clusters, we analyze the dynamics traced by stars belonging to the Brightest Cluster Galaxies (BCGs) and their surrounding diffuse component, forming the intracluster light (ICL), and compare it to the dynamics traced by dark matter and galaxies identified in the simulations. We compute scaling relations between the BCG and cluster… ▽ More Using the DIANOGA hydrodynamical zoom-in simulation set of galaxy clusters, we analyze the dynamics traced by stars belonging to the Brightest Cluster Galaxies (BCGs) and their surrounding diffuse component, forming the intracluster light (ICL), and compare it to the dynamics traced by dark matter and galaxies identified in the simulations. We compute scaling relations between the BCG and cluster velocity dispersions and their corresponding masses (i.e. $M_\mathrm{BCG}^{\star}$- $σ_\mathrm{BCG}^{\star}$, $M_{200}$- $σ_{200}$, $M_\mathrm{BCG}^{\star}$- $M_{200}$, $σ_\mathrm{BCG}^{\star}$- $σ_{200}$), we find in general a good agreement with observational results. Our simulations also predict $σ_\mathrm{BCG}^{\star}$- $σ_{200}$ relation to not change significantly up to redshift $z=1$, in line with a relatively slow accretion of the BCG stellar mass at late times. We analyze the main features of the velocity dispersion profiles, as traced by stars, dark matter, and galaxies. As a result, we discuss that observed stellar velocity dispersion profiles in the inner cluster regions are in excellent agreement with simulations. We also report that the slopes of the BCG velocity dispersion profile from simulations agree with what is measured in observations, confirming the existence of a robust correlation between the stellar velocity dispersion slope and the cluster velocity dispersion (thus, cluster mass) when the former is computed within $0.1 R_{500}$. Our results demonstrate that simulations can correctly describe the dynamics of BCGs and their surrounding stellar envelope, as determined by the past star-formation and assembly histories of the most massive galaxies of the Universe. △ Less

Submitted 1 September, 2021; originally announced September 2021.

Comments: 16 pages, 14 figures

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

Halo cluster shapes: Insights from simulated galaxies and ICL with prospects for weak lensing applications

Authors: Elizabeth J. Gonzalez, Cinthia Ragone-Figueroa, Carlos J. Donzelli, Martín Makler, Diego García Lambas, Gian Luigi Granato

Abstract: We present a detailed study of the shapes and alignments of different galaxy cluster components using hydrodynamical simulations. We compute shape parameters from the Dark Matter (DM) distribution, the galaxy members and the intra-cluster light (ICL). We assess how well the DM cluster shape can be constrained by means of the identified galaxy member positions and the ICL. Further, we address the d… ▽ More We present a detailed study of the shapes and alignments of different galaxy cluster components using hydrodynamical simulations. We compute shape parameters from the Dark Matter (DM) distribution, the galaxy members and the intra-cluster light (ICL). We assess how well the DM cluster shape can be constrained by means of the identified galaxy member positions and the ICL. Further, we address the dilution factor introduced when estimating the cluster elongation using weak-lensing stacking techniques, which arises due to the misalignment between the total surface mass distribution and the distribution of luminous tracers. The dilution is computed considering the alignment between the DM and the Brightest Cluster Galaxy, the galaxy members and the ICL. Our study shows that distributions of galaxy members and ICL are less spherical than the DM component, although both are well aligned with the semi-major axis of the later. We find that the distribution of galaxy members hosted in more concentrated subhalos is more elongated than the distribution of the DM. Moreover, these galaxies are better aligned with the dark matter component compared to the distribution of galaxies hosted in less concentrated subhalos. We conclude that the positions of galaxy members can be used as suitable tracers to estimate the cluster surface density orientation, even when a low number of members is considered. Our results provide useful information for interpreting the constraints on the shapes of galaxy clusters in observational studies. △ Less

Submitted 7 July, 2021; originally announced July 2021.

Comments: Submitted MNRAS. 12 pages, 7 figures (+ appendix)

Brightest Cluster Galaxies: the centre can(not?) hold

Authors: Roberto De Propris, Michael J. West, Felipe Andrade-Santos, Cinthia Ragone-Figueroa, Elena Rasia, William Forman, Christine Jones, Rain Kipper, Stefano Borgani, Diego Garcia Lambas, Elena A. Romashkova, Kishore C. Patra

Abstract: We explore the persistence of the alignment of brightest cluster galaxies (BCGs) with their local environment. We find that a significant fraction of BCGs do not coincide with the centroid of the X-ray gas distribution and/or show peculiar velocities (they are not at rest with respect to the cluster mean). Despite this, we find that BCGs are generally aligned with the cluster mass distribution eve… ▽ More We explore the persistence of the alignment of brightest cluster galaxies (BCGs) with their local environment. We find that a significant fraction of BCGs do not coincide with the centroid of the X-ray gas distribution and/or show peculiar velocities (they are not at rest with respect to the cluster mean). Despite this, we find that BCGs are generally aligned with the cluster mass distribution even when they have significant offsets from the X-ray centre and significant peculiar velocities. The large offsets are not consistent with simple theoretical models. To account for these observations BCGs must undergo mergers preferentially along their major axis, the main infall direction. Such BCGs may be oscillating within the cluster potential after having been displaced by mergers or collisions, or the dark matter halo itself may not yet be relaxed. △ Less

Submitted 8 November, 2020; v1 submitted 19 October, 2020; originally announced October 2020.

Comments: In press, MNRAS

Dust evolution in zoom-in cosmological simulations of galaxy formation

Authors: Gian Luigi Granato, Cinthia Ragone-Figueroa, Antonela Taverna, Laura Silva, Milena Valentini, Stefano Borgani, Pierluigi Monaco, Giuseppe Murante, Luca Tornatore

Abstract: We present cosmological zoom-in hydro-dynamical simulations for the formation of disc galaxies, implementing dust evolution and dust promoted cooling of hot gas. We couple an improved version of our previous treatment of dust evolution, which adopts the two-size approximation to estimate the grain size distribution, with the MUPPI star formation and feedback sub-resolution model. Our dust evolutio… ▽ More We present cosmological zoom-in hydro-dynamical simulations for the formation of disc galaxies, implementing dust evolution and dust promoted cooling of hot gas. We couple an improved version of our previous treatment of dust evolution, which adopts the two-size approximation to estimate the grain size distribution, with the MUPPI star formation and feedback sub-resolution model. Our dust evolution model follows carbon and silicate dust separately. To distinguish differences induced by the chaotic behaviour of simulations from those genuinely due to different simulation set-up, we run each model six times, after introducing tiny perturbations in the initial conditions. With this method, we discuss the role of various dust-related physical processes and the effect of a few possible approximations adopted in the literature. Metal depletion and dust cooling affect the evolution of the system, causing substantial variations in its stellar, gas and dust content. We discuss possible effects on the Spectral Energy Distribution of the significant variations of the size distribution and chemical composition of grains, as predicted by our simulations during the evolution of the galaxy. We compare dust surface density, dust-to-gas ratio and small-to-big grain mass ratio as a function of galaxy radius and gas metallicity predicted by our fiducial run with recent observational estimates for three disc galaxies of different masses. The general agreement is good, in particular taking into account that we have not adjusted our model for this purpose. △ Less

Submitted 7 February, 2021; v1 submitted 12 October, 2020; originally announced October 2020.

Comments: 21 pages, 15 figures, MNRAS, accepted

The DIANOGA simulations of galaxy clusters: characterizing star formation in proto-clusters

Authors: L. Bassini, E. Rasia, S. Borgani, G. L. Granato, C. Ragone-Figueroa, V. Biffi, A. Ragagnin, K. Dolag, W. Lin, G. Murante, N. R. Napolitano, G. Taffoni, L. Tornatore, Y. Wang

Abstract: We studied the star formation rate (SFR) in cosmological hydrodynamical simulations of galaxy (proto-)clusters in the redshift range $0<z<4$, comparing them to recent observational studies; we also investigated the effect of varying the parameters of the star formation model on galaxy properties such as SFR, star-formation efficiency, and gas fraction. We analyze a set of zoom-in cosmological hydr… ▽ More We studied the star formation rate (SFR) in cosmological hydrodynamical simulations of galaxy (proto-)clusters in the redshift range $0<z<4$, comparing them to recent observational studies; we also investigated the effect of varying the parameters of the star formation model on galaxy properties such as SFR, star-formation efficiency, and gas fraction. We analyze a set of zoom-in cosmological hydrodynamical simulations centred on twelve clusters. The simulations are carried out with the GADGET-3 TreePM/SPH code which includes various subgrid models to treat unresolved baryonic physics, including AGN feedback. Simulations do not reproduce the high values of SFR observed within protoclusters cores, where the values of SFR are underpredicted by a factor $\gtrsim 4$ both at $z\sim2$ and $z\sim 4$. The difference arises as simulations are unable to reproduce the observed starburst population and is worsened at $z\sim 2$ because simulations underpredict the normalization of the main sequence of star forming galaxies (i.e., the correlation between stellar mass and SFR) by a factor of $\sim 3$. As the low normalization of the main sequence seems to be driven by an underestimated gas fraction, it remains unclear whether numerical simulations miss starburst galaxies due to a too low predicted gas fractions or too low star formation efficiencies. Our results are stable against varying several parameters of the star formation subgrid model and do not depend on the details of the AGN feedback. △ Less

Submitted 20 July, 2020; v1 submitted 24 June, 2020; originally announced June 2020.

Comments: 20 pages, 19 figures. Accepted for publication in A&A

Journal ref: A&A 642, A37 (2020)

Evolution and Role of Mergers in the BCG-Cluster Alignment. A View from Cosmological Hydro-Simulations

Authors: Cinthia Ragone-Figueroa, Gian Luigi Granato, Stefano Borgani, Roberto De Propris, Diego Garcia Lambas, Giuseppe Murante, Elena Rasia, Michael West

Abstract: Contradictory results have been reported on the time evolution of the alignment between clusters and their Brightest Cluster Galaxy (BCG). We study this topic by analyzing cosmological hydro-simulations of 24 massive clusters with $M_{200}|_{z=0} \gtrsim 10^{15}\, M_\odot$, plus 5 less massive with $1 \times 10^{14} \lesssim M_{200}|_{z=0} \lesssim 7 \times 10^{14}\, M_\odot$, which have already p… ▽ More Contradictory results have been reported on the time evolution of the alignment between clusters and their Brightest Cluster Galaxy (BCG). We study this topic by analyzing cosmological hydro-simulations of 24 massive clusters with $M_{200}|_{z=0} \gtrsim 10^{15}\, M_\odot$, plus 5 less massive with $1 \times 10^{14} \lesssim M_{200}|_{z=0} \lesssim 7 \times 10^{14}\, M_\odot$, which have already proven to produce realistic BCG masses. We compute the BCG alignment with both the distribution of cluster galaxies and the dark matter (DM) halo. At redshift $z=0$, the major axes of the simulated BCGs and their host cluster galaxy distributions are aligned on average within 20$^\circ$. The BCG alignment with the DM halo is even tighter. The alignment persists up to $z\lesssim2$ with no evident evolution. This result continues, although with a weaker signal, when considering the projected alignment. The cluster alignment with the surrounding distribution of matter ($3R_{200}$) is already in place at $z\sim4$ with a typical angle of $35^\circ$, before the BCG-Cluster alignment develops. The BCG turns out to be also aligned with the same matter distribution, albeit always to a lesser extent. These results taken together might imply that the BCG-Cluster alignment occurs in an outside-in fashion. Depending on their frequency and geometry, mergers can promote, destroy or weaken the alignments. Clusters that do not experience recent major mergers are typically more relaxed and aligned with their BCG. In turn, accretions closer to the cluster elongation axis tend to improve the alignment as opposed to accretions closer to the cluster minor axis. △ Less

Submitted 19 May, 2020; v1 submitted 27 April, 2020; originally announced April 2020.

Comments: 10 pages, 8 figures. MNRAS accepted. Comments welcome

Black hole mass of central galaxies and cluster mass correlation in cosmological hydro-dynamical simulations

Authors: Luigi Bassini, Elena Rasia, Stefano Borgani, Cinthia Ragone-Figueroa, Veronica Biffi, Klaus Dolag, Massimo Gaspari, Gian Luigi Granato, Giuseppe Murante, Giuliano Taffoni, Luca Tornatore

Abstract: Recently, relations connecting the SMBH mass of central galaxies and global properties of the hosting cluster, such as temperature and mass, were observed. We investigate the correlation between SMBH mass and cluster mass and temperature, their establishment and evolution. We compare their scatter to that of the classical $M_{\rm BH}-M_{\rm BCG}$ relation. We study how gas accretion and BH-BH merg… ▽ More Recently, relations connecting the SMBH mass of central galaxies and global properties of the hosting cluster, such as temperature and mass, were observed. We investigate the correlation between SMBH mass and cluster mass and temperature, their establishment and evolution. We compare their scatter to that of the classical $M_{\rm BH}-M_{\rm BCG}$ relation. We study how gas accretion and BH-BH mergers contribute to SMBH growth across cosmic time. We employed 135 groups and clusters with a mass range $1.4\times 10^{13}M_{\odot}-2.5\times 10^{15} M_{\odot}$ extracted from a set of 29 zoom-in cosmological hydro-dynamical simulations where the baryonic physics is treated with various sub-grid models, including feedback by AGN. In our simulations we find that $M_{\rm BH}$ correlates well with $M_{500}$ and $T_{500}$, with the scatter around these relations compatible within $2σ$ with the scatter around $M_{\rm BH}-M_{\rm BCG}$ at $z=0$. The $M_{\rm BH}-M_{500}$ relation evolves with time, becoming shallower at lower redshift as a direct consequence of hierarchical structure formation. On average, in our simulations the contribution of gas accretion to the total SMBH mass dominates for the majority of the cosmic time ($z>0.4$), while in the last 2 Gyr the BH-BH mergers become a larger contributor. During this last process, substructures hosting SMBHs are disrupted in the merger process with the BCG and the unbound stars enrich the diffuse stellar component rather than increase BCG mass. From the results obtained in our simulations with simple sub-grid models we conclude that the scatter around the $M_{\rm BH}-T_{500}$ relation is comparable to the scatter around the $M_{\rm BH}-M_{\rm BCG}$ relation and that, given the observational difficulties related to the estimation of the BCG mass, clusters temperature and mass can be a useful proxy for the SMBHs mass, especially at high redshift. △ Less

Submitted 27 October, 2019; v1 submitted 7 March, 2019; originally announced March 2019.

Comments: 16 pages, 17 figures, 3 tables

Journal ref: A&A 630, A144 (2019)

Dust Evolution in Galaxy Cluster Simulations

Authors: Eda Gjergo, Gian Luigi Granato, Giuseppe Murante, Cinthia Ragone-Figueroa, Luca Tornatore, Stefano Borgani

Abstract: We implement a state-of-the-art treatment of the processes affecting the production and Interstellar Medium (ISM) evolution of carbonaceous and silicate dust grains within SPH simulations. We trace the dust grain size distribution by means of a two-size approximation. We test our method on zoom-in simulations of four massive ($M_{200} \geq 3 \times 10^{14} M_{\odot}$) galaxy clusters. We predict t… ▽ More We implement a state-of-the-art treatment of the processes affecting the production and Interstellar Medium (ISM) evolution of carbonaceous and silicate dust grains within SPH simulations. We trace the dust grain size distribution by means of a two-size approximation. We test our method on zoom-in simulations of four massive ($M_{200} \geq 3 \times 10^{14} M_{\odot}$) galaxy clusters. We predict that during the early stages of assembly of the cluster at $z \gtrsim 3$, where the star formation activity is at its maximum in our simulations, the proto-cluster regions are rich of dusty gas. Compared to the case in which only dust production in stellar ejecta is active, if we include processes occurring in the cold ISM,the dust content is enhanced by a factor $2-3$. However, the dust properties in this stage turn out to be significantly different than those observationally derived for the {\it average} Milky Way dust, and commonly adopted in calculations of dust reprocessing. We show that these differences may have a strong impact on the predicted spectral energy distributions. At low redshift in star forming regions our model reproduces reasonably well the trend of dust abundances over metallicity as observed in local galaxies. However we under-produce by a factor of 2 to 3 the total dust content of clusters estimated observationally at low redshift, $z \lesssim 0.5$ using IRAS, Planck and Herschel satellites data. This discrepancy can be solved by decreasing the efficiency of sputtering which erodes dust grains in the hot Intracluster Medium (ICM). △ Less

Submitted 12 June, 2018; v1 submitted 18 April, 2018; originally announced April 2018.

Comments: 19 pages, 15 figures, accepted on June 8th, 2018, by MNRAS

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

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

Kinetic AGN Feedback Effects on Cluster Cool Cores Simulated using SPH

Authors: Paramita Barai, Giuseppe Murante, Stefano Borgani, Massimo Gaspari, Gian Luigi Granato, Pierluigi Monaco, Cinthia Ragone-Figueroa

Abstract: We implement novel numerical models of AGN feedback in the SPH code GADGET-3, where the energy from a supermassive black hole (BH) is coupled to the surrounding gas in the kinetic form. Gas particles lying inside a bi-conical volume around the BH are imparted a one-time velocity (10,000 km/s) increment. We perform hydrodynamical simulations of isolated cluster (total mass 10^14 /h M_sun), which is… ▽ More We implement novel numerical models of AGN feedback in the SPH code GADGET-3, where the energy from a supermassive black hole (BH) is coupled to the surrounding gas in the kinetic form. Gas particles lying inside a bi-conical volume around the BH are imparted a one-time velocity (10,000 km/s) increment. We perform hydrodynamical simulations of isolated cluster (total mass 10^14 /h M_sun), which is initially evolved to form a dense cool core, having central T<10^6 K. A BH resides at the cluster center, and ejects energy. The feedback-driven fast wind undergoes shock with the slower-moving gas, which causes the imparted kinetic energy to be thermalized. Bipolar bubble-like outflows form propagating radially outward to a distance of a few 100 kpc. The radial profiles of median gas properties are influenced by BH feedback in the inner regions (r<20-50 kpc). BH kinetic feedback, with a large value of the feedback efficiency, depletes the inner cool gas and reduces the hot gas content, such that the initial cool core of the cluster is heated up within a time 1.9 Gyr, whereby the core median temperature rises to above 10^7 K, and the central entropy flattens. Our implementation of BH thermal feedback (using the same efficiency as kinetic), within the star-formation model, cannot do this heating, where the cool core remains. The inclusion of cold gas accretion in the simulations produces naturally a duty cycle of the AGN with a periodicity of 100 Myr. △ Less

Submitted 9 June, 2016; v1 submitted 25 May, 2016; originally announced May 2016.

Comments: 22 pages, 11 figures, version accepted for publication in MNRAS, references and minor revisions added

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 early phases of galaxy clusters formation in IR: coupling hydrodynamical simulations with GRASIL3D

Authors: Gian Luigi Granato, Cinthia Ragone-Figueroa, Rosa Dominguez-Tenreiro, Aura Obreja, Stefano Borgani, Gabriella De Lucia, Giuseppe Murante

Abstract: We compute and study the infrared and sub-mm properties of high redshift ($z \gtrsim 1$) simulated clusters and proto-clusters. The results of a large set of hydro-dynamical zoom-in simulations including active galactic nuclei (AGN) feedback, have been treated with the recently developed radiative transfer code GRASIL-3D, which accounts for the effect of dust reprocessing in an arbitrary geometry.… ▽ More We compute and study the infrared and sub-mm properties of high redshift ($z \gtrsim 1$) simulated clusters and proto-clusters. The results of a large set of hydro-dynamical zoom-in simulations including active galactic nuclei (AGN) feedback, have been treated with the recently developed radiative transfer code GRASIL-3D, which accounts for the effect of dust reprocessing in an arbitrary geometry. Here, we have slightly generalized the code to adapt it to the present purpose. Then we have post-processed boxes of physical size 2 Mpc encompassing each of the 24 most massive clusters identified at z=0, at several redshifts between 0.5 and 3, producing IR and sub-mm mock images of these regions and SEDs of the radiation coming out from them. While this field is in its infancy from the observational point of view, rapid development is expected in the near future thanks to observations performed in the far IR and sub-mm bands. Notably, we find that in this spectral regime our prediction are little affected by the assumption required by this post-processing, and the emission is mostly powered by star formation rather than accretion onto super massive black hole (SMBH). The comparison with the little observational information currently available, highlights that the simulated cluster regions never attain the impressive star formation rates suggested by these observations. This problem becomes more intriguing taking into account that the brightest cluster galaxies (BCGs) in the same simulations turn out to be too massive. It seems that the interplay between the feedback schemes and the star formation model should be revised, possibly incorporating a positive feedback mode. △ Less

Submitted 25 March, 2015; v1 submitted 18 December, 2014; originally announced December 2014.

Comments: 14 pages, 11 figures, MNRAS in press. Minor editorial improvements

Temperature Structure of the Intra-Cluster Medium from SPH and AMR simulations

Authors: Elena Rasia, Erwin T. Lau, Stefano Borgani, Daisuke Nagai, Klaus Dolag, Camille Avestruz, Gian Luigi Granato, Pasquale Mazzotta, Giuseppe Murante, Kaylea Nelson, Cinthia Ragone-Figueroa

Abstract: Analyses of cosmological hydrodynamic simulations of galaxy clusters suggest that X-ray masses can be underestimated by 10% to 30%. The largest bias originates by both violation of hydrostatic equilibrium and an additional temperature bias caused by inhomogeneities in the X-ray emitting intra-cluster medium (ICM). To elucidate on this large dispersion among theoretical predictions, we evaluate the… ▽ More Analyses of cosmological hydrodynamic simulations of galaxy clusters suggest that X-ray masses can be underestimated by 10% to 30%. The largest bias originates by both violation of hydrostatic equilibrium and an additional temperature bias caused by inhomogeneities in the X-ray emitting intra-cluster medium (ICM). To elucidate on this large dispersion among theoretical predictions, we evaluate the degree of temperature structures in cluster sets simulated either with smoothed-particle-hydrodynamics (SPH) and adaptive-mesh-refinement (AMR) codes. We find that the SPH simulations produce larger temperature variations connected to the persistence of both substructures and their stripped cold gas. This difference is more evident in no-radiative simulations, while it is reduced in the presence of radiative cooling. We also find that the temperature variation in radiative cluster simulations is generally in agreement with the observed one in the central regions of clusters. Around R_500 the temperature inhomogeneities of the SPH simulations can generate twice the typical hydrostatic-equilibrium mass bias of the AMR sample. We emphasize that a detailed understanding of the physical processes responsible for the complex thermal structure in ICM requires improved resolution and high sensitivity observations in order to extend the analysis to higher temperature systems and larger cluster-centric radii. △ Less

Submitted 7 August, 2014; v1 submitted 17 June, 2014; originally announced June 2014.

Comments: 13 pages, 12 figures, 4 tables

Journal ref: 2014, ApJ, 791, 96

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

Brightest cluster galaxies in cosmological simulations: achievements and limitations of AGN feedback models

Authors: Cinthia Ragone-Figueroa, Gian Luigi Granato, Giuseppe Murante, Stefano Borgani, Weiguang Cui

Abstract: We analyze the basic properties of Brightest Cluster Galaxies (BCGs) produced by state of the art cosmological zoom-in hydrodynamical simulations. These simulations have been run with different sub-grid physics included. Here we focus on the results obtained with and without the inclusion of the prescriptions for supermassive black hole (SMBH) growth and of the ensuing Active Galactic Nuclei (AGN)… ▽ More We analyze the basic properties of Brightest Cluster Galaxies (BCGs) produced by state of the art cosmological zoom-in hydrodynamical simulations. These simulations have been run with different sub-grid physics included. Here we focus on the results obtained with and without the inclusion of the prescriptions for supermassive black hole (SMBH) growth and of the ensuing Active Galactic Nuclei (AGN) feedback. The latter process goes in the right direction of decreasing significantly the overall formation of stars. However, BCGs end up still containing too much stellar mass, a problem that increases with halo mass, and having an unsatisfactory structure. This is in the sense that their effective radii are too large, and that their density profiles feature a flattening on scales much larger than observed. We also find that our model of thermal AGN feedback has very little effect on the stellar velocity dispersions, which turn out to be very large. Taken together, these problems, which to some extent can be recognized also in other numerical studies typically dealing with smaller halo masses, indicate that on one hand present day sub-resolution models of AGN feedback are not effective enough in diminishing the global formation of stars in the most massive galaxies, but on the other hand they are relatively too effective in their centers. It is likely that a form of feedback generating large scale gas outflows from BCGs precursors, and a more widespread effect over the galaxy volume, can alleviate these difficulties. △ Less

Submitted 22 September, 2013; v1 submitted 14 August, 2013; originally announced August 2013.

Comments: 17 pages, 14 figures, accepted for publication on MNRAS, comments welcome

Effects of Baryon Mass Loss on Profiles of Large Galactic Dark Matter Haloes

Authors: Cinthia Ragone-Figueroa, Gian Luigi Granato, Mario G. Abadi

Abstract: We perform controlled numerical experiments to asses the effect of baryon mass loss on the inner structure of large galactic dark matter haloes. This mass expulsion is intended to mimic both the supernovae and AGN feedbacks, as well as the evolution of stellar populations. This study is meant in particular for precursors of massive Early Type Galaxies, wherein strong AGN feedback (often dubbed "QS… ▽ More We perform controlled numerical experiments to asses the effect of baryon mass loss on the inner structure of large galactic dark matter haloes. This mass expulsion is intended to mimic both the supernovae and AGN feedbacks, as well as the evolution of stellar populations. This study is meant in particular for precursors of massive Early Type Galaxies, wherein strong AGN feedback (often dubbed "QSO mode" in galaxy formation models) has been proposed to remove on a short timescale, of the order of a few dynamical times, a substantial fraction of their baryons. In a previous paper we evaluated the observational consequences (size increase) of this process on the galactic structure (Ragone-Figueroa & Granato 2011). Here we focus on the distribution of dark matter in the galactic region. It is shown that the inner region of the DM halo expands and its density profile flattens by a sizeable amount, with little dependence on the expulsion timescale. We also evaluate the effect of the commonly made approximation of treating the baryonic component as a potential that changes in intensity without any variation in shape. This approximation leads to some underestimates of the halo expansion and its slope flattening. We conclude that cuspy density profiles in ETGs could be difficult to reconcile with an effective AGN (or stellar) feedback during the evolution of these systems. △ Less

Submitted 17 April, 2012; v1 submitted 7 February, 2012; originally announced February 2012.

Comments: 9 Pages, 5 figures, MNRAS accepted. Final version with minor revision

A Consistent Comparison of Bias Models using Observational Data

Authors: A. Papageorgiou, M. Plionis, S. Basilakos, C. Ragone-Figueroa

Abstract: We investigate five different models for the dark matter halo bias, ie., the ratio of the fluctuations of mass tracers to those of the underlying mass, by comparing their cosmological evolution using optical QSO and galaxy bias data at different redshifts, consistently scaled to the WMAP7 cosmology. Under the assumption that each halo hosts one extragalactic mass tracer, we use a $χ^2$ minimizatio… ▽ More We investigate five different models for the dark matter halo bias, ie., the ratio of the fluctuations of mass tracers to those of the underlying mass, by comparing their cosmological evolution using optical QSO and galaxy bias data at different redshifts, consistently scaled to the WMAP7 cosmology. Under the assumption that each halo hosts one extragalactic mass tracer, we use a $χ^2$ minimization procedure to determine the free parameters of the bias models as well as to statistically quantify their ability to represent the observational data. Using the Akaike information criterion we find that the model that represents best the observational data is the Basilakos & Plionis (2001; 2003) model with the tracer merger extension of Basilakos, Plionis & Ragone-Figueroa (2008) model. The only other statistically equivalent model, as indicated by the same criterion, is the Tinker et al. (2010) model. Finally, we find an average, over the different models, dark matter halo mass that hosts optical QSOs of: $M_h\simeq 2.7 (\pm 0.6) \times 10^{12} h^{-1} M_{\odot}$, while the corresponding value for optical galaxies is: $M_h\simeq 6.3 (\pm 2.1) \times 10^{11} h^{-1} M_{\odot}$. △ Less

Submitted 23 January, 2012; originally announced January 2012.

Comments: MNRAS in press, 12 pages, 6 color figures, 4 tables

Puffing up early-type galaxies by baryonic mass loss: numerical experiments

Authors: Cinthia Ragone-Figueroa, Gian Luigi Granato

Abstract: Observations performed in the last few years indicate that most massive early-type galaxies (ETGs) observed at redshift z>1 exhibit sizes smaller by a factor of a few than local ETGs of analogous stellar mass. We present numerical simulations of the effect of baryonic mass loss on the structure of a spheroidal stellar system, embedded in a dark matter halo. This process, invoked as a possible expl… ▽ More Observations performed in the last few years indicate that most massive early-type galaxies (ETGs) observed at redshift z>1 exhibit sizes smaller by a factor of a few than local ETGs of analogous stellar mass. We present numerical simulations of the effect of baryonic mass loss on the structure of a spheroidal stellar system, embedded in a dark matter halo. This process, invoked as a possible explanation of the observed size increase of ETGs since z \sim 2, could be caused either by QSO/starburst driven galactic winds, promptly ejecting from Early Type Galaxies (ETGs) the residual gas and halting star formation (galactic winds), or by stellar mass returned to the ISM in the final stages of stellar evolution. Indeed, we find that a conceivable loss of about 50% of the baryonic mass can produce a significant size increase. However, the puffing up due to galactic winds occurs when the stellar populations are much younger than the estimated ages >0.5 Gyr of compact high-z ETGs. Therefore, while it may have had a role in deciding the final structure of ETGs, it cannot explain the evolution observed so far of their size-mass relation; its signature should be searched for in much younger systems. Conversely, the mass loss due to stellar evolution could cause a relatively modest expansion of passively evolving stellar systems later on, contributing to, without dominating, the observed evolution of their mass-size relationship. △ Less

Submitted 7 March, 2011; v1 submitted 25 January, 2011; originally announced January 2011.

Comments: MNRAS, accepted 2011 March 07. Received 2011 February 25; in original form 2011 January 20

The Relation Between Halo Shape, Velocity Dispersion and Formation Time

Authors: C. Ragone-Figueroa, M. Plionis, M. Merchan, S. Gottlober, G. Yepes

Abstract: We use dark matter haloes identified in the MareNostrum Universe and galaxy groups identified in the Sloan Data Release 7 galaxy catalogue, to study the relation between halo shape and halo dynamics, parametrizing out the mass of the systems. A strong shape-dynamics, independent of mass, correlation is present in the simulation data, which we find it to be due to different halo formation times. Ea… ▽ More We use dark matter haloes identified in the MareNostrum Universe and galaxy groups identified in the Sloan Data Release 7 galaxy catalogue, to study the relation between halo shape and halo dynamics, parametrizing out the mass of the systems. A strong shape-dynamics, independent of mass, correlation is present in the simulation data, which we find it to be due to different halo formation times. Early formation time haloes are, at the present epoch, more spherical and have higher velocity dispersions than late forming-time haloes. The halo shape-dynamics correlation, albeit weaker, survives the projection in 2D (ie., among projected shape and 1-D velocity dispersion). A similar shape-dynamics correlation, independent of mass, is also found in the SDSS DR7 groups of galaxies and in order to investigate its cause we have tested and used, as a proxy of the group formation time, a concentration parameter. We have found, as in the case of the simulated haloes, that less concentrated groups, corresponding to late formation times, have lower velocity dispersions and higher elongations than groups with higher values of concentration, corresponding to early formation times. △ Less

Submitted 11 May, 2010; originally announced May 2010.

Comments: MNRAS in press (10 pages, 10 figures)

The Halo Mass-Bias Redshift Evolution in the $Λ$CDM Cosmology

Authors: S. Basilakos, M. Plionis, C. Ragone-Figueroa

Abstract: We derive an analytic model for the redshift evolution of linear-bias, allowing for interactions and merging of the mass-tracers, by solving a second order differential equation based on linear perturbation theory and the Friedmann-Lemaitre solutions of the cosmological field equations. We then study the halo-mass dependence of the bias evolution, using the dark matter halo distribution in a… ▽ More We derive an analytic model for the redshift evolution of linear-bias, allowing for interactions and merging of the mass-tracers, by solving a second order differential equation based on linear perturbation theory and the Friedmann-Lemaitre solutions of the cosmological field equations. We then study the halo-mass dependence of the bias evolution, using the dark matter halo distribution in a $Λ$CDM simulation in order to calibrate the free parameters of the model. Finally, we compare our theoretical predictions with available observational data and find a good agreement. In particular, we find that the bias of optical QSO's evolve differently than those selected in X-rays and that their corresponding typical dark matter halo mass is $\sim 10^{13} h^{-1} M_{\odot}$ and $\magcir 5 \times 10^{13} h^{-1} M_{\odot}$, respectively. △ Less

Submitted 25 January, 2008; originally announced January 2008.

Comments: 8 pages, 5 figures, accepted for publication in ApJ

Environmental Effects of Dark Matter Haloes: The Clustering-Substructure relation of Group-size Haloes

Authors: N. Espino-Briones, M. Plionis, C. Ragone-Figueroa

Abstract: We estimate the two-point correlation function of dark matter haloes, with masses >10^{13} h^{-1} Mo, that have or not significant substructure. The haloes are identified with a friends of friends algorithm in a large LCDM simulation at two redshift snapshots (z=0 and 1), while halo substructure is determined using an observationally driven method. We find in both epochs a clear and significant… ▽ More We estimate the two-point correlation function of dark matter haloes, with masses >10^{13} h^{-1} Mo, that have or not significant substructure. The haloes are identified with a friends of friends algorithm in a large LCDM simulation at two redshift snapshots (z=0 and 1), while halo substructure is determined using an observationally driven method. We find in both epochs a clear and significant signal by which haloes with substructure are more clustered than those with no-substructure. This is true for all the considered halo mass ranges, although for the highest halo masses the signal is noisy and present only out to ~20 h^{-1} Mpc. There is also a smooth increase of the halo correlation length with increasing amplitude of the halo substructure. We also find that substructured haloes are typically located in high-density large-scale environments, while the opposite is true for non-substructured haloes. If the haloes found in high-density regions have a relatively earlier formation time, as suggested by recent works, then they do indeed have more time to cluster than haloes, of a similar mass, which form later in the low-density regions. In such a case one would have naively expected that the former (earlier formed) haloes would typically be dynamically more relaxed than the latter (later formed). However, the higher merging and interaction rate,expected in high-density regions, could disrupt their relatively relaxed dynamical state and thus be the cause for the higher fraction of haloes with substructure found in such regions. △ Less

Submitted 17 July, 2007; originally announced July 2007.

Comments: ApJ. Letters (accepted)

Environmental Influences on the Morphology and Dynamics of Group Size Haloes

Authors: Cinthia Ragone-Figueroa, Manolis Plionis

Abstract: We use group size haloes identified with a ``friends of friends'' (FOF) algorithm in a concordance $Λ\rm{CDM}$ GADGET2 (dark matter only) simulation to investigate the dependence of halo properties on the environment at $z=0$. The study is carried out using samples of haloes at different distances from their nearest massive {\em cluster} halo. We find that the fraction of haloes with substructur… ▽ More We use group size haloes identified with a ``friends of friends'' (FOF) algorithm in a concordance $Λ\rm{CDM}$ GADGET2 (dark matter only) simulation to investigate the dependence of halo properties on the environment at $z=0$. The study is carried out using samples of haloes at different distances from their nearest massive {\em cluster} halo. We find that the fraction of haloes with substructure typically increases in high density regions. The halo mean axial ratio $<c/a>$ also increases in overdense regions, a fact which is true for the whole range of halo mass studied. This can be explained as a reflection of an earlier halo formation time in high-density regions, which gives haloes more time to evolve and become more spherical. Moreover, this interpretation is supported by the fact that, at a given halo-cluster distance, haloes with substructure are more elongated than their equal mass counterparts with no substructure, reflecting that the virialization (and thus sphericalization) process is interrupted by merger events. The velocity dispersion of low mass haloes with strong substructure shows a significant increase near massive clusters with respect to equal mass haloes with low-levels of substructure or with haloes found in low-density environments. The alignment signal between the shape and the velocity ellipsoid principal axes decreases going from lower to higher density regions, while such an alignment is stronger for haloes without substructure. We also find, in agreement with other studies, a tendency of halo major axes to be aligned and of minor axes to lie roughly perpendicular with the orientation of the filament within which the halo is embedded, an effect which is stronger in the proximity of the massive clusters. △ Less

Submitted 16 March, 2007; originally announced March 2007.

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

Journal ref: Mon.Not.Roy.Astron.Soc.377:1785-1794,2007

Morphological & Dynamical Properties of low redshift 2dFGRS Groups

Authors: M. Plionis, S. Basilakos, C. Ragone-Figueroa

Abstract: We estimate the average group morphological and dynamical characteristics of the Percolation-Inferred Galaxy Group (2PIGG) catalogue within z~0.08, for which the group space density is roughly constant. We quantify the different biases that enter in the determination of these characteristics and we devise statistical correction procedures to recover their bias free values. We find that the onl… ▽ More We estimate the average group morphological and dynamical characteristics of the Percolation-Inferred Galaxy Group (2PIGG) catalogue within z~0.08, for which the group space density is roughly constant. We quantify the different biases that enter in the determination of these characteristics and we devise statistical correction procedures to recover their bias free values. We find that the only acceptable morphological model is that of prolate, or triaxial with pronounced prolatness, group shapes having a roughly Gaussian intrinsic axial ratio distribution with mean ~0.46 and dispersion of ~0.16. After correcting for various biases, the most important of which is a redshift dependant bias, the median values of the virial mass and virial radius of groups with 4 to 30 galaxy members, is: Mv ~6 x 10^12 h_{72}^{-1} M_solar, Rv~ 0.4 h^{-1}_{72} Mpc, which are significantly smaller than recent literature values that do not take into account the previously mentioned biases. The group mean crossing time is ~1.5 Gyr's, independent of the group galaxy membership. We also find that there is a correlation of the group size, velocity dispersion and virial mass with the number of group member galaxies, a manifestation of the hierarchy of cosmic structures. △ Less

Submitted 22 August, 2006; originally announced August 2006.

Comments: ApJ in press

Journal ref: Astrophys.J.650:770-776,2006