-
Is GN-z11 powered by a super-Eddington massive black hole?
Authors:
Maulik Bhatt,
Simona Gallerani,
Andrea Ferrara,
Chiara Mazzucchelli,
Valentina D'Odorico,
Milena Valentini,
Tommaso Zana,
Emanuele Paolo Farina,
Srija Chakraborty
Abstract:
Observations of $z \sim 6$ quasars powered by supermassive black holes (SMBHs; $M_{\rm BH} \sim 10^{8-10}\, M_\odot$) challenge our current understanding of early black hole (BH) formation and evolution. The advent of the James Webb Space Telescope (JWST) has enabled the study of massive BHs (MBHs; $M_{\rm BH}\sim 10^{6-7} \ \mathrm{M}_\odot$) up to $z\sim 11$, thus bridging the properties of…
▽ More
Observations of $z \sim 6$ quasars powered by supermassive black holes (SMBHs; $M_{\rm BH} \sim 10^{8-10}\, M_\odot$) challenge our current understanding of early black hole (BH) formation and evolution. The advent of the James Webb Space Telescope (JWST) has enabled the study of massive BHs (MBHs; $M_{\rm BH}\sim 10^{6-7} \ \mathrm{M}_\odot$) up to $z\sim 11$, thus bridging the properties of $z\sim 6$ quasars to their ancestors. The JWST spectroscopic observations of GN-z11, a well-known $z=10.6$ star-forming galaxy, have been interpreted with the presence of a super-Eddington (Eddington ratio $\equiv \,λ_{\rm Edd}\sim 5.5$) accreting MBH. To test this hypothesis, we used a zoom-in cosmological simulation of galaxy formation and BH co-evolution. We first tested the simulation results against the observed probability distribution function (PDF) of $λ_{\rm Edd}$ found in $z\sim 6$ quasars. Then, in the simulation we selected the BHs that satisfy the following criteria: (a) $10 < z < 11 $, (b) $M_{\rm BH} > 10^6 \ \mathrm{M}_\odot$. Next, we apply the extreme value statistics to the PDF of $λ_{\rm Edd}$ resulting from the simulation and we find that the probability of observing a $z\sim 10-11$ MBH accreting with $λ_{\rm Edd} \sim 5.5$ in the volume surveyed by JWST is very low ($<0.2\%$). We compared our predictions with those in the literature, and discuss the main limitations of our work. Our simulation cannot explain the JWST observations of GN-z11. This might be due to (i) poor resolution and statistics in simulations, (ii) simplistic sub-grid models (e.g. BH accretion and seeding), (iii) uncertainties in the data analysis and interpretation.
△ Less
Submitted 8 April, 2024; v1 submitted 24 January, 2024;
originally announced January 2024.
-
Are there more galaxies than we see around high-$z$ quasars?
Authors:
Tommaso Zana,
Stefano Carniani,
David Prelogović,
Fabio Vito,
Viola Allevato,
Andrea Ferrara,
Simona Gallerani,
Eleonora Parlanti
Abstract:
Whether or not $z \gtrsim 6$ quasars lie in the most massive dark-matter halos of the Universe is still a subject of dispute. While most theoretical studies support this scenario, current observations yield discordant results when they probe the halo mass through the detection rate of quasar companion galaxies. Feedback processes from supermassive black holes and dust obscuration have been blamed…
▽ More
Whether or not $z \gtrsim 6$ quasars lie in the most massive dark-matter halos of the Universe is still a subject of dispute. While most theoretical studies support this scenario, current observations yield discordant results when they probe the halo mass through the detection rate of quasar companion galaxies. Feedback processes from supermassive black holes and dust obscuration have been blamed for this discrepancy, but the impact of these effects is complex and far from being clearly understood. This paper aims to improve the interpretation of current far-infrared observations by taking into account the cosmological volume probed by the Atacama Large Millimeter/submillimeter Array Telescope and to explain the observational discrepancies. We statistically investigate the detection rate of quasar companions in current observations and verify if they match the expected distribution from various theoretical models, once convolved with the ALMA field-of-view, through the use of Monte Carlo simulations. We demonstrate that the telescope geometrical bias is fundamental and can alone explain the scatter in the number of detected satellite galaxies in different observations. We conclude that the resulting companion densities depend on the chosen galaxy distributions. According to our fiducial models, current data favour a density scenario where quasars lie in dark-matter halos of viral mass $M_{\rm vir} \gtrsim 10^{12}~{\rm M_{\odot}}$, in agreement with most theoretical studies. According to our analysis, each quasar has about 2 companion galaxies, with a [CII] luminosity $L_{\rm [CII]} \gtrsim 10^8~{\rm L}_{\odot}$, within a distance of about 1~Mpc from the quasar.
△ Less
Submitted 7 September, 2023;
originally announced September 2023.
-
Probing $z \gtrsim 6$ massive black holes with gravitational waves
Authors:
Srija Chakraborty,
Simona Gallerani,
Tommaso Zana,
Alberto Sesana,
Milena Valentini,
David Izquierdo-Villalba,
Fabio Di Mascia,
Fabio Vito,
Paramita Barai
Abstract:
We investigate the coalescence of massive black hole ($M_{\rm BH}\gtrsim 10^{6}~\rm M_{\odot}$) binaries (MBHBs) at $6<z<10$ by adopting a suite of cosmological hydrodynamical simulations of galaxy formation, zoomed-in on biased ($ >3 σ$) overdense regions ($M_h\sim 10^{12}~\rm M_{\odot}$ dark matter halos at $z = 6$) of the Universe. We first analyse the impact of different resolutions and AGN fe…
▽ More
We investigate the coalescence of massive black hole ($M_{\rm BH}\gtrsim 10^{6}~\rm M_{\odot}$) binaries (MBHBs) at $6<z<10$ by adopting a suite of cosmological hydrodynamical simulations of galaxy formation, zoomed-in on biased ($ >3 σ$) overdense regions ($M_h\sim 10^{12}~\rm M_{\odot}$ dark matter halos at $z = 6$) of the Universe. We first analyse the impact of different resolutions and AGN feedback prescriptions on the merger rate, assuming instantaneous mergers. Then, we compute the halo bias correction factor due to the overdense simulated region. Our simulations predict merger rates that range between 3 - 15 $\rm yr^{-1}$ at $z\sim 6$, depending on the run considered, and after correcting for a bias factor of $\sim 20-30$.
For our fiducial model, we further consider the effect of delay in the MBHB coalescence due to dynamical friction. We find that 83 per cent of MBHBs will merge within the Hubble time, and 21 per cent within 1 Gyr, namely the age of the Universe at $z > 6$. We finally compute the expected properties of the gravitational wave (GW) signals and find the fraction of LISA detectable events with high signal-to-noise ratio (SNR $>$ 5) to range between 66-69 per cent. However, identifying the electro-magnetic counterpart of these events remains challenging due to the poor LISA sky localization that, for the loudest signals ($\mathcal M_c\sim 10^6~\rm M_{\odot}$ at $z=6$), is around 10 $\rm deg^2$.
△ Less
Submitted 16 May, 2023;
originally announced May 2023.
-
Bubbles and outflows: the novel JWST/NIRSpec view of the z=1.59 obscured quasar XID2028
Authors:
Giovanni Cresci,
Giulia Tozzi,
Michele Perna,
Marcella Brusa,
Cosimo Marconcini,
Alessandro Marconi,
Stefano Carniani,
Marisa Brienza,
Marcello Giroletti,
Francesco Belfiore,
Michele Ginolfi,
Filippo Mannucci,
Lorenzo Ulivi,
Jan Scholtz,
Giacomo Venturi,
Santiago Arribas,
Hanna Übler,
Francesco D'Eugenio,
Matilde Mingozzi,
Barbara Balmaverde,
Alessandro Capetti,
Eleonora Parlanti,
Tommaso Zana
Abstract:
Quasar feedback in the form of powerful outflows is invoked as a key mechanism to quench star formation in galaxies, although direct observational evidence is still scarce and debated. Here we present Early Release Science JWST NIRSpec IFU observations of the z=1.59 prototypical obscured Active Galactic Nucleus (AGN) XID2028: This target represents a unique test case for studying quasar feedback a…
▽ More
Quasar feedback in the form of powerful outflows is invoked as a key mechanism to quench star formation in galaxies, although direct observational evidence is still scarce and debated. Here we present Early Release Science JWST NIRSpec IFU observations of the z=1.59 prototypical obscured Active Galactic Nucleus (AGN) XID2028: This target represents a unique test case for studying quasar feedback at the peak epoch of AGN-galaxy co-evolution because extensive multi-wavelength coverage is available and a massive and extended outflow is detected in the ionised and molecular components. With the unprecedented sensitivity and spatial resolution of the JWST, the NIRSpec dataset reveals a wealth of structures in the ionised gas kinematics and morphology that were previously hidden in the seeing-limited ground-based data. In particular, we find evidence of an interaction between the interstellar medium of the galaxy and the quasar-driven outflow and radio jet that produces an expanding bubble from which the fast and extended wind detected in previous observations emerges. The new observations confirm the complex interplay between the AGN jet, wind and the interstellar medium of the host galaxy, highlighting the role of low-luminosity radio jets in AGN feedback. They also clearly show the new window that NIRSpec opens for detailed studies of feedback at high redshift.
△ Less
Submitted 20 March, 2023; v1 submitted 26 January, 2023;
originally announced January 2023.
-
The feedback of quasars on their galactic environment
Authors:
Andrea Ferrara,
Tommaso Zana,
Simona Gallerani,
Laura Sommovigo
Abstract:
Quasar outflows might either quench (negative) or enhance (positive feedback) star formation in galaxies located in the quasar environment. The possible outcome depends on 4 parameters: the quasar ($σ$) and satellite ($σ_*$) halo velocity dispersion, their relative distance, $d$, and satellite disk radius, $r_d$. We find that: (i) small satellites with $σ_* < 164\ σ_{200}^{2/3}\, \rm km\ s^{-1}$ h…
▽ More
Quasar outflows might either quench (negative) or enhance (positive feedback) star formation in galaxies located in the quasar environment. The possible outcome depends on 4 parameters: the quasar ($σ$) and satellite ($σ_*$) halo velocity dispersion, their relative distance, $d$, and satellite disk radius, $r_d$. We find that: (i) small satellites with $σ_* < 164\ σ_{200}^{2/3}\, \rm km\ s^{-1}$ have their star formation quenched; (ii) in larger satellites, star formation, and hence UV/FIR luminosity, is instead boosted by $>80$\% in a burst with a typical duration of $5-10$ Myr, if the following positive feedback criterion is met: ${d}/{r_d} < 15 (Q/η)^{1/2} σ_{200}$, where $Q \approx 1$ is the satellite disk Toomre parameter; the disruption parameter (see eq. 17) must be $η>1$ to prevent complete satellite gas removal. We compare our predictions with ALMA data finding that observed satellites of $z\simeq 6$ QSOs on average form stars at a $3\times$ higher rate with respect to field galaxies at the same redshift. Further tests of the model are suggested.
△ Less
Submitted 25 January, 2023;
originally announced January 2023.
-
The black hole and host galaxy growth in an isolated $z\sim 6$ QSO observed with ALMA
Authors:
R. Tripodi,
C. Feruglio,
F. Fiore,
M. Bischetti,
V. D'Odorico,
S. Carniani,
S. Cristiani,
S. Gallerani,
R. Maiolino,
A. Marconi,
A. Pallottini,
E. Piconcelli,
L. Vallini,
T. Zana
Abstract:
The outstanding mass growth of supermassive black holes (SMBHs) at the Reionisation Epoch and how it is related to the concurrent growth of their host galaxies, poses challenges to theoretical models aimed at explaining how these systems formed in short timescales (<1 Gyr). To trace the average evolutionary paths of quasi-stellar objects (QSOs) and their host galaxies in the BH mass-host mass (…
▽ More
The outstanding mass growth of supermassive black holes (SMBHs) at the Reionisation Epoch and how it is related to the concurrent growth of their host galaxies, poses challenges to theoretical models aimed at explaining how these systems formed in short timescales (<1 Gyr). To trace the average evolutionary paths of quasi-stellar objects (QSOs) and their host galaxies in the BH mass-host mass ($M_{\rm dyn}$) plane, we compare the star formation rate (SFR), derived from the accurate estimate of the dust temperature and the dust mass ($T_{\rm dust}, M_{\rm dust}$), with the BH accretion rate. To this aim, we analysed a deep, $900$ pc resolution ALMA observation of the sub-mm continuum, [CII] and H$_2$O of the $z\sim 6$ QSO J2310+1855, enabling a detailed study of dust properties and cold gas kinematics. We performed an accurate SED analysis obtaining a dust temperature of $T_{\rm dust} = 71$ K and a dust mass of $M_{\rm dust}= 4.4 \times 10^8\ \rm M_{\odot}$. The implied AGN-corrected SFR is $1240 \ \rm M_{\odot}yr^{-1}$, a factor of 2 smaller than previously reported for this QSO. We derived the best estimate of the dynamical mass $M_{\rm dyn} = 5.2\times 10^{10}\ \rm M_{\odot}$ within $r = 1.7$ kpc, based on a dynamical model of the system. We found that ${\rm SFR}/M_{\rm dyn}>\dot M_{\rm BH}/M_{\rm BH}$, suggesting that AGN feedback might be efficiently acting to slow down the SMBH accretion, while the stellar mass assembly is still vigorously taking place in the host galaxy. In addition, we were also able to detect high-velocity emission on the red and blue sides of the [CII] emission line, that traces a gaseous outflow, and for the first time, we mapped a spatially-resolved water vapour disk through the H$_2$O v=0 $3_{(2,2)}-3_{(1,3)}$ emission line detected at $ν_{\rm obs} = 274.074$ GHz, whose kinematic properties and size are broadly consistent with those of the [CII] disk.
△ Less
Submitted 30 July, 2022; v1 submitted 7 July, 2022;
originally announced July 2022.
-
Morphological decomposition of TNG50 galaxies: methodology and catalogue
Authors:
Tommaso Zana,
Alessandro Lupi,
Matteo Bonetti,
Massimo Dotti,
Yetli Rosas-Guevara,
David Izquierdo-Villalba,
Silvia Bonoli,
Lars Hernquist,
Dylan Nelson
Abstract:
We present MORDOR (MORphological DecOmposeR, a new algorithm for structural decomposition of simulated galaxies based on stellar kinematics. The code measures the properties of up to five structural components (a thin/cold and a thick/warm disc, a classical and a secular bulge, and a spherical stellar halo), and determines the properties of a stellar bar (if present). A comparison with other algor…
▽ More
We present MORDOR (MORphological DecOmposeR, a new algorithm for structural decomposition of simulated galaxies based on stellar kinematics. The code measures the properties of up to five structural components (a thin/cold and a thick/warm disc, a classical and a secular bulge, and a spherical stellar halo), and determines the properties of a stellar bar (if present). A comparison with other algorithms presented in the literature yields overall good agreement, with MORDOR displaying a higher flexibility in correctly decomposing systems and identifying bars in crowded environments (e.g. with ongoing fly-bys, often observable in cosmological simulations). We use MORDOR to analyse galaxies in the TNG50 simulation and find the following: ($i$) the thick disc component undergoes the strongest evolution in the binding energy-circularity plane, as expected when disc galaxies decrease their turbulent-rotational support with cosmic time; ($ii$) smaller galaxies (with stellar mass, $10^{9} \lesssim M_{*} / {\rm M_{\odot}} \leq 5 \times 10^{9}$) undergo a major growth in their disc components after $z\sim 1$, whereas ($iii$) the most massive galaxies ($5 \times 10^{10} < M_{*} / {\rm M_{\odot}} \leq 5\times10^{11}$) evolve toward more spheroidal dominated objects down to $z=0$ due to frequent gravitational interactions with satellites; ($iv$) the fraction of barred galaxies grows rapidly at high redshift and stabilizes below $z\sim 2$, except for the most massive galaxies that show a decrease in the bar occupation fraction at low redshift; ($v$) galaxies with $M_{*} \sim 10^{11}~{\rm M_{\odot}}$ exhibit the highest relative occurrence of bars at $z=0$, in agreement with observational studies. We publicly release MORDOR and the morphological catalogue of TNG50 galaxies.
△ Less
Submitted 8 July, 2022; v1 submitted 9 June, 2022;
originally announced June 2022.
-
Feedback effect on the observable properties of $z>6$ AGN
Authors:
Fabio Vito,
Fabio Di Mascia,
Simona Gallerani,
Tommaso Zana,
Andrea Ferrara,
Stefano Carniani,
Roberto Gilli
Abstract:
Active galactic nuclei (AGN) feedback has a major impact onto the supermassive black-hole (SMBH) growth, the properties of the host galaxies, and their cosmic evolution. We investigate the effects of different kinetic feedback prescriptions on the observable properties of AGN and their host galaxies at $z>6$ in a suite of zoom-in cosmological simulations. We find that kinetic feedback decreases th…
▽ More
Active galactic nuclei (AGN) feedback has a major impact onto the supermassive black-hole (SMBH) growth, the properties of the host galaxies, and their cosmic evolution. We investigate the effects of different kinetic feedback prescriptions on the observable properties of AGN and their host galaxies at $z>6$ in a suite of zoom-in cosmological simulations. We find that kinetic feedback decreases the column density of the interstellar medium (ISM) in the host galaxy by up to a factor of $\approx10$, especially when the SMBHs reach high accretion rates ($\approx10-30\,\mathrm{M_\odot\,yr^{-1}}$). In particular, kinetic feedback is required to extend the ISM size to $>1$ kpc and match the observed sizes of the gas reservoirs in $z>6$ AGN host galaxies. Moreover, it produces unobscured lines of sight along which the AGN can be detected in the rest-frame UV band with magnitudes consistent with observed values of $z>6$ AGN. The assumed geometry of the outflow plays an important role in shaping the observed properties of high-redshift AGN. We find that a biconical geometry is favored over a spherical one to reproduce the observed properties, but it overestimates the number of multiple AGN systems detectable in X-ray observations. This result suggests that simplistic BH seeding recipes widely employed in cosmological simulations produce too many X-ray detectable multiple AGN at $z=6-7$, thus soliciting the adoption of more physically motivated seeding prescriptions.
△ Less
Submitted 18 May, 2022;
originally announced May 2022.
-
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.
-
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.
-
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.
-
Infrared emission of $z \sim 6$ galaxies: AGN imprints
Authors:
F. Di Mascia,
S. Gallerani,
C. Behrens,
A. Pallottini,
S. Carniani,
A. Ferrara,
P. Barai,
F. Vito,
T. Zana
Abstract:
We investigate the infrared (IR) emission of high-redshift ($z\sim 6$), highly star-forming (${ {\rm SFR} > 100}$ $M_{\rm \odot} {\rm yr}^{-1}$) galaxies, with/without Active Galactic Nuclei (AGN), using a suite of cosmological simulations featuring dust radiative transfer. Synthetic Spectral Energy Distributions (SEDs) are used to quantify the relative contribution of stars/AGN to dust heating. I…
▽ More
We investigate the infrared (IR) emission of high-redshift ($z\sim 6$), highly star-forming (${ {\rm SFR} > 100}$ $M_{\rm \odot} {\rm yr}^{-1}$) galaxies, with/without Active Galactic Nuclei (AGN), using a suite of cosmological simulations featuring dust radiative transfer. Synthetic Spectral Energy Distributions (SEDs) are used to quantify the relative contribution of stars/AGN to dust heating. In dusty (${M_{\rm d}\gtrsim 3\times 10^7 M_{\rm \odot}}$) galaxies, $\gtrsim 50-90 \%$ of the UV radiation is obscured by dust inhomogeneities on scales ${\gtrsim 100}$ pc. In runs with AGN, a clumpy, warm ($\approx 250$ K) dust component co-exists with a colder ($\approx 60$ K) and more diffuse one, heated by stars. Warm dust provides up to ${50 \%}$ of the total IR luminosity, but only $\lesssim 0.1 \%$ of the total mass content. The AGN boosts the MIR flux by ${10-100 \times}$ with respect to star forming galaxies, without significantly affecting the FIR. Our simulations successfully reproduce the observed SED of bright (${M_{\rm UV}\sim -26}$) ${z\sim 6}$ quasars, and show that these objects are part of complex, dust-rich merging systems, containing multiple sources (accreting BHs and/or star forming galaxies) in agreement with recent HST and ALMA observations. Our results show that the proposed ORIGINS missions will be able to investigate the MIR properties of dusty star forming galaxies and to obtain good quality spectra of bright quasars at $z\sim 6$. Finally, the MIR-to-FIR flux ratio of faint (${M_{\rm UV}\sim -24}$) AGN is ${>10\times}$ higher than for normal star forming galaxies. This implies that combined JWST/ORIGINS/ALMA observations will be crucial to identify faint and/or dust-obscured AGN in the distant Universe.
△ Less
Submitted 17 February, 2021;
originally announced February 2021.
-
Global torques and stochasticity as the drivers of massive black hole pairing in the young Universe
Authors:
Elisa Bortolas,
Pedro R. Capelo,
Tommaso Zana,
Lucio Mayer,
Matteo Bonetti,
Massimo Dotti,
Melvyn B. Davies,
Piero Madau
Abstract:
The forthcoming Laser Interferometer Space Antenna (LISA) will probe the population of coalescing massive black hole (MBH) binaries up to the onset of structure formation. Here we simulate the galactic-scale pairing of $\sim10^6 M_\odot$ MBHs in a typical, non-clumpy main-sequence galaxy embedded in a cosmological environment at $z = 7-6$. In order to increase our statistical sample, we adopt a st…
▽ More
The forthcoming Laser Interferometer Space Antenna (LISA) will probe the population of coalescing massive black hole (MBH) binaries up to the onset of structure formation. Here we simulate the galactic-scale pairing of $\sim10^6 M_\odot$ MBHs in a typical, non-clumpy main-sequence galaxy embedded in a cosmological environment at $z = 7-6$. In order to increase our statistical sample, we adopt a strategy that allows us to follow the evolution of six secondary MBHs concomitantly. We find that the magnitude of the dynamical-friction induced torques is significantly smaller than that of the large-scale, stochastic gravitational torques arising from the perturbed and morphologically evolving galactic disc, suggesting that the standard dynamical friction treatment is inadequate for realistic galaxies at high redshift. The dynamical evolution of MBHs is very stochastic, and a variation in the initial orbital phase can lead to a drastically different time-scale for the inspiral. Most remarkably, the development of a galactic bar in the host system either significantly accelerates the inspiral by dragging a secondary MBH into the centre, or ultimately hinders the orbital decay by scattering the MBH in the galaxy outskirts. The latter occurs more rarely, suggesting that galactic bars overall promote MBH inspiral and binary coalescence. The orbital decay time can be an order of magnitude shorter than what would be predicted relying on dynamical friction alone. The stochasticity, and the important role of global torques, have crucial implications for the rates of MBH coalescences in the early Universe: both have to be accounted for when making predictions for the upcoming LISA observatory.
△ Less
Submitted 25 August, 2020; v1 submitted 5 May, 2020;
originally announced May 2020.
-
The buildup of strongly barred galaxies in the TNG100 simulation
Authors:
Yetli Rosas-Guevara,
Silvia Bonoli,
Massimo Dotti,
Tommaso Zana,
Dylan Nelson,
Annalisa Pillepich,
Luis C. Ho,
David Izquierdo-Villalba,
Lars Hernquist,
Ruediger Pakmor
Abstract:
We analyse the properties of strongly barred disc galaxies using the TNG100 simulation, a cosmological hydrodynamical realisation of the IllustrisTNG suite. We identify 270 disc galaxies at $z=0$ in the stellar mass range $M_{*}=10^{10.4-11}M_{\odot}$, of which 40 per cent are barred. Of the detected bars, more than half are strong. We find that the fraction of barred galaxies increases with stell…
▽ More
We analyse the properties of strongly barred disc galaxies using the TNG100 simulation, a cosmological hydrodynamical realisation of the IllustrisTNG suite. We identify 270 disc galaxies at $z=0$ in the stellar mass range $M_{*}=10^{10.4-11}M_{\odot}$, of which 40 per cent are barred. Of the detected bars, more than half are strong. We find that the fraction of barred galaxies increases with stellar mass, in agreement with observational results. Strongly barred galaxies exhibit, overall, lower gas-to-stellar mass ratios compared to unbarred galaxies. The majority of barred galaxies are quenched (sSFR $\sim10^{-11.7} $yr$^{-1}$), whereas unbarred galaxies continue to be active (sSFR $\sim10^{-10.3}$yr$^{-1}$) on the main sequence of star-forming galaxies. We explore the evolution of strongly barred and unbarred galaxies to investigate their formation and quenching histories. We find that strong bars form between $0.5< z< 1.5$, with more massive galaxies hosting older bars. Strong bars form in galaxies with an early-established prominent disc component, undergoing periods of enhanced star formation and black hole accretion, possibly assisted by cosmological inflows. Unbarred galaxies, on the other hand, assemble most of their mass and disc component at late times. The nuclear region of strongly barred galaxies quenches shortly after bar formation, while unbarred galaxies remain active across time. Our findings are indicative of bar quenching, possibly assisted by nuclear feedback processes.
We conclude that the cosmological environment, together with small scale feedback processes, determine the chances of a galaxy to form a bar and to rapidly quench its central region.
△ Less
Submitted 30 March, 2022; v1 submitted 1 August, 2019;
originally announced August 2019.
-
Barred galaxies in cosmological zoom-in simulations: the importance of feedback
Authors:
Tommaso Zana,
Pedro R. Capelo,
Massimo Dotti,
Lucio Mayer,
Alessandro Lupi,
Francesco Haardt,
Silvia Bonoli,
Sijing Shen
Abstract:
Bars are a key factor in the long-term evolution of spiral galaxies, in their unique role in redistributing angular momentum and transporting gas and stars on large scales. The Eris-suite simulations are cosmological zoom-in, N-body, smoothed-particle hydrodynamic simulations built to follow the formation and evolution of a Milky Way-sized galaxy across the build-up of the large scale structure. H…
▽ More
Bars are a key factor in the long-term evolution of spiral galaxies, in their unique role in redistributing angular momentum and transporting gas and stars on large scales. The Eris-suite simulations are cosmological zoom-in, N-body, smoothed-particle hydrodynamic simulations built to follow the formation and evolution of a Milky Way-sized galaxy across the build-up of the large scale structure. Here we analyse and describe the outcome of two particular simulations taken from the Eris suite - ErisBH and Eris2k - which mainly differ in the prescriptions employed for gas cooling, star formation, and feedback from supernovae and black holes. Our study shows that the enhanced effective feedback in Eris2k, due to the collective effect of the different micro-physics implementations, results in a galaxy which is less massive than its ErisBH counterpart till z~1. However, when the stellar content is large enough so that global dynamical instabilities can be triggered, the galaxy in Eris2k develops a stronger and more extended bar with respect to ErisBH. We demonstrate that he structural properties and time evolution of the two bars are very different. Our results highlight the importance of accurate sub-grid prescriptions in cosmological zoom-in simulations of the process of galaxy formation and evolution, and the possible use of a statistical sample of barred galaxies to assess the strength of the stellar feedback.
△ Less
Submitted 2 July, 2019; v1 submitted 17 October, 2018;
originally announced October 2018.
-
Bar resilience to flybys in a cosmological framework
Authors:
Tommaso Zana,
Massimo Dotti,
Pedro R. Capelo,
Lucio Mayer,
Francesco Haardt,
Sijing Shen,
Silvia Bonoli
Abstract:
It has been proposed that close interactions with satellite galaxies can significantly perturb the morphology of the main galaxy. However, the dynamics of an already formed bar following the interaction with the external environment has not been studied in detail in a fully cosmological context. In this work, analysing the cosmological zoom-in simulation Eris2k, we study the effects that a very un…
▽ More
It has been proposed that close interactions with satellite galaxies can significantly perturb the morphology of the main galaxy. However, the dynamics of an already formed bar following the interaction with the external environment has not been studied in detail in a fully cosmological context. In this work, analysing the cosmological zoom-in simulation Eris2k, we study the effects that a very unequal-mass flyby crossing the stellar disc has on the stability of the pre-existing bar. We characterize the evolution of the bar strength and length showing that the perturbation exerted by the flyby shuffles the orbits of stars for less than one Gyr. After this time, the bar shows a remarkable resilience, reforming with properties comparable to those it had before the interaction. Our work shows that close unequal-mass encounters, the most frequent interactions occurring during the evolution of cosmic structures, have (i) an overall minor impact on the global evolution of the bar in the long term, still (ii) the effect is destructive and (iii) a very weak interaction is sufficient to dismantle a strong bar leading to its "apparent death". As a consequence, due to the non-negligible duration of the bar-less period, a fraction of observed spiral galaxies classified as non-barred could be prone to bar formation.
△ Less
Submitted 20 July, 2018; v1 submitted 9 May, 2018;
originally announced May 2018.
-
External versus internal triggers of bar formation in cosmological zoom-in simulations
Authors:
Tommaso Zana,
Massimo Dotti,
Pedro R. Capelo,
Silvia Bonoli,
Francesco Haardt,
Lucio Mayer,
Daniele Spinoso
Abstract:
The emergence of a large-scale stellar bar is one of the most striking features in disc galaxies. By means of state-of-the-art cosmological zoom-in simulations, we study the formation and evolution of bars in Milky Way-like galaxies in a fully cosmological context, including the physics of gas dissipation, star formation, and supernova feedback. Our goal is to characterise the actual trigger of th…
▽ More
The emergence of a large-scale stellar bar is one of the most striking features in disc galaxies. By means of state-of-the-art cosmological zoom-in simulations, we study the formation and evolution of bars in Milky Way-like galaxies in a fully cosmological context, including the physics of gas dissipation, star formation, and supernova feedback. Our goal is to characterise the actual trigger of the non-axisymmetric perturbation that leads to the strong bar observable in the simulations at z=0, discriminating between an internal/secular versus an external/tidal origin. To this aim, we run a suite of cosmological zoom-in simulations altering the original history of galaxy-satellite interactions at a time when the main galaxy, though already bar-unstable, does not feature any non-axisymmetric structure yet. We find that the main effect of a late minor merger and of a close fly-by is to delay the time of bar formation and those two dynamical events are not directly responsible for the development of the bar and do not alter significantly its global properties (e.g. its final extension). We conclude that, once the disc has grown to a mass large enough to sustain global non-axisymmetric modes, then bar formation is inevitable.
△ Less
Submitted 20 July, 2018; v1 submitted 5 May, 2017;
originally announced May 2017.
