Researcher profile

Stefano Borgani

Stefano Borgani contributes to research discovery and scholarly infrastructure.

ResearcherAffiliation not importedOpen to collaborate

Trust snapshot

Quick read

Trust 21 - EmergingVerification L1Unclaimed author
11works
0followers
6topics
4close collaborators

Actions

Decide how to stay connected

Follow researcher0

Identity and collaboration

How to connect with this researcher

Claiming links this public author record to a researcher profile and unlocks direct collaboration workflows.

Log in to claim

Direct collaboration

Open a focused conversation when the fit is right

Claim this author entity first to unlock direct invitations.

Research graph

See the researcher in context

Open full explorer

Inspect adjacent work, topics, institutions and collaborators without jumping out to a separate graph page.

Building this graph slice

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

Published work

11 published item(s)

preprint2025arXiv

Dynamical friction and massive black hole orbits: analytical predictions and numerical solutions

We investigate the orbital decay of a massive BH embedded in a dark matter halo and a stellar bulge, using both analytical and numerical simulations with the aim of developing and validating a reliable dynamical friction (DF) correction across simulation resolutions. We develop a Python-based library to solve the equations of motion of the BH and provide an analytical framework for the numerical results. Then, we carry out simulations at different resolutions and for different softening choices using the Tree-PM code OpenGADGET3, where we implement an improved DF correction based on a kernel-weighted local density estimation. Our results demonstrate that the DF correction significantly accelerates BH sinking and ensures convergence with increasing resolution, closely matching analytical predictions. We find that in low-resolution regimes - particularly when the BH mass is smaller than that of the background particles - our DF model still effectively controls BH dynamics. Contrary to expectations, the inclusion of a stellar bulge can delay sinking due to numerical heating, an effect partially mitigated by the DF correction. We conclude that our refined DF implementation provides a robust framework for modeling BH dynamics both in controlled simulation setups of galaxies and in large-scale cosmological simulations. This will be crucial for future simulation campaigns, to enable more accurate predictions of AGN accretion and feedback, and to estimate gravitational-wave event rates.

preprint2022arXiv

\textsc{The Three Hundred} project: The \textsc{Gizmo-Simba} run

We introduce \textsc{Gizmo-Simba}, a new suite of galaxy cluster simulations within \textsc{The Three Hundred} project. \textsc{The Three Hundred} consists of zoom re-simulations of 324 clusters with $M_{200}\gtrsim 10^{14.8}M_\odot$ drawn from the MultiDark-Planck $N$-body simulation, run using several hydrodynamic and semi-analytic codes. The \textsc{Gizmo-Simba} suite adds a state-of-the-art galaxy formation model based on the highly successful {\sc Simba} simulation, mildly re-calibrated to match $z=0$ cluster stellar properties. Comparing to \textsc{The Three Hundred} zooms run with \textsc{Gadget-X}, we find intrinsic differences in the evolution of the stellar and gas mass fractions, BCG ages, and galaxy colour-magnitude diagrams, with \textsc{Gizmo-Simba} generally providing a good match to available data at $z \approx 0$. \textsc{Gizmo-Simba}'s unique black hole growth and feedback model yields agreement with the observed BH scaling relations at the intermediate-mass range and predicts a slightly different slope at high masses where few observations currently lie. \textsc{Gizmo-Simba} provides a new and novel platform to elucidate the co-evolution of galaxies, gas, and black holes within the densest cosmic environments.

preprint2022arXiv

Dust evolution with MUPPI in cosmological volumes

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.

preprint2022arXiv

Galaxies in the central regions of simulated galaxy clusters

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.

preprint2022arXiv

Weighing Cosmic Structures with Clusters of Galaxies and the Intergalactic Medium

We present an analysis aimed at combining cosmological constraints from number counts of galaxy clusters identified through the Sunyaev-Zeldovich effect, obtained with the South Pole Telescope (SPT), and from Lyman-$α$ spectra obtained with the MIKE/HIRES and X-shooter spectrographs. The SPT cluster analysis relies on mass calibration based on weak lensing measurements, while the Lyman-$α$ analysis is built over a suite of hydrodynamical simulations for the extraction of mock spectra. The resulting constraints exhibit a tension ($\sim 3.3σ$) between the low $σ_8$ values preferred by the low-redshift cluster data, $σ_8=0.74 ^{+0.03}_{-0.04}$, and the higher one preferred by the high-redshift Lyman-$α$ data, $σ_8=0.91 ^{+0.03}_{-0.03}$. We present a detailed analysis in order to understand the origin of this tension and, in particular, to establish whether it arises from systematic uncertainties related to the assumptions underlying the analyses of cluster counts and/or Lyman-$α$ forest. We found this tension to be robust with respect to the choice of modeling of the IGM, even when including possible systematics from unaccounted sub-Damped Lyman-$α$ (DLA) and Lyman-limit systems (LLS) in the Lyman-$α$ data. We conclude that to solve this tension from the SPT side would require a large bias on the cluster mass estimate, or from the Lyman-$α$ side large unaccounted errors on the Lyman-$α$ mean fluxes, respectively. Our results have important implications for future analyses based on cluster number counts from future large photometric surveys (e.g. Euclid and LSST) and on larger samples of high-redshift quasar spectra (e.g. DESI and WEAVE surveys). If confirmed at the much higher statistical significance reachable by such surveys, this tension could represent a significant challenge for the standard $Λ$CDM paradigm.

preprint2021arXiv

Dust evolution in zoom-in cosmological simulations of galaxy formation

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.

preprint2020arXiv

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

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.

preprint2020arXiv

Iron in X-COP: tracing enrichment in cluster outskirts with high accuracy abundance profiles

We present the first metal abundance profiles for a representative sample of massive clusters. Our measures extend to $R_{500}$ and are corrected for a systematic error plaguing previous outskirt estimates. Our profiles flatten out at large radii, admittedly not a new result, however the radial range and representative nature of our sample extends its import well beyond previous findings. We find no evidence of segregation between cool-core and non-cool-core systems beyond $\sim 0.3 R_{500}$, implying that, as was found for thermodynamic properties (Ghirardini et al, 2019), the physical state of the core does not affect global cluster properties. Our mean abundance within $R_{500}$ shows a very modest scatter, $< $15%, suggesting the enrichment process must be quite similar in all these massive systems. This is a new finding and has significant implications on feedback processes. Together with results from thermodynamic properties presented in a previous X-COP paper, it affords a coherent picture where feedback effects do not vary significantly from one system to another. By combing ICM with stellar measurements we have found the amount of Fe diffused in the ICM to be about ten times higher than that locked in stars. Although our estimates suggest, with some strength, that the measured iron mass in clusters is well in excess of the predicted one, systematic errors prevent us from making a definitive statement. Further advancements will only be possible when systematic uncertainties, principally those associated to stellar masses, both within and beyond $R_{500}$, can be reduced.

preprint2019arXiv

Impact of AGN feedback on galaxies and their multiphase ISM across cosmic time

We present simulations of galaxy formation, based on the GADGET-3 code, in which a sub-resolution model for star formation and stellar feedback is interfaced with a new model for AGN feedback. Our sub-resolution model describes a multiphase ISM, accounting for hot and cold gas within the same resolution element: we exploit this feature to investigate the impact of coupling AGN feedback energy to the different phases of the ISM over cosmic time. Our fiducial model considers that AGN feedback energy coupling is driven by the covering factors of the hot and cold phases. We perform a suite of cosmological hydrodynamical simulations of disc galaxies ($M_{\rm halo, \, DM} \simeq 2 \cdot 10^{12}$ M$_{\odot}$, at $z=0$), to investigate: $(i)$ the effect of different ways of coupling AGN feedback energy to the multiphase ISM; $(ii)$ the impact of different prescriptions for gas accretion (i.e. only cold gas, both cold and hot gas, with the additional possibility of limiting gas accretion from cold gas with high angular momentum); $(iii)$ how different models of gas accretion and coupling of AGN feedback energy affect the coevolution of supermassive BHs and their host galaxy. We find that at least a share of the AGN feedback energy has to couple with the diffuse gas, in order to avoid an excessive growth of the BH mass. When the BH only accretes cold gas, it experiences a growth that is faster than in the case in which both cold and hot gas are accreted. If the accretion of cold gas with high angular momentum is reduced, the BH mass growth is delayed, the BH mass at $z=0$ is reduced by up to an order of magnitude, and the BH is prevented from accreting below $z \lesssim 2$, when the galaxy disc forms.

preprint2017arXiv

Simulation-based marginal likelihood for cluster strong lensing cosmology

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

preprint2010arXiv

A decline and fall in the future of Italian Astronomy?

On May 27th 2010, the Italian astronomical community learned with concern that the National Institute for Astrophysics (INAF) was going to be suppressed, and that its employees were going to be transferred to the National Research Council (CNR). It was not clear if this applied to all employees (i.e. also to researchers hired on short-term contracts), and how this was going to happen in practice. In this letter, we give a brief historical overview of INAF and present a short chronicle of the few eventful days that followed. Starting from this example, we then comment on the current situation and prospects of astronomical research in Italy.