Source author record

Gian Luigi Granato

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

ResearcherUnclaimed source record

astro-ph.CO astro-ph.GA astro-ph Applications astro-ph.HE astro-ph.IM physics.data-an

Catalog footprint

What is connected

18works

7topics

4close collaborators

Actions

Connect this record

Open graph Browse works

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

Building this map preview

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

preprint2025arXiv

High$-z$ [OI] emission lines: ColdSIM simulations and ALMA observations

Neutral-oxygen [OI] far-infrared emission lines at $63μ$m and $145μ$m are powerful probes of the physical conditions in the interstellar medium, although they have not been fully exploited in high-redshift studies. We investigate the connection between [OI] emission lines and key galaxy properties, such as star formation rate (SFR) and H$_2$ content. Our predictions are compared with existing observations and new data analysed in this work. We post-process the outputs of the ColdSIM cosmological simulations with the DESPOTIC model, taking into account [OI]$63μ$m self-absorption by cold foreground material. A Random Forest algorithm is employed to accelerate computations and new observational ALMA data for galaxies at redshift $z\simeq 5-7$ are used to validate our model. Our predictions show significant [OI]$63μ$m luminosities ($\approx 10^8\,\rm L_\odot$) for galaxies with SFRs of $\approx 10^2\,\rm M_\odot\,{\rm yr}^{-1}$. The $145μ$m line luminosity is typically $15 \%$ the [OI]$63μ$m one and is a factor $\approx 2-20$ below high-$z$ observations. Both [OI] lines correlate with SFR and molecular mass, but exhibit flattening in scaling relations with metallicity and stellar mass. Foreground self-absorption reduces the [OI] flux by a factor of $2-4$, consistent with empirical corrections in observational studies. We find typical line ratios of [OI]$63μ$m / [CII]$158μ$m $\approx 1$ and [OI]$145μ$m / [CII]$158μ$m $ \approx 0.2 \, -$ consistent with $z\gtrsim 6$ observations, but only when [OI]$63μ$m self-absorption is included. Both [OI]$63μ$m and [OI]$145μ$m lines serve as tracers of star formation and molecular gas at high redshift. Their joint detection can provide constraints on the properties of the early interstellar medium and self-absorption of the [OI]$63μ$m line.

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.

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.

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.

preprint2016arXiv

Kinetic AGN Feedback Effects on Cluster Cool Cores Simulated using SPH

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.

preprint2015arXiv

The early phases of galaxy clusters formation in IR: coupling hydrodynamical simulations with GRASIL3D

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.

preprint2014arXiv

Black hole accretion preferentially occurs in gas rich galaxies

We have investigated the gas content of a sample of several hundred AGN host galaxies at z$<$1 and compared it with a sample of inactive galaxies, matched in bins of stellar mass and redshift. Gas masses have been inferred from the dust masses, obtained by stacked Herschel far-IR and sub-mm data in the GOODS and COSMOS fields, under reasonable assumptions and metallicity scaling relations for the dust-to-gas ratio. We find that AGNs are on average hosted in galaxies much more gas rich than inactive galaxies. In the vast majority of stellar mass bins, the average gas content of AGN hosts is higher than in inactive galaxies. The difference is up to a factor of ten higher in low stellar mass galaxies, with a significance of 6.5$σ$. In almost half of the AGN sample the gas content is three times higher than in the control sample of inactive galaxies. Our result strongly suggests that the probability of having an AGN activated is simply driven by the amount of gas in the host galaxy; this can be explained in simple terms of statistical probability of having a gas cloud falling into the gravitational potential of the black hole. The increased probability of an AGN being hosted by a star-forming galaxy, identified by previous works, may be a consequence of the relationship between gas content and AGN activity, found in this paper, combined with the Schmidt-Kennicutt law for star formation.

preprint2013arXiv

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

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.

preprint2012arXiv

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

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.

preprint2012arXiv

Fast radiative transfer of dust reprocessing in semi-analytic models with artificial neural networks

A serious concern for semi-analytical galaxy formation models, aiming to simulate multi-wavelength surveys and to thoroughly explore the model parameter space, is the extremely time consuming numerical solution of the radiative transfer of stellar radiation through dusty media. To overcome this problem, we have implemented an artificial neural network algorithm in the radiative transfer code GRASIL, in order to significantly speed up the computation of the infrared SED. The ANN we have implemented is of general use, in that its input neurons are defined as those quantities effectively determining the shape of the IR SED. Therefore, the training of the ANN can be performed with any model and then applied to other models. We made a blind test to check the algorithm, by applying a net trained with a standard chemical evolution model (i.e. CHE_EVO) to a mock catalogue extracted from the SAM MORGANA, and compared galaxy counts and evolution of the luminosity functions in several near-IR to sub-mm bands, and also the spectral differences for a large subset of randomly extracted models. The ANN is able to excellently approximate the full computation, but with a gain in CPU time by $\sim 2$ orders of magnitude. It is only advisable that the training covers reasonably well the range of values of the input neurons in the application. Indeed in the sub-mm at high redshift, a tiny fraction of models with some sensible input neurons out of the range of the trained net cause wrong answer by the ANN. These are extreme starbursting models with high optical depths, favorably selected by sub-mm observations, and difficult to predict a priori.

preprint2011arXiv

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

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.

preprint2008arXiv

Evolution of asymptotic giant branch stars II. Optical to far-infrared isochrones with improved TP-AGB models

We present a large set of theoretical isochrones, whose distinctive features mostly reside on the greatly improved treatment of the thermally pulsing asymptotic giant branch (TP-AGB) phase. Essentially, we have coupled the TP-AGB tracks described in Paper I, at their stages of pre-flash quiescent H-shell burning, with the evolutionary tracks for the previous evolutionary phases from Girardi et al. (2000). Theoretical isochrones for any intermediate value of age and metallicity are then derived by interpolation in the grids. We take care that the isochrones keep, to a good level of detail, the several peculiarities present in these TP-AGB tracks. Theoretical isochrones are then converted to about 20 different photometric systems -- including traditional ground-based systems, and those of recent major wide-field surveys such as SDSS, OGLE, DENIS, 2MASS, UKIDSS, etc., -- by means of synthetic photometry applied to an updated library of stellar spectra, suitably extended to include C-type stars. Finally, we correct the predicted photometry by the effect of circumstellar dust during the mass-losing stages of the AGB evolution, which allows us to improve the results for the optical-to-infrared systems, and to simulate mid- and far-IR systems such as those of Spitzer and AKARI. Access to the data is provided both via a web repository of static tables (http://stev.oapd.inaf.it/dustyAGB07 and CDS), and via an interactive web interface (http://stev.oapd.inaf.it/cmd) that provides tables for any intermediate value of age and metallicity, for several photometric systems, and for different choices of dust properties.

preprint2007arXiv

Exploring Supermassive Black Hole Growth with ALMA

Massive tori with $\approx 10^{8-9}M_{\odot}$ are predicted to extend on $\sim $100 pc scale around the centre of elliptical galaxy progenitors by a model of a supermassive black hole (SMBH) growth coeval to the spheroidal population of the host galaxy. Direct detection of such massive tori would cast light on a key physical condition that allows the rapid growth of SMBHs and the appearance of QSOs at high redshift. For this reason, we examine the detectability of such structures at substantial redshift with the Atacama Large Millimeter Array (ALMA). We propose that submillimeter galaxies (SMGs) are the best targets to test our predictions. In order to assess the observational feasibility, we estimate the expected number counts of SMGs with massive tori and check the detectability with the ALMA instrument, the unique facility which can resolve the central region of high redshift objects. Our work shows that ALMA will be able to resolve and detect high-$J$ ($J >$ 4) CO emissions from $\sim$100 pc scale extended massive tori up to $z\approx2$. Observations of lensed SMGs will yield excellent spatial resolution, allowing even to resolve their massive tori at higher redshift. We discuss further the detectability of the HCN molecule, as a better tracer of the high density gas expected in such tori. The final goal of these kind of observations is to pinpoint possible physical mechanisms that storage in the very central galactic regions very large amount gas on timescale of several 10$^{8}$ yr.

preprint2001arXiv

Dust and Nebular Emission in Star Forming Galaxies

Star forming galaxies exhibit a variety of physical conditions, from quiescent normal spirals to the most powerful dusty starbursts. In order to study these complex systems, we need a suitable tool to analyze the information coming from observations at all wavelengths. We present a new spectro-photometric model which considers in a consistent way starlight as reprocessed by gas and dust. We discuss preliminary results to interpret some observed properties of VLIRGs.