Showing 1–50 of 139 results for author: Dotti, M

An effective model for the tidal disruption of satellites undergoing minor mergers with axisymmetric primaries

Authors: Ludovica Varisco, Massimo Dotti, Matteo Bonetti, Elisa Bortolas, Alessandro Lupi

Abstract: According to the hierarchical formation paradigm, galaxies form through mergers of smaller entities and massive black holes (MBHs), if lurking at their centers, migrate to the nucleus of the newly formed galaxy, where they form binary systems. The formation and evolution of MBH binaries, and in particular their coalescence timescale, is very relevant for current and future facilities aimed at dete… ▽ More According to the hierarchical formation paradigm, galaxies form through mergers of smaller entities and massive black holes (MBHs), if lurking at their centers, migrate to the nucleus of the newly formed galaxy, where they form binary systems. The formation and evolution of MBH binaries, and in particular their coalescence timescale, is very relevant for current and future facilities aimed at detecting the gravitational-wave signal produced by the MBH close to coalescence. While most of the studies targeting this process are based on hydrodynamic simulations, the high computational cost makes a complete parameter space exploration prohibitive. Semi-analytic approaches represent a valid alternative, but they require ad-hoc prescriptions for the mass loss of the merging galaxies in minor mergers due to tidal stripping, which is not commonly considered or at most modelled assuming very idealised geometries. In this work, we propose a novel, effective model for the tidal stripping in axisymmetric potentials, to be implemented in semi-analytic models. We validate our semi-analytic approach against N-body simulations considering different galaxy sizes, inclinations, and eccentricities, finding only a moderate dependence on the orbit eccentricity. In particular, we find that, for almost circular orbits, our model mildly overestimates the mass loss, and this is due to the adjustment of the stellar distribution after the mass is removed. Nonetheless, the model exhibits a very good agreement with simulations in all the considered conditions, and thus represents an extremely powerful addition to semi-analytic calculations. △ Less

Submitted 26 June, 2024; originally announced June 2024.

Comments: 15 pages, 8 figures, 2 tables. Accepted for publication on A&A

Size matters: are we witnessing super-Eddington accretion in high-redshift black holes from JWST?

Authors: Alessandro Lupi, Alessandro Trinca, Marta Volonteri, Massimo Dotti, Chiara Mazzucchelli

Abstract: Observations by the James Webb Space Telescope of the Universe at $z\gtrsim 4$ have shown that massive black holes (MBHs) appear extremely overmassive compared to the local correlation for active galactic nuclei. In some cases, these objects might even reach half the stellar mass inferred for the galaxy. Understanding how such objects formed and grew to this masses has then become a big challenge… ▽ More Observations by the James Webb Space Telescope of the Universe at $z\gtrsim 4$ have shown that massive black holes (MBHs) appear extremely overmassive compared to the local correlation for active galactic nuclei. In some cases, these objects might even reach half the stellar mass inferred for the galaxy. Understanding how such objects formed and grew to this masses has then become a big challenge for theoretical models, with different ideas ranging from heavy seed to super-Eddington accretion phases. Here, we take a different approach, and try to infer how accurate these MBH mass estimates are and whether we really need to revise our physical models. By considering how the emerging spectrum (both the continuum and the broad lines) of an accreting MBH changes close to and above the Eddington limit, we infer a much larger uncertainty in the MBH mass estimates relative to that of local counterparts, up to an order of magnitude, and a potential preference for lower masses and higher accretion rates, which i) move them closer to the local correlations, and ii) might indicate that we are witnessing for the first time a widespread phase of very rapid accretion. △ Less

Submitted 25 June, 2024; originally announced June 2024.

Comments: 8 pages, 5 figures, submitted to A&A

Ticking away: the long-term X-ray timing and spectral evolution of eRO-QPE2

Authors: R. Arcodia, I. Linial, G. Miniutti, A. Franchini, M. Giustini, M. Bonetti, A. Sesana, R. Soria, J. Chakraborty, M. Dotti, E. Kara, A. Merloni, G. Ponti, F. Vincentelli

Abstract: Quasi-periodic eruptions (QPEs) are repeated X-ray flares from galactic nuclei. Despite some diversity in the recurrence and amplitude of eruptions, their striking regularity has motivated theorists to associate QPEs with orbital systems. Among the known QPE sources, eRO-QPE2 has shown the most regular flare timing and luminosity since its discovery. We report here on its long-term evolution over… ▽ More Quasi-periodic eruptions (QPEs) are repeated X-ray flares from galactic nuclei. Despite some diversity in the recurrence and amplitude of eruptions, their striking regularity has motivated theorists to associate QPEs with orbital systems. Among the known QPE sources, eRO-QPE2 has shown the most regular flare timing and luminosity since its discovery. We report here on its long-term evolution over $\sim3.3\,$yr from discovery and find that: i) the average QPE recurrence time per epoch has decreased over time, albeit not at a uniform rate; ii) the distinct alternation between consecutive long and short recurrence times found at discovery has not been significant since; iii) the spectral properties, namely flux and temperature of both eruptions and quiescence components, have remained remarkably consistent within uncertainties. We attempted to interpret these results as orbital period and eccentricity decay coupled with orbital and disk precession. However, since gaps between observations are too long, we are not able to distinguish between an evolution dominated by just a decreasing trend, or by large modulations (e.g. due to the precession frequencies at play). In the former case, the observed period decrease is roughly consistent with that of a star losing orbital energy due to hydrodynamic gas drag from disk collisions, although the related eccentricity decay is too fast and additional modulations have to contribute too. In the latter case, no conclusive remarks are possible on the orbital evolution and the nature of the orbiter due to the many effects at play. However, these two cases come with distinctive predictions for future X-ray data: in the former, we expect all future observations to show a shorter recurrence time than the latest epoch, while in the latter we expect some future observations to be found with a larger recurrence, hence an apparent temporary period increase. △ Less

Submitted 24 June, 2024; originally announced June 2024.

Comments: Recently submitted to A&A. Comments welcome

A quasar-galaxy merger at $z\sim 6.2$: rapid host growth via accretion of two massive satellite galaxies

Authors: Roberto Decarli, Federica Loiacono, Emanuele Paolo Farina, Massimo Dotti, Alessandro Lupi, Romain A. Meyer, Marco Mignoli, Antonio Pensabene, Michael A. Strauss, Bram Venemans, Jinyi Yang, Fabian Walter, Julien Wolf, Eduardo Bañados, Laura Blecha, Sarah Bosman, Chris L. Carilli, Andrea Comastri, Thomas Connor, Tiago Costa, Anna-Christina Eilers, Xiaohui Fan, Roberto Gilli, Hyunsung D. Jun, Weizhe Liu , et al. (16 additional authors not shown)

Abstract: We present JWST/NIRSpec Integral Field Spectroscopy in the rest-frame optical bands of the system PJ308-21, a quasar at $z=6.2342$ caught as its host galaxy interacts with companion galaxies. We detect spatially extended emission of several emission lines (H$α$, H$β$, [OIII], [NII], [SII], HeII), which we use to study the properties of the ionized phase of the interstellar medium: the source and h… ▽ More We present JWST/NIRSpec Integral Field Spectroscopy in the rest-frame optical bands of the system PJ308-21, a quasar at $z=6.2342$ caught as its host galaxy interacts with companion galaxies. We detect spatially extended emission of several emission lines (H$α$, H$β$, [OIII], [NII], [SII], HeII), which we use to study the properties of the ionized phase of the interstellar medium: the source and hardness of the photoionizing radiation field, metallicity, dust reddening, electron density and temperature, and star formation. We also marginally detect continuum starlight emission associated with the companion sources. We find that at least two independent satellite galaxies are part of the system. While the quasar host appears highly enriched and obscured, with AGN-like photoionization conditions, the western companion shows minimal dust extinction, low metallicity ($Z\sim0.4$ Z$_\odot$), and star-formation driven photoionization. The eastern companion shows higher extinction and metallicity ($Z\sim0.8$ Z$_\odot$) compared to the western companion, and it is at least partially photoionized by the nearby quasar. We do not find any indication of AGN in the companion sources. Our study shows that while the quasar host galaxy is already very massive ($M_{\rm dyn}>10^{11}$ M$_\odot$), it is still rapidly building up by accreting two relatively massive ($M_{\rm star}\sim 10^{10}$ M$_\odot$) companion sources. This dataset showcases the power of JWST in exposing the build-up of massive galaxies in the first Gyr of the Universe. △ Less

Submitted 10 June, 2024; originally announced June 2024.

Comments: 15 pages, 16 figures. Accepted for publication in A&A

Repeating partial disruptions and two-body relaxation

Authors: Luca Broggi, Nicholas C. Stone, Taeho Ryu, Elisa Bortolas, Massimo Dotti, Matteo Bonetti, Alberto Sesana

Abstract: Two-body relaxation may drive stars onto near-radial orbits around a massive black hole, resulting in a tidal disruption event (TDE). In some circumstances, stars are unlikely to undergo a single terminal disruption, but rather to have a sequence of many grazing encounters with the black hole. It has long been unclear what is the physical outcome of this sequence: each of these encounters can only… ▽ More Two-body relaxation may drive stars onto near-radial orbits around a massive black hole, resulting in a tidal disruption event (TDE). In some circumstances, stars are unlikely to undergo a single terminal disruption, but rather to have a sequence of many grazing encounters with the black hole. It has long been unclear what is the physical outcome of this sequence: each of these encounters can only liberate a small amount of stellar mass, but may significantly alter the orbit of the star. We study the phenomenon of repeating partial tidal disruptions (pTDEs) by building a semi-analytical model that accounts for mass loss and tidal excitation. In the empty loss cone regime, where two-body relaxation is weak, we estimate the number of consecutive partial disruptions that a star can undergo, on average, before being significantly affected by two-body encounters. We find that in this empty loss cone regime, a star will be destroyed in a sequence of weak pTDEs, possibly explaining the tension between the low observed TDE rate and its higher theoretical estimates. △ Less

Submitted 19 June, 2024; v1 submitted 8 April, 2024; originally announced April 2024.

Comments: Open Journal of Astrophysics, accepted for publication

Partial alignment between jets and megamasers: coherent or selective accretion?

Authors: Massimo Dotti, Riccardo Buscicchio, Francesco Bollati, Roberto Decarli, Walter Del Pozzo, Alessia Franchini

Abstract: Spins play a crucial role in the appearance, evolution and occupation fraction of massive black holes. To date observational estimates of massive black hole spins are scarce, and the assumptions commonly made in such estimates have been recently questioned. Similarly, theoretical models for massive black hole spin evolution, while reproducing the few observational constraints, are based on possibl… ▽ More Spins play a crucial role in the appearance, evolution and occupation fraction of massive black holes. To date observational estimates of massive black hole spins are scarce, and the assumptions commonly made in such estimates have been recently questioned. Similarly, theoretical models for massive black hole spin evolution, while reproducing the few observational constraints, are based on possibly oversimplified assumptions. New independent constraints on massive black hole spins are therefore of primary importance. We present a rigorous statistical analysis of the relative orientation of radio jets and mega-maser discs in ten low-redshift galaxies. We find a strong preference for (partial) alignment between jets and mega-masers, that can be attributed to two different causes: "coherent" and "selective" accretion. In the first case the partial alignment is due to an anisotropy in the gas reservoir fueling the growth of MBHs. In the second case the spin-dependent anisotropic feedback allows for long-lived accretion only if the orbits of the gas inflows are almost aligned to the MBH equatorial plane. A discussion of the implications of the two accretion scenarios on the evolution of MBHs is presented, together with an outlook on future observational tests aiming at discriminating between the two scenarios and at checking if any of the two apply to different redshifts and black hole mass regimes. △ Less

Submitted 26 March, 2024; originally announced March 2024.

Comments: 8 pages, 3 figures, 2 tables

BANG-MaNGA: A census of kinematic discs and bulges across mass and star formation in the local Universe

Authors: Fabio Rigamonti, Luca Cortese, Francesco Bollati, Stefano Covino, Massimo Dotti, A. Fraser-McKelvie, Francesco Haardt

Abstract: We investigate the relevance of kinematically identified bulges, discs and their role relative to galaxy quenching. We utilize an analysis of the SDSS-MaNGA survey conducted with the GPU-based code BANG which simultaneously models galaxy photometry and kinematics to decompose galaxies into their structural components. Below M~1011 Msun, galaxies exhibit a wide range of dynamical properties, determ… ▽ More We investigate the relevance of kinematically identified bulges, discs and their role relative to galaxy quenching. We utilize an analysis of the SDSS-MaNGA survey conducted with the GPU-based code BANG which simultaneously models galaxy photometry and kinematics to decompose galaxies into their structural components. Below M~1011 Msun, galaxies exhibit a wide range of dynamical properties, determined by the relative prominence of a dispersion-supported inner region and a rotationally-supported disc. Our analysis reveals a natural separation between these classes, with only a minor fraction of stellar mass retained by structures exhibiting intermediate dynamical support. When examining galaxies in terms of their star formation activity, an apparent decrease in rotational support is observed as they move below the star-forming main sequence. This behaviour is evident with luminosity-weighted tracers of kinematics, while it almost vanishes with mass-weighted tracers. Luminosity-weighted quantities not only capture differences in kinematics but also in the stellar population, potentially leading to biased interpretations of galaxy dynamical properties and quenching. Our findings suggest that quenching does not imply almost any structural transformation in galaxies below M~10^11 Msun. Processes as disc fading more likely account for observed differences in mass-weighted and luminosity-weighted galaxy properties; when the galactic disc ceases star formation, its mass-to-light ratio grows without any significant morphological transformation. The picture is remarkably different above M~10^11 Msun. Regardless of the tracer used, a substantial increase in galaxy dispersion support is observed along with a significant structural change. A different quenching mechanism, most likely associated with mergers, dominates. Notably, this mechanism is confined to a very limited range of high masses. △ Less

Submitted 14 March, 2024; originally announced March 2024.

Comments: Accepted for publication in A&A. 16 pages, 10 figures

A quasar-galaxy merger at $z\sim6.2$: black hole mass and quasar properties from the NIRSpec spectrum

Authors: Federica Loiacono, Roberto Decarli, Marco Mignoli, Emanuele Paolo Farina, Eduardo Bañados, Sarah Bosman, Anna-Christina Eilers, Jan-Torge Schindler, Michael A. Strauss, Marianne Vestergaard, Feige Wang, Laura Blecha, Chris L. Carilli, Andrea Comastri, Thomas Connor, Tiago Costa, Massimo Dotti, Xiaohui Fan, Roberto Gilli, Hyunsung D. Jun, Weizhe Liu, Alessandro Lupi, Madeline A. Marshall, Chiara Mazzucchelli, Romain A. Meyer , et al. (9 additional authors not shown)

Abstract: We present JWST/NIRSpec integral field data of the quasar PJ308-21 at $z=6.2342$. As shown by previous ALMA and HST imaging, the quasar has two companion sources, interacting with the quasar host galaxy. The high-resolution G395H/290LP NIRSpec spectrum covers the $2.87-5.27\ \rm μm$ wavelength range and shows the rest-frame optical emission of the quasar with exquisite quality ($S/N\sim 100-400$ p… ▽ More We present JWST/NIRSpec integral field data of the quasar PJ308-21 at $z=6.2342$. As shown by previous ALMA and HST imaging, the quasar has two companion sources, interacting with the quasar host galaxy. The high-resolution G395H/290LP NIRSpec spectrum covers the $2.87-5.27\ \rm μm$ wavelength range and shows the rest-frame optical emission of the quasar with exquisite quality ($S/N\sim 100-400$ per spectral element). Based on the H$β$ line from the broad line region, we obtain an estimate of the black hole mass $M_{\rm BH,Hβ}\sim 2.7\times 10^{9}\ \rm M_{\odot}$. This value is within a factor $\lesssim 1.5$ of the H$α$-based black hole mass from the same spectrum ($M_{\rm BH, Hα}\sim 1.93\times 10^{9}\ \rm M_{\odot}$) and is consistent with a previous estimate relying on the MgII $λ2799$ ($M_{\rm BH, MgII}\sim 2.65\times 10^{9}\ \rm M_{\odot}$). All these $M_{\rm BH}$ are within the $\sim 0.5$ dex intrinsic scatter of the adopted mass calibrations. The high Eddington ratio of PJ308-21 $λ_{\rm Edd,Hβ}\sim 0.67$ ($λ_{\rm Edd,Hα}\sim 0.96$) is in line with the overall quasar population at $z \gtrsim 6$. The relative strengths of the [OIII], FeII and H$β$ lines are consistent with the empirical "Eigenvector 1" correlations as observed for low redshift quasars. We find evidence for blueshifted [OIII] $λ5007$ emission with a velocity offset $Δv_{\rm [OIII]}=-1922\pm 39$ km s$^{-1}$ from the systemic velocity and a $\rm FWHM([OIII])=2776^{+75}_{-74}$ km s$^{-1}$. This may be the signature of an outflow from the nuclear region, despite the true values of $Δv_{\rm [OIII]}$ and $\rm FWHM([OIII])$ are likely more uncertain due to the blending with H$β$ and FeII lines. Our study demonstrates the unique capabilities of NIRSpec in capturing quasar spectra at cosmic dawn and studying their properties in unprecedented detail. △ Less

Submitted 20 February, 2024; originally announced February 2024.

Comments: 17 pages, 9 figures, accepted for publication in A&A

Journal ref: A&A 685, A121 (2024)

Astrophysical and relativistic modeling of the recoiling black-hole candidate in quasar 3C 186

Authors: Matteo Boschini, Davide Gerosa, Om Sharan Salafia, Massimo Dotti

Abstract: The compact object in quasar 3C 186 is one of the most promising recoiling black-hole candidates, exhibiting both an astrometric displacement between the quasar and the host galaxy as well as a spectroscopic shift between broad and narrow lines. 3C 186 also presents a radio jet which, when projected onto the plane of the sky, appears to be perpendicular to the quasar/galaxy displacement. Assuming… ▽ More The compact object in quasar 3C 186 is one of the most promising recoiling black-hole candidates, exhibiting both an astrometric displacement between the quasar and the host galaxy as well as a spectroscopic shift between broad and narrow lines. 3C 186 also presents a radio jet which, when projected onto the plane of the sky, appears to be perpendicular to the quasar/galaxy displacement. Assuming a gravitational-wave kick is indeed responsible for the properties of 3C 186 and using state-of-the-art relativistic modeling, we show that current observations allow for exquisite modeling of the recoiling black hole. Most notably, we find that the kick velocity, the black-hole spin, and the line of sight are almost collinear and the former appear perpendicular to each other only because of a strong projection effect. The targeted configuration requires substantial fine-tuning: while there exists a region in the black-hole binary parameter space that is compatible with 3C 186, the observed system appears to be a rare occurrence. Using archival radio observations, we explore different strategies that could potentially confirm or rule out our interpretation. In particular, we develop two observational tests that rely on the brightness ratio between the approaching and receding jet as well as the asymmetry of the jet lobes. While the available radio data provide loose constraints, deeper observations have the unique potential of unveiling the nature of 3C 186. △ Less

Submitted 21 June, 2024; v1 submitted 13 February, 2024; originally announced February 2024.

Comments: 12 pages, 8 figures

Journal ref: A&A 686, A245 (2024)

On the connection between AGN radiative feedback and massive black hole spin

Authors: Francesco Bollati, Alessandro Lupi, Massimo Dotti, Francesco Haardt

Abstract: We present a novel implementation for active galactic nucleus (AGN) feedback through ultra-fast winds in the code gizmo. Our feedback recipe accounts for the angular dependence of radiative feedback upon black hole spin. We self-consistently evolve in time i) the gas accretion process from resolved scales to an unresolved AGN disc, ii) the evolution of the spin of the massive black hole (MBH), iii… ▽ More We present a novel implementation for active galactic nucleus (AGN) feedback through ultra-fast winds in the code gizmo. Our feedback recipe accounts for the angular dependence of radiative feedback upon black hole spin. We self-consistently evolve in time i) the gas accretion process from resolved scales to an unresolved AGN disc, ii) the evolution of the spin of the massive black hole (MBH), iii) the injection of AGN-driven winds into the resolved scales, and iv) the spin-induced anisotropy of the overall feedback process. We test our implementation by following the propagation of the wind-driven outflow into a homogeneous medium, and we compare the results against simple analytical models. Then, we consider an isolated galaxy setup and there we study the impact of the AGN feedback on the evolution of the MBH and the of the host galaxy. We find that: i) AGN feedback limits the gas inflow that powers the MBH, with a consequent weak impact on the host galaxy characterized by a star formation (SF) suppression of about a factor of two in the nuclear region; ii) the impact of AGN feedback on the host galaxy and on MBH growth is primarily determined by the AGN luminosity, rather than by its angular pattern set by the MBH spin; iii) the imprint of the angular pattern of the AGN radiation emission manifest in a more clear way at high accretion rates. At such high rates the more isotropic angular patterns, proper to higher spin values, sweep away gas in the nuclear region more easily, hence causing a slower MBH mass and spin growths and a higher quenching of the SF. We argue that the influence of spin-dependent anisotropy of AGN feedback on MBH and galaxy evolution is likely to be relevant in those scenarios characterized by high and prolonged MBH accretion episodes and by high AGN wind-galaxy coupling. Such conditions are more frequently met in galaxy mergers and/or high redshift galaxies. △ Less

Submitted 13 November, 2023; originally announced November 2023.

A fast test for the identification and confirmation of massive black hole binaries

Authors: Massimo Dotti, Fabio Rigamonti, Stefano Rinaldi, Walter Del Pozzo, Roberto Decarli, Riccardo Buscicchio

Abstract: We present a new observational test to identify massive black hole binaries in large multi-epoch spectroscopical catalogues and to probe the real nature of already proposed binary candidates. The test is tailored for binaries with separations large enough to allow each black hole to retain its own broad line region. In this limit the fast AGN variability typically observed over months cannot be as… ▽ More We present a new observational test to identify massive black hole binaries in large multi-epoch spectroscopical catalogues and to probe the real nature of already proposed binary candidates. The test is tailored for binaries with separations large enough to allow each black hole to retain its own broad line region. In this limit the fast AGN variability typically observed over months cannot be associated to the much longer binary period, and it is assumed (as for the case of single black holes) to be the consequence of the evolution of the innermost regions of the two accretion discs. A simple analysis of the cross-correlation between different parts of individual broad emission lines can therefore identify the presence of two massive black holes whose continua vary independently of each other. Our analysis indicates that to be less affected by the noise in the spectra the broad lines should be divided in two close-to-equal-flux parts. This ensures that in the single massive black hole scenario the cross-correlation will always be high. With monitoring campaigns similar to those performed for reverberation mapping studies, on the other way, a binary can show any value of the cross-correlation and can therefore be distinguished from a standard AGN. The new test can be performed over timescales orders of magnitude shorter than the alternative tests already discussed in literature, and can be a powerful complement to the massive black hole binary search strategies already in place. △ Less

Submitted 9 January, 2024; v1 submitted 10 October, 2023; originally announced October 2023.

Comments: 9 pages, 6 figures, accepted for publication on Astronomy & Astrophysics

Journal ref: A&A 680, A69 (2023)

The fragility of thin discs in galaxies -- II. Thin discs as tracers of the assembly history of galaxies

Authors: Pablo M. Galán-de Anta, Pedro R. Capelo, Eugene Vasiliev, Massimo Dotti, Marc Sarzi, Enrico Maria Corsini, Lorenzo Morelli

Abstract: Thin galactic discs and nuclear stellar discs (NSDs) are fragile structures that can be easily disturbed by merger events. By studying the age of the stellar populations in present-day discs, we can learn about the assembly history of galaxies and place constraints on their past merger events. Following on the steps of our initial work, we explore the fragility of such disc structures in intermedi… ▽ More Thin galactic discs and nuclear stellar discs (NSDs) are fragile structures that can be easily disturbed by merger events. By studying the age of the stellar populations in present-day discs, we can learn about the assembly history of galaxies and place constraints on their past merger events. Following on the steps of our initial work, we explore the fragility of such disc structures in intermediate-mass-ratio dry encounters using the previously constructed $N$-body model of the Fornax galaxy NGC 1381 (FCC 170), which hosts both a thin galactic disc and a NSD. We dismiss major and minor encounters, as the former were previously shown to easily destroy thin-disc structures, whereas the latter take several Hubble times to complete in the specific case of FCC 170. The kinematics and structure of the thin galactic disc are dramatically altered by the mergers, whereas the NSD shows a remarkable resilience, exhibiting only a smooth increase of its size when compared to the model evolved in isolation. Our results suggest that thin galactic discs are better tracers for intermediate-mass-ratio mergers, while NSDs may be more useful for major encounters. Based on our simulations and previous analysis of the stellar populations, we concluded that FCC 170 has not experienced any intermediate-mass-ratio dry encounters for at least $\sim$10 Gyr, as indicated by the age of its thin-disc stellar populations. △ Less

Submitted 22 May, 2023; originally announced May 2023.

Comments: 11 pages, 10 figures

Decomposing galaxies with BANG: an automated morpho-kinematical decomposition of the SDSS-DR17 MaNGA survey

Authors: Fabio Rigamonti, Massimo Dotti, Stefano Covino, Francesco Haardt, Luca Cortese, Marco Landoni, Ludovica Varisco

Abstract: From a purely photometric perspective galaxies are generally decomposed into a bulge+disc system, with bulges being dispersion-dominated and discs rotationally-supported. However, recent observations have demonstrated that such a framework oversimplifies complexity, especially if one considers galaxy kinematics. To address this issue we introduced with the GPU-based code \textsc{bang} a novel appr… ▽ More From a purely photometric perspective galaxies are generally decomposed into a bulge+disc system, with bulges being dispersion-dominated and discs rotationally-supported. However, recent observations have demonstrated that such a framework oversimplifies complexity, especially if one considers galaxy kinematics. To address this issue we introduced with the GPU-based code \textsc{bang} a novel approach that employs analytical potential-density pairs as galactic components, allowing for a computationally fast, still reliable fit of the morphological and kinematic properties of galaxies. Here we apply \textsc{bang} to the SDSS-MaNGA survey, estimating key parameters such as mass, radial extensions, and dynamics, for both bulges and discs of +10,000 objects. We test our methodology against a smaller subsample of galaxies independently analysed with an orbit-based algorithm, finding agreement in the recovered total stellar mass. We also manage to reproduce well-established scaling relations, demonstrating how proper dynamical modelling can result in tighter correlations and provide corrections to standard approaches. Finally, we propose a more general way of decomposing galaxies into "hot" and "cold" components, showing a correlation with orbit-based approaches and visually determined morphological type. Unexpected tails in the "hot-to-total" mass-ratio distribution are present for galaxies of all morphologies, possibly due to visual morphology misclassifications. △ Less

Submitted 9 January, 2024; v1 submitted 5 May, 2023; originally announced May 2023.

Comments: 13 pages, 10 figures, accepted for pubblication on MNRAS

Quasi-periodic eruptions from impacts between the secondary and a rigidly precessing accretion disc in an extreme mass-ratio inspiral system

Authors: Alessia Franchini, Matteo Bonetti, Alessandro Lupi, Giovanni Miniutti, Elisa Bortolas, Margherita Giustini, Massimo Dotti, Alberto Sesana, Riccardo Arcodia, Taeho Ryu

Abstract: X-ray quasi-periodic eruptions (QPEs) represent a recently discovered example of extreme X-ray variability associated with supermassive black holes. These are high-amplitude bursts recurring every few hours that are detected in the soft X-ray band from the nuclei of nearby galaxies whose optical spectra lack the broad emission lines typically observed in unobscured active galaxies. The physical or… ▽ More X-ray quasi-periodic eruptions (QPEs) represent a recently discovered example of extreme X-ray variability associated with supermassive black holes. These are high-amplitude bursts recurring every few hours that are detected in the soft X-ray band from the nuclei of nearby galaxies whose optical spectra lack the broad emission lines typically observed in unobscured active galaxies. The physical origin of this new X-ray variability phenomenon is still unknown and several theoretical models have been presented. However, no attempt has been made so far to account for the varying QPE recurrence time and luminosity in individual sources, nor for the diversity of the QPE phenomenology in the different known erupters. We present a semi-analytical model based on an extreme mass-ratio inspiral (EMRI) system where the secondary intersects, along its orbit, a rigidly precessing accretion disc surrounding the primary. We assume that QPEs result from emission from an adiabatically expanding, initially optically thick gas cloud expelled from the disc plane at each impact. We produced synthetic X-ray light curves, which we then compared with X-ray data from four QPE sources: GSN 069, eRO-QPE1, eRO-QPE2, and RX J1301.9+2747. Our model aptly reproduces the diversity of QPE properties between the considered objects and it is also able to naturally account for the varying QPE amplitudes and recurrence times in individual sources. Future implementations will enable us to refine the match with the data and to estimate the system parameters precisely, making additional use of multi-epoch QPE data. We briefly discuss the nature of the secondary object, as well as the possible implications of our findings for the EMRI population at large. △ Less

Submitted 14 June, 2023; v1 submitted 3 April, 2023; originally announced April 2023.

Comments: 11 pages, 7 figures, accepted for publication in A&A

Journal ref: A&A 675, A100 (2023)

Gravitational waves from an eccentric population of primordial black holes orbiting Sgr A$^{\star}$

Authors: Stefano Bondani, Matteo Bonetti, Luca Broggi, Francesco Haardt, Alberto Sesana, Massimo Dotti

Abstract: Primordial black holes (PBH), supposedly formed in the very early Universe, have been proposed as a possible viable dark matter candidate. In this work we characterize the expected gravitational wave (GW) losses from a population of PBHs orbiting Sgr A$^{\star}$, the super-massive black hole at the Galactic center (GC), and assess the signal detectability by the planned space-borne interferometer… ▽ More Primordial black holes (PBH), supposedly formed in the very early Universe, have been proposed as a possible viable dark matter candidate. In this work we characterize the expected gravitational wave (GW) losses from a population of PBHs orbiting Sgr A$^{\star}$, the super-massive black hole at the Galactic center (GC), and assess the signal detectability by the planned space-borne interferometer LISA and by the proposed next generation space-borne interferometer $μ$Ares. Assuming that PBHs indeed form the entire diffuse mass allowed to reside within the orbit of the S2 star, we compute an upper limit to the expected GW signal both from resolved and non-resolved sources, under the further assumptions of monochromatic mass function and thermally distributed eccentricities. By comparing with our previous work where PBHs on circular orbits were assumed, we show for 1 M$_{\odot}$ PBHs how the GW signal from high harmonics over a 10 year data stream increases by a factor of six the chances of LISA detectability, from the $\approx 10\%$ of the circular case, to $\approx 60\%$, whereas multiple sources can be identified in $20\%$ of our mock populations. The background signal, made by summing up all non resolved sources, should be certainly detectable thanks to the PBHs with higher eccentricity evolving under two body relaxation. In the case of $μ$Ares, because of its improved sensitivity in the $μ$Hz band, one third of the entire population of PBHs orbiting Sgr A$^{\star}$ would be resolved. The background noise from the remaining non resolved sources should be detectable as well. Finally we present the results for different PBH masses. △ Less

Submitted 19 February, 2024; v1 submitted 22 March, 2023; originally announced March 2023.

Comments: 11 pages, 9 figures

Journal ref: Phys. Rev. D 109, 043005 (2024)

Another X-ray UFO without a momentum-boosted molecular outflow. ALMA CO(1-0) observations of the galaxy pair IRAS 05054+1718

Authors: F. Bonanomi, C. Cicone, P. Severgnini, V. Braito, C. Vignali, J. N. Reeves, M. Sirressi, I. Montoya Arroyave, R. Della Ceca, L. Ballo, M. Dotti

Abstract: We present ALMA CO(1-0) observations of the nearby LIRG galaxy pair IRAS05054+1718 with a new analysis of X-ray data collected between 2012 and 2021 using NuSTAR, Swift, and XMM-Newton. The western component of the pair, NED01, hosts a Seyfert 1.9 nucleus launching a powerful X-ray UFO. Our X-ray spectral analysis suggests the UFO could be variable or multi-component in velocity and constrains its… ▽ More We present ALMA CO(1-0) observations of the nearby LIRG galaxy pair IRAS05054+1718 with a new analysis of X-ray data collected between 2012 and 2021 using NuSTAR, Swift, and XMM-Newton. The western component of the pair, NED01, hosts a Seyfert 1.9 nucleus launching a powerful X-ray UFO. Our X-ray spectral analysis suggests the UFO could be variable or multi-component in velocity and constrains its momentum flux to $\dot p^{X-ray}_{out} \sim (4\pm2)\times 10^{34}$ gcms$^{-2}$. ALMA CO(1-0) observations include also the eastern component of the pair, a LIRG with no clear evidence for an AGN. We study the CO(1-0) kinematics in the two galaxies using the 3D-BAROLO code. In both sources, we can model the bulk of the CO(1-0) emission with rotating disks and, after subtracting the best-fit models, we detect compact residual emission at S/N=15 within $\sim3$kpc from the centre. A molecular outflow in NED01, if present, cannot be brighter than such residuals, implying an upper limit on its outflow rate of $\dot{M}^{mol}_{out} \lesssim 19\pm14~M_{\odot}~yr^{-1}$ and on its momentum rate of $\dot p^{mol}_{out} \lesssim (2.7\pm2.4) \times 10^{34}$gcms$^{-1}$. Combined with the revised energetics of the X-ray wind, we derive an upper limit on the momentum rate ratio of $\dot{p}^{mol}_{out}/\dot{p}^{X-ray}_{out}<0.67$. We discuss these results in the context of the expectations of AGN feedback models, and we propose the X-ray disk wind in NED01 has not significantly impacted the molecular gas reservoir (yet), and we can constrain its effect to be much smaller than expectations of AGN ''energy-driven'' feedback models. We also consider and discuss the hypothesis of asymmetries of the molecular disk not properly captured by the 3D-BAROLO code. Our results highlight the challenges in testing the predictions of popular AGN disk-wind feedback theories, even with good quality multi-wavelength observations. △ Less

Submitted 16 May, 2023; v1 submitted 1 March, 2023; originally announced March 2023.

Comments: Accepted by A&A

Journal ref: A&A 673, A46 (2023)

The fragility of thin discs in galaxies -- I. Building tailored N-body galaxy models

Authors: Pablo M. Galán-de Anta, Eugene Vasiliev, Marc Sarzi, Massimo Dotti, Pedro R. Capelo, Andrea Incatasciato, Lorenzo Posti, Lorenzo Morelli, Enrico Maria Corsini

Abstract: Thin stellar discs on both galactic and nuclear, sub-kpc scales are believed to be fragile structures that would be easily destroyed in major mergers. In turn, this makes the age-dating of their stellar populations a useful diagnostics for the assembly history of galaxies. We aim at carefully exploring the fragility of such stellar discs in intermediate- and low- mass encounters, using high-resolu… ▽ More Thin stellar discs on both galactic and nuclear, sub-kpc scales are believed to be fragile structures that would be easily destroyed in major mergers. In turn, this makes the age-dating of their stellar populations a useful diagnostics for the assembly history of galaxies. We aim at carefully exploring the fragility of such stellar discs in intermediate- and low- mass encounters, using high-resolution N-body simulations of galaxy models with structural and kinematic properties tailored to actually observed galaxies. As a first but challenging step, we create a dynamical model of FCC 170, a nearly edge-on galaxy in the Fornax cluster with multiple galactic components and including both a galactic scale and nuclear stellar disc (NSD), using detailed kinematic data from the Multi Unit Spectroscopic Explorer and a novel method for constructing distribution function-based self-consistent galaxy models. We then create N-body realisations of this model and demonstrate that it remains in equilibrium and preserves its properties over many Gyr, when evolved with a sufficiently high particle number. However, the NSD is more prone to numerical heating, which gradually increases its thickness by up to 22 per cent in 10 Gyr even in our highest-resolution runs. Nevertheless, these N-body models can serve as realistic representations of actual galaxies in merger simulations. △ Less

Submitted 7 February, 2023; originally announced February 2023.

Comments: 13 pages, 9 figures

Dynamical evolution of massive black hole pairs in the presence of spin-dependent radiative feedback

Authors: Francesco Bollati, Alessandro Lupi, Massimo Dotti, Francesco Haardt

Abstract: The putative ubiquity of massive black holes (MBH) at the center of galaxies, and the hierarchical progress of structure formation along the cosmic history, together necessarily imply the existence of a large population of cosmic MBH binaries. Such systems are understood to be the loudest sources of gravitational waves at mHz frequencies, the regime that will be probed by the next Laser Interferom… ▽ More The putative ubiquity of massive black holes (MBH) at the center of galaxies, and the hierarchical progress of structure formation along the cosmic history, together necessarily imply the existence of a large population of cosmic MBH binaries. Such systems are understood to be the loudest sources of gravitational waves at mHz frequencies, the regime that will be probed by the next Laser Interferometer Space Antenna (LISA). It has been proposed that the rate at which MBHs pair and then bind to form binaries is critically dependent upon the feedback exerted by the MBHs on the surrounding gaseous environment. Using the publicly available code GIZMO, we perform a suite of simulations aimed at studying the dynamics of a MBH pair embedded in a gaseous disk on 100 pc scale. By means of dedicated modules, we follow the dynamics of MBHs in the presence of different spin-dependent radiative feedback models, and compare the results to a benchmark case with no feedback at all. Our main finding is that feedback causes the secondary MBH to shrink its orbit at a reduced pace, when compared to models where feedback is absent. Moreover, such slower inspiral occurs on eccentric orbits, as feedback has the net effect of hampering the circularization process. Though idealized in many aspects, our study highlights and quantifies the importance of including spin-dependent feedback recipes in hydrodynamic simulations of MBH pairs, and ultimately in assessing the cosmological coalescence rate of such systems in view of their detection through gravitational waves. △ Less

Submitted 16 December, 2022; originally announced December 2022.

Morphological decomposition of TNG50 galaxies: methodology and catalogue

Authors: Tommaso Zana, Alessandro Lupi, Matteo Bonetti, Massimo Dotti, Yetli Rosas-Guevara, David Izquierdo-Villalba, Silvia Bonoli, Lars Hernquist, Dylan Nelson

Abstract: 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 algor… ▽ More 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. △ Less

Submitted 8 July, 2022; v1 submitted 9 June, 2022; originally announced June 2022.

Comments: 20 pages, 20 figures. Submitted to MNRAS. MORDOR is available at https://github.com/thanatom/mordor. The catalogue is available at www.tng-project.org/zana22

Maximally informed Bayesian modelling of disc galaxies

Authors: Fabio Rigamonti, Massimo Dotti, Stefano Covino, Francesco Haardt, Marco Landoni, Walter Del Pozzo, Alessandro Lupi, Stefano Zibetti

Abstract: 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 full… ▽ More 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. △ Less

Submitted 17 May, 2022; v1 submitted 13 May, 2022; originally announced May 2022.

Comments: 14 pages, 7 figures, accepted for publication in MNRAS

Extreme mass ratio inspirals and tidal disruption events in nuclear clusters. I. Time dependent rates

Authors: Luca Broggi, Elisa Bortolas, Matteo Bonetti, Alberto Sesana, Massimo Dotti

Abstract: In this paper we develop a computationally efficient, two-population, time-dependent Fokker-Plank approach in the two dimensions of energy and angular momentum to study the rates of tidal disruption events (TDEs), extreme mass ratio inspirals (EMRIs) and direct plunges occurring around massive black holes (MBHs) in galactic nuclei. We test our code by exploring a wide range of the astrophysically… ▽ More In this paper we develop a computationally efficient, two-population, time-dependent Fokker-Plank approach in the two dimensions of energy and angular momentum to study the rates of tidal disruption events (TDEs), extreme mass ratio inspirals (EMRIs) and direct plunges occurring around massive black holes (MBHs) in galactic nuclei. We test our code by exploring a wide range of the astrophysically relevant parameter space, including MBH masses, galaxy central densities and inner density slopes. We find that mass segregation and, more in general, the time dependency of the distribution function regulate the event rate: TDEs always decline with time, whereas EMRIs and plunges reach a maximum and undergo a subsequent nearly exponential decay. Once suitably normalized, the rates associated to different choices of MBH mass and galaxy density overlap nearly perfectly. Based on this, we provide a simple scaling that allows to reproduce the time-dependent event rates for any MBH mass and underlying galactic nucleus. Although our peak rates are in general agreement with the literature relying on the steady-state (non-time dependent) assumption, those can be sustained on a timescale that strongly depends on the properties of the system. In particular this can be much shorter than a Gyr for relatively light MBHs residing in dense systems. This warns against using steady state models to compute global TDE, EMRI and plunge rates and calls for a more sophisticated, time dependent treatment of the problem. △ Less

Submitted 12 May, 2022; originally announced May 2022.

Comments: 14 pages, 10 figures, submitted to MNRAS

Optical follow-up of the tick-tock massive black hole binary candidate

Authors: Massimo Dotti, Matteo Bonetti, Fabio Rigamonti, Elisa Bortolas, Matteo Fossati, Roberto Decarli, Stefano Covino, Alessandro Lupi, Alessia Franchini, Alberto Sesana, Giorgio Calderone

Abstract: The observation of a population of massive black hole binaries (MBHBs) is key for our complete understanding of galaxy mergers and for the characterization of the expected gravitational waves (GWs) signal. However, MBHBs still remain elusive with only a few candidates proposed to date. Among these, SDSSJ143016.05+230344.4 ('tick-tock' hereafter) is the only candidate with a remarkably well sampled… ▽ More The observation of a population of massive black hole binaries (MBHBs) is key for our complete understanding of galaxy mergers and for the characterization of the expected gravitational waves (GWs) signal. However, MBHBs still remain elusive with only a few candidates proposed to date. Among these, SDSSJ143016.05+230344.4 ('tick-tock' hereafter) is the only candidate with a remarkably well sampled light curve showing a clear reduction of the modulation period and amplitude over three years of observations. This particular feature has been recently claimed to be the signature of a MBHB that is about to merge. In this paper, we provide an optical follow-up of the tick-tock source using the Rapid Eye Mount (REM) telescope. The decreasing luminosity observed in our follow up is hardly explained within the binary scenario. We speculate about an alternative scenario that might explain the observed light curve through relativistic Lense-Thirring precession of an accretion disc around a single massive black hole. △ Less

Submitted 14 November, 2022; v1 submitted 12 May, 2022; originally announced May 2022.

Comments: 9 pages, 2 figure, accepted for publication in MNRAS

NOEMA observations support a recoiling black hole in 3C186

Authors: G. Castignani, E. Meyer, M. Chiaberge, F. Combes, T. Morishita, R. Decarli, A. Capetti, M. Dotti, G. R. Tremblay, C. A. Norman

Abstract: 3C186 is a powerful radio loud quasar (QSO) at the center of a cool-core cluster at z=1.06. Previous studies reported evidence for a projected spatial offset of ~1'' between the isophotal center of the galaxy and the point-source QSO as well as a spectral shift of ~2000 km/s between the narrow and broad line region of the system. In this work we report high-resolution molecular gas CO(4-3) observa… ▽ More 3C186 is a powerful radio loud quasar (QSO) at the center of a cool-core cluster at z=1.06. Previous studies reported evidence for a projected spatial offset of ~1'' between the isophotal center of the galaxy and the point-source QSO as well as a spectral shift of ~2000 km/s between the narrow and broad line region of the system. In this work we report high-resolution molecular gas CO(4-3) observations of the system taken with NOEMA interferometer. We clearly detect a large reservoir of molecular gas, $M_{H_2}\sim8\times10^{10}~M_\odot$, that is cospatial with the host galaxy and likely associated with a rotating disk-like structure. We firmly confirm both the spatial offset of the galaxy's gas reservoir with respect to the continuum emission of the QSO and the spectral offset with respect to the redshift of the broad line region. Our morphological and kinematical analysis confirms that the most likely scenario to explain the 3C186 system is that the QSO is a kicked super-massive black hole (SMBH), which we believe may have resulted from a strong gravitational wave recoil as two SMBHs coalesced after the merger of their host galaxies. △ Less

Submitted 12 April, 2022; originally announced April 2022.

Comments: 9 pages, Astronomy & Astrophysics Letters in press

Extreme Mass Ratio Inspirals triggered by Massive Black Hole Binaries: from Relativistic Dynamics to Cosmological Rates

Authors: Giovanni Mazzolari, Matteo Bonetti, Alberto Sesana, Riccardo M. Colombo, Massimo Dotti, Giuseppe Lodato, David Izquierdo-Villalba

Abstract: Extreme mass ratio inspirals (EMRIs) are compact binary systems characterized by a mass-ratio $q=m/M$ in the range $~10^{-9}-10^{-4}$ and represent primary gravitational wave (GW) sources for the forthcoming Laser Interferometer Space Antenna (LISA). While their standard formation channel involves relaxation processes deflecting compact objects on very low angular momentum orbits around the centra… ▽ More Extreme mass ratio inspirals (EMRIs) are compact binary systems characterized by a mass-ratio $q=m/M$ in the range $~10^{-9}-10^{-4}$ and represent primary gravitational wave (GW) sources for the forthcoming Laser Interferometer Space Antenna (LISA). While their standard formation channel involves relaxation processes deflecting compact objects on very low angular momentum orbits around the central massive black hole, a number of alternative formation channels has been proposed, including binary tidal break-up, migration in accretion disks and secular and chaotic dynamics around a massive black hole binary (MBHB). In this work, we take an extensive closer look at this latter scenario, investigating how EMRIs can be triggered by a MBHBs, formed in the aftermath of galaxy mergers. By employing a suite of relativistic three-body simulations, we evaluate the efficiency of EMRI formation for different parameters of the MBHB, assessing the importance of both secular and chaotic dynamics. By modelling the distribution of compact objects in galaxy nuclei, we estimate the resulting EMRI formation rate, finding that EMRI are produced in a sharp burst, with peak rates that are 10-100 times higher than the standard two-body relaxation channel, lasting for 10$^6$--10$^8$ years. By coupling our results with an estimate of the cosmic MBHB merger rate, we finally forecast that LISA could observe ${\cal O}(10)$ EMRIs per year formed by this channel. △ Less

Submitted 8 August, 2022; v1 submitted 11 April, 2022; originally announced April 2022.

Comments: 18 pages, 15 figures, accepted for publication in MNRAS

Re-weighting dark matter in disc galaxies: a new general relativistic observational test

Authors: Davide Astesiano, Sergio L. Cacciatori, Massimo Dotti, Francesco Haardt, Federico Re

Abstract: A recent analysis of data from the ESA Gaia mission demonstrated that the kinematics of stars in the Milky Way can be modelled without invoking the presence of dark matter whatsoever. Indeed, the higher-than-Keplerian velocities observed in outer stars can be ascribed to the properties of the general relativistic (GR) metric assumed to describe the Galaxy. Here we generalize the concept, and deriv… ▽ More A recent analysis of data from the ESA Gaia mission demonstrated that the kinematics of stars in the Milky Way can be modelled without invoking the presence of dark matter whatsoever. Indeed, the higher-than-Keplerian velocities observed in outer stars can be ascribed to the properties of the general relativistic (GR) metric assumed to describe the Galaxy. Here we generalize the concept, and derive the most general exact GR model for a stationary axi-symmetric dynamically cold dust structure. We explicitly show how deviations from the commonly adopted Newtonian dynamics are indeed manifest even at low velocities and low densities. We provide for the first time a detailed description of the frequency shift experienced by photons travelling from any emission site within an external disc galaxy to the observer, relating the outcome of the shift measurement to the gravitational properties of the galaxy. Finally, we devise a novel, groundbreaking observational test potentially able to fully characterize the GR metric under the minimal set of assumptions mentioned above. The proposed experiment exploits the effects non-diagonal GR terms have on the frequency shift of photons, ultimately providing a test to evaluate whether dark matter is actually required by disc galaxy kinematics. △ Less

Submitted 11 April, 2022; originally announced April 2022.

Comments: 7 pages, 0 figures

Disc instability and bar formation: view from the IllustrisTNG simulations

Authors: David Izquierdo-Villalba, Silvia Bonoli, Yetli Rosas-Guevara, Volker Springel, Simon D. M. White, Tommaso Zana, Massimo Dotti, Daniele Spinoso, Matteo Bonetti, Alessandro Lupi

Abstract: 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 m… ▽ More 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. △ Less

Submitted 19 May, 2022; v1 submitted 15 March, 2022; originally announced March 2022.

Comments: 15 pages, 12 Figures, Accepted by MNRAS

Astrophysics with the Laser Interferometer Space Antenna

Authors: Pau Amaro Seoane, Jeff Andrews, Manuel Arca Sedda, Abbas Askar, Quentin Baghi, Razvan Balasov, Imre Bartos, Simone S. Bavera, Jillian Bellovary, Christopher P. L. Berry, Emanuele Berti, Stefano Bianchi, Laura Blecha, Stephane Blondin, Tamara Bogdanović, Samuel Boissier, Matteo Bonetti, Silvia Bonoli, Elisa Bortolas, Katelyn Breivik, Pedro R. Capelo, Laurentiu Caramete, Federico Cattorini, Maria Charisi, Sylvain Chaty , et al. (134 additional authors not shown)

Abstract: The Laser Interferometer Space Antenna (LISA) will be a transformative experiment for gravitational wave astronomy, and, as such, it will offer unique opportunities to address many key astrophysical questions in a completely novel way. The synergy with ground-based and space-born instruments in the electromagnetic domain, by enabling multi-messenger observations, will add further to the discovery… ▽ More The Laser Interferometer Space Antenna (LISA) will be a transformative experiment for gravitational wave astronomy, and, as such, it will offer unique opportunities to address many key astrophysical questions in a completely novel way. The synergy with ground-based and space-born instruments in the electromagnetic domain, by enabling multi-messenger observations, will add further to the discovery potential of LISA. The next decade is crucial to prepare the astrophysical community for LISA's first observations. This review outlines the extensive landscape of astrophysical theory, numerical simulations, and astronomical observations that are instrumental for modeling and interpreting the upcoming LISA datastream. To this aim, the current knowledge in three main source classes for LISA is reviewed; ultracompact stellar-mass binaries, massive black hole binaries, and extreme or intermediate mass ratio inspirals. The relevant astrophysical processes and the established modeling techniques are summarized. Likewise, open issues and gaps in our understanding of these sources are highlighted, along with an indication of how LISA could help making progress in the different areas. New research avenues that LISA itself, or its joint exploitation with upcoming studies in the electromagnetic domain, will enable, are also illustrated. Improvements in modeling and analysis approaches, such as the combination of numerical simulations and modern data science techniques, are discussed. This review is intended to be a starting point for using LISA as a new discovery tool for understanding our Universe. △ Less

Submitted 25 May, 2023; v1 submitted 11 March, 2022; originally announced March 2022.

Journal ref: Living Reviews in Relativity, Volume 26, Article number: 2 (2023)

On the detectability of gravitational waves from primordial black holes orbiting Sgr A*

Authors: Stefano Bondani, Francesco Haardt, Alberto Sesana, Enrico Barausse, Massimo Dotti

Abstract: In this work we characterize the expected gravitational wave signal detectable by the planned space-borne interferometer LISA and the proposed next generation space-borne interferometer $μ$Ares arising from a population of primordial black holes orbiting Sgr A*, the super-massive black hole at the Galactic center. Assuming that such objects indeed form the entire diffuse mass allowed by the observ… ▽ More In this work we characterize the expected gravitational wave signal detectable by the planned space-borne interferometer LISA and the proposed next generation space-borne interferometer $μ$Ares arising from a population of primordial black holes orbiting Sgr A*, the super-massive black hole at the Galactic center. Assuming that such objects indeed form the entire diffuse mass allowed by the observed orbit of S2 in the Galactic center, under the simplified assumption of circular orbits and monochromatic mass function, we assess the expected signal in gravitational waves, either from resolved and non-resolved sources. We estimate a small but non negligible chance of $\simeq$ 10% of detecting one single 1 M$_{\odot}$ primordial black hole with LISA in a 10-year-long data stream, while the background signal due to unresolved sources would essentially elude any reasonable chance of detection. On the contrary, $μ$Ares, with a $\simeq$ 3 orders-of-magnitude better sensitivity at $\simeq$ 10$^{-5}$ Hz, would be able to resolve $\simeq$ 140 solar mass primordial black holes in the same amount of time, while the unresolved background should be observable with an integrated signal-to-noise ratio $\gtrsim$ 100. Allowing the typical PBH mass to be in the range 0.01-10 M$_{\odot}$ would increase LISA chance of detection to $\simeq$ 40% towards the lower limit of the mass spectrum. In the case of $μ$Ares, instead, we find a "sweet spot" just about 1 M$_{\odot}$, a mass for which the number of resolvable events is indeed maximized. △ Less

Submitted 2 August, 2022; v1 submitted 10 March, 2022; originally announced March 2022.

Comments: 10 pages, 11 figures; Referee reviewed

Journal ref: Phys. Rev. D 106, 043015, 2022

Dual AGN in the Horizon-AGN simulation and their link to galaxy and massive black hole mergers, with an excursus on multiple AGN

Authors: Marta Volonteri, Hugo Pfister, Ricarda Beckmann, Massimo Dotti, Yohan Dubois, Warren Massonneau, Gibwa Musoke, Michael Tremmel

Abstract: The occurrence of dual active galactic nuclei (AGN) on scales of a few tens of kpc can be used to study merger-induced accretion on massive black holes (MBHs) and to derive clues on MBH mergers, using dual AGN as a parent population of precursors. We investigate the properties of dual AGN in the cosmological simulation Horizon-AGN. We create catalogs of dual AGN selected with distance and luminosi… ▽ More The occurrence of dual active galactic nuclei (AGN) on scales of a few tens of kpc can be used to study merger-induced accretion on massive black holes (MBHs) and to derive clues on MBH mergers, using dual AGN as a parent population of precursors. We investigate the properties of dual AGN in the cosmological simulation Horizon-AGN. We create catalogs of dual AGN selected with distance and luminosity criteria, plus sub-catalogs where further mass cuts are applied. We divide the sample into dual AGN hosted in different galaxies, on the way to a merger, and into those hosted in one galaxy, after the galaxy merger has happened. We find that the relation between MBH and galaxy mass is similar to that of general AGN population and we compare the properties of dual AGN also with a control sample, discussing differences and similarities in masses and Eddington ratios. The typical mass ratio of galaxy mergers associated to dual AGN is 0.2, with mass loss in the smaller galaxy decreasing the mass ratio as the merger progresses. Between 30 and 80 per cent of dual AGN with separations between 4 and 30~kpc can be matched to an ensuing MBH merger. The dual AGN fraction increases with redshift and with separation threshold, although above 50~kpc the increase of multiple AGN limits that of duals. Multiple AGN are generally associated with massive halos, and mass loss of satellites shapes the galaxy-halo relation. △ Less

Submitted 17 June, 2022; v1 submitted 14 December, 2021; originally announced December 2021.

Comments: Published in MNRAS

The evolution of the barred galaxy population in the TNG50 simulation

Authors: Yetli Rosas-Guevara, Silvia Bonoli, Massimo Dotti, David Izquierdo-Villalba, Alessandro Lupi, Tommaso Zana, Matteo Bonetti, Dylan Nelson, Volker Springel, Lars Hernquist, Mark Vogelsberger

Abstract: 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… ▽ More 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$. △ Less

Submitted 30 March, 2022; v1 submitted 9 October, 2021; originally announced October 2021.

Comments: 19 pages and 16 figures. Accepted in MNRAS and published catalogues in https://www.tng-project.org/data/docs/specifications/\#sec5j

Circumbinary disc self-gravity governing supermassive black hole binary mergers

Authors: Alessia Franchini, Alberto Sesana, Massimo Dotti

Abstract: Understanding the interaction of massive black hole binaries with their gaseous environment is crucial since at sub-parsec scales the binary is too wide for gravitational wave emission to take over and to drive the two black holes to merge. We here investigate the interaction between a massive black hole binary and a self-gravitating circumbinary disc using 3D smoothed particle hydrodynamics simul… ▽ More Understanding the interaction of massive black hole binaries with their gaseous environment is crucial since at sub-parsec scales the binary is too wide for gravitational wave emission to take over and to drive the two black holes to merge. We here investigate the interaction between a massive black hole binary and a self-gravitating circumbinary disc using 3D smoothed particle hydrodynamics simulations. We find that, when the disc self-gravity regulates the angular momentum transport, the binary semi-major axis decreases regardless the choice of disc masses and temperatures, within the range we explored. In particular, we find that the disc initial temperature (hence the disc aspect ratio) has little effect on the evolution of the binary since discs with the same mass self-regulate towards the same temperature. Initially warmer discs cause the binary to shrink on a slightly shorter timescale until the disc has reached the self-regulated equilibrium temperature. More massive discs drive the binary semi-major axis to decrease at a faster pace compared to less massive discs and result in faster binary eccentricity growth even after the initial-condition-dependent transient evolution. Finally we investigate the effect that the initial cavity size has on the binary-disc interaction and we find that, in the self-gravitating regime, an initially smaller cavity leads to a much faster binary shrinking, as expected. Our results are especially important for very massive black hole binaries such as those in the PTA band, for which gas self gravity cannot be neglected. △ Less

Submitted 24 June, 2021; originally announced June 2021.

Comments: 9 pages, submitted to MNRAS on June 1st 2021

Stellar hardening of massive black hole binaries: the impact of the host rotation

Authors: Ludovica Varisco, Elisa Bortolas, Massimo Dotti, Alberto Sesana

Abstract: Massive black hole binaries (MBHBs) with masses of ~ 10^4 to ~ 10^10 of solar masses are one of the main targets for currently operating and forthcoming space-borne gravitational wave observatories. In this paper, we explore the effect of the stellar host rotation on the bound binary hardening efficiency, driven by three-body stellar interactions. As seen in previous studies, we find that the cent… ▽ More Massive black hole binaries (MBHBs) with masses of ~ 10^4 to ~ 10^10 of solar masses are one of the main targets for currently operating and forthcoming space-borne gravitational wave observatories. In this paper, we explore the effect of the stellar host rotation on the bound binary hardening efficiency, driven by three-body stellar interactions. As seen in previous studies, we find that the centre of mass (CoM) of a prograde MBHB embedded in a rotating environment starts moving on a nearly circular orbit about the centre of the system shortly after the MBHB binding. In our runs, the oscillation radius is approximately 0.25 ( approximately 0.1) times the binary influence radius for equal mass MBHBs (MBHBs with mass ratio 1:4). Conversely, retrograde binaries remain anchored about the centre of the host. The binary shrinking rate is twice as fast when the binary CoM exhibits a net orbital motion, owing to a more efficient loss cone repopulation even in our spherical stellar systems. We develop a model that captures the CoM oscillations of prograde binaries; we argue that the CoM angular momentum gain per time unit scales with the internal binary angular momentum, so that most of the displacement is induced by stellar interactions occurring around the time of MBHB binding, while the subsequent angular momentum enhancement gets eventually quashed by the effect of dynamical friction. The effect of the background rotation on the MBHB evolution may be relevant for LISA sources, that are expected to form in significantly rotating stellar systems. △ Less

Submitted 13 September, 2021; v1 submitted 29 April, 2021; originally announced April 2021.

Comments: 10 pages, 9 figures, 2 tables. Accepted for publication on MNRAS

Binary black hole signatures in polarized light curves

Authors: Massimo Dotti, Matteo Bonetti, Daniel J. D'Orazio, Zoltan Haiman, Luis C. Ho

Abstract: Variable active galactic nuclei showing periodic light curves have been proposed as massive black hole binary (MBHB) candidates. In such scenarios the periodicity can be due to relativistic Doppler-boosting of the emitted light. This hypothesis can be tested through the timing of scattered polarized light. Following the results of polarization studies in type I nuclei and of dynamical studies of M… ▽ More Variable active galactic nuclei showing periodic light curves have been proposed as massive black hole binary (MBHB) candidates. In such scenarios the periodicity can be due to relativistic Doppler-boosting of the emitted light. This hypothesis can be tested through the timing of scattered polarized light. Following the results of polarization studies in type I nuclei and of dynamical studies of MBHBs with circumbinary discs, we assume a coplanar equatorial scattering ring, whose elements contribute differently to the total polarized flux, due to different scattering angles, levels of Doppler boost, and line-of-sight time delays. We find that in the presence of a MBHB, both the degree of polarization and the polarization angle have periodic modulations. The minimum of the polarization degree approximately coincides with the peak of the light curve, regardless of the scattering ring size. The polarization angle oscillates around the semi-minor axis of the projected MBHB orbital ellipse, with a frequency equal either to the binary's orbital frequency (for large scattering screen radii), or twice this value (for smaller scattering structures). These distinctive features can be used to probe the nature of periodic MBHB candidates and to compile catalogs of the most promising sub-pc MBHBs. The identification of such polarization features in gravitational-wave detected MBHBs would enormously increase the amount of physical information about the sources, allowing the measurement of the individual masses of the binary components, and the orientation of the line of nodes on the sky, even for monochromatic gravitational wave signals. △ Less

Submitted 26 March, 2021; originally announced March 2021.

Comments: 11 pages, 6 figures, submitted to MNRAS

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

Authors: Elisa Bortolas, Matteo Bonetti, Massimo Dotti, Alessandro Lupi, Pedro R. Capelo, Lucio Mayer, Alberto Sesana

Abstract: 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-l… ▽ More 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. △ Less

Submitted 7 March, 2022; v1 submitted 12 March, 2021; originally announced March 2021.

Comments: 19 pages, 12 figures, 2 tables. MNRAS, accepted

MSC Class: 85A05

A possible sub-kiloparsec dual AGN buried behind the galaxy curtain

Authors: P. Severgnini, V. Braito, C. Cicone, P. Saracco, C. Vignali, R. Serafinelli, R. Della Ceca, M. Dotti, F. Cusano, D. Paris, G. Pruto, A. Zaino, L. Ballo, M. Landoni

Abstract: Although thousands of galaxy mergers are known, only a handful of sub-kiloparsec-scale supermassive black hole (SMBH) pairs have been confirmed so far, leaving a huge gap between the observed and predicted numbers of such objects. In this work, we present a detailed analysis of the Sloan Digital Sky Survey optical spectrum and of near-infrared (NIR) diffraction limited imaging of SDSS~J1431+4358.… ▽ More Although thousands of galaxy mergers are known, only a handful of sub-kiloparsec-scale supermassive black hole (SMBH) pairs have been confirmed so far, leaving a huge gap between the observed and predicted numbers of such objects. In this work, we present a detailed analysis of the Sloan Digital Sky Survey optical spectrum and of near-infrared (NIR) diffraction limited imaging of SDSS~J1431+4358. This object is a local radio-quiet type 2 active galactic nucleus (AGN) previously selected as a double AGN candidate on the basis of the double-peaked [OIII] emission line. The NIR adaptive optics-assisted observations were obtained at the Large Binocular Telescope with the LUCI+FLAO camera. We found that most of the prominent optical emission lines are characterized by a double-peaked profile, mainly produced by AGN photoionization. Our spectroscopical analysis disfavors the hypothesis that the double-peaked emission lines in the source are the signatures of outflow kinematics, leaving open the possibility that we are detecting either the rotation of a single narrow-line region or the presence of two SMBHs orbiting around a common central potential. The latter scenario is further supported by the high-spatial resolution NIR imaging: after subtracting the dominant contribution of the stellar bulge component in the host galaxy, we detect two faint nuclear sources at r<0.5 kpc projected separation. Interestingly, the two sources have a position angle consistent with that defined by the two regions where the [OIII] double peaks most likely originate. Aside from the discovery of a promising sub-kiloparsec scale dual AGN, our analysis shows the importance of an appropriate host galaxy subtraction in order to achieve a reliable estimate of the incidence of dual AGNs at small projected separations. △ Less

Submitted 24 February, 2021; v1 submitted 16 December, 2020; originally announced December 2020.

Comments: 11 pages, 8 figures, 1 table, accepted by A&A (revised to match the A&A version)

Journal ref: A&A 646, A153 (2021)

Non-isotropic feedback from accreting spinning black holes

Authors: Luca Sala, Elia Cenci, Pedro R. Capelo, Alessandro Lupi, Massimo Dotti

Abstract: 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 ne… ▽ More 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. △ Less

Submitted 13 April, 2022; v1 submitted 12 November, 2020; originally announced November 2020.

Comments: 14 pages, 6 figures, accepted for publication in MNRAS. Updated with published version, minimal changes

Journal ref: MNRAS, 500, 4788 (2021)

Black hole spin evolution in warped accretion discs

Authors: Elia Cenci, Luca Sala, Alessandro Lupi, Pedro R. Capelo, Massimo Dotti

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 m… ▽ More 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. △ Less

Submitted 13 April, 2022; v1 submitted 12 November, 2020; originally announced November 2020.

Comments: 10 pages, 5 figures, accepted for publication in MNRAS. Updated with published version, minimal changes

Journal ref: MNRAS, 500, 3719 (2021)

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

Authors: Matteo Bonetti, Elisa Bortolas, Alessandro Lupi, Massimo Dotti

Abstract: 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 i… ▽ More 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. △ Less

Submitted 23 January, 2021; v1 submitted 16 October, 2020; originally announced October 2020.

Comments: 17 pages, 9 figures; accepted for publication in MNRAS

Unveiling sub-parsec supermassive black hole binary candidates in active galactic nuclei

Authors: Roberto Serafinelli, Paola Severgnini, Valentina Braito, Roberto Della Ceca, Cristian Vignali, Filippo Ambrosino, Claudia Cicone, Alessandra Zaino, Massimo Dotti, Alberto Sesana, Vittoria E. Gianolli, Lucia Ballo, Valentina La Parola, Gabriele A. Matzeu

Abstract: Elusive supermassive black hole binaries (SMBHBs) are thought to be the penultimate stage of galaxy mergers, preceding a final coalescence phase. SMBHBs are sources of continuous gravitational waves, possibly detectable by pulsar timing arrays; the identification of candidates could help in performing targeted gravitational wave searches. Due to their origin in the innermost parts of active galact… ▽ More Elusive supermassive black hole binaries (SMBHBs) are thought to be the penultimate stage of galaxy mergers, preceding a final coalescence phase. SMBHBs are sources of continuous gravitational waves, possibly detectable by pulsar timing arrays; the identification of candidates could help in performing targeted gravitational wave searches. Due to their origin in the innermost parts of active galactic nuclei (AGN), X-rays are a promising tool to unveil the presence of SMBHBs, by means of either double Fe K$α$ emission lines or periodicity in their light curve. Here we report on a new method to select SMBHBs by means of the presence of a periodic signal in their Swift-BAT 105-months light curves. Our technique is based on the Fisher's exact g-test and takes into account the possible presence of colored noise. Among the 553 AGN selected for our investigation, only the Seyfert 1.5 Mrk 915 emerged as possible candidate for a SMBHB; from the subsequent analysis of its light curve we find a period $P_0=35\pm2$ months, and the null hypothesis is rejected at the $3.7σ$ confidence level. We also present a detailed analysis of the BAT light curve of the only previously X-ray-selected binary candidate source in the literature, the Seyfert 2 galaxy MCG+11-11-032. We find $P_0=26.3\pm0.6$ months, consistent with the one inferred from previously reported double Fe K$α$ emission lines. △ Less

Submitted 14 September, 2020; originally announced September 2020.

Comments: 11 pages, 10 figures; accepted for publication in The Astrophysical Journal

Black hole mergers from dwarf to massive galaxies with the NewHorizon and Horizon-AGN simulations

Authors: Marta Volonteri, Hugo Pfister, Ricarda S. Beckman, Yohan Dubois, Monica Colpi, Christopher J. Conselice, Massimo Dotti, Garreth Martin, Ryan Jackson, Katarina Kraljic, Christophe Pichon, Maxime Trebitsch, Sukyoung. K. Yi, Julien Devriendt, Sebastien Peirani

Abstract: Massive black hole (MBH) coalescences are powerful sources of low-frequency gravitational waves. To study these events in the cosmological context we need to trace the large-scale structure and cosmic evolution of a statistical population of galaxies, from dim dwarfs to bright galaxies. To cover such a large range of galaxy masses, we analyse two complementary simulations: Horizon-AGN with a large… ▽ More Massive black hole (MBH) coalescences are powerful sources of low-frequency gravitational waves. To study these events in the cosmological context we need to trace the large-scale structure and cosmic evolution of a statistical population of galaxies, from dim dwarfs to bright galaxies. To cover such a large range of galaxy masses, we analyse two complementary simulations: Horizon-AGN with a large volume and low resolution which tracks the high-mass (> 1e7 Msun) MBH population, and NewHorizon with a smaller volume but higher resolution that traces the low-mass (< 1e7 Msun) MBH population. While Horizon-AGN can be used to estimate the rate of inspirals for Pulsar Timing Arrays, NewHorizon can investigate MBH mergers in a statistical sample of dwarf galaxies for LISA, which is sensitive to low-mass MBHs. We use the same method to analyse the two simulations, post-processing MBH dynamics to account for time delays mostly determined by dynamical friction and stellar hardening. In both simulations, MBHs typically merge long after the galaxies do, so that the galaxy morphology at the time of the MBH merger is no longer determined by the galaxy merger from which the MBH merger originated. These time delays cause a loss of high-z MBH coalescences, shifting the peak of the MBH merger rate to z~1-2. This study shows how tracking MBH mergers in low-mass galaxies is crucial to probing the MBH merger rate for LISA and investigate the properties of the host galaxies. △ Less

Submitted 15 September, 2020; v1 submitted 11 May, 2020; originally announced May 2020.

Comments: MNRAS in press

Global torques and stochasticity as the drivers of massive black hole pairing in the young Universe

Authors: Elisa Bortolas, Pedro R. Capelo, Tommaso Zana, Lucio Mayer, Matteo Bonetti, Massimo Dotti, Melvyn B. Davies, Piero Madau

Abstract: The forthcoming Laser Interferometer Space Antenna (LISA) will probe the population of coalescing massive black hole (MBH) binaries up to the onset of structure formation. Here we simulate the galactic-scale pairing of $\sim10^6 M_\odot$ MBHs in a typical, non-clumpy main-sequence galaxy embedded in a cosmological environment at $z = 7-6$. In order to increase our statistical sample, we adopt a st… ▽ More The forthcoming Laser Interferometer Space Antenna (LISA) will probe the population of coalescing massive black hole (MBH) binaries up to the onset of structure formation. Here we simulate the galactic-scale pairing of $\sim10^6 M_\odot$ MBHs in a typical, non-clumpy main-sequence galaxy embedded in a cosmological environment at $z = 7-6$. In order to increase our statistical sample, we adopt a strategy that allows us to follow the evolution of six secondary MBHs concomitantly. We find that the magnitude of the dynamical-friction induced torques is significantly smaller than that of the large-scale, stochastic gravitational torques arising from the perturbed and morphologically evolving galactic disc, suggesting that the standard dynamical friction treatment is inadequate for realistic galaxies at high redshift. The dynamical evolution of MBHs is very stochastic, and a variation in the initial orbital phase can lead to a drastically different time-scale for the inspiral. Most remarkably, the development of a galactic bar in the host system either significantly accelerates the inspiral by dragging a secondary MBH into the centre, or ultimately hinders the orbital decay by scattering the MBH in the galaxy outskirts. The latter occurs more rarely, suggesting that galactic bars overall promote MBH inspiral and binary coalescence. The orbital decay time can be an order of magnitude shorter than what would be predicted relying on dynamical friction alone. The stochasticity, and the important role of global torques, have crucial implications for the rates of MBH coalescences in the early Universe: both have to be accounted for when making predictions for the upcoming LISA observatory. △ Less

Submitted 25 August, 2020; v1 submitted 5 May, 2020; originally announced May 2020.

Comments: Accepted for publication in MNRAS; 15 pages, 10 Figures, 2 Tables

MSC Class: 85A05

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

Authors: Matteo Bonetti, Elisa Bortolas, Alessandro Lupi, Massimo Dotti, Sandra I. Raimundo

Abstract: 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 circulari… ▽ More 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. △ Less

Submitted 4 April, 2020; v1 submitted 11 February, 2020; originally announced February 2020.

Comments: 7 pages, 3 figures, accepted for publication in MNRAS

From galactic nuclei to the halo outskirts: tracing supermassive black holes across cosmic history and environments

Authors: David Izquierdo-Villalba, Silvia Bonoli, Massimo Dotti, Alberto Sesana, Yetli Rosas-Guevara, Daniele Spinoso

Abstract: We study the mass assembly and spin evolution of supermassive black holes (BHs) across cosmic time as well as the impact of gravitational recoil on the population of nuclear and wandering black holes (wBHs) by using the semi-analytical model LGalaxies run on top of Millennium merger trees. We track spin changes that BHs experience during both coalescence events and gas accretion phases. For the la… ▽ More We study the mass assembly and spin evolution of supermassive black holes (BHs) across cosmic time as well as the impact of gravitational recoil on the population of nuclear and wandering black holes (wBHs) by using the semi-analytical model LGalaxies run on top of Millennium merger trees. We track spin changes that BHs experience during both coalescence events and gas accretion phases. For the latter, we assume that spin changes are coupled with the bulge assembly. This assumption leads to predictions for the median spin values of z=0 BHs that depend on whether they are hosted by pseudobulges, classical bulges or ellipticals, being a ~ 0.9, 0.7 and 0.4, respectively. The outcomes of the model display a good consistency with z<4 quasar luminosity functions and the z=0 BH mass function, spin values and black hole-bulge mass correlation. Regarding the wBHs, we assume that they can originate from both the disruption of satellite galaxies (orphan wBH) and ejections due to gravitational recoils (ejected wBH). The model points to a number density of wBHs that increases with decreasing redshift, although this population is always ~2 dex smaller than the one of nuclear black holes. At all redshifts, wBHs are typically hosted in M_{halo}>10^{13} M_{\odot} and M_{stellar}>10^{10} M_{\odot}, being orphan wBHs the dominant type. Besides, independently of redshift and halo mass, \textit{ejected} wBHs inhabit the central regions (<0.3R_{200}) of the host DM halo, while orphan wBHs linger at larger scales (>0.5R_{200}). Finally, we find that gravitational recoils cause a progressive depletion of nuclear BHs with decreasing redshift and stellar mass. Moreover, ejection events lead to changes in the predicted local black hole-bulge relation, in particular for BHs in pseudobulges, for which the relation is flattened at M_{bulge}>10^{10.2} M{\odot} and the scatter increase up to ~3dex. △ Less

Submitted 3 June, 2020; v1 submitted 28 January, 2020; originally announced January 2020.

Comments: 30 pages, 17 figures, accepted by MNRAS

The Quest for Dual and Binary Supermassive Black Holes: A Multi-Messenger View

Authors: Alessandra De Rosa, Cristian Vignali, Tamara Bogdanović, Pedro R. Capelo, Maria Charisi, Massimo Dotti, Bernd Husemann, Elisabeta Lusso, Lucio Mayer, Zsolt Paragi, Jessie Runnoe, Alberto Sesana, Lisa Steinborn, Stefano Bianchi, Monica Colpi, Luciano Del Valle, Sándor Frey, Krisztina É. Gabányi, Margherita Giustini, Matteo Guainazzi, Zoltan Haiman, Noelia Herrera Ruiz, Rubén Herrero-Illana, Kazushi Iwasawa, S. Komossa , et al. (5 additional authors not shown)

Abstract: The quest for binary and dual supermassive black holes (SMBHs) at the dawn of the multi-messenger era is compelling. Detecting dual active galactic nuclei (AGN) -- active SMBHs at projected separations larger than several parsecs -- and binary AGN -- probing the scale where SMBHs are bound in a Keplerian binary -- is an observational challenge. The study of AGN pairs (either dual or binary) also r… ▽ More The quest for binary and dual supermassive black holes (SMBHs) at the dawn of the multi-messenger era is compelling. Detecting dual active galactic nuclei (AGN) -- active SMBHs at projected separations larger than several parsecs -- and binary AGN -- probing the scale where SMBHs are bound in a Keplerian binary -- is an observational challenge. The study of AGN pairs (either dual or binary) also represents an overarching theoretical problem in cosmology and astrophysics. The AGN triggering calls for detailed knowledge of the hydrodynamical conditions of gas in the imminent surroundings of the SMBHs and, at the same time, their duality calls for detailed knowledge on how galaxies assemble through major and minor mergers and grow fed by matter along the filaments of the cosmic web. This review describes the techniques used across the electromagnetic spectrum to detect dual and binary AGN candidates and proposes new avenues for their search. The current observational status is compared with the state-of-the-art numerical simulations and models for formation of dual and binary AGN. Binary SMBHs are among the loudest sources of gravitational waves (GWs) in the Universe. The search for a background of GWs at nHz frequencies from inspiralling SMBHs at low redshifts, and the direct detection of signals from their coalescence by the Laser Interferometer Space Antenna in the next decade, make this a theme of major interest for multi-messenger astrophysics. This review discusses the future facilities and observational strategies that are likely to significantly advance this fascinating field. △ Less

Submitted 17 January, 2020; originally announced January 2020.

Comments: 69 pages, 21 figures. Accepted for publication in New Astronomy Reviews

Journal ref: New Astronomy Reviews, Vol. 86, id. 101525 (2019)

Real galaxy mergers from galaxy pair catalogs

Authors: Hugo Pfister, Massimo Dotti, Clothilde Laigle, Yohan Dubois, Marta Volonteri

Abstract: Mergers of galaxies are extremely violent events shaping their evolution. Such events are thought to trigger starbursts and, possibly, black hole accretion. Nonetheless, it is still not clear how to know the fate of a galaxy pair from the data available at a given time, limiting our ability to constrain the exact role of mergers. In this paper we use the lightcone of the Horizon-AGN simulation, fo… ▽ More Mergers of galaxies are extremely violent events shaping their evolution. Such events are thought to trigger starbursts and, possibly, black hole accretion. Nonetheless, it is still not clear how to know the fate of a galaxy pair from the data available at a given time, limiting our ability to constrain the exact role of mergers. In this paper we use the lightcone of the Horizon-AGN simulation, for which we know the fate of each pair, to test three selection processes aiming at identifying true merging pairs. We find that the simplest one (selecting objects within two thresholds on projected distance $d$ and redshift difference $Δz$) gives similar results than the most complex one (based on a neural network analyzing $d$, $Δz$, redshift of the primary, masses/star formation rates/aspect ratio of both galaxies). Our best thresholds are $d_\mathrm{th}\sim100\mathrm{\, kpc}$ and $Δz_\mathrm{th} \sim 10^{-3}$, in agreement with recent results. △ Less

Submitted 8 January, 2020; originally announced January 2020.

Comments: submitted to MNRAS

On the eccentricity evolution of massive black hole binaries in stellar backgrounds

Authors: Matteo Bonetti, Alexander Rasskazov, Alberto Sesana, Massimo Dotti, Francesco Haardt, Nathan W. C. Leigh, Manuel Arca Sedda, Giacomo Fragione, Elena Rossi

Abstract: We study the dynamical evolution of eccentric massive black hole binaries (MBHBs) interacting with unbound stars by means of an extensive set of three body scattering experiments. Compared to previous studies, we extend the investigation down to a MBHB mass ratio of $q=m_2/m_1=10^{-4}$, where $m_1$ and $m_2$ are the masses of the primary and secondary hole respectively. Contrary to a simple extrap… ▽ More We study the dynamical evolution of eccentric massive black hole binaries (MBHBs) interacting with unbound stars by means of an extensive set of three body scattering experiments. Compared to previous studies, we extend the investigation down to a MBHB mass ratio of $q=m_2/m_1=10^{-4}$, where $m_1$ and $m_2$ are the masses of the primary and secondary hole respectively. Contrary to a simple extrapolation from higher mass ratios, we find that for $q\lesssim 10^{-3}$ the eccentricity growth rate becomes negative, i.e., the binary {\it circularises} as it shrinks. This behaviour is due to the subset of interacting stars captured in metastable counter-rotating orbits; those stars tend to extract angular momentum from the binary, promoting eccentricity growth for $q>10^{-3}$, but tend to inject angular momentum into the binary driving it towards circularisation for $q<10^{-3}$. The physical origin of this behaviour requires a detailed study of the orbits of this subset of stars and is currently under investigation. Our findings might have important consequences for intermediate MBHs (IMBHs) inspiralling onto MBHs (e.g. a putative $10^3\rm M_{\odot}$ black hole inspiralling onto SgrA$^*$). △ Less

Submitted 4 February, 2020; v1 submitted 7 January, 2020; originally announced January 2020.

Comments: 6 pages, 5 figures; accepted for publication in MNRAS

The buildup of strongly barred galaxies in the TNG100 simulation

Authors: Yetli Rosas-Guevara, Silvia Bonoli, Massimo Dotti, Tommaso Zana, Dylan Nelson, Annalisa Pillepich, Luis C. Ho, David Izquierdo-Villalba, Lars Hernquist, Ruediger Pakmor

Abstract: We analyse the properties of strongly barred disc galaxies using the TNG100 simulation, a cosmological hydrodynamical realisation of the IllustrisTNG suite. We identify 270 disc galaxies at $z=0$ in the stellar mass range $M_{*}=10^{10.4-11}M_{\odot}$, of which 40 per cent are barred. Of the detected bars, more than half are strong. We find that the fraction of barred galaxies increases with stell… ▽ More We analyse the properties of strongly barred disc galaxies using the TNG100 simulation, a cosmological hydrodynamical realisation of the IllustrisTNG suite. We identify 270 disc galaxies at $z=0$ in the stellar mass range $M_{*}=10^{10.4-11}M_{\odot}$, of which 40 per cent are barred. Of the detected bars, more than half are strong. We find that the fraction of barred galaxies increases with stellar mass, in agreement with observational results. Strongly barred galaxies exhibit, overall, lower gas-to-stellar mass ratios compared to unbarred galaxies. The majority of barred galaxies are quenched (sSFR $\sim10^{-11.7} $yr$^{-1}$), whereas unbarred galaxies continue to be active (sSFR $\sim10^{-10.3}$yr$^{-1}$) on the main sequence of star-forming galaxies. We explore the evolution of strongly barred and unbarred galaxies to investigate their formation and quenching histories. We find that strong bars form between $0.5< z< 1.5$, with more massive galaxies hosting older bars. Strong bars form in galaxies with an early-established prominent disc component, undergoing periods of enhanced star formation and black hole accretion, possibly assisted by cosmological inflows. Unbarred galaxies, on the other hand, assemble most of their mass and disc component at late times. The nuclear region of strongly barred galaxies quenches shortly after bar formation, while unbarred galaxies remain active across time. Our findings are indicative of bar quenching, possibly assisted by nuclear feedback processes. We conclude that the cosmological environment, together with small scale feedback processes, determine the chances of a galaxy to form a bar and to rapidly quench its central region. △ Less

Submitted 30 March, 2022; v1 submitted 1 August, 2019; originally announced August 2019.

Comments: 18 pages, 15 figures, and 1 table. Accepted in MNRAS. Corrected pdf compilation

ALMA and HST kiloparsec-scale imaging of a quasar-galaxy merger at $z\approx 6.2$

Authors: Roberto Decarli, Massimo Dotti, Eduardo Bañados, Emanuele Paolo Farina, Fabian Walter, Chris Carilli, Xiaohui Fan, Chiara Mazzucchelli, Marcel Neeleman, Mladen Novak, Dominik Riechers, Michael A. Strauss, Bram P. Venemans, Yujin Yang, Ran Wang

Abstract: We present kpc-scale ALMA and HST imaging of the quasar PJ308-21 at $z$=$6.2342$, tracing dust, gas (via the [CII] 158 $μ$m line) and young stars. At a resolution of $\sim0.3''$ ($\approx1.7$ kpc), the system is resolved over $>4''$ ($>$20 kpc). In particular, it features a main component, identified to be the quasar host galaxy, centered on the accreting supermassive black hole; and two other ext… ▽ More We present kpc-scale ALMA and HST imaging of the quasar PJ308-21 at $z$=$6.2342$, tracing dust, gas (via the [CII] 158 $μ$m line) and young stars. At a resolution of $\sim0.3''$ ($\approx1.7$ kpc), the system is resolved over $>4''$ ($>$20 kpc). In particular, it features a main component, identified to be the quasar host galaxy, centered on the accreting supermassive black hole; and two other extended components on the West and East side, one redshifted and the other blueshifted relative to the quasar. The [CII] emission of the entire system stretches over $>$1500 km/s along the line of sight. All the components of the system are observed in dust, [CII], and rest-frame UV emission. The inferred [CII] luminosities [(0.9-4.6)$\times 10^9$ L$_\odot$], dust luminosities [(0.15-2.6)$\times10^{12}$ L$_\odot$], and rest-frame UV luminosities [(6.6-15)$\times10^{10}$ L$_\odot$], their ratios, and the implied gas/dust masses and star formation rates [11-290 M$_\odot$ yr$^{-1}$] are typical of high-redshift star-forming galaxies. A toy model of a single satellite galaxy that is tidally stripped by the interaction with the quasar host galaxy can account for the observed velocity and spatial extent of the two extended components. An outflow interpretation of the unique features in PJ308-21 is not supported by the data. PJ308-21 is thus one of the earliest galaxy mergers imaged at cosmic dawn. △ Less

Submitted 12 June, 2019; originally announced June 2019.

Comments: 9 pages, 6 figures. Accepted for publication in ApJ

Testing the blast-wave AGN feedback scenario in MCG-03-58-007

Authors: Mattia Sirressi, Claudia Cicone, Paola Severgnini, Valentina Braito, Massimo Dotti, Roberto Della Ceca, James Reeves, Gabriele Matzeu, Cristian Vignali, Lucia Ballo

Abstract: We report the first Atacama large millimeter/submillimeter array observations of MCG-03-58-007, a local ($z=0.03236\pm0.00002$, this work) AGN ($L_{AGN}\sim10^{45}~\rm erg~s^{-1}$), hosting a powerful X-ray ultra-fast ($v=0.1c$) outflow (UFO). The CO(1-0) line emission is observed across $\sim18\,$kpc scales with a resolution of $\sim 1\,\rm kpc$. About 78\% of the CO(1-0) luminosity traces a gala… ▽ More We report the first Atacama large millimeter/submillimeter array observations of MCG-03-58-007, a local ($z=0.03236\pm0.00002$, this work) AGN ($L_{AGN}\sim10^{45}~\rm erg~s^{-1}$), hosting a powerful X-ray ultra-fast ($v=0.1c$) outflow (UFO). The CO(1-0) line emission is observed across $\sim18\,$kpc scales with a resolution of $\sim 1\,\rm kpc$. About 78\% of the CO(1-0) luminosity traces a galaxy-size rotating disk. However, after subtracting the emission due to such rotating disk, we detect with a S/N=20 a residual emission in the central $\sim 4\,$kpc. Such residuals may trace a low velocity ($v_{LOS}=170\,\rm km\,s^{-1}$) outflow. We compare the momentum rate and kinetic power of such putative molecular outflow with that of the X-ray UFO and find $\dot{P}_{out}/\dot{P}_{UFO}=0.3\pm0.2$ and $\dot{E}_{mol}/\dot{E}_{UFO}\sim4\cdot10^{-3}$. This result is at odds with the energy-conserving scenario suggested by the large momentum boosts measured in some other molecular outflows. An alternative interpretation of the residual CO emission would be a compact rotating structure, distinct from the main disk, which would be a factor of $\sim10-100$ more extended and massive than typical circumnuclear disks revealed in Seyferts. However, in both scenarios, our results rule out the hypothesis of a momentum-boosted molecular outflow in this AGN, despite the presence of a powerful X-ray UFO. [Abridged] △ Less

Submitted 10 August, 2019; v1 submitted 3 June, 2019; originally announced June 2019.

Comments: Accepted for publication in MNRAS. 13 pages, 11 figures

Neutron star binary orbits in their host potential: effect on early r-process enrichment

Authors: Matteo Bonetti, Albino Perego, Massimo Dotti, Gabriele Cescutti

Abstract: Coalescing neutron star binary (NSB) systems are primary candidates for $r$-process enrichment of galaxies. The recent detection of $r$-process elements in ultra-faint dwarf (UFD) galaxies and the abundances measured in classical dwarfs challenges the NSB merger scenario both in terms of coalescence time scales and merger locations. In this paper, we focus on the dynamics of NSBs in the gravitatio… ▽ More Coalescing neutron star binary (NSB) systems are primary candidates for $r$-process enrichment of galaxies. The recent detection of $r$-process elements in ultra-faint dwarf (UFD) galaxies and the abundances measured in classical dwarfs challenges the NSB merger scenario both in terms of coalescence time scales and merger locations. In this paper, we focus on the dynamics of NSBs in the gravitational potentials of different types of host galaxies and on its impact on the subsequent galactic enrichment. We find that, for a $\sim t^{-1}$ delay time distribution, even when receiving a low kick ($\sim 10~{\rm km~s^{-1}}$) from the second supernova explosion, in shallow dwarf galaxy potentials NSBs tend to merge with a large off-set from the host galaxy. This results in a significant geometrical dilution of the amount of produced $r-$process elements that fall back and pollute the host galaxy gas reservoir. The combination of dilution and small number statistics produces a large scatter in the expected $r$-process enrichment within a single UFD or classical dwarf galaxy. Comparison between our results and observed europium abundances reveals a tension that even a systematic choice of optimistic parameters in our models cannot release. Such a discrepancy could point to the need of additional $r$-process production sites that suffer less severe dilution or to a population of extremely fast merging binaries. △ Less

Submitted 10 September, 2019; v1 submitted 28 May, 2019; originally announced May 2019.

Comments: 18 pages, 11 figures, MNRAS accepted