Hydrodynamical N-body Simulations Research Articles

Simulating the Arrival of the Southern Substructure in the Galaxy Cluster Abell 1758

Abell 1758 (z ∼ 0.278) is a galaxy cluster composed of two structures, A1758N and A1758S, separated by ∼2.2 Mpc. The northern cluster is itself a dissociative merging cluster that has already been modeled by dedicated simulations. Recent radio observations revealed the existence of a previously undetected bridge connecting A1758N and A1758S. New simulations are now needed to take into account the presence of A1758S. We wish to evaluate which orbital configuration would be compatible with a bridge between the clusters. Using N-body hydrodynamical simulations that build upon the previous model, we explore different scenarios that could have led to the current observed configuration. Five types of orbital approaches were tested: radial, tangential, vertical, postapocentric, and outgoing. We found that the incoming simulated scenarios are generally consistent with mild enhancements of gas density between the approaching clusters. The mock X-ray images exhibit a detectable bridge in all cases. Compared to measurements of Chandra data, the amplitude of the X-ray excess is overestimated by a factor of ∼2–3 in the best simulations. The scenario of tangential approach proved to be the one that best matches the properties of the profiles of X-ray surface brightness. The scenarios of radial approach of vertical approach are also marginally compatible.

Beyond the surface: hydrodynamical N-body simulations of the interacting dwarf galaxies NGC 5238 and UGC 8760

From deep imaging data obtained with the Large Binocular Telescope as part of the Smallest Scale of Hierarchy Survey (SSH), we have discovered low-surface brightness tidal features around NGC 5238 and UGC 8760, two nearby and relatively isolated dwarf galaxies with stellar masses of approximately $10^8 and $2 respectively. In this study, we present detailed hydrodynamical $N$-body simulations that explain the observed faint substructures as the outcome of interactions between the dwarf galaxies and smaller satellite systems. We show that the asymmetric stellar distribution of NGC 5238 and the low-luminosity substructures observed to the northeast of UGC 8760 can be well attributed to recent interactions with smaller galaxies, each with a stellar mass roughly a few $10^5 50 times less massive than their respective hosts. In the simulations, these satellites have stellar and dark-matter masses consistent with the ones predicted by Lambda CDM cosmology and share properties similar to those of local dwarf galaxies with similar stellar masses. The satellite-to-main galaxy mass ratio is approximately 1:10 in both cases. This satellite population aligns closely with predictions from cosmological simulations in terms of the number and mass relative to the host galaxy mass.

The Rotation of Classical Bulges in Barred Galaxies in the Presence of Gas

Barred galaxies often develop a box/peanut pseudobulge, but they can also host a nearly spherical classical bulge, which is known to gain rotation due to the bar. We aim to explore how the presence of gas impacts the rotation of classical bulges. We carried out a comprehensive set of hydrodynamical N-body simulations with different combinations of bulge masses and gas fractions. In these models, both massive bulges and high gas content tend to inhibit the formation of strong bars. For low-mass bulges, the resulting bar is stronger in cases of low gas content. In the stronger bar models, bulges acquire more angular momentum and thus display considerable rotational velocity. Such bulges also develop anisotropic velocity dispersions and become triaxial in shape. We found that the rotation of the bulge becomes less pronounced as the gas fraction is increased from 0 to 30%. These results indicate that the gas content has a significant effect on the dynamics of the classical bulge, because it influences bar strength. Particularly in the case of the low-mass bulges (10% bulge mass fraction), all of the measured rotational and structural properties of the classical bulge depend strongly and systematically on the gas content of the galaxy.

Unravelling the collision scenario of the dissociative galaxy cluster Abell 56 through hydrodynamic simulations

ABSTRACT In galaxy cluster collisions, the gas can be separated from dark matter haloes. Abell 56 displays signatures of a dissociative bullet-like merger with a possible high-inclination angle between the plane of orbit and the sky. Our objective is to provide a comprehensive description of the features observed in the collision scenario of Abell 56. Additionally, we aim to apply a potential weak lensing mass bias correction attributed to the merger to evaluate its impact on our findings. To investigate this, we perform tailored hydrodynamical N-body simulations, varying the impact parameter. We initially identified an early scenario at 0.12 Gyr after the central passage that reproduces some observational features. However, the mean temperature of 9.7 keV exceeded the observed value. Our best model corresponds to the late scenario at 0.52 Gyr after the pericentre, reproducing observed features of Abell 56, with an inclination of 58°. These features include the offset of 103 kpc between the main gas density peak and the south dark matter density peak, gas morphology, a line-of-sight relative velocity of 184 km s−1, and a mean temperature of 6.7 keV. This late model provides a plausible scenario to describe the dynamics of Abell 56. The weak lensing mass bias did not significantly impact the overall dynamics of this cluster merger.

Modelling post-reionization H i distributions in fuzzy dark matter cosmologies using conditional normalizing flows

ABSTRACT Upcoming 21-cm intensity mapping experiments like the Square Kilometer Array (SKA) hold significant potential to constrain the properties of dark matter. In this work, we model neutral hydrogen (H i) distributions using high-resolution hydrodynamical N-body simulations of both cold dark matter (CDM) and fuzzy dark matter (FDM) cosmologies in the post-reionization redshift range of z = 3.42−4.94. We show that the H i abundance decreases in FDM-like cosmologies. Extreme FDM models with m ∼ 10−22 eV are at odds with a range of measurements. Due to the increased halo bias, the H i bias increases, paralleled by the damped Lyman-α (DLA) bias which we infer from the cross-section of DLAs. The distribution of the latter in extreme FDM models has a high median at the low-mass end, which can be traced to the high column density of cosmic filaments. FDM models exhibit a very similar abundance of DLAs compared to CDM while sub-DLAs are already less abundant. We study the prospects of detecting the brightest H i peaks with SKA1-Low at z = 4.94, indicating moderate signal-to-noise ratios (SNRs) at angular resolution θA = 2 arcmin with a rapidly declining SNR for lower values of θA. After training the conditional normalizing flow network HIGlow on 2D H i maps, we interpolate its latent space of axion masses to predict the peak flux for a new, synthetic FDM cosmology, finding good agreement with expectations. This work thus underscores the potential of normalizing flows in capturing complex, non-linear structures within H i maps, offering a versatile tool for conditional sample generation and prediction tasks.

Galactic coronae in Milky Way-like galaxies: the role of stellar feedback in gas accretion

ABSTRACT Star-forming galaxies like the Milky Way are surrounded by a hot gaseous halo at the virial temperature – the so-called galactic corona – that plays a fundamental role in their evolution. The interaction between the disc and the corona has been shown to have a direct impact on accretion of coronal gas onto the disc with major implications for galaxy evolution. In this work, we study the gas circulation between the disc and the corona of star-forming galaxies like the Milky Way. We use high-resolution hydrodynamical N-body simulations of a Milky Way-like galaxy with the inclusion of an observationally motivated galactic corona. In doing so, we use SMUGGLE, an explicit interstellar medium (ISM), and stellar feedback model coupled with the moving-mesh code arepo. We find that the reservoir of gas in the galactic corona is sustaining star formation: the gas accreted from the corona is the primary fuel for the formation of new stars, helping in maintaining a nearly constant level of cold gas mass in the galactic disc. Stellar feedback generates a gas circulation between the disc and the corona (the so-called galactic fountain) by ejecting different gas phases that are eventually re-accreted onto the disc. The accretion of coronal gas is promoted by its mixing with the galactic fountains at the disc–corona interface, causing the formation of intermediate temperature gas that enhances the cooling of the hot corona. We find that this process acts as a positive feedback mechanism, increasing the accretion rate of coronal gas onto the galaxy.

Simulating nearly edge-on sloshing in the galaxy cluster Abell 2199

ABSTRACT Off-axis collisions between galaxy clusters may induce the phenomenon of sloshing, causing dense gas to be dragged from the cool core of a cluster, resulting in a spiral of enhanced X-ray emission. Abell 2199 displays signatures of sloshing in its core and it is possible that the orbital plane of the collision is seen nearly edge-on. We aim to evaluate whether the features of Abell 2199 can be explained by a sloshing spiral seen under a large inclination angle. To address this, we perform tailored hydrodynamical N-body simulations of a non-frontal collision with a galaxy group of $M_{200}=1.6\times 10^{13}\, {\rm M_{\odot }}$. We obtain a suitable scenario in which the group passed by the main cluster core 0.8 Gyr ago, with a pericentric separation of 292 kpc. A good agreement is obtained from the temperature maps as well as the residuals from a β-model fit to the simulated X-ray emission. We find that under an inclination of i = 70° the simulation results remain consistent with the observations.

Kinematical Analysis of Substructure in the Southern Periphery of the Large Magellanic Cloud

Abstract We report the first 3D kinematical measurements of 88 stars in the direction of several recently discovered substructures in the southern periphery of the Large Magellanic Cloud (LMC) using a combination of Gaia proper motions and radial velocities from the APOGEE-2 survey. More specifically, we explore stars in assorted APOGEE-2 pointings in a region of the LMC periphery where various overdensities of stars have previously been identified in maps of stars from Gaia and DECam. By using a model of the LMC disk rotation, we find that a sizable fraction of the APOGEE-2 stars have extreme space velocities that are distinct from, and not a simple extension of, the LMC disk. Using N-body hydrodynamical simulations of the past dynamical evolution and interaction of the LMC and Small Magellanic Cloud (SMC), we explore whether the extreme-velocity stars may be accounted for as tidal debris created in the course of that interaction. We conclude that the combination of LMC and SMC debris produced from their interaction is a promising explanation, although we cannot rule out other possible origins, and that these new data should be used to constrain future simulations of the LMC–SMC interaction. We also conclude that many of the stars in the southern periphery of the LMC lie outside of the LMC plane by several kiloparsecs. Given that the metallicity of these stars suggests that they are likely of Magellanic origin, our results suggest that a wider exploration of the past interaction history of the Magellanic Clouds is needed.

Clustering of Gravitational Wave and Supernovae events: a multitracer analysis in Luminosity Distance Space

We study the clustering of Gravitational Wave (GW) merger events and Supernovae IA (SN), as cosmic tracers in Luminosity Distance Space. We modify the publicly available CAMB code to numerically evaluate auto- and cross- power spectra for the different sources, including Luminosity Distance Space distortion effects generated by peculiar velocities and lensing convergence. We perform a multitracer Fisher analysis to forecast expected constraints on cosmological and GW bias coefficients, using outputs from hydrodynamical N-body simulations to determine the bias fiducial model and considering future observations from the Vera Rubin Observatory and Einstein Telescope (ET), both single and in a 3 detector network configuration. We find that adding SN to the GW merger dataset considerably improves the forecast, mostly by breaking significant parameter degeneracies, with final constraints comparable to those obtainable from a Euclid-like survey. GW merger bias is forecasted to be detectable with good significance even in the single ET case.

Dancing in the void: hydrodynamical N-body simulations of the extremely metal-poor galaxy DDO 68

ABSTRACT Using hydrodynamical N-body simulations, we show that the observed structure and kinematics of the extremely metal-poor dwarf irregular galaxy DDO 68 is compatible with a merger event with at least two smaller satellite galaxies. We were able to obtain a self-consistent model that simultaneously reproduces several of its observed features, including the very asymmetric and disturbed shape of the stellar component, the overall ${\rm H\, \small {\rm I}}$ distribution and its velocity field, the arc-like stellar structure to the west, and the low surface brightness stellar stream to the north. The model implies the interaction of the main progenitor of DDO 68 with two systems with dynamical masses $7\times 10^8\, \mathrm{ M}_{\odot }$ and almost $10^8\, \mathrm{ M}_{\odot }$ – 1/20 and 1/150 times the dynamical mass of DDO 68, respectively. We show that the merger between DDO 68 and the most massive of its satellites offers a route to explain the large offset of DDO 68 from the mass–metallicity relation. Assuming that the interacting galaxies have metallicities prior to the merger compatible with those of galaxies with similar stellar masses, we provide quantitative evidence that gas mixing alone does not suffice at diluting the gas of the two components; according to our simulations, the ${\rm H\, \small {\rm II}}$ regions observed along the cometary tail trace the low metallicity of the accreted satellite rather than that of DDO 68’s main body. In this case, the mass corresponding to the low metallicity is that of the secondary body and DDO 68 becomes consistent with the mass–metallicity relation.

Black hole spin evolution in warped accretion discs

ABSTRACT Massive black holes (BHs) inhabiting galactic nuclei can be described by two parameters only, i.e. mass and spin, that change through cosmic time in response to accretion and merger events. While most numerical simulations accurately track the BH mass, spin evolution is rarely taken into account. In this work, we implement and validate a self-consistent sub-grid model for the evolution of the BH mass and spin via gas accretion in the hydrodynamics code gizmo. The model assumes that accretion from resolved scales does not occur instantaneously but is mediated by a sub-grid geometrically thin α-disc. After validating our model semi-analytically, we test it in an idealized environment consisting of a circumnuclear disc, where gas accretion on to the accretion disc is consistently determined by gizmo. In the absence of any accretion-related feedback, the spin evolution closely traces that observed in the semi-analytical models, and depends on the free parameters of our implementation, such as the initial BH spin, angular momentum of the accretion disc, and radius at which the gas inflow circularizes. In gizmo, we also couple our model with the biconical-outflow model presented in a companion paper, wherein the feedback axis is always aligned with the BH spin. In this last case, the evolution of the central BH differs significantly from the previous cases, since the feedback process modifies the gas dynamics and its inflow rates from resolved scales. Such an interaction cannot be modelled by simple semi-analytical models and should be treated using full N-body hydrodynamical simulations.

Simulations of the merging galaxy cluster Abell 2034: what determines the level of separation between gas and dark matter

ABSTRACT Cluster mergers are an important laboratory for studying the behaviour of dark matter (DM) and intracluster gas. There are dissociative collisions that can separate the intracluster gas from the DM. Abell 2034 presents clear dissociative features observed by X-rays and gravitational lensing. The cluster, at z = 0.114, consists of two substructures with mass ratio of about 1:2.2, separated by ∼720 kpc. The X-ray emission peak is offcentred from the south DM peak by ∼350 kpc. Using N-body hydrodynamical simulations, we aim to reconstruct the dynamic history of the collision, reproducing the observed features, and also to explore the conditions that led to the dissociation. Our best model assuming that the collision is close to the plane of the sky, with a small impact parameter, observed 0.26 Gyr after central passage, reproduces the observed features of this cluster, such as the offset between X-ray and DM peaks, X-ray morphology, and temperatures. We explored several variations using different gas and DM concentrations for each cluster. The level of dissociation was quantified by the distances between X-ray and DM peaks, and also by the gas retention in the cluster cores. We found that the ratio of central gas densities is more important than the ratio of central DM densities in determining the level of dissociation.

Extreme close encounters between proto-Mercury and proto-Venus in terrestrial planet formation

ABSTRACT Modern models of terrestrial planet formation require solids depletion interior to 0.5–0.7 au in the planetesimal disc to explain the small mass of Mercury. The Earth and Venus analogues emerge after ∼100 Myr collisional growth, while Mercury forms in the diffusive tails of the planetesimal disc. We carried out 250 N-body simulations of planetesimal discs with mass confined to 0.7–1.0 au to study the statistics of close encounters that were recently proposed as an explanation for the high iron mass fraction in Mercury. We formed 39 Mercury analogues in total and all proto-Mercury analogues were scattered inwards by proto-Venus. Proto-Mercury typically experiences six extreme close encounters (closest approach smaller than six Venus radii) with Proto-Venus after Proto-Venus acquires 0.7 Venus Mass. At such close separation, the tidal interaction can already affect the orbital motion significantly such that the N-body treatment itself is invalid. More and closer encounters are expected should tidal dissipation of orbital angular momentum accounted. Hybrid N-body hydrodynamic simulations, treating orbital and encounter dynamics self-consistently, are desirable to evaluate the degree of tidal mantle stripping of proto-Mercury.

Dark matter cores and cusps in spiral galaxies and their explanations

We compare proposed solutions to the core vs cusp issue of spiral galaxies, which has also been framed as a diversity problem, and demonstrate that the cuspiness of dark matter halos is correlated with the stellar surface brightness. We compare the rotation curve fits to the SPARC sample from a self-interacting dark matter (SIDM) model, which self-consistently includes the impact of baryons on the halo profile, and hydrodynamical N-body simulations with cold dark matter (CDM) . The SIDM model predicts a strong correlation between the core size and the stellar surface density, and it provides the best global fit to the data. The CDM simulations without strong baryonic feedback effects fail to explain the large dark matter cores seen in low surface brightness galaxies. On the other hand, with strong feedback, CDM simulations do not produce galaxy analogs with high stellar and dark matter densities, and therefore they have trouble in explaining the rotation curves of high surface brightness galaxies. This implies that current feedback implementations need to be modified. We also explicitly show how the concentration-mass and stellar-to-halo mass relations together lead to a radial acceleration relation (RAR) in an averaged sense, and reiterate the point that the RAR does not capture the diversity of galaxy rotation curves in the inner regions. These results make a strong case for SIDM as the explanation for the cores and cusps of field galaxies.

Simulations of gas sloshing induced by a newly discovered gas poor substructure in galaxy cluster Abell 1644

ABSTRACT Collision events lead to peculiar morphologies in the intracluster gas of galaxies clusters. That seems to be the case of Abell 1644, a nearby galaxy cluster, composed of three main structures: the southern cluster that exhibits a spiral-like morphology, A1644S; the northern cluster seen in X-ray observations, A1644N1; and the recently discovered substructure, A1644N2. By means of N-body hydrodynamical simulations, we attempt to reconstruct the dynamical history of this system. These simulations resulted in two specific scenarios: (i) The collision between A1644S and A1644N2. Our best model has an inclination between the merger plane and the plane of the sky of 30°, and reaches the best morphology 1.6 Gyr after the pericentric passage. At this instant A1644N2 is gas poor, becoming nearly undetectable in X-ray emission. This model shows a good agreement with observations; (ii) The collision between A1644S and A1644N1. This approach did not give rise to results as satisfactory as the first scenario, due to great disturbances in density and mismatching temperature maps. As a complementary study, we perform a three-cluster simulation using as base the best-fitting model to reproduce the current state of A1644 with the three main structures. This scenario presented a good agreement to the global morphology of the observations. Thus, we find that the more likely scenario is a collision between A1644S and the newly discovered A1644N2, where A1644N1 may be present as long as it does not greatly interfere in the formation of the spiral feature.

Confirmation of NIKA2 investigation of the Sunyaev-Zel’dovich effect by using synthetic clusters of galaxies

The NIKA2 Sunyaev-Zel’dovich Large Program (SZLP) is focused on mapping the thermal SZ signal of a representative sample of selected Planck and ACT clusters spanning the redshift range 0.5 < z < 0.9. Hydrodynamical N-body simulations prove to be a powerful tool to endorse NIKA2 capabilities for estimating the impact of IntraCluster Medium (ICM) disturbances when re- covering the pressure radial profiles. For this goal we employ a subsample of objects, carefully extracted from the catalog Marenostrum MUltidark SImulations of galaxy Clusters (MUSIC), spanning equivalent redshift and mass ranges as the LPSZ. The joint analysis of real observations of the tSZ with NIKA2 and Planck enables to validate the NIKA2 pipeline and to estimate the ICM pressure profiles. Moreover, the possibility to identify a priori the dynamical state of the selected synthetic clusters allows us to verify the impact on the recovered ICM profile shapes and their scatters. Morphological analysis of maps of the Compton parameter seems to be a way to observationally segregate the sample based on the dynamical state in relaxed and disturbed synthetic clusters.

Metal enrichment signatures of the first stars on high-z DLAs

We use numerical N-body hydrodynamical simulations with varying PopIII stellar models to investigate the possibility of detecting first star signatures with observations of high-redshift damped Ly$\alpha$ absorbers (DLAs). The simulations include atomic and molecular cooling, star formation, energy feedback and metal spreading due to the evolution of stars with a range of masses and metallicities. Different initial mass functions (IMFs) and corresponding metal-dependent yields and lifetimes are adopted to model primordial stellar populations. The DLAs in the simulations are selected according to either the local gas temperature (temperature selected) or the host mass (mass selected). We find that 3\% (40\%) of mass (temperature) selected high-$z$ ($z\ge5.5$) DLAs retain signatures of pollution from PopIII stars, independently from the first star model. Such DLAs have low halo mass ($<10^{9.6}\,\rm M_{\odot}$), metallicity ($<10^{-3}\,\rm Z_{\odot}$) and star formation rate ($<10^{-1.5}\,\rm M_{\odot}\,yr^{-1}$). { Metal abundance ratios of DLAs imprinted in the spectra of QSO} can be useful tools to infer the properties of the polluting stellar generation and to constrain the first star mass ranges. Comparing the abundance ratios derived from our simulations to those observed in DLAs at $z\ge5$, we find that most of these DLAs are consistent within errors with PopII stars dominated enrichment and strongly disfavor the pollution pattern of very massive first stars (i.e. 100~$\rm M_{\odot}$-500~$\rm M_{\odot}$). However, some of them could still result from the pollution of first stars in the mass range [0.1, 100]~$\rm M_{\odot}$. In particular, we find that the abundance ratios from SDSS J1202+3235 are consistent with those expected from PopIII enrichment dominated by massive (but not extreme) first stars.

Phase-space Analysis in the Group and Cluster Environment: Time Since Infall and Tidal Mass Loss

Abstract Using the latest cosmological hydrodynamic N-body simulations of groups and clusters, we study how location in phase-space coordinates at z = 0 can provide information on environmental effects acting in clusters. We confirm the results of previous authors showing that galaxies tend to follow a typical path in phase-space as they settle into the cluster potential. As such, different regions of phase-space can be associated with different times since first infalling into the cluster. However, in addition, we see a clear trend between total mass loss due to cluster tides and time since infall. Thus, we find location in phase-space provides information on both infall time and tidal mass loss. We find the predictive power of phase-space diagrams remains even when projected quantities are used (i.e., line of sight velocities, and projected distances from the cluster). We provide figures that can be directly compared with observed samples of cluster galaxies and we also provide the data used to make them as supplementary data to encourage the use of phase-space diagrams as a tool to understand cluster environmental effects. We find that our results depend very weakly on galaxy mass or host mass, so the predictions in our phase-space diagrams can be applied to groups or clusters alike, or to galaxy populations from dwarfs up to giants.

Effect of primordial non-Gaussianities on galaxy clusters scaling relations

Galaxy clusters are a valuable source of cosmological information. Their formation and evolution depends on the underlying cosmology and on the statistical nature of the primordial density fluctuations. In this work we investigate the impact of primordial non-gaussianities (PNG) on the scaling properties of galaxy clusters. We performed a series of cosmological hydrodynamic N-body simulations featuring adiabatic gas physics and different levels of non-Gaussian initial conditions within the $\Lambda$CDM framework. We focus on the T-M, S-M, Y-M and Yx-M scalings relating the total cluster mass with temperature, entropy and SZ cluster integrated pressure that reflect the thermodynamical state of the intra-cluster medium. Our results show that PNG have an impact on cluster scalings laws. The mass power-law indexes of the scalings are almost unaffected by the existence of PNG but the amplitude and redshift evolution of their normalizations are clearly affected. The effect is stronger for the evolution of the Y-M and Yx-M normalizations, which change by as much as 22% and 16% when $f_{NL}$ varies from -500 to 500, respectively. These results are consistent with the view that positive/negative $f_{NL}$ affect cluster profiles due to an increase/decrease of cluster concentrations. At low values of $f_{NL}$, as suggested by present Planck constraints on a scale invariant $f_{NL}$, the impact on the scalings normalizations is only a few percent, which is small when compared with the effect of additional gas physics and other cosmological effects such as dark energy. However if $f_{NL}$ is in fact a scale dependent parameter, PNG may have larger positive/negative amplitudes at clusters scales and therefore our results suggest that PNG should be taken into account when galaxy cluster data is used to infer cosmological parameters or to asses the constraining power of future cluster surveys.

Multiscale mass transport inz ∼6 galactic discs: fuelling black holes

By using AMR cosmological hydrodynamic N-body zoom-in simulations, with the RAMSES code, we studied the mass transport processes onto galactic nuclei from high redshift up to $z\sim6$. Due to the large dynamical range of the simulations we were able to study the mass accretion process on scales from $\sim50[kpc]$ to $\sim$ few $1[pc]$. We studied the BH growth on to the galactic center in relation with the mass transport processes associated to both the Reynolds stress and the gravitational stress on the disc. Such methodology allowed us to identify the main mass transport process as a function of the scales of the problem. We found that in simulations that include radiative cooling and SNe feedback, the SMBH grows at the Eddington limit for some periods of time presenting $\langle f_{EDD}\rangle\approx 0.5$ throughout its evolution. The $\alpha$ parameter is dominated by the Reynolds term, $\alpha_R$, with $\alpha_R\gg 1$. The gravitational part of the $\alpha$ parameter, $\alpha_G$, has an increasing trend toward the galactic center at higher redshifts, with values $\alpha_G\sim 1$ at radii <$\sim$ few $ 10^1[pc]$ contributing to the BH fueling. In terms of torques, we also found that gravity has an increasing contribution toward the galactic center at earlier epochs with a mixed contribution above $\sim 100 [pc]$. This complementary work between pressure gradients and gravitational potential gradients allows an efficient mass transport on the disc with average mass accretion rates of the order $\sim$ few $1 [M_{\odot}/yr]$. These level of SMBH accretion rates found in our cosmological simulations are needed in all models of SMBH growth that attempt to explain the formation of redshift $6-7$ quasars.