BZPEERPapers, people, evidence and research paths
Menu

Discover

SearchFind papers, people, topics, institutions and opportunities from one place.GraphSee how papers, authors, topics and evidence connect.

Research tools

HomeContinue from the papers, people and topics that matter now.CollectionsTurn saved papers into reading lists, evidence maps and shareable context.ResearchReconstruct the path that led you from search to paper to question.CollaborationMove from a paper, topic or expert request into a focused thread.PublishCreate full-text research articles with metadata, formulas and versions.

Network

ExpertsFind specialists whose work and reviews explain why they can help.CommunitiesJoin research groups organized around papers, topics and shared evidence.

Opportunities

ListingsBrowse roles, calls and collaborations connected to real research context.

Account

Log inReturn to your saved papers, graph views and active threads.Create accountStart saving papers, building a research trail and joining teams.
Log inCreate account

Source author record

Alessandro Lupi

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

ResearcherUnclaimed source record
astro-ph.GAastro-ph.COastro-ph.HEastro-ph.IMastro-ph.SRgr-qcphysics.comp-ph

Catalog footprint

What is connected

25works
7topics
4close collaborators

Actions

Connect this record

Log in to claim
Open graphBrowse works

Research graph

See the researcher in context

Open full explorer

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.

Published work

25 published item(s)

preprint2025arXiv

Super-Eddington accretion in high-redshift quasar hosts: Black-hole driven outflows, galaxy quenching, and the nature of little red dots

The advent of the James Webb Space Telescope has revolutionised our understanding of the high-redshift Universe through its detection of bright, massive galaxies up to $z\gtrsim 10$ and its identification of peculiar sources called `little red dots' (LRDs). The origin of both classes of objects remains uncertain but is likely linked to the formation and early growth of the first massive black holes (MBHs), which may be more easily explained by invoking phases of super-Eddington accretion. In this study, we used a state-of-the-art zoom-in cosmological simulation of a quasar host to investigate whether these objects could resemble any of the peculiar sources observed with JWST during their assembly. We find that the impact of MBH feedback on star formation is typically moderate, with outflows preferentially escaping perpendicular to the galactic disc. However, for approximately ten percent of the galaxy's lifetime, the system enters a distinct quenched phase following rapid MBH growth driven by super-Eddington accretion. This phase culminates in a powerful feedback event, during which the MBH jet and disc-driven winds interact directly with the galactic disc and carve out a central cavity. We also find that, during the history of the quasar host progenitor, the spectral properties of the system can resemble both LRDs and quenched galaxies, depending on the specific evolutionary stage considered. These findings suggest that both conditions may represent transient phases in the life cycle of high-redshift galaxies.

0 citations0 reviewsNo code linkedOpen access
preprint2023arXiv

A general relativistic extension to mesh-free methods for hydrodynamics

The detection of gravitational waves has opened a new era for astronomy, allowing for the combined use of gravitational wave and electromagnetic emissions to directly probe the physics of compact objects, still poorly understood. So far, the theoretical modelling of these sources has mainly relied on standard numerical techniques as grid-based methods or smoothed particle hydrodynamics, with only a few recent attempts at using new techniques as moving-mesh schemes. Here, we introduce a general relativistic extension to the mesh-less hydrodynamic schemes in the code GIZMO, which benefits from the use of Riemann solvers and at the same time perfectly conserves angular momentum thanks to a generalised leap-frog integration scheme. We benchmark our implementation against many standard tests for relativistic hydrodynamics, either in one or three dimensions, and also test the ability to preserve the equilibrium solution of a Tolman-Oppenheimer-Volkoff compact star. In all the presented tests, the code performs extremely well, at a level at least comparable to other numerical techniques.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

Black hole spin evolution in warped accretion discs

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 onto 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 the radius at which the gas inflow circularises. 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 modeled by simple semi-analytical models and should be treated using full N-body hydrodynamical simulations.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

Disc instability and bar formation: view from the IllustrisTNG simulations

We make use of z = 0 samples of strongly barred and unbarred disc galaxies from the TNG100 and TNG50 cosmological hydrodynamical simulations to assess the performance of the simple disc instability criterion proposed by Efstathiou, Lake & Negroponte (1982) (ELN-criterion). We find that strongly barred galaxies generally assemble earlier, are more star-dominated in their central regions, and have more massive and more compact discs than unbarred galaxies. The ELN-criterion successfully identifies ~75% and ~80% of the strongly barred and the unbarred galaxies, respectively. Strongly barred galaxies that the criterion fails to identify tend to have more extended discs, higher spin values and bars that assembled later than is typical for the bulk of the barred population. The bars in many of these cases appear to be produced by an interaction with a close neighbour (i.e. to be externally triggered) rather than to result from secular growth in the disc. On the other hand, we find that unbarred galaxies misclassified as barred by the ELN-criterion typically have stellar discs similar to those of barred galaxies, although more extended in the vertical direction and less star-dominated in their central regions, possibly reflecting later formation times. In addition, the bulge component of these galaxies is significantly more prominent at early times than in the strongly barred sample. Thus, the ELN-criterion robustly identifies secular bar instabilities in most simulated disc galaxies, but additional environmental criteria are needed to account for interaction-induced bar formation.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

Enhanced star formation in $z\sim6$ quasar companions

Quasars powered by supermassive black holes (MBH, $>10^8~M_{\odot}$) at $z\sim 6$ are predicted to reside in cosmic over-dense regions. However, observations so far could not confirm this expectation due to limited statistics. The picture is further complicated by the possible effects of quasar outflows (i.e. feedback) that could either suppress or stimulate the star formation rate (SFR) of companion galaxies, thus modifying the expected bias. Here we quantify feedback effects on the properties and detectability of companions by comparing cosmological zoom-in simulations of a quasar in which feedback is either included or turned-off. With respect to the no-feedback case, companions (a) directly impacted by the outflow have their SFR increased by a factor $2-3$, and (b) tend to be more massive. Both effects shift the [CII]158$μ$m and UV luminosity functions toward brighter magnitudes. This leads us to conclude that quasar feedback slightly increases the effective quasar bias, boosting the number density of observable quasar companions, in agreement with what has been found around the brightest quasars of recent ALMA [CII] surveys. Deeper observations performed with JWST and/or ALMA will improve the statistical significance of this result by detecting a larger number of fainter quasar companions.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

Forming massive seed black holes in high-redshift quasar host progenitors

The presence of massive black holes (BHs) with masses of order $10^9\rm\, M_\odot$, powering bright quasars when the Universe was less than 1 Gyr old, poses strong constraints on their formation mechanism. Several scenarios have been proposed to date to explain massive BH formation, from the low-mass seed BH remnants of the first generation of stars to the massive seed BHs resulting from the rapid collapse of massive gas clouds. However, the plausibility of some of these scenarios to occur within the progenitors of high-z quasars has not yet been thoroughly explored. In this work, we investigate, by combining dark-matter only N-body simulations with a semi-analytic framework, whether the conditions for the formation of massive seed BHs from synchronised atomic-cooling halo pairs and/or dynamically-heated mini-haloes are fulfilled in the overdense regions where the progenitors of a typical high-redshift quasar host form and evolve. Our analysis shows that the peculiar conditions in such regions, i.e. strong halo clustering and high star formation rates, are crucial to produce a non-negligible number of massive seed BH host candidates: we find $\approx1400$ dynamically heated metal-free mini-haloes, including one of these which evolves to a synchronised pair and ends up in the massive quasar-host halo by $z=6$. This demonstrates that the progenitors of high-redshift quasar host haloes can harbour early massive seed BHs. Our results further suggest that multiple massive seed BHs may form in or near the quasar host's progenitors, potentially merging at lower redshifts and yielding gravitational wave events.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

High-redshift quasars and their host galaxies II: multiphase gas and stellar kinematics

Observations of $z \gtrsim 6$ quasars provide information on the early phases of the most massive black holes (MBHs) and galaxies. Current observations at sub-mm wavelengths trace cold and warm gas, and future observations will extend information to other gas phases and the stellar properties. The goal of this study is to examine the gas life cycle in a $z \gtrsim 6$ quasar: from accretion from the halo to the galaxy and all the way into the MBH, to how star formation and the MBH itself affect the gas properties. Using a very-high resolution cosmological zoom-in simulation of a $z=7$ quasar including state-of-the-art non-equilibrium chemistry, MBH formation, growth and feedback, we investigate the distribution of the different gas phases in the interstellar medium across cosmic time. We assess the morphological evolution of the quasar host using different tracers (star- or gas-based) and the thermodynamic distribution of the MBH accretion-driven outflows, finding that obscuration in the disc is mainly due to molecular gas, with the atomic component contributing at larger scales and/or above/below the disc plane. Moreover, our results also show that molecular outflows, if present, are more likely the result of gas being lifted near the MBH than production within the wind because of thermal instabilities. Finally, we also discuss how different gas phases can be employed to dynamically constrain the MBH mass, and argue that resolutions below $\sim 100$ pc yield unreliable estimates because of the strong contribution of the nuclear stellar component to the potential at larger scales.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

Maximally informed Bayesian modelling of disc galaxies

Dissecting the underlying structure of galaxies is of main importance in the framework of galaxy formation and evolution theories. While a classical bulge+disc decomposition of disc galaxies is usually taken as granted, this is only rarely solidly founded upon the full exploitation of the richness of data arising from spectroscopic studies with integral field units. In this work we describe a fully Bayesian estimation method of the global structure of disc galaxies which makes use of the wealth of photometric, kinematic, and mass-to-light ratio data, and that can be seen as a first step towards a machine-learning approach, certainly needed when dealing with larger samples of galaxies. Ours is a novel, hybrid line of action in tackling the problem of galactic parameter estimation, neither purely photometric nor orbit-based. Being rooted on a nested sampler, our code, which is available publicly as an online repository, allows for a statistical assessment of the need for multiple components in the dissecting process. As a first case-study the GPU-optimized code is applied to the S0 galaxy NGC-7683, finding that in this galaxy a pseudo-bulge, possibly the remnant of a bar-like structure, does exist in the center of the system. These results are then tested against the publicly available, orbit-based code DYNAMITE, finding substantial agreement.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

Morphological decomposition of TNG50 galaxies: methodology and catalogue

We present MORDOR (MORphological DecOmposeR, a new algorithm for structural decomposition of simulated galaxies based on stellar kinematics. The code measures the properties of up to five structural components (a thin/cold and a thick/warm disc, a classical and a secular bulge, and a spherical stellar halo), and determines the properties of a stellar bar (if present). A comparison with other algorithms presented in the literature yields overall good agreement, with MORDOR displaying a higher flexibility in correctly decomposing systems and identifying bars in crowded environments (e.g. with ongoing fly-bys, often observable in cosmological simulations). We use MORDOR to analyse galaxies in the TNG50 simulation and find the following: ($i$) the thick disc component undergoes the strongest evolution in the binding energy-circularity plane, as expected when disc galaxies decrease their turbulent-rotational support with cosmic time; ($ii$) smaller galaxies (with stellar mass, $10^{9} \lesssim M_{*} / {\rm M_{\odot}} \leq 5 \times 10^{9}$) undergo a major growth in their disc components after $z\sim 1$, whereas ($iii$) the most massive galaxies ($5 \times 10^{10} < M_{*} / {\rm M_{\odot}} \leq 5\times10^{11}$) evolve toward more spheroidal dominated objects down to $z=0$ due to frequent gravitational interactions with satellites; ($iv$) the fraction of barred galaxies grows rapidly at high redshift and stabilizes below $z\sim 2$, except for the most massive galaxies that show a decrease in the bar occupation fraction at low redshift; ($v$) galaxies with $M_{*} \sim 10^{11}~{\rm M_{\odot}}$ exhibit the highest relative occurrence of bars at $z=0$, in agreement with observational studies. We publicly release MORDOR and the morphological catalogue of TNG50 galaxies.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

Non-isotropic feedback from accreting spinning black holes

Active galactic nuclei (AGN) are massive black holes (BHs) caught in the act of accreting gas at the centre of their host galaxies. Part of the accreting mass is converted to energy and released into the surrounding medium, in a process loosely referred to as AGN feedback. Most numerical simulations include AGN feedback as a sub-grid model, wherein energy or momentum (or both) is coupled to the nearby gas. In this work, we implement a new momentum-driven model in the hydrodynamics code GIZMO, in which accretion from large scales is mediated by a sub-grid accretion disc model, and gas particles are stochastically kicked over a bi-conical region, to mimic observed kinetic winds. The feedback cone's axis can be set either parallel to the angular momentum of the gas surrounding the BH or to the BH spin direction, which is self-consistently evolved within the accretion-disc model. Using a circumnuclear disc (CND) as a test bed, we find that (i) the conical shape of the outflow is always visible and is weakly dependent on the launching orientation and aperture, resulting in comparable mass inflows and outflows; (ii) the cone's orientation is also similar amongst our tests, and it is not always the same as the initial value, due to the interaction with the CND playing a crucial role in shaping the outflow; and (iii) the velocity of the outflow, instead, differs and strongly depends on the interplay with the CND.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

On the low ortho-to-para H2 ratio in star-forming filaments

The formation of stars and planetary systems is a complex phenomenon, which relies on the interplay of multiple physical processes. Nonetheless, it represents a crucial stage for our understanding of the Universe, and in particular of the conditions leading to the formation of key molecules (e.g. water) on comets and planets. Herschel observations demonstrated that stars form out of gaseous filamentary structures in which the main constituent is molecular hydrogen (H$_2$). Depending on its nuclear spin H$_2$ can be found in two forms: `ortho' with parallel spins and `para' where the spins are anti-parallel. The relative ratio among these isomers, i.e. the ortho-to-para ratio (OPR), plays a crucial role in a variety of processes related to the thermodynamics of star-forming gas and to the fundamental chemistry affecting the deuteration of water in molecular clouds, commonly used to determine the origin of water in Solar System's bodies. Here, for the first time, we assess the evolution of the OPR starting from the warm neutral medium, by means of state-of-the-art three-dimensional magneto-hydrodynamic simulations of turbulent molecular clouds. Our results show that star-forming clouds exhibit a low OPR ($\ll 0.1$) already at moderate densities ($\sim$1000 cm$^{-3}$). We also constrain the cosmic rays ionisation rate, finding that $10^{-16}\,\rm s^{-1}$ is the lower limit required to explain the observations of diffuse clouds. Our results represent a step forward in the understanding of the star and planet formation process providing a robust determination of the chemical initial conditions for both theoretical and observational studies.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

The evolution of the barred galaxy population in the TNG50 simulation

We use the magnetic-hydrodynamical simulation TNG50 to study the evolution of barred massive disc galaxies. Massive spiral galaxies are already present as early as $z=4$, and bar formation takes place already at those early times. The bars grow longer and stronger as the host galaxies evolve, with the bar sizes increasing at a pace similar to that of the disc scale lengths. The bar fraction mildly evolves with redshift for galaxies with $M_{*}\geq10^{10}M\odot$, being greater than $\sim40\%$ at $0.5<z<3$ and $\sim30\%$ at $z=0$. When bars larger than a given physical size ($\geq 2\,\rm kpc$) or the angular resolution limit of twice the I-band angular PSF FWHM of the HST are considered, the bar fraction dramatically decreases with increasing redshift, reconciling the theoretical predictions with observational data. We find that barred galaxies have an older stellar population, lower gas fractions and star formation rates than unbarred galaxies. In most cases, the discs of barred galaxies assembled earlier and faster than the discs of unbarred galaxies. We also find that barred galaxies are typical in haloes with larger concentrations and smaller spin parameters than unbarred galaxies. Furthermore, the inner regions of barred galaxies are more baryon-dominated than those of unbarred galaxies but have comparable global stellar mass fractions. Our findings suggest that the bar population could be used as a potential tracer of the buildup of disc galaxies and their host haloes. With this paper, we release a catalogue of barred galaxies in TNG50 at $6$ redshifts between $z=4$ and $z=0$.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

The role of bars on the dynamical-friction driven inspiral of massive perturbers

In this paper, we systematically explore the impact of a galactic bar on the inspiral time-scale of a massive object (MO) within a Milky Way-like galaxy. We integrate the orbit of MOs in a multi-component galaxy model via a semi-analytical approach that accounts for dynamical friction generalized to rotationally supported backgrounds. We compare the MO evolution in a galaxy featuring a Milky Way-like rotating bar to the evolution within an analogous axisymmetric galaxy without the bar. In agreement with previous studies, we find that the bar presence may significantly affect the inspiral, sometimes making it shorter by a factor of a few, sometimes hindering it for a Hubble time. The erratic behaviour is mainly impacted by the relative phase at which the MO encounters the stronger bar-induced resonances. In particular, the effect of the bar is more prominent for initially in-plane, prograde MOs, especially those crossing the bar co-rotation radius or outer Lindblad resonance. In the barred galaxy, we find the sinking of the most massive MOs (>~10^7.5 Msun) approaching the galaxy from large separations (>~8 kpc) to be most efficiently hampered. Neglecting the effect of global torques associated with the non-symmetric mass distribution is thus not advisable even within an idealized, smooth galaxy model; we further note that spiral patterns are unlikely to affect the inspiral due to their transient and fluctuating nature. We speculate that the sinking efficiency of massive black holes involved in minor galaxy mergers may be hampered in barred galaxies, making them less likely to host a gravitational wave signal accessible to low-frequency detectors.

0 citations0 reviewsNo code linkedOpen access
preprint2021arXiv

Dynamical evolution of massive perturbers in realistic multi-component galaxy models I: implementation and validation

Galaxies are self-gravitating structures composed by several components encompassing spherical, axial and triaxial symmetry. Although real systems feature heterogeneous components whose properties are intimately connected, semi-analytical approaches often exploit the linearity of the Poisson's equation to represent the potential and mass distribution of a multi-component galaxy as the sum of the individual components. In this work, we expand the semi-analytical framework developed in Bonetti et al. (2020) by including both a detailed implementation of the gravitational potential of exponential disc (modelled with a ${\rm sech}^2$ and an exponential vertical profile) and an accurate prescription for the dynamical friction experienced by massive perturbers in composite galaxy models featuring rotating disc structures. Such improvements allow us to evolve arbitrary orbits either within or outside the galactic disc plane. We validate the results obtained by our numerical model against public semi-analytical codes as well as full N-body simulations, finding that our model is in excellent agreement to the codes it is compared with. The ability to reproduce the relevant physical processes responsible for the evolution of massive perturber orbits and its computational efficiency make our framework perfectly suited for large parameter-space exploration studies.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

A new proxy to estimate the cosmic-ray ionisation rate in dense cores

Cosmic rays are a global source of ionisation, and the ionisation fraction represents a fundamental parameter in the interstellar medium. Ions couple to magnetic fields, affect the chemistry, and the dynamics of star-forming regions as well as planetary atmospheres. However, the cosmic-ray ionisation rate represents one of the bottlenecks for astrochemical models, and its determination is one of the most puzzling problems in astrophysics. While for diffuse clouds reasonable values have been provided from H$_3^+$ observations, for dense clouds, due to the lack of rotational transitions, this is not possible, and estimates are strongly biased by the employed model. We present here an analytical expression, obtained from first principles, to estimate the cosmic-ray ionisation rate from observational quantities. The theoretical predictions are validated with high-resolution three-dimensional numerical simulations and applied to the well known core L1544; we obtained an estimate of $ζ_2 \sim 2-3 \times 10^{-17}$ s$^{-1}$. Our results and the analytical formulae provided represent the first model-independent, robust tool to probe the cosmic-ray ionisation rate in the densest part of star-forming regions (on spatial scales of $R \leq 0.05$ pc). An error analysis is presented to give statistical relevance to our study.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Difficulties in Mid-Infrared selection of AGN in dwarf galaxies

While massive black holes (MBHs) are known to inhabit all massive galaxies, their ubiquitous presence in dwarf galaxies has not been confirmed yet, with only a limited number of sources detected so far. Recently, some studies proposed infrared emission as an alternative way to identify MBHs in dwarfs, based on a similar approach usually applied to quasars. In this study, by accurately combining optical and infrared data taking into account resolution effects and source overlapping, we investigate in detail the possible limitations of this approach with current ground-based facilities, finding a quite low ($\sim$0.4 per cent) fraction of active MBH in dwarfs that are luminous in mid-infrared, consistent with several previous results. Our results suggest that the infrared selection is strongly affected by several limitations that make the identification of MBHs in dwarf galaxies currently prohibitive, especially because of the very poor resolution compared to optical surveys, and the likely contamination by nearby sources, although we find a few good candidates worth further follow-ups. Optical, X-ray and radio observations, therefore, still represent the most secure way to search for MBH in dwarfs.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Dynamical friction-driven orbital circularisation in rotating discs: a semi-analytical description

We present and validate a novel semi-analytical approach to study the effect of dynamical friction on the orbits of massive perturbers in rotating stellar discs. We find that dynamical friction efficiently circularises the orbit of co-rotating perturbers, while it constantly increases the eccentricity of counter-rotating ones until their angular momenta reverse, then once again promoting circularisation. Such "drag toward circular corotation" could shape the distribution of orientations of kinematically decoupled cores in disc galaxies, naturally leading to the observed larger fraction of co-rotating cores.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Globular Cluster Formation from Colliding Substructure

We investigate a scenario where the formation of Globular Clusters (GCs) is triggered by high-speed collisions between infalling atomic-cooling subhalos during the assembly of the main galaxy host, a special dynamical mode of star formation that operates at high gas pressures and is intimately tied to LCDM hierarchical galaxy assembly. The proposed mechanism would give origin to "naked" globulars, as colliding dark matter subhalos and their stars will simply pass through one another while the warm gas within them clashes at highly supersonic speed and decouples from the collisionless component, in a process reminiscent of the Bullet galaxy cluster. We find that the resulting shock-compressed layer cools on a timescale that is typically shorter than the crossing time, first by atomic line emission and then via fine-structure metal-line emission, and is subject to gravitational instability and fragmentation. Through a combination of kinetic theory approximation and high-resolution $N$-body simulations, we show that this model may produce: (a) a GC number-halo mass relation that is linear down to dwarf galaxy scales and agrees with the trend observed over five orders of magnitude in galaxy mass; (b) a population of old globulars with a median age of 12 Gyr and an age spread similar to that observed; (c) a spatial distribution that is biased relative to the overall mass profile of the host; and (d) a bimodal metallicity distribution with a spread similar to that observed in massive galaxies.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

The [CII]--SFR correlation in dwarf galaxies across cosmic time

Current galaxy observations suggest that a roughly linear correlation exists between the [CII] emission and the star formation rate, either as spatially-resolved or integrated quantities. Observationally, this correlation seems to be independent of metallicity, but the very large scatter does not allow to properly assess whether this is true. On the other hand, theoretical models tend to suggest a metallicity dependence of the correlation. In this study, we investigate the metallicity evolution of the correlation via a high-resolution zoom-in cosmological simulation of a dwarf galaxy employing state-of-the-art sub-grid modelling for gas cooling, star formation, and stellar feedback, and that self-consistently evolves the abundances of metal elements out of equilibrium. Our results suggest that the correlation should evolve with metallicity, in agreement with theoretical predictions, but also that this evolution can be hardly detected in observations, because of the large scatter. We also find that most of the [CII] emission is associated with neutral gas at low-intermediate densities, whereas the highest emissivity is produced by the densest regions around star-forming regions.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

The AGORA high-resolution galaxy simulations comparison project: Public data release

As part of the AGORA High-resolution Galaxy Simulations Comparison Project (Kim et al. 2014, 2016) we have generated a suite of isolated Milky Way-mass galaxy simulations using 9 state-of-the-art gravito-hydrodynamics codes widely used in the numerical galaxy formation community. In these simulations we adopted identical galactic disk initial conditions, and common physics models (e.g., radiative cooling and ultraviolet background by a standardized package). Subgrid physics models such as Jeans pressure floor, star formation, supernova feedback energy, and metal production were carefully constrained. Here we release the simulation data to be freely used by the community. In this release we include the disk snapshots at 0 and 500Myr of evolution per each code as used in Kim et al. (2016), from simulations with and without star formation and feedback. We encourage any member of the numerical galaxy formation community to make use of these resources for their research - for example, compare their own simulations with the AGORA galaxies, with the common analysis yt scripts used to obtain the plots shown in our papers, also available in this release.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Clumpy high-z galaxies as a testbed for feedback-regulated galaxy formation

We study the dependence of fragmentation in massive gas-rich galaxy disks at $z > 1$ on feedback model and hydrodynamical method, employing the GASOLINE2 SPH code and the lagrangian mesh-less code GIZMO in finite mass mode. We compare non-cosmological galaxy disk runs with standard blastwave supernovae (SN)feedback, which introduces delayed cooling in order to drive winds, and runs with the new superbubble SN feedback, which produces winds naturally by modelling the detailed physics of SN-driven bubbles and leads to efficient self-regulation of star formation. We find that, with blastwave feedback, massive star forming clumps form in comparable number and with very similar masses in GASOLINE2 and GIZMO. The typical masses are in the range $10^7-10^8 M_{\odot}$, lower than in most previous works, while giant clumps with masses above $10^9 M_{\odot}$ are exceedingly rare. With superbubble feedback, instead, massive bound star forming clumps do not form because galaxies never undergo a phase of violent disk instability. Only sporadic, unbound star forming overdensities lasting only a few tens of Myr can arise that are triggered by perturbations of massive satellite companions. We conclude that there is a severe tension between explaining massive star forming clumps observed at $z > 1$ primarily as the result of disk fragmentation driven by gravitational instability and the prevailing view of feedback-regulated galaxy formation. The link between disk stability and star formation efficiency should thus be regarded as a key testing ground for galaxy formation theory.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Hydrodynamical simulations of the tidal stripping of binary stars by massive black holes

In a galactic nucleus, a star on a low angular momentum orbit around the central massive black hole can be fully or partially disrupted by the black hole tidal field, lighting up the compact object via gas accretion. This phenomenon can repeat if the star, not fully disrupted, is on a closed orbit. Because of the multiplicity of stars in binary systems, also binary stars may experience in pairs such a fate, immediately after being tidally separated. The consumption of both the binary components by the black hole is expected to power a double-peaked flare. In this paper, we perform for the first time, with GADGET2, a suite of smoothed particle hydrodynamics simulations of binary stars around a galactic central black hole in the Newtonian regime. We show that accretion luminosity light curves from double tidal disruptions reveal a more prominent knee, rather than a double peak, when decreasing the impact parameter of the encounter and when elevating the difference between the mass of the star which leaves the system after binary separation and the mass of the companion. The detection of a knee can anticipate the onset of periodic accretion luminosity flares if one of the stars, only partially disrupted, remains bound to the black hole after binary separation. Thus knees could be precursors of periodic flares, which can then be predicted, followed up and better modelled. Analytical estimates in the black hole mass range $10^5-10^8 \rm M_{\rm \odot}$ show that the knee signature is enhanced in the case of black holes of mass $10^6-10^7 \rm M_{\rm \odot}$.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Massive black hole and gas dynamics in mergers of galaxy nuclei - II. Black hole sinking in star-forming nuclear discs

Mergers of gas-rich galaxies are key events in the hierarchical built-up of cosmic structures, and can lead to the formation of massive black hole binaries. By means of high-resolution hydrodynamical simulations we consider the late stages of a gas-rich major merger, detailing the dynamics of two circumnuclear discs, and of the hosted massive black holes during their pairing phase. During the merger gas clumps with masses of a fraction of the black hole mass form because of fragmentation. Such high-density gas is very effective in forming stars, and the most massive clumps can substantially perturb the black hole orbits. After $\sim 10$ Myr from the start of the merger a gravitationally bound black hole binary forms at a separation of a few parsecs, and soon after, the separation falls below our resolution limit of $0.39$ pc. At the time of binary formation the original discs are almost completely disrupted because of SNa feedback, while on pc scales the residual gas settles in a circumbinary disc with mass $\sim 10^5 M_\odot$. We also test that binary dynamics is robust against the details of the SNa feedback employed in the simulations, while gas dynamics is not. We finally highlight the importance of the SNa time-scale on our results.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Constraining the high redshift formation of black hole seeds in nuclear star clusters with gas inflows

In this paper we explore a possible route of black hole seed formation that appeal to a model by Davies, Miller & Bellovary who considered the case of the dynamical collapse of a dense cluster of stellar black holes subjected to an inflow of gas. Here, we explore this case in a broad cosmological context. The working hypotheses are that (i) nuclear star clusters form at high redshifts in pre-galactic discs hosted in dark matter halos, providing a suitable environment for the formation of stellar black holes in their cores, (ii) major central inflows of gas occur onto these clusters due to instabilities seeded in the growing discs and/or to mergers with other gas-rich halos, and that (iii) following the inflow, stellar black holes in the core avoid ejection due to the steepening to the potential well, leading to core collapse and the formation of a massive seed of $<~ 1000\, \rm M_\odot$. We simulate a cosmological box tracing the build up of the dark matter halos and there embedded baryons, and explore cluster evolution with a semi-analytical model. We show that this route is feasible, peaks at redshifts $z <~ 10$ and occurs in concomitance with the formation of seeds from other channels. The channel is competitive relative to others, and is independent of the metal content of the parent cluster. This mechanism of gas driven core collapse requires inflows with masses at least ten times larger than the mass of the parent star cluster, occurring on timescales shorter than the evaporation/ejection time of the stellar black holes from the core. In this respect, the results provide upper limit to the frequency of this process.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Massive black hole and gas dynamics in galaxy nuclei mergers. I. Numerical implementation

Numerical effects are known to plague adaptive mesh refinement (AMR) codes when treating massive particles, e.g. representing massive black holes (MBHs). In an evolving background, they can experience strong, spurious perturbations and then follow unphysical orbits. We study by means of numerical simulations the dynamical evolution of a pair MBHs in the rapidly and violently evolving gaseous and stellar background that follows a galaxy major merger. We confirm that spurious numerical effects alter the MBH orbits in AMR simulations, and show that numerical issues are ultimately due to a drop in the spatial resolution during the simulation, drastically reducing the accuracy in the gravitational force computation. We therefore propose a new refinement criterion suited for massive particles, able to solve in a fast and precise way for their orbits in highly dynamical backgrounds. The new refinement criterion we designed enforces the region around each massive particle to remain at the maximum resolution allowed, independently upon the local gas density. Such maximally-resolved regions then follow the MBHs along their orbits, and effectively avoids all spurious effects caused by resolution changes. Our suite of high resolution, adaptive mesh-refinement hydrodynamic simulations, including different prescriptions for the sub-grid gas physics, shows that the new refinement implementation has the advantage of not altering the physical evolution of the MBHs, accounting for all the non trivial physical processes taking place in violent dynamical scenarios, such as the final stages of a galaxy major merger.

0 citations0 reviewsNo code linkedOpen access