Showing 1–50 of 79 results for author: Lupi, A

Exploring Active Galactic Nuclei and Little Red Dots with the Obelisk simulation

Authors: M. Volonteri, M. Trebitsch, Y. Dubois, J. E. Greene, C. -A. Dong-Paez, M. Habouzit, A. Lupi, Y. Ma, R. S. Beckmann, P. Dayal

Abstract: The James Webb Space telescope has discovered an abundant population of broad line emitters, typical signposts for Active Galactic Nuclei (AGN). Many of these sources have red colors and a compact appearance that has led to naming them `Little Red Dots'. In this paper we develop a detailed framework to estimate the photometry of AGN embedded in galaxies extracted from the Obelisk cosmological simu… ▽ More The James Webb Space telescope has discovered an abundant population of broad line emitters, typical signposts for Active Galactic Nuclei (AGN). Many of these sources have red colors and a compact appearance that has led to naming them `Little Red Dots'. In this paper we develop a detailed framework to estimate the photometry of AGN embedded in galaxies extracted from the Obelisk cosmological simulation to understand the properties of color-selected Little Red Dots (cLRDs) in the context of the full AGN and massive black hole population. We find that using realistic spectral energy distributions (SEDs) and attenuation for AGN we can explain the shape of the cLRD SED as long as galaxies host a sufficiently luminous AGN that is not too much or too little attenuated. When attenuation is too low or too high, AGN do not enter the cLRD selection, because the AGN dominates over the host galaxy too much in blue filters, or it does not contribute to photometry anywhere, respectively. cLRDs are also characterized by high Eddington ratios, possibility super-Eddington, and/or high ratios between black hole and stellar mass. △ Less

Submitted 23 August, 2024; originally announced August 2024.

Comments: Submitted to A&A

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

Authors: Santi Roca-Fàbrega, Ji-hoon Kim, Joel R. Primack, Anna Genina, Minyong Jung, Alessandro Lupi, Kentaro Nagamine, Johnny W. Powell, Thomas R. Quinn, Yves Revaz, Ikkoh Shimizu, Héctor Velázquez, the AGORA Collaboration

Abstract: The AGORA Cosmorun (arXiv:2106.09738) is a set of hydrodynamical cosmological zoom-in simulations carried out within the AGORA High-resolution Galaxy Simulations Comparison Project (arXiv:1308.2669,arXiv:1610.03066). These simulations show the formation and evolution of a Milky Way-sized galaxy using eight of the most widely used numerical codes in the community (Art-I, Enzo, Ramses, Changa, Gadge… ▽ More The AGORA Cosmorun (arXiv:2106.09738) is a set of hydrodynamical cosmological zoom-in simulations carried out within the AGORA High-resolution Galaxy Simulations Comparison Project (arXiv:1308.2669,arXiv:1610.03066). These simulations show the formation and evolution of a Milky Way-sized galaxy using eight of the most widely used numerical codes in the community (Art-I, Enzo, Ramses, Changa, Gadget-3, Gear, Gizmo, and Arepo). In this short report, we describe the public release of the raw output data from all of these simulations at z = 8, 7, 6, 5, 4, 3, 2 (plus at z=1, 0 when available), and several metadata files containing the halo centers, virial quantities, and merger trees. The data from even thinner timesteps will be released as soon as the upcoming collaboration papers (VII-IX) are submitted and accepted. △ Less

Submitted 1 August, 2024; originally announced August 2024.

Comments: ArXiv Data Release from the AGORA Collaboration. arXiv admin note: text overlap with arXiv:2001.04354

A SPectroscopic survey of biased halos In the Reionization Era (ASPIRE): Broad-line AGN at $z=4-5$ revealed by JWST/NIRCam WFSS

Authors: Xiaojing Lin, Feige Wang, Xiaohui Fan, Zheng Cai, Jaclyn B. Champagne, Fengwu Sun, Marta Volonteri, Jinyi Yang, Joseph F. Hennawi, Eduardo Bañados, Aaron Barth, Anna-Christina Eilers, Emanuele Paolo Farina, Weizhe Liu, Xiangyu Jin, Hyunsung D. Jun, Alessandro Lupi, Koki Kakiichi, Chiara Mazzucchelli, Masafusa Onoue, Zhiwei Pan, Elia Pizzati, Sofía Rojas-Ruiz, Jan-Torge Schindler, Benny Trakhtenbrot , et al. (11 additional authors not shown)

Abstract: Low-luminosity AGNs with low-mass black holes (BHs) in the early universe are fundamental to understanding the BH growth and their co-evolution with the host galaxies. Utilizing JWST NIRCam Wide Field Slitless Spectroscopy (WFSS), we perform a systematic search for broad-line ${\rm Hα}$ emitters (BHAEs) at $z\approx 4-5$ in 25 fields of the ASPIRE (A SPectroscopic survey of biased halos In the Rei… ▽ More Low-luminosity AGNs with low-mass black holes (BHs) in the early universe are fundamental to understanding the BH growth and their co-evolution with the host galaxies. Utilizing JWST NIRCam Wide Field Slitless Spectroscopy (WFSS), we perform a systematic search for broad-line ${\rm Hα}$ emitters (BHAEs) at $z\approx 4-5$ in 25 fields of the ASPIRE (A SPectroscopic survey of biased halos In the Reionization Era) project, covering a total area of 275 arcmin$^2$. We identify 16 BHAEs with FWHM of the broad components spanning from $\sim$ 1000 km s$^{-1}$ to 3000 km s$^{-1}$. Assuming the broad linewidths arise due to Doppler broadening around BHs, the implied BH masses range from $10^7$ to $10^{8}~M_\odot$, with broad ${\rm Hα}$-converted bolometric luminosity of $10^{44.5}-10^{45.5}$ erg s$^{-1}$ and Eddington ratios of $0.07-0.47$. The spatially extended structure of the F200W stacked image may trace the stellar light from the host galaxies. The ${\rm Hα}$ luminosity function indicates an increasing AGN fraction towards the higher ${\rm Hα}$ luminosities. We find possible evidence for clustering of BHAEs: two sources are at the same redshift with a projected separation of 519 kpc; one BHAE appears as a composite system residing in an overdense region with three close companion ${\rm Hα}$ emitters. Three BHAEs exhibit blueshifted absorption troughs indicative of the presence of high-column-density gas. We find the broad-line and photometrically selected BHAE samples exhibit different distributions in the optical continuum slopes, which can be attributed to their different selection methods. The ASPIRE broad-line ${\rm Hα}$ sample provides a good database for future studies of faint AGN populations at high redshift. △ Less

Submitted 24 July, 2024; originally announced July 2024.

Comments: 19 pages, 13 figures, 4 tables. Accepted by the ApJ

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

Supermassive black hole formation via collisions in black hole clusters

Authors: Benjamin Gaete, Dominik R. G. Schleicher, Alessandro Lupi, Bastian Reinoso, Michael Fellhauer, Marcelo C. Vergara

Abstract: More than 300 supermassive black holes have been detected at redshifts larger than 6, and they are abundant in the centers of local galaxies. Their formation mechanisms, however, are still rather unconstrained. A possible origin of these supermassive black holes could be through mergers in dense black hole clusters, forming as a result of mass segregation within nuclear star clusters in the center… ▽ More More than 300 supermassive black holes have been detected at redshifts larger than 6, and they are abundant in the centers of local galaxies. Their formation mechanisms, however, are still rather unconstrained. A possible origin of these supermassive black holes could be through mergers in dense black hole clusters, forming as a result of mass segregation within nuclear star clusters in the center of galaxies. In this study, we present the first systematic investigation of the evolution of such black hole clusters where the effect of an external potential is taken into account. Such a potential could be the result of gas inflows into the central region, for example as a result of galaxy mergers. We show here that the efficiency for the formation of a massive central object is mostly regulated by the ratio of cluster velocity dispersion divided by the speed of light, potentially reaching efficiencies of 0.05-0.08 in realistic systems. Our results show that this scenario is potentially feasible and may provide seeds black hole of at least 10^3 solar masses. We conclude that the formation of seed black holes via this channel should be taken into account in statistical assessments of the black hole population. △ Less

Submitted 18 June, 2024; originally announced June 2024.

Comments: 15 pages, 11 figures, 1 table. Accepted for publication with A&A

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

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)

The Santa Barbara Binary-Disk Code Comparison

Authors: Paul C. Duffell, Alexander J. Dittmann, Daniel J. D'Orazio, Alessia Franchini, Kaitlin M. Kratter, Anna B. T. Penzlin, Enrico Ragusa, Magdalena Siwek, Christopher Tiede, Haiyang Wang, Jonathan Zrake, Adam M. Dempsey, Zoltan Haiman, Alessandro Lupi, Michal Pirog, Geoffrey Ryan

Abstract: We have performed numerical calculations of a binary interacting with a gas disk, using eleven different numerical methods and a standard binary-disk setup. The goal of this study is to determine whether all codes agree on a numerically converged solution, and to determine the necessary resolution for convergence and the number of binary orbits that must be computed to reach an agreed-upon relaxed… ▽ More We have performed numerical calculations of a binary interacting with a gas disk, using eleven different numerical methods and a standard binary-disk setup. The goal of this study is to determine whether all codes agree on a numerically converged solution, and to determine the necessary resolution for convergence and the number of binary orbits that must be computed to reach an agreed-upon relaxed state of the binary-disk system. We find that all codes can agree on a converged solution (depending on the diagnostic being measured). The zone spacing required for most codes to reach a converged measurement of the torques applied to the binary by the disk is roughly 1% of the binary separation in the vicinity of the binary components. For our disk model to reach a relaxed state, codes must be run for at least 200 binary orbits, corresponding to about a viscous time for our parameters, $0.2 (a^2 Ω_B /ν)$ binary orbits, where $ν$ is the kinematic viscosity. We did not investigate dependence on binary mass ratio, eccentricity, disk temperature, or disk viscosity; therefore, these benchmarks may act as guides towards expanding converged solutions to the wider parameter space but might need to be updated in a future study that investigates dependence on system parameters. We find the most major discrepancies between codes resulted from the dimensionality of the setup (3D vs 2D disks). Beyond this, we find good agreement in the total torque on the binary between codes, although the partition of this torque between the gravitational torque, orbital accretion torque, and spin accretion torque depends sensitively on the sink prescriptions employed. In agreement with previous studies, we find a modest difference in torques and accretion variability between 2D and 3D disk models. We find cavity precession rates to be appreciably faster in 3D than in 2D. △ Less

Submitted 20 June, 2024; v1 submitted 20 February, 2024; originally announced February 2024.

Comments: ApJ Accepted

Testing analytical methods to derive the cosmic-ray ionisation rate in cold regions via synthetic observations

Authors: E. Redaelli, S. Bovino, A. Lupi, T. Grassi, D. Gaete-Espinoza, G. Sabatini, P. Caselli

Abstract: Cosmic rays (CRs) heavily impact the chemistry and physics of cold and dense star-forming regions. However, characterising their ionisation rate is still challenging from an observational point of view. In the past, a few analytical formulas have been proposed to infer the cosmic-ray ionization rate $ζ_2$ from molecular line observations. These have been derived from the chemical kinetics of the i… ▽ More Cosmic rays (CRs) heavily impact the chemistry and physics of cold and dense star-forming regions. However, characterising their ionisation rate is still challenging from an observational point of view. In the past, a few analytical formulas have been proposed to infer the cosmic-ray ionization rate $ζ_2$ from molecular line observations. These have been derived from the chemical kinetics of the involved species, but they have not been validated using synthetic data processed with a standard observative pipeline. We aim to bridge this gap. We perform the radiative transfer on a set of three-dimensional magneto-hydrodynamical simulations of prestellar cores, exploring different initial $ζ_2$, evolutionary stages, types of radiative transfer (e.g. assuming local-thermodynamic-equilibrium conditions), and telescope responses. We then compute the column densities of the involved tracers to determine $ζ_2$, using, in particular, the equation proposed by Bovino et. al (2020) and by Caselli et al. (1998) both used nowadays. Our results confirm that the method of Bovino et al. (2020) accurately retrieves the actual $ζ_2$ within a factor of $2-3$, in the physical conditions explored in our tests. Since we also explore a non-local thermodynamic equilibrium radiative transfer, this work indirectly offers insights into the excitation temperatures of common transitions at moderate volume densities ($n\approx 10^5 \, \rm cm^{-3}$). We have also performed a few tests using the formula proposed by Caselli et al. (1998), which overestimates the actual $ζ_2$ by at least two orders of magnitudes. We also consider a new derivation of this method, which, however, still leads to large overestimates. △ Less

Submitted 16 February, 2024; originally announced February 2024.

Comments: Accepted for publication in A&A (15th Feb.)

The AGORA High-resolution Galaxy Simulations Comparison Project IV: Halo and Galaxy Mass Assembly in a Cosmological Zoom-in Simulation at $z\le2$

Authors: Santi Roca-Fàbrega, Ji-hoon Kim, Joel R. Primack, Minyong Jung, Anna Genina, Loic Hausammann, Hyeonyong Kim, Alessandro Lupi, Kentaro Nagamine, Johnny W. Powell, Yves Revaz, Ikkoh Shimizu, Clayton Strawn, Héctor Velázquez, Tom Abel, Daniel Ceverino, Bili Dong, Thomas R. Quinn, Eun-jin Shin, Alvaro Segovia-Otero, Oscar Agertz, Kirk S. S. Barrow, Corentin Cadiou, Avishai Dekel, Cameron Hummels , et al. (3 additional authors not shown)

Abstract: In this fourth paper from the AGORA Collaboration, we study the evolution down to redshift $z=2$ and below of a set of cosmological zoom-in simulations of a Milky Way mass galaxy by eight of the leading hydrodynamic simulation codes. We also compare this CosmoRun suite of simulations with dark matter-only simulations by the same eight codes. We analyze general properties of the halo and galaxy at… ▽ More In this fourth paper from the AGORA Collaboration, we study the evolution down to redshift $z=2$ and below of a set of cosmological zoom-in simulations of a Milky Way mass galaxy by eight of the leading hydrodynamic simulation codes. We also compare this CosmoRun suite of simulations with dark matter-only simulations by the same eight codes. We analyze general properties of the halo and galaxy at $z=4$ and 3, and before the last major merger, focusing on the formation of well-defined rotationally-supported disks, the mass-metallicity relation, the specific star formation rate, the gas metallicity gradients, and the non-axisymmetric structures in the stellar disks. Codes generally converge well to the stellar-to-halo mass ratios predicted by semi-analytic models at $z\sim$2. We see that almost all the hydro codes develop rotationally-supported structures at low redshifts. Most agree within 0.5 dex with the observed MZR at high and intermediate redshifts, and reproduce the gas metallicity gradients obtained from analytical models and low-redshift observations. We confirm that the inter-code differences in the halo assembly history reported in the first paper of the collaboration also exist in CosmoRun, making the code-to-code comparison more difficult. We show that such differences are mainly due to variations in code-dependent parameters that control the time-stepping strategy of the gravity solver. We find that variations in the early stellar feedback can also result in differences in the timing of the low-redshift mergers. All the simulation data down to $z=2$ and the auxiliary data will be made publicly available. △ Less

Submitted 9 February, 2024; originally announced February 2024.

Comments: Submitted to ApJ. Visit the AGORA Collaboration website (www.agorasimulations.org <http://www.agorasimulations.org/>) for more information

The AGORA High-resolution Galaxy Simulations Comparison Project. V: Satellite Galaxy Populations In A Cosmological Zoom-in Simulation of A Milky Way-mass Halo

Authors: Minyong Jung, Santi Roca-Fàbrega, Ji-hoon Kim, Anna Genina, Loic Hausammann, Hyeonyong Kim, Alessandro Lupi, Kentaro Nagamine, Johnny W. Powell, Yves Revaz, Ikkoh Shimizu, Héctor Velázquez, Daniel Ceverino, Joel R. Primack, Thomas R. Quinn, Clayton Strawn, Tom Abel, Avishai Dekel, Bili Dong, Boon Kiat Oh, Romain Teyssier

Abstract: We analyze and compare the satellite halo populations at $z\sim2$ in the high-resolution cosmological zoom-in simulations of a $10^{12}\,{\rm M}_{\odot}$ target halo ($z=0$ mass) carried out on eight widely-used astrophysical simulation codes ({\sc Art-I}, {\sc Enzo}, {\sc Ramses}, {\sc Changa}, {\sc Gadget-3}, {\sc Gear}, {\sc Arepo-t}, and {\sc Gizmo}) for the {\it AGORA} High-resolution Galaxy… ▽ More We analyze and compare the satellite halo populations at $z\sim2$ in the high-resolution cosmological zoom-in simulations of a $10^{12}\,{\rm M}_{\odot}$ target halo ($z=0$ mass) carried out on eight widely-used astrophysical simulation codes ({\sc Art-I}, {\sc Enzo}, {\sc Ramses}, {\sc Changa}, {\sc Gadget-3}, {\sc Gear}, {\sc Arepo-t}, and {\sc Gizmo}) for the {\it AGORA} High-resolution Galaxy Simulations Comparison Project. We use slightly different redshift epochs near $z=2$ for each code (hereafter ``$z\sim2$') at which the eight simulations are in the same stage in the target halo's merger history. After identifying the matched pairs of halos between the {\it CosmoRun} simulations and the DMO simulations, we discover that each {\it CosmoRun} halo tends to be less massive than its DMO counterpart. When we consider only the halos containing stellar particles at $z\sim2$, the number of satellite {\it galaxies} is significantly fewer than that of dark matter halos in all participating {\it AGORA} simulations, and is comparable to the number of present-day satellites near the Milky Way or M31. The so-called ``missing satellite problem' is fully resolved across all participating codes simply by implementing the common baryonic physics adopted in {\it AGORA} and the stellar feedback prescription commonly used in each code, with sufficient numerical resolution ($\lesssim100$ proper pc at $z=2$). We also compare other properties such as the stellar mass$-$halo mass relation and the mass$-$metallicity relation. Our work highlights the value of comparison studies such as {\it AGORA}, where outstanding problems in galaxy formation theory are studied simultaneously on multiple numerical platforms. △ Less

Submitted 7 February, 2024; originally announced February 2024.

Comments: Accepted for publication to the ApJ, 19 pages, 10 figures. Visit the AGORA Collaboration website (www.agorasimulations.org) for more information. For summary video, please see https://www.youtube.com/watch?v=HShRJYvPe1k

Journal ref: ApJ 964 123 (2024)

The AGORA High-resolution Galaxy Simulations Comparison Project. VI. Similarities and Differences in the Circumgalactic Medium

Authors: Clayton Strawn, Santi Roca-Fàbrega, Joel R. Primack, Ji-hoon Kim, Anna Genina, Loic Hausammann, Hyeonyong Kim, Alessandro Lupi, Kentaro Nagamine, Johnny W. Powell, Yves Revaz, Ikkoh Shimizu, Héctor Velázquez, Tom Abel, Daniel Ceverino, Bili Dong, Minyong Jung, Thomas R. Quinn, Eun-jin Shin, Kirk S. S. Barrow, Avishai Dekel, Boon Kiat Oh, Nir Mandelker, Romain Teyssier, Cameron Hummels , et al. (4 additional authors not shown)

Abstract: We analyze the circumgalactic medium (CGM) for eight commonly-used cosmological codes in the AGORA collaboration. The codes are calibrated to use identical initial conditions, cosmology, heating and cooling, and star formation thresholds, but each evolves with its own unique code architecture and stellar feedback implementation. Here, we analyze the results of these simulations in terms of the str… ▽ More We analyze the circumgalactic medium (CGM) for eight commonly-used cosmological codes in the AGORA collaboration. The codes are calibrated to use identical initial conditions, cosmology, heating and cooling, and star formation thresholds, but each evolves with its own unique code architecture and stellar feedback implementation. Here, we analyze the results of these simulations in terms of the structure, composition, and phase dynamics of the CGM. We show properties such as metal distribution, ionization levels, and kinematics are effective tracers of the effects of the different code feedback and implementation methods, and as such they can be highly divergent between simulations. This is merely a fiducial set of models, against which we will in the future compare multiple feedback recipes for each code. Nevertheless, we find that the large parameter space these simulations establish can help disentangle the different variables that affect observable quantities in the CGM, e.g., showing that abundances for ions with higher ionization energy are more strongly determined by the simulation's metallicity, while abundances for ions with lower ionization energy are more strongly determined by the gas density and temperature. △ Less

Submitted 7 February, 2024; originally announced February 2024.

Comments: 32 pages, 17 figures. Published in ApJ

Journal ref: ApJ 962 29 (2024)

Electromagnetic signatures from accreting massive black hole binaries in time domain photometric surveys

Authors: Fabiola Cocchiararo, Alessia Franchini, Alessandro Lupi, Alberto Sesana

Abstract: We study spectral and time variability of accreting massive black hole binaries (MBHBs) at milli-pc separations surrounded by a geometrically thin circumbinary disc. We present the first computed spectral energy distribution (SED) and light curves (LCs) from 3D hyper-Lagrangian resolution hydrodynamic simulations of these systems. We model binaries with mass of $10^6$ M$_\odot$, eccentricities e=0… ▽ More We study spectral and time variability of accreting massive black hole binaries (MBHBs) at milli-pc separations surrounded by a geometrically thin circumbinary disc. We present the first computed spectral energy distribution (SED) and light curves (LCs) from 3D hyper-Lagrangian resolution hydrodynamic simulations of these systems. We model binaries with mass of $10^6$ M$_\odot$, eccentricities e=0, 0.9 and mass ratio q=0.1,1. The circumbinary disc has initial aspect ratio of 0.1, features an adiabatic equation of state, and evolves under the effect of viscous heating, black body cooling and self gravity. To compute the SED, we consider black body emission from each disc element and we add an X-ray corona with luminosity proportional to that of the mini-discs around each black hole. We find significant variability of the SED, especially at high energies, which translates into LCs displaying modulations of a factor of ~ 2 in optical and of ~ 10 in UV and X-rays. We focus on flux variability in the optical band which will be probed by the Vera Rubin Observatory (VRO). Modulations on the orbital period and half of the orbital period are evident in all systems. In equal mass binaries, we find another longer timescale modulation, linked to an over-density forming at the inner edge of the disc. Considering the VRO properties, we find that equal mass, circular binaries are unlikely to be identified, due to the lack of prominent peaks in their Fourier spectra. Conversely, unequal mass and/or eccentric binaries can be singled out up to z ~ 0.5 (for systems with $L_{\rm bol}\approx10^{42}$ erg s$^{-1}$) and z ~ 2 (for systems with $L_{\rm bol}\approx10^{44}$ erg s$^{-1}$). Identifying electromagnetic signatures of MBHBs at separations $\sim 10^{-4}-10^{-2}$ pc is crucial for understanding the physics of the future Laser Interferometer Space Antenna sources and the origin of the GW background. △ Less

Submitted 7 February, 2024; originally announced February 2024.

Comments: 14 pages, 11 figures, A&A submitted

Emission signatures from sub-pc Post-Newtonian binaries embedded in circumbinary discs

Authors: Alessia Franchini, Matteo Bonetti, Alessandro Lupi, Alberto Sesana

Abstract: We study the dynamical evolution of quasi-circular equal mass massive black hole binaries embedded in circumbinary discs from separations of $\sim 100R_{\rm g}$ down to the merger, following the post merger evolution. The binary orbit evolves owing to the presence of the gaseous disc and the addition of Post-Newtonian (PN) corrections up to the 2.5 PN order, therefore including the dissipative gra… ▽ More We study the dynamical evolution of quasi-circular equal mass massive black hole binaries embedded in circumbinary discs from separations of $\sim 100R_{\rm g}$ down to the merger, following the post merger evolution. The binary orbit evolves owing to the presence of the gaseous disc and the addition of Post-Newtonian (PN) corrections up to the 2.5 PN order, therefore including the dissipative gravitational wave back-reaction. We investigate two cases of a relatively cold and warm circumbinary discs, with aspect ratios $H/R=0.03,\,0.1$ respectively, employing 3D hyper-Lagrangian resolution simulations with the {\sc gizmo}-MFM code. We extract spectral energy distributions and light curves in different frequency bands (i.e. X-ray, optical and UV) from the simulations. We find a clear two orders of magnitude drop in the X-ray flux right before merger if the disc is warm while we identify a significant increase in the UV flux regardless of the disc temperature. The optical flux shows clear distinctive modulations on the binary orbital period and on the cavity edge period, regardless of the disc temperature. We find that the presence of a cold disc can accelerate the coalescence of the binary by up to 130 seconds over the last five days of inspiral, implying a phase shift accumulation of about $π\,$radians compared to the binary evolution in vacuum. These differences are triggered by the presence of the gaseous disc and might have implications on the waveforms that can be in principle detected. We discuss the implications that these distinctive signatures might have for existing and upcoming time domain surveys and for multimessenger astronomy. △ Less

Submitted 18 January, 2024; originally announced January 2024.

Comments: 10 pages, A&A submitted

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

Sustained super-Eddington accretion in high-redshift quasars

Authors: Alessandro Lupi, Giada Quadri, Marta Volonteri, Monica Colpi, John A. Regan

Abstract: Observations of $z \gtrsim 6$ quasars provide information on the early evolution of the most massive black holes (MBHs) and galaxies. Current observations, able to trace both gas and stellar properties, reveal a population of MBHs that is significantly more massive than expected from the local MBH-stellar mass relation. The population lies on, but mostly above, the relation observed in the nearby… ▽ More Observations of $z \gtrsim 6$ quasars provide information on the early evolution of the most massive black holes (MBHs) and galaxies. Current observations, able to trace both gas and stellar properties, reveal a population of MBHs that is significantly more massive than expected from the local MBH-stellar mass relation. The population lies on, but mostly above, the relation observed in the nearby Universe. This suggests that these objects grew very rapidly. To explain their presence when the Universe was less than 1~Gyr old and to assess the physical conditions for their rapid growth, we explored whether episodes of accretion above the Eddington limit can occur across cosmic epochs. By employing state-of-the-art high-resolution cosmological zoom-in simulations of a $z\sim 7$ quasar, where different accretion regimes are included consistently, together with their associated radiative and kinetic feedback, we show that super-Eddington phases can be sustained for relatively long timescales (tens of millions of years). This allows the MBH to rapidly grow by up to three orders of magnitude, depending on the strength of the kinetic feedback. We also show by means of a semianalytic calculation that the MBH spin remains moderate and does not take on extremely high values during the super-Eddington phases. This results in a lower feedback efficiency, which may allow the rapid growth required to explain over-massive high-redshift MBHs. △ Less

Submitted 24 May, 2024; v1 submitted 13 December, 2023; originally announced December 2023.

Comments: 8 pages, 5 figures, matches the published version

Journal ref: A&A 686, A256 (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.

Massive Black Holes in Galactic Nuclei

Authors: David Izquierdo-Villalba, Alessandro Lupi, John Regan, Matteo Bonetti, Alessia Franchini

Abstract: Massive black holes are key inhabitants of the nuclei of galaxies. Moreover, their astrophysical relevance has gained significant traction in recent years, thanks especially to the amazing results that are being (or will be) delivered by instruments such as the James Webb Space Telescope, Pulsar Timing Array projects and LISA. In this Chapter, we aim to detail a broad set of aspects related to the… ▽ More Massive black holes are key inhabitants of the nuclei of galaxies. Moreover, their astrophysical relevance has gained significant traction in recent years, thanks especially to the amazing results that are being (or will be) delivered by instruments such as the James Webb Space Telescope, Pulsar Timing Array projects and LISA. In this Chapter, we aim to detail a broad set of aspects related to the astrophysical nature of massive black holes embedded in galactic nuclei, with a particular focus on recent and upcoming advances in the field. In particular, we will address questions such as: What shapes the relations connecting the mass of massive black holes with the properties of their host galaxies? How do massive black holes form in the early Universe? What mechanisms keep on feeding them so that they can attain very large masses at z = 0? How do binaries composed of two massive black holes form and coalesce into a single, larger black hole? Here we present these topics from a mainly theoretical viewpoint and discuss how present and upcoming facilities may enhance our understanding of massive black holes in the near future. △ Less

Submitted 6 November, 2023; originally announced November 2023.

Comments: All authors equally contributed to the Chapter writing. To appear as Chapter 3 in the book "Black Holes in the Era of Gravitational Wave Astronomy", ed. Arca Sedda, Bortolas, Spera, pub. Elsevier

Shaping the unseen: the influence of baryons and environment on low-mass, high-redshift dark matter haloes in the SIEGE simulations

Authors: R. Pascale, F. Calura, A. Lupi, J. Rosdahl, E. Lacchin, M. Meneghetti, C. Nipoti, E. Vanzella, E. Vesperini, A. Zanella

Abstract: We use zoom-in, hydrodynamical, cosmological $N$-body simulations tracing the formation of the first stellar clumps from the SImulating the Environments where Globular clusters Emerged (SIEGE) project, to study key structural properties of dark matter haloes when the Universe was only $0.92$ Gyr old. The very high-resolution (maximum physical resolution 0.3 h$^{-1}$ pc at $z=6.14$, smallest dark-m… ▽ More We use zoom-in, hydrodynamical, cosmological $N$-body simulations tracing the formation of the first stellar clumps from the SImulating the Environments where Globular clusters Emerged (SIEGE) project, to study key structural properties of dark matter haloes when the Universe was only $0.92$ Gyr old. The very high-resolution (maximum physical resolution 0.3 h$^{-1}$ pc at $z=6.14$, smallest dark-matter particle mass $164\,M_{\odot}$) allows us to reach the very low mass end of the stellar-to-halo mass relation ($M_{\rm vir}=10^{7.5-9.5}\,M_{\odot}$) to study the processes that mould dark matter haloes during the first stages of structure formation. We investigate the role of baryonic cooling and stellar feedback, modeled from individual stars, in shaping haloes, and of environmental effects as accretion of dark matter along cosmic filaments and mergers. We find that the onset of star formation (typically for $\log M_{\rm vir}/M_{\odot}\simeq7.6$) causes the inner cusp in the haloes density profile to flatten into a core with constant density and size proportionally to the halo virial mass. Even at these mass scales, we confirm that baryons make haloes that have formed stars rounder in the central regions than haloes that have not formed stars yet, with median minor-to-major $\langle q \rangle$ and intermediate-to-major $\langle s \rangle$ axes 0.66 and 0.84, respectively. Our morphological analysis shows that, at $z=6.14$, haloes are largely prolate in the outer parts, with the major axis aligned along filaments of the cosmic web or towards smaller sub-haloes, with the degree of elongation having no significant dependence on the halo mass. △ Less

Submitted 4 September, 2023; originally announced September 2023.

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

A SPectroscopic survey of biased halos In the Reionization Era (ASPIRE): JWST Reveals a Filamentary Structure around a z=6.61 Quasar

Authors: Feige Wang, Jinyi Yang, Joseph F. Hennawi, Xiaohui Fan, Fengwu Sun, Jaclyn B. Champagne, Tiago Costa, Melanie Habouzit, Ryan Endsley, Zihao Li, Xiaojing Lin, Romain A. Meyer, Jan-Torge Schindler, Yunjing Wu, Eduardo Bañados, Aaron J. Barth, Aklant K. Bhowmick, Rebekka Bieri, Laura Blecha, Sarah Bosman, Zheng Cai, Luis Colina, Thomas Connor, Frederick B. Davies, Roberto Decarli , et al. (34 additional authors not shown)

Abstract: We present the first results from the JWST ASPIRE program (A SPectroscopic survey of biased halos In the Reionization Era). This program represents an imaging and spectroscopic survey of 25 reionization-era quasars and their environments by utilizing the unprecedented capabilities of NIRCam Wide Field Slitless Spectroscopy (WFSS) mode. ASPIRE will deliver the largest ($\sim280~{\rm arcmin}^2$) gal… ▽ More We present the first results from the JWST ASPIRE program (A SPectroscopic survey of biased halos In the Reionization Era). This program represents an imaging and spectroscopic survey of 25 reionization-era quasars and their environments by utilizing the unprecedented capabilities of NIRCam Wide Field Slitless Spectroscopy (WFSS) mode. ASPIRE will deliver the largest ($\sim280~{\rm arcmin}^2$) galaxy redshift survey at 3-4 $μ$m among JWST Cycle-1 programs and provide extensive legacy values for studying the formation of the earliest supermassive black holes (SMBHs), the assembly of galaxies, early metal enrichment, and cosmic reionization. In this first ASPIRE paper, we report the discovery of a filamentary structure traced by the luminous quasar J0305-3150 and ten [OIII] emitters at $z=6.6$. This structure has a 3D galaxy overdensity of $δ_{\rm gal}=12.6$ over 637 cMpc$^3$, one of the most overdense structures known in the early universe, and could eventually evolve into a massive galaxy cluster. Together with existing VLT/MUSE and ALMA observations of this field, our JWST observations reveal that J0305-3150 traces a complex environment where both UV-bright and dusty galaxies are present, and indicate that the early evolution of galaxies around the quasar is not simultaneous. In addition, we discovered 31 [OIII] emitters in this field at other redshifts, $5.3<z<6.7$, with half of them situated at $z\sim5.4$ and $z\sim6.2$. This indicates that star-forming galaxies, such as [OIII] emitters, are generally clustered at high redshifts. These discoveries demonstrate the unparalleled redshift survey capabilities of NIRCam WFSS and the potential of the full ASPIRE survey dataset. △ Less

Submitted 19 April, 2023; originally announced April 2023.

Comments: accepted for publication in ApJL

A SPectroscopic survey of biased halos In the Reionization Era (ASPIRE): A First Look at the Rest-frame Optical Spectra of $z > 6.5$ Quasars Using JWST

Authors: Jinyi Yang, Feige Wang, Xiaohui Fan, Joseph F. Hennawi, Aaron J. Barth, Eduardo Bañados, Fengwu Sun, Weizhe Liu, Zheng Cai, Linhua Jiang, Zihao Li, Masafusa Onoue, Jan-Torge Schindler, Yue Shen, Yunjing Wu, Aklant K. Bhowmick, Rebekka Bieri, Laura Blecha, Sarah Bosman, Jaclyn B. Champagne, Luis Colina, Thomas Connor, Tiago Costa, Frederick B. Davies, Roberto Decarli , et al. (31 additional authors not shown)

Abstract: Studies of rest-frame optical emission in quasars at $z>6$ have historically been limited by the wavelengths accessible by ground-based telescopes. The James Webb Space Telescope (JWST) now offers the opportunity to probe this emission deep into the reionization epoch. We report the observations of eight quasars at $z>6.5$ using the JWST/NIRCam Wide Field Slitless Spectroscopy, as a part of the ''… ▽ More Studies of rest-frame optical emission in quasars at $z>6$ have historically been limited by the wavelengths accessible by ground-based telescopes. The James Webb Space Telescope (JWST) now offers the opportunity to probe this emission deep into the reionization epoch. We report the observations of eight quasars at $z>6.5$ using the JWST/NIRCam Wide Field Slitless Spectroscopy, as a part of the ''A SPectroscopic survey of biased halos In the Reionization Era (ASPIRE)" program. Our JWST spectra cover the quasars' emission between rest frame $\sim$ 4100 and 5100 Å. The profiles of these quasars' broad H$β$ emission lines span a FWHM from 3000 to 6000 $\rm{km~s^{-1}}$. The H$β$-based virial black hole (BH) masses, ranging from 0.6 to 2.1 billion solar masses, are generally consistent with their MgII-based BH masses. The new measurements based on the more reliable H$β$ tracer thus confirm the existence of billion solar-mass BHs in the reionization epoch. In the observed [OIII] $λλ$4960,5008 doublets of these luminous quasars, broad components are more common than narrow core components ($\le~1200~\rm{km~s^{-1}}$), and only one quasar shows stronger narrow components than broad. Two quasars exhibit significantly broad and blueshifted [OIII] emission, thought to trace galactic-scale outflows, with median velocities of $-610~\rm{km~s^{-1}}$ and $-1430~\rm{km~s^{-1}}$ relative to the [CII] $158\,μ$m line. All eight quasars show strong optical FeII emission, and follow the Eigenvector 1 relations defined by low-redshift quasars. The entire ASPIRE program will eventually cover 25 quasars and provide a statistical sample for the studies of the BHs and quasar spectral properties. △ Less

Submitted 19 April, 2023; originally announced April 2023.

Comments: 13 pages, 4 figures, accepted for publication in ApJL

The importance of live binary evolution in numerical simulations of binaries embedded in circumbinary discs

Authors: Alessia Franchini, Alessandro Lupi, Alberto Sesana, Zoltan Haiman

Abstract: The shrinking of a binary orbit driven by the interaction with a gaseous circumbinary disc, initially advocated as a potential way to catalyze the binary merger, has been recently debated in the case of geometrically thick (i.e. with $H/R\gtrsim 0.1$) discs. However, a clear consensus is still missing mainly owing to numerical limitations, such as fixed orbit binaries or lack of resolution inside… ▽ More The shrinking of a binary orbit driven by the interaction with a gaseous circumbinary disc, initially advocated as a potential way to catalyze the binary merger, has been recently debated in the case of geometrically thick (i.e. with $H/R\gtrsim 0.1$) discs. However, a clear consensus is still missing mainly owing to numerical limitations, such as fixed orbit binaries or lack of resolution inside the cavity carved by the binary in its circumbinary disc. In this work, we asses the importance of evolving the binary orbit by means of hydrodynamic simulations performed with the code {\sc gizmo} in meshless-finite-mass mode. In order to model the interaction between equal mass circular binaries and their locally isothermal circumbinary discs, we enforce hyper-Lagrangian resolution inside the cavity. We find that fixing the binary orbit ultimately leads to an overestimate of the gravitational torque that the gas exerts on the binary, and in an underestimate of the torque due to the accretion of material onto the binary components. Furthermore, we find that the modulation of the accretion rate on the binary orbital period is strongly suppressed in the fixed orbit simulation, while it is clearly present in the live binary simulations. This has potential implications for the prediction of the observable periodicities in massive black hole binary candidates. △ Less

Submitted 7 April, 2023; originally announced April 2023.

Comments: 7 pages, 8 figures, accepted for publication in 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)

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.

Non-ideal magnetohydrodynamics of self-gravitating filaments

Authors: Nicol Gutiérrez-Vera, Tommaso Grassi, Stefano Bovino, Alessandro Lupi, Daniele Galli, Dominik R. G. Schleicher

Abstract: Filaments have been studied in detail through observations and simulations. A range of numerical works have separately investigated how chemistry and diffusion effects, as well as magnetic fields and their structure impact the gas dynamics of the filament. However, non-ideal effects have hardly been explored thus far. We investigate how non-ideal magnetohydrodynamic (MHD) effects, combined with a… ▽ More Filaments have been studied in detail through observations and simulations. A range of numerical works have separately investigated how chemistry and diffusion effects, as well as magnetic fields and their structure impact the gas dynamics of the filament. However, non-ideal effects have hardly been explored thus far. We investigate how non-ideal magnetohydrodynamic (MHD) effects, combined with a simplified chemical model affect the evolution and accretion of a star-forming filament. We modeled an accreting self-gravitating turbulent filament using lemongrab, a one-dimensional (1D) non-ideal MHD code that includes chemistry. We explore the influence of non-ideal MHD, the orientation and strength of the magnetic field, and the cosmic ray ionization rate, on the evolution of the filament, with particular focus on the width and accretion rate. We find that the filament width and the accretion rate are determined by the magnetic field properties, including the initial strength, the coupling with the gas controlled by the cosmic ray ionization rate, and the orientation of the magnetic field with respect to the accretion flow direction. Increasing the cosmic-ray ionization rate leads to a behavior closer to that of ideal MHD, reducing the magnetic pressure support and, hence, damping the accretion efficiency with a consequent broadening of the filament width. For the same reason, we obtained a narrower width and a larger accretion rate when we reduced the initial magnetic field strength. Overall, while these factors affect the final results by approximately a factor of~2, removing the non-ideal MHD effects results in a much greater variation (up to a factor of~7). The inclusion of non-ideal MHD effects and the cosmic-ray ionization is crucial for the study of self-gravitating filaments and in determining critical observable quantities, such as the filament width and accretion rate. △ Less

Submitted 28 November, 2022; originally announced November 2022.

Comments: 11 pages, 9 figures, to be published in Astronomy & Astrophysics

Journal ref: A&A 670, A38 (2023)

Detection of stellar light from quasar host galaxies at redshifts above 6

Authors: Xuheng Ding, Masafusa Onoue, John D. Silverman, Yoshiki Matsuoka, Takuma Izumi, Michael A. Strauss, Knud Jahnke, Camryn L. Phillips, Junyao Li, Marta Volonteri, Zoltan Haiman, Irham Taufik Andika, Kentaro Aoki, Shunsuke Baba, Rebekka Bieri, Sarah E. I. Bosman, Connor Bottrell, Anna-Christina Eilers, Seiji Fujimoto, Melanie Habouzit, Masatoshi Imanishi, Kohei Inayoshi, Kazushi Iwasawa, Nobunari Kashikawa, Toshihiro Kawaguchi , et al. (19 additional authors not shown)

Abstract: The detection of starlight from the host galaxies of quasars during the reionization epoch ($z>6$) has been elusive, even with deep HST observations. The current highest redshift quasar host detected, at $z=4.5$, required the magnifying effect of a foreground lensing galaxy. Low-luminosity quasars from the Hyper Suprime-Cam Subaru Strategic Program (HSC-SSP) mitigate the challenge of detecting the… ▽ More The detection of starlight from the host galaxies of quasars during the reionization epoch ($z>6$) has been elusive, even with deep HST observations. The current highest redshift quasar host detected, at $z=4.5$, required the magnifying effect of a foreground lensing galaxy. Low-luminosity quasars from the Hyper Suprime-Cam Subaru Strategic Program (HSC-SSP) mitigate the challenge of detecting their underlying, previously-undetected host galaxies. Here we report rest-frame optical images and spectroscopy of two HSC-SSP quasars at $z>6$ with JWST. Using NIRCam imaging at 3.6$μ$m and 1.5$μ$m and subtracting the light from the unresolved quasars, we find that the host galaxies are massive (stellar masses of $13\times$ and $3.4\times$ $10^{10}$ M$_{\odot}$, respectively), compact, and disk-like. NIRSpec medium-resolution spectroscopy shows stellar absorption lines in the more massive quasar, confirming the detection of the host. Velocity-broadened gas in the vicinity of these quasars enables measurements of their black hole masses ($1.4\times 10^9$ and $2.0\times$ $10^{8}$ M$_{\odot}$, respectively). Their location in the black hole mass - stellar mass plane is consistent with the distribution at low redshift, suggesting that the relation between black holes and their host galaxies was already in place less than a billion years after the Big Bang. △ Less

Submitted 23 June, 2023; v1 submitted 25 November, 2022; originally announced November 2022.

Comments: Matched to the published Nature version of the article. 27 pages, 4 main figures, 1 table, 6 supplementary figures, 2 supplementary table

A general relativistic extension to mesh-free methods for hydrodynamics

Authors: Alessandro Lupi

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

Submitted 4 January, 2023; v1 submitted 11 October, 2022; originally announced October 2022.

Comments: 17 pages, 15 figures, matches the published version

Sub-parsec resolution cosmological simulations of star-forming clumps at high redshift with feedback of individual stars

Authors: F. Calura, A. Lupi, J. Rosdahl, E. Vanzella, M. Meneghetti, P. Rosati, E. Vesperini, E. Lacchin, R. Pascale, R. Gilli

Abstract: We introduce a new set of zoom-in cosmological simulations with sub-pc resolution, intended to model extremely faint, highly magnified star-forming stellar clumps, detected at z=6.14 thanks to gravitational lensing. The simulations include feedback from individual massive stars (in both the pre-supernova and supernova phases), generated via stochastic, direct sampling of the stellar initial mass f… ▽ More We introduce a new set of zoom-in cosmological simulations with sub-pc resolution, intended to model extremely faint, highly magnified star-forming stellar clumps, detected at z=6.14 thanks to gravitational lensing. The simulations include feedback from individual massive stars (in both the pre-supernova and supernova phases), generated via stochastic, direct sampling of the stellar initial mass function. We adopt a modified 'delayed cooling' feedback scheme, specifically created to prevent artificial radiative loss of the energy injected by individual stars in very dense gas (n~10^3-10^5 cm^{-3}). The sites where star formation ignites are characterised by maximum densities of the order of 10^5 cm^{-3} and gravitational pressures P/k>10^7 K/cm^3, corresponding to the values of the local, turbulent regions where the densest stellar aggregates form. The total stellar mass at z=6.14 is 3.4x10^7 M_sun, in satisfactory agreement with the observed stellar mass of the observed systems. The most massive clumps have masses of ~10^6 M_sun and half-mass sizes of ~100 pc. These sizes are larger than the observed ones, including also other samples of lensed high-redshift clumps, and imply an average density one order of magnitude lower than the observed one. In the size-mass plane, our clumps populate a sequence that is intermediate between the ones of observed high-redshift clumps and local dSph galaxies. △ Less

Submitted 16 August, 2022; v1 submitted 27 June, 2022; originally announced June 2022.

Comments: MNRAS, accepted for publication. 24 pages, 17 figures

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

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

Enhanced star formation in $z\sim6$ quasar companions

Authors: Tommaso Zana, Simona Gallerani, Stefano Carniani, Fabio Vito, Andrea Ferrara, Alessandro Lupi, Fabio Di Mascia, Paramita Barai

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

Submitted 7 April, 2022; originally announced April 2022.

Comments: 20 pages, 13 figures; Accepted for publications in MNRAS

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)

Resolving massive black hole binaries evolution via adaptive particle-splitting

Authors: Alessia Franchini, Alessandro Lupi, Alberto Sesana

Abstract: The study of the interaction of a massive black hole binary with its gaseous environment is crucial in order to be able to predict merger rates and possible electromagnetic counterparts of gravitational wave signals. The evolution of the binary semi-major axis resulting from this interaction has been recently debated, and a clear consensus is still missing, also because of several numerical limita… ▽ More The study of the interaction of a massive black hole binary with its gaseous environment is crucial in order to be able to predict merger rates and possible electromagnetic counterparts of gravitational wave signals. The evolution of the binary semi-major axis resulting from this interaction has been recently debated, and a clear consensus is still missing, also because of several numerical limitations, i.e. fixed orbit binaries or lack of resolution inside the cavity carved by the binary in its circumbinary disc. Using on-the-fly particle-splitting in the 3D meshless code gizmo, we achieve hyper-Lagrangian resolution, which allows us to properly resolve the dynamics inside the cavity, and in particular for the first time the discs that form around the two components of a live binary surrounded by a locally isothermal gaseous circumbinary disc. We show that the binary orbit decays with time for very cold and very warm discs and that the result of the interaction in the intermediate regime is strongly in uenced by the disc viscosity as this essentially regulates the fraction of mass contained in the discs around the binary components as well as the fraction that is accreted by the binary. We find the balance between these two quantities to determine whether the binary semi-major axis decreases with time. △ Less

Submitted 3 May, 2022; v1 submitted 14 January, 2022; originally announced January 2022.

Comments: 8 pages, published in ApJL

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

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

Authors: Alessandro Lupi, Stefano Bovino, Tommaso Grassi

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

Submitted 4 January, 2022; v1 submitted 6 September, 2021; originally announced September 2021.

Comments: 12 pages, 9 figures, 2 tables, matches the published version

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

Authors: Alessandro Lupi, Marta Volonteri, Roberto Decarli, Stefano Bovino, Joseph Silk

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

Submitted 4 January, 2022; v1 submitted 3 September, 2021; originally announced September 2021.

Comments: 21 pages, 16 figures, 1 table, matches the accepted version

The AGORA High-resolution Galaxy Simulations Comparison Project. III: Cosmological zoom-in simulation of a Milky Way-mass halo

Authors: Santi Roca-Fàbrega, Ji-hoon Kim, Loic Hausammann, Kentaro Nagamine, Johnny W. Powell, Ikkoh Shimizu, Daniel Ceverino, Alessandro Lupi, Joel R. Primack, Thomas Quinn, Yves Revaz, Héctor Velázquez, Tom Abel, Michael Buehlmann, Avishai Dekel, Bili Dong, Oliver Hahn, Cameron B. Hummels, Ki-won Kim, Britton D. Smith, Clayton J. Strawn, Romain Teyssier, Matthew Turk

Abstract: We present a suite of high-resolution cosmological zoom-in simulations to $z=4$ of a $10^{12}\,{\rm M}_{\odot}$ halo at $z=0$, obtained using seven contemporary astrophysical simulation codes widely used in the numerical galaxy formation community. Physics prescriptions for gas cooling, heating, and star formation, are similar to the ones used in our previous {\it AGORA} disk comparison but now ac… ▽ More We present a suite of high-resolution cosmological zoom-in simulations to $z=4$ of a $10^{12}\,{\rm M}_{\odot}$ halo at $z=0$, obtained using seven contemporary astrophysical simulation codes widely used in the numerical galaxy formation community. Physics prescriptions for gas cooling, heating, and star formation, are similar to the ones used in our previous {\it AGORA} disk comparison but now account for the effects of cosmological processes. In this work, we introduce the most careful comparison yet of galaxy formation simulations run by different code groups, together with a series of four calibration steps each of which is designed to reduce the number of tunable simulation parameters adopted in the final run. After all the participating code groups successfully completed the calibration steps, we reach a suite of cosmological simulations with similar mass assembly histories down to $z=4$. With numerical accuracy that resolves the internal structure of a target halo, we find that the codes overall agree well with one another in e.g., gas and stellar properties, but also show differences in e.g., circumgalactic medium properties. We argue that, if adequately tested in accordance with our proposed calibration steps and common parameters, the results of high-resolution cosmological zoom-in simulations can be robust and reproducible. New code groups are invited to join this comparison by generating equivalent models by adopting the common initial conditions, the common easy-to-implement physics package, and the proposed calibration steps. Further analyses of the simulations presented here will be in forthcoming reports from our Collaboration. △ Less

Submitted 17 June, 2021; originally announced June 2021.

Comments: Accepted for publication to the ApJ. The simulation snapshots used in this paper will be made publicly available in a user-friendly interface. The announcement of the data release will be posted on the AGORA Project webpage (www.agorasimulations.org). Early access to the data can be provided upon request to the authors, or to the project coordinator ([email protected], [email protected])

Chemical analysis of prestellar cores in Ophiuchus yields short timescales and rapid collapse

Authors: Stefano Bovino, Alessandro Lupi, Andrea Giannetti, Giovanni Sabatini, Dominik R. G. Schleicher, Friedrich Wyrowski, Karl M. Menten

Abstract: Sun-like stars form from the contraction of cold and dense interstellar clouds. How the collapse proceeds and what are the main physical processes driving it, however, is still under debate and a final consensus on the timescale of the process has not been reached. Does this contraction proceed slowly, sustained by strong magnetic fields and ambipolar diffusion, or is it driven by fast collapse wi… ▽ More Sun-like stars form from the contraction of cold and dense interstellar clouds. How the collapse proceeds and what are the main physical processes driving it, however, is still under debate and a final consensus on the timescale of the process has not been reached. Does this contraction proceed slowly, sustained by strong magnetic fields and ambipolar diffusion, or is it driven by fast collapse with gravity dominating the entire process? One way to answer this question is to measure the age of prestellar cores through statistical methods based on observations or via reliable chemical chronometers, which should better reflect the physical conditions of the cores. Here we report APEX observations of ortho-H$_2$D$^+$ and para-D$_2$H$^+$ for six cores in the Ophiuchus complex and combine them with detailed three-dimensional magneto-hydrodynamical simulations including chemistry, providing a range of ages for the observed cores up to 200 kyr. The outcome of our simulations and subsequent analysis provides a good match with the observational results in terms of physical (core masses and volume densities) and dynamical parameters such as the Mach number and the virial parameter. We show that models of fast collapse successfully reproduce the observed range of chemical abundance ratios as the timescales to reach the observed stages is comparable to the dynamical time of the cores (i.e. the free-fall time) and much shorter than the ambipolar diffusion time, measured from the electron fraction in the simulations. To confirm that this ratio can be used to distinguish between different star-formation scenarios a larger (statistically relevant) sample of star-forming cores should be explored. △ Less

Submitted 17 June, 2021; v1 submitted 5 May, 2021; originally announced May 2021.

Comments: submitted to A&A, 12 pages, 5 figures, 4 tables

Chemical post-processing of magneto-hydrodynamical simulations of star-forming regions: robustness and pitfalls

Authors: Simón Ferrada-Chamorro, Alessandro Lupi, Stefano Bovino

Abstract: A common approach to model complex chemistry in numerical simulations is via post-processing of existing magneto-hydrodynamic simulations, relying on computing the evolution of chemistry over the dynamic history of a subset of particles from within the raw simulation. Here, we validate such a technique, assessing its ability to recover the abundances of chemical species, using the chemistry packag… ▽ More A common approach to model complex chemistry in numerical simulations is via post-processing of existing magneto-hydrodynamic simulations, relying on computing the evolution of chemistry over the dynamic history of a subset of particles from within the raw simulation. Here, we validate such a technique, assessing its ability to recover the abundances of chemical species, using the chemistry package KROME. We also assess, for the first time, the importance of the main free input parameters, by means of a direct comparison with a self-consistent state-of-the-art simulation in which chemistry was directly coupled to hydrodynamics. We have found that the post-processing is highly reliable, with an accuracy at the percent level, even when the most relaxed input parameters are employed. In particular, our results show that the number of particles used does not affect significantly the average properties, although it suppresses the appearance of possibly important spatial features. On the other hand, the choice of the integration time-step plays a crucial role. Longer integration time-steps can produce large errors, as the post-processing solution will be forced towards chemical equilibrium, a condition that does not always necessarily apply. When the interpolation-based reconstruction of chemical properties is performed, the errors further increase up to a factor of $\sim2$. Concluding, our results suggest that this technique is extremely useful when exploring the relative quantitative effect of different chemical parameters and/or networks, without the need of re-running simulations multiple times, but some care should be taken in the choice of particles sub-sample and integration time-step. △ Less

Submitted 4 May, 2021; originally announced May 2021.

Comments: 11 pages, 6 figures, 3 tables

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

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

Authors: Alessandro Lupi, Zoltan Haiman, Marta Volonteri

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

Submitted 4 January, 2022; v1 submitted 9 February, 2021; originally announced February 2021.

Comments: 16 pages, 9 figures, 3 tables, matches the published version

A Virgo Environmental Survey Tracing Ionised Gas Emission (VESTIGE).IX. The effects of ram pressure stripping down to the scale of individual HII regions in the dwarf galaxy IC 3476

Authors: A. Boselli, A. Lupi, B. Epinat, P. Amram, M. Fossati, J. P. Anderson, S. Boissier, M. Boquien, G. Consolandi, P. Cote, J. C. Cuillandre, L. Ferrarese, L. Galbany, G. Gavazzi, J. A. Gomez-Lopez, S. Gwyn, G. Hensler, J. Hutchings, H. Kuncarayakti, A. Longobardi, E. W. Peng, H. Plana, J. Postma, J. Roediger, Y. Roehlly , et al. (3 additional authors not shown)

Abstract: We study the IB(s)m galaxy IC 3476 observed in the context of VESTIGE, a blind narrow-band Halpha+[NII] imaging survey of the Virgo cluster. The deep narrow-band (NB) image reveals a very pertubed ionised gas distribution, characterised by a prominent banana-shaped structure in the front of the galaxy formed of giant HII regions crossing the stellar disc, with star forming structures at ~8 kpc fro… ▽ More We study the IB(s)m galaxy IC 3476 observed in the context of VESTIGE, a blind narrow-band Halpha+[NII] imaging survey of the Virgo cluster. The deep narrow-band (NB) image reveals a very pertubed ionised gas distribution, characterised by a prominent banana-shaped structure in the front of the galaxy formed of giant HII regions crossing the stellar disc, with star forming structures at ~8 kpc from the edges of the stellar disc, detected also in a deep FUV ASTROSAT/UVIT image. This particular morphology indicates that the galaxy is undergoing an almost edge-on ram pressure stripping event. The NB image also shows that the star formation activity is totally quenched in the leading edge of the disc, where the gas has been removed during the interaction. The SED fitting analysis indicates that this quenching episode is very recent (~50 Myr), and roughly corresponds to an increase of the star formation activity in the inner regions with respect to what expected for secular evolution. The analysis of these data, whose angular resolution allows the study of the induced effects of the perturbation down to the scale of individual HII regions, also suggests that the increase of the star formation activity is due to the compression of the gas along the stellar disc of the galaxy, which is able to increase its mean electron density and boost the star formation process producing bright HII regions. The hydrodynamic interaction has deeply perturbed the velocity field of the ionised gas component while leaving unaffected that of the stellar disc. The comparison of the data with hydrodynamic simulations accounting for the different gas phases (atomic, molecular, ionised) consistently indicates that the perturbing event is very recent, once again confirming that ram pressure stripping is a violent phenomenon able to perturb on short timescales the evolution of galaxies in rich environments. △ Less

Submitted 14 December, 2020; originally announced December 2020.

Comments: Accepted for publication on Astronomy & Astrophysics

Journal ref: A&A 646, A139 (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

Constraints on the [CII] luminosity of a proto-globular cluster at z~6 obtained with ALMA

Authors: F. Calura, E. Vanzella, S. Carniani, R. Gilli, P. Rosati, M. Meneghetti, R. Paladino, R. Decarli, M. Brusa, A. Lupi, Q. D'Amato, P. Bergamini, . G. B. Caminha, -

Abstract: We report on ALMA observations of D1, a system at z~6.15 with stellar mass M_* ~ 10^7 M_sun containing globular cluster (GC) precursors, strongly magnified by the galaxy cluster MACS J0416.1-2403. Since the discovery of GC progenitors at high redshift, ours is the first attempt to probe directly the physical properties of their neutral gas through infrared observations. A careful analysis of our d… ▽ More We report on ALMA observations of D1, a system at z~6.15 with stellar mass M_* ~ 10^7 M_sun containing globular cluster (GC) precursors, strongly magnified by the galaxy cluster MACS J0416.1-2403. Since the discovery of GC progenitors at high redshift, ours is the first attempt to probe directly the physical properties of their neutral gas through infrared observations. A careful analysis of our dataset, performed with a suitable procedure designed to identify faint narrow lines and which can test various possible values for the unknown linewidth value, allowed us to identify a 4-sigma tentative detection of [CII] emission with intrinsic luminosity L_CII=(2.9 +/- 1.4) 10^6 L_sun, one of the lowest values ever detected at high redshift. This study offers a first insight on previously uncharted regions of the L_CII-SFR relation. Despite large uncertainties affecting our measure of the star formation rate, if taken at face value our estimate lies more than 1 dex below the values observed in local and high redshift systems. Our weak detection indicates a deficiency of [CII] emission, possibly ascribed to various explanations, such as a low-density gas and/or a strong radiation field caused by intense stellar feedback, and a low metal content. From the non-detection in the continuum we derive constraints on the dust mass, with 3-sigma upper limit values as low as a few 10^4 M_sun, consistent with the values measured in local metal-poor galaxies. △ Less

Submitted 14 October, 2020; originally announced October 2020.

Comments: MNRAS, accepted. 15 pages, 5 figures

Shaping the structure of a GMC with radiation and winds

Authors: D. Decataldo, A. Lupi, A. Ferrara, A. Pallottini, M. Fumagalli

Abstract: We study the effect of stellar feedback (photodissociation/ionization, radiation pressure and winds) on the evolution of a Giant Molecular Cloud (GMC), by means of a 3D radiative transfer, hydro-simulation implementing a complex chemical network featuring ${\rm H}_2$ formation and destruction. We track the formation of individual stars with mass $M>1\,{\rm M}_\odot$ with a stochastic recipe. Each… ▽ More We study the effect of stellar feedback (photodissociation/ionization, radiation pressure and winds) on the evolution of a Giant Molecular Cloud (GMC), by means of a 3D radiative transfer, hydro-simulation implementing a complex chemical network featuring ${\rm H}_2$ formation and destruction. We track the formation of individual stars with mass $M>1\,{\rm M}_\odot$ with a stochastic recipe. Each star emits radiation according to its spectrum, sampled with 10 photon bins from near-infrared to extreme ultra-violet bands; winds are implemented by energy injection in the neighbouring cells. We run a simulation of a GMC with mass $M=10^5\,{\rm M}_\odot$, following the evolution of different gas phases. Thanks to the simultaneous inclusion of different stellar feedback mechanisms, we identify two stages in the cloud evolution: (1) radiation and winds carve ionized, low-density bubbles around massive stars, while FUV radiation dissociates most ${\rm H}_2$ in the cloud, apart from dense, self-shielded clumps; (2) rapid star formation (SFR$\simeq 0.1\,{\rm M}_\odot\,{\rm yr}^{-1}$) consumes molecular gas in the dense clumps, so that UV radiation escapes and ionizes the remaining HI gas in the GMC. ${\rm H}_2$ is exhausted in $1.6$ Myr, yielding a final star formation efficiency of 36 per cent. The average intensity of FUV and ionizing fields increases almost steadily with time; by the end of the simulation ($t=2.5$ Myr) we find $\langle G_0 \rangle \simeq 10^3$ (in Habing units), and a ionization parameter $\langle U_{\rm ion} \rangle \simeq 10^2$, respectively. The ionization field has also a more patchy distribution than the FUV one within the GMC. Throughout the evolution, the escape fraction of ionizing photons from the cloud is $f_{\rm ion, esc} < 0.03$. △ Less

Submitted 16 September, 2020; originally announced September 2020.

Comments: 16 pages, 12 figures

Journal ref: Monthly Notices of the Royal Astronomical Society, Volume 497, Issue 4, October 2020, Pages 4718-4732

Missing [CII] emission from early galaxies

Authors: S. Carniani, A. Ferrara, R. Maiolino, M. Castellano, S. Gallerani, A. Fontana, M. Kohandel, A. Lupi, A. Pallottini, L. Pentericci, L. Vallini, E. Vanzella

Abstract: ALMA observations have revealed that [CII] 158$μ$m line emission in high-z galaxies is ~2-3$\times$ more extended than the UV continuum emission. Here we explore whether surface brightness dimming (SBD) of the [CII] line is responsible for the reported [CII] deficit, and the large $L_{\rm [OIII]}/L_{\rm [CII]}$ luminosity ratio measured in early galaxies. We first analyse archival ALMA images of n… ▽ More ALMA observations have revealed that [CII] 158$μ$m line emission in high-z galaxies is ~2-3$\times$ more extended than the UV continuum emission. Here we explore whether surface brightness dimming (SBD) of the [CII] line is responsible for the reported [CII] deficit, and the large $L_{\rm [OIII]}/L_{\rm [CII]}$ luminosity ratio measured in early galaxies. We first analyse archival ALMA images of nine z>6 galaxies observed in both [CII] and [OIII]. After performing several uv-tapering experiments to optimize the identification of extended line emission, we detect [CII] emission in the whole sample, with an extent systematically larger than the [CII] emission. Next, we use interferometric simulations to study the effect of SBD on the line luminosity estimate. About 40% of the extended [CII] component might be missed at an angular resolution of 0.8$^{\prime\prime}$, implying that $L_{\rm [CII]}$ is underestimated by a factor $\approx2$ in data at low (<7) signal-to-noise ratio . By combining these results, we conclude that $L_{\rm [CII]}$ of z>6 galaxies lies, on average, slightly below the local $L_{\rm [CII]}-SFR$ relation ($Δ^{z=6-9}=-0.07\pm0.3$), but within the intrinsic dispersion of the relation. SBD correction also yields $L_{\rm [OIII]}/L_{\rm [CII]}<10$, i.e. more in line with current hydrodynamical simulations. △ Less

Submitted 7 December, 2020; v1 submitted 16 June, 2020; originally announced June 2020.