-
The complex effect of gas cooling and turbulence on AGN-driven outflow properties
Authors:
K. Zubovas,
M. Tartėnas,
M. A. Bourne
Abstract:
(abridged) Accretion onto supermassive black holes (SMBHs) at close to the Eddington rate can influence the host galaxy via powerful winds. Theoretical models of such winds can explain observational correlations between SMBHs and their host galaxies and the powerful multi-phase outflows observed in a number of active galaxies. Analytic models usually assume spherical symmetry and a smooth gas dist…
▽ More
(abridged) Accretion onto supermassive black holes (SMBHs) at close to the Eddington rate can influence the host galaxy via powerful winds. Theoretical models of such winds can explain observational correlations between SMBHs and their host galaxies and the powerful multi-phase outflows observed in a number of active galaxies. Analytic models usually assume spherical symmetry and a smooth gas distribution with an adiabatic equation of state. However, the interstellar medium in real galaxies is clumpy and cooling is important, complicating the analysis. We used a suite of idealised hydrodynamical simulations to isolate the effects of turbulence and cooling on the development and global properties of AGN wind-driven outflows on kiloparsec scales. We measured the outflow velocity, mass outflow rate and momentum and energy loading factors as the system evolved over 1.2 Myr and estimated plausible observationally derived values.
We find that adiabatic simulations approximately reproduce the analytical estimates of outflow properties independently of turbulence or clumpiness. However, cooling reduces the outflow energy rate by 1-2 orders of magnitude in the smooth simulations and by up to one order of magnitude in the turbulent ones. The interplay between cooling and turbulence depends on AGN luminosity: in Eddington-limited AGN, turbulence enhances the coupling between the AGN wind and the gas, while the opposite happens in lower-luminosity simulations. This occurs because dense gas clumps are resilient to low-luminosity AGN feedback but get driven away by high-luminosity AGN feedback. The overall properties of multi-phase outflowing gas in our simulations qualitatively agree with observations of multi-phase outflows. We also find that using `observable' outflow properties leads to their parameters being underestimated by a factor of a few compared with real values.
△ Less
Submitted 26 September, 2024;
originally announced September 2024.
-
Stirring the cosmic pot: how black hole feedback shapes the matter power spectrum in the Fable simulations
Authors:
Sergio Martin-Alvarez,
Vid Iršič,
Sophie Koudmani,
Martin Bourne,
Leah Bigwood,
Debora Sijacki
Abstract:
Understanding the impact of baryonic physics on cosmic structure formation is crucial for accurate cosmological predictions, especially as we usher in the era of large galaxy surveys with the Rubin Observatory as well as the Euclid and Roman Space Telescopes. A key process that can redistribute matter across a large range of scales is feedback from accreting supermassive black holes. How exactly t…
▽ More
Understanding the impact of baryonic physics on cosmic structure formation is crucial for accurate cosmological predictions, especially as we usher in the era of large galaxy surveys with the Rubin Observatory as well as the Euclid and Roman Space Telescopes. A key process that can redistribute matter across a large range of scales is feedback from accreting supermassive black holes. How exactly these active galactic nuclei (AGN) operate from sub-parsec to Mega-parsec scales however remains largely unknown. To understand this, we investigate how different AGN feedback models in the Fable simulation suite affect the cosmic evolution of the matter power spectrum (MPS).
Our analysis reveals that AGN feedback significantly suppresses clustering at scales $k \sim 10\,h\,cMpc^{-1}$, with the strongest effect at redshift $z = 0$ causing a reduction of $\sim 10\%$ with respect to the dark matter-only simulation. This is due to the efficient feedback in both radio (low Eddington ratio) and quasar (high Eddington ratio) modes in our fiducial Fable model. We find that variations of the quasar and radio mode feedback with respect to the fiducial Fable model have distinct effects on the MPS redshift evolution, with the radio mode being more effective on larger scales and later epochs. Furthermore, MPS suppression is dominated by AGN feedback effects inside haloes at $z = 0$, while for $z \gtrsim 1$ the matter distribution both inside and outside of haloes shapes the MPS suppression. Hence, future observations probing earlier cosmic times beyond $z \sim 1$ will be instrumental in constraining the nature of AGN feedback.
△ Less
Submitted 25 July, 2024;
originally announced July 2024.
-
A dormant, overmassive black hole in the early Universe
Authors:
Ignas Juodžbalis,
Roberto Maiolino,
William M. Baker,
Sandro Tacchella,
Jan Scholtz,
Francesco D'Eugenio,
Raffaella Schneider,
Alessandro Trinca,
Rosa Valiante,
Christa DeCoursey,
Mirko Curti,
Stefano Carniani,
Jacopo Chevallard,
Anna de Graaff,
Santiago Arribas,
Jake S. Bennett,
Martin A. Bourne,
Andrew J. Bunker,
Stéphane Charlot,
Brian Jiang,
Sophie Koudmani,
Michele Perna,
Brant Robertson,
Debora Sijacki,
Hannah Übler
, et al. (3 additional authors not shown)
Abstract:
Recent observations have found a large number of supermassive black holes already in place in the first few hundred million years after Big Bang. The channels of formation and growth of these early, massive black holes are not clear, with scenarios ranging from heavy seeds to light seeds experiencing bursts of high accretion rate. Here we present the detection, from the JADES survey, of broad Halp…
▽ More
Recent observations have found a large number of supermassive black holes already in place in the first few hundred million years after Big Bang. The channels of formation and growth of these early, massive black holes are not clear, with scenarios ranging from heavy seeds to light seeds experiencing bursts of high accretion rate. Here we present the detection, from the JADES survey, of broad Halpha emission in a galaxy at z=6.68, which traces a black hole with mass of ~ 4 * 10^8 Msun and accreting at a rate of only 0.02 times the Eddington limit. The host galaxy has low star formation rate (~ 1 Msun/yr, a factor of 3 below the star forming main sequence). The black hole to stellar mass ratio is ~ 0.4, i.e. about 1,000 times above the local relation, while the system is closer to the local relations in terms of dynamical mass and velocity dispersion of the host galaxy. This object is most likely the tip of the iceberg of a much larger population of dormant black holes around the epoch of reionisation. Its properties are consistent with scenarios in which short bursts of super-Eddington accretion have resulted in black hole overgrowth and massive gas expulsion from the accretion disk; in between bursts, black holes spend most of their life in a dormant state.
△ Less
Submitted 6 March, 2024;
originally announced March 2024.
-
A unified accretion disc model for supermassive black holes in galaxy formation simulations: method and implementation
Authors:
Sophie Koudmani,
Rachel S. Somerville,
Debora Sijacki,
Martin A. Bourne,
Yan-Fei Jiang,
Kasar Profit
Abstract:
It is well established that supermassive black hole (SMBH) feedback is crucial for regulating the evolution of massive, if not all, galaxies. However, modelling the interplay between SMBHs and their host galaxies is challenging due to the vast dynamic range. Previous simulations have utilized simple subgrid models for SMBH accretion, while recent advancements track the properties of the unresolved…
▽ More
It is well established that supermassive black hole (SMBH) feedback is crucial for regulating the evolution of massive, if not all, galaxies. However, modelling the interplay between SMBHs and their host galaxies is challenging due to the vast dynamic range. Previous simulations have utilized simple subgrid models for SMBH accretion, while recent advancements track the properties of the unresolved accretion disc, usually based on the thin $α$-disc model. However, this neglects accretion in the radiatively inefficient regime, expected to occur through a thick disc for a significant portion of an SMBH's lifetime. To address this, we present a novel 'unified' accretion disc model for SMBHs, harnessing results from the analytical advection-dominated inflow-outflow solution (ADIOS) model and state-of-the-art GR(R)MHD simulations. Going from low to high Eddington ratios, our model transitions from an ADIOS flow to a thin $α$-disc via a truncated disc, incorporating self-consistently SMBH spin evolution due to Lense-Thirring precession. Utilizing the moving mesh code AREPO, we perform simulations of single and binary SMBHs within gaseous discs to validate our model and assess its impact. The disc state significantly affects observable luminosities, and we predict markedly different electromagnetic counterparts in SMBH binaries. Crucially, the assumed disc model shapes SMBH spin magnitudes and orientations, parameters that gravitational wave observatories like LISA and IPTA are poised to constrain. Our simulations emphasize the importance of accurately modelling SMBH accretion discs and spin evolution, as they modulate the available accretion power, profoundly shaping the interaction between SMBHs and their host galaxies.
△ Less
Submitted 3 July, 2024; v1 submitted 13 December, 2023;
originally announced December 2023.
-
Dynamics and spin alignment in massive, gravito-turbulent circumbinary discs around supermassive black hole binaries
Authors:
Martin A. Bourne,
Davide Fiacconi,
Debora Sijacki,
Joanna M. Piotrowska,
Sophie Koudmani
Abstract:
Parsec-scale separation supermassive black hole binaries in the centre of gas-rich galaxy merger remnants could be surrounded by massive circumbinary discs (CBDs). Black hole mass and spin evolution during the gas-rich binary inspiral are crucial in determining the direction and power of relativistic jets that radio observations with LOFAR and SKAO will probe, and for predicting gravitational wave…
▽ More
Parsec-scale separation supermassive black hole binaries in the centre of gas-rich galaxy merger remnants could be surrounded by massive circumbinary discs (CBDs). Black hole mass and spin evolution during the gas-rich binary inspiral are crucial in determining the direction and power of relativistic jets that radio observations with LOFAR and SKAO will probe, and for predicting gravitational wave (GW) emission that IPTA and LISA will measure. We present 3D hydrodynamic simulations capturing gas-rich, self-gravitating CBDs around a $2\times 10^6$ M$_{\odot}$ supermassive black hole binary, that probe different mass ratios, eccentricities and inclinations. We employ a sub-grid Shakura-Sunyaev accretion disc to self-consistently model black hole mass and spin evolution together with super-Lagrangian refinement techniques to resolve gas flows, streams and mini-discs within the cavity, which play a fundamental role in torquing and feeding the binary. We find that higher mass ratio and eccentric binaries result in larger cavities, while retrograde binaries result in smaller cavities. All of the simulated binaries are expected to shrink with net gravitational torques being negative. Unlike previous simulations, we do not find preferential accretion onto the secondary black hole. This implies smaller chirp masses at coalescence and hence a weaker GW background. Critically this means that spin-alignment is faster than the binary inspiral timescale even for low mass ratios. When considering initially misaligned systems, the orientation of the mini-discs around each black hole can vary significantly. We discuss the implications of this behaviour for black hole spin alignment and highlight the need for broader parameter space studies of misaligned systems to understand the impact on black hole recoil velocities.
△ Less
Submitted 12 September, 2024; v1 submitted 28 November, 2023;
originally announced November 2023.
-
Simulations of spin-driven AGN jets in gas-rich galaxy mergers
Authors:
Rosie Y. Talbot,
Debora Sijacki,
Martin A. Bourne
Abstract:
In this work, we use hydrodynamical simulations to explore the effects of kinetic AGN jet feedback on the progression and outcome of the major merger of two isolated, gas-rich galaxies. We present simulations that use the moving-mesh code AREPO to follow the progression of the merger through first passage and up to the final coalescence, modelling the black holes at the centres of both of the merg…
▽ More
In this work, we use hydrodynamical simulations to explore the effects of kinetic AGN jet feedback on the progression and outcome of the major merger of two isolated, gas-rich galaxies. We present simulations that use the moving-mesh code AREPO to follow the progression of the merger through first passage and up to the final coalescence, modelling the black holes at the centres of both of the merging galaxies using our prescription for black hole accretion via an $α$-disc and feedback in the form of a spin-driven jet. We find that the jets drive large-scale, multiphase outflows which launch large quantities of cold gas out to distances greater than 100 kpc and with velocities that reach $\sim 2500 \, {\rm km \, s^{-1}}$. Gas in the outflows that decelerates, cools and falls back on the galaxies can provide a rich source of fuel for the black hole, leading to intense episodes of jet activity in which the jet can become significantly misaligned. The presence of AGN jets affects the growth of the stellar component: star formation is moderately suppressed at all times during the merger and the peak of the star formation rate, attained during the final coalescence of the galaxies, is reduced by a factor of $\sim 2$. Analysis of simulations such as these will play a central role in making precise predictions for multimessenger investigations of dual radio-AGN, which next-generation observational facilities such as LISA, Athena and SKA will make possible.
△ Less
Submitted 16 February, 2024; v1 submitted 12 June, 2023;
originally announced June 2023.
-
Recent Progress in Modeling the Macro- and Micro-Physics of Radio Jet Feedback in Galaxy Clusters
Authors:
Martin A. Bourne,
Hsiang-Yi Karen Yang
Abstract:
Radio jets and the lobes they inflate are common in cool-core clusters and are known to play a critical role in regulating the heating and cooling of the intracluster medium (ICM). This is an inherently multi-scale problem, and much effort has been made to understand the processes governing the inflation of lobes and their impact on the cluster, as well as the impact of the environment on the jet-…
▽ More
Radio jets and the lobes they inflate are common in cool-core clusters and are known to play a critical role in regulating the heating and cooling of the intracluster medium (ICM). This is an inherently multi-scale problem, and much effort has been made to understand the processes governing the inflation of lobes and their impact on the cluster, as well as the impact of the environment on the jet-ICM interaction, on both macro- and microphysical scales. Developments of new numerical techniques and improving computational resources have seen simulations of jet feedback in galaxy clusters become ever more sophisticated. This ranges from modelling ICM plasma physics processes such as the effects of magnetic fields, cosmic rays and viscosity to including jet feedback in cosmologically evolved cluster environments in which the ICM thermal and dynamic properties are shaped by large-scale structure formation. In this review, we discuss the progress made over the last ~decade in capturing both the macro- and microphysical processes in numerical simulations, highlighting both the current state of the field as well as open questions and potential ways in which these questions can be addressed in the future.
△ Less
Submitted 13 June, 2023; v1 submitted 28 April, 2023;
originally announced May 2023.
-
Blandford-Znajek jets in galaxy formation simulations: exploring the diversity of outflows produced by spin-driven AGN jets in Seyfert galaxies
Authors:
Rosie Y. Talbot,
Debora Sijacki,
Martin A. Bourne
Abstract:
Recent observations of Seyfert galaxies indicate that low power, misaligned jets can undergo significant interaction with the gas in the galactic disc and may be able to drive large-scale, multiphase outflows. We apply our novel sub-grid model for Blandford-Znajek jets to simulations of the central regions of Seyferts, in which a black hole is embedded in a dense, sub-kpc circumnuclear disc (CND)…
▽ More
Recent observations of Seyfert galaxies indicate that low power, misaligned jets can undergo significant interaction with the gas in the galactic disc and may be able to drive large-scale, multiphase outflows. We apply our novel sub-grid model for Blandford-Znajek jets to simulations of the central regions of Seyferts, in which a black hole is embedded in a dense, sub-kpc circumnuclear disc (CND) and surrounded by a dilute circumgalactic medium (CGM). We find that the variability of the accretion flow is highly sensitive both to the jet power and to the CND thermodynamics and, ultimately, is determined by the complex interplay between jet-driven outflows and backflows. Even at moderate Eddington ratios, AGN jets are able to significantly alter the thermodynamics and kinematics of CNDs and entrain up to 10% of their mass in the outflow. Mass outflow rates and kinetic powers of the warm outflowing component are in agreement with recent observations for black holes with similar bolometric luminosities, with outflow velocities that are able to reach 500 km/s. Depending on their power and direction, jets are able to drive a wide variety of large-scale outflows, ranging from light, hot and collimated structures to highly mass-loaded, multiphase, bipolar winds. This diversity of jet-driven outflows highlights the importance of applying physically motivated models of AGN feedback to realistic galaxy formation contexts. Such simulations will play a crucial role in accurately interpreting the wealth of data that next generation facilities such as JWST, SKA and Athena will provide.
△ Less
Submitted 7 June, 2022; v1 submitted 2 November, 2021;
originally announced November 2021.
-
Blandford-Znajek jets in galaxy formation simulations: method and implementation
Authors:
Rosie Y. Talbot,
Martin A. Bourne,
Debora Sijacki
Abstract:
Jets launched by active galactic nuclei (AGN) are believed to play a significant role in shaping the properties of galaxies and provide an energetically viable mechanism through which galaxies can become quenched. Here we present a novel AGN feedback model, which we have incorporated into the AREPO code, that evolves the black hole mass and spin as the accretion flow proceeds through a thin $α$-di…
▽ More
Jets launched by active galactic nuclei (AGN) are believed to play a significant role in shaping the properties of galaxies and provide an energetically viable mechanism through which galaxies can become quenched. Here we present a novel AGN feedback model, which we have incorporated into the AREPO code, that evolves the black hole mass and spin as the accretion flow proceeds through a thin $α$-disc which we self-consistently couple to a Blandford-Znajek jet. We apply our model to the central region of a typical radio-loud Seyfert galaxy embedded in a hot circumgalactic medium (CGM). We find that jets launched into high pressure environments thermalise efficiently due to the formation of recollimation shocks and the vigorous instabilities that these shocks excite increase the efficiency of the mixing of CGM and jet material. The beams of more overpressured jets, however, are not as readily disrupted by instabilities so the majority of the momentum flux at the jet base is retained out to the head, where the jet terminates in a reverse shock. All jets entrain a significant amount of cold circumnuclear disc material which, while energetically insignificant, dominates the lobe mass together with the hot, entrained CGM material. The jet power evolves significantly due to effective self-regulation by the black hole, fed by secularly-driven, intermittent mass flows. The direction of jets launched directly into the circumnuclear disc changes considerably due to effective Bardeen-Petterson torquing. Interestingly, these jets obliterate the innermost regions of the disc and drive large-scale, multi-phase, turbulent, bipolar outflows.
△ Less
Submitted 29 July, 2021; v1 submitted 20 November, 2020;
originally announced November 2020.
-
AGN jet feedback on a moving mesh: gentle cluster heating by weak shocks and lobe disruption
Authors:
Martin A. Bourne,
Debora Sijacki
Abstract:
While there is overwhelming observational evidence of AGN-driven jets in galaxy clusters and groups, if and how the jet energy is delivered to the ambient medium remains unanswered. Here we perform very high resolution AGN jet simulations within a live, cosmologically evolved cluster with the moving mesh code AREPO. We find that mock X-ray and radio lobe properties are in good agreement with obser…
▽ More
While there is overwhelming observational evidence of AGN-driven jets in galaxy clusters and groups, if and how the jet energy is delivered to the ambient medium remains unanswered. Here we perform very high resolution AGN jet simulations within a live, cosmologically evolved cluster with the moving mesh code AREPO. We find that mock X-ray and radio lobe properties are in good agreement with observations with different power jets transitioning from FR-I to FR-II-like morphologies. During the lobe inflation phase, heating by both internal and bow shocks contributes to lobe energetics, and $\sim 40$ per cent of the feedback energy goes into the $PdV$ work done by the expanding lobes. Low power jets are more likely to simply displace gas during lobe inflation, but higher power jets become more effective at driving shocks and heating the intracluster medium (ICM), although shocks rarely exceed $\mathcal{M}\sim 2-3$. Once the lobe inflation phase ceases, cluster weather significantly impacts the lobe evolution. Lower power jet lobes are more readily disrupted and mixed with the ICM, depositing up to $\sim 70$ per cent of the injected energy, however, ultimately the equivalent of $\gtrsim 50$ per cent of the feedback energy ends up as potential energy of the system. Even though the mean ICM entropy is increased up to $80$~Myr after the jets switch off, AGN heating is gentle, inducing no large variations in cluster radial profiles in accord with observations.
△ Less
Submitted 14 June, 2021; v1 submitted 28 August, 2020;
originally announced August 2020.
-
An ultra-wide bandwidth (704 to 4032 MHz) receiver for the Parkes radio telescope
Authors:
G. Hobbs,
R. N. Manchester,
A. Dunning,
A. Jameson,
P. Roberts,
D. George,
J. A. Green,
J. Tuthill,
L. Toomey,
J. F. Kaczmarek,
S. Mader,
M. Marquarding,
A. Ahmed,
S. W. Amy,
M. Bailes,
R. Beresford,
N. D. R. Bhat,
D. C. -J. Bock,
M. Bourne,
M. Bowen,
M. Brothers,
A. D. Cameron,
E. Carretti,
N. Carter,
S. Castillo
, et al. (47 additional authors not shown)
Abstract:
We describe an ultra-wide-bandwidth, low-frequency receiver ("UWL") recently installed on the Parkes radio telescope. The receiver system provides continuous frequency coverage from 704 to 4032 MHz. For much of the band (~60%) the system temperature is approximately 22K and the receiver system remains in a linear regime even in the presence of strong mobile phone transmissions. We discuss the scie…
▽ More
We describe an ultra-wide-bandwidth, low-frequency receiver ("UWL") recently installed on the Parkes radio telescope. The receiver system provides continuous frequency coverage from 704 to 4032 MHz. For much of the band (~60%) the system temperature is approximately 22K and the receiver system remains in a linear regime even in the presence of strong mobile phone transmissions. We discuss the scientific and technical aspects of the new receiver including its astronomical objectives, as well as the feed, receiver, digitiser and signal-processor design. We describe the pipeline routines that form the archive-ready data products and how those data files can be accessed from the archives. The system performance is quantified including the system noise and linearity, beam shape, antenna efficiency, polarisation calibration and timing stability.
△ Less
Submitted 2 November, 2019;
originally announced November 2019.
-
AGN jet feedback on a moving mesh: lobe energetics and X-ray properties in a realistic cluster environment
Authors:
Martin A. Bourne,
Debora Sijacki,
Ewald Puchwein
Abstract:
Jet feedback from active galactic nuclei (AGN) harboured by brightest cluster galaxies is expected to play a fundamental role in regulating cooling in the intracluster medium (ICM). While observations and theory suggest energy within jet lobes balances ICM radiative losses, the modus operandi of energy communication with the ICM remains unclear. We present simulations of very high resolution AGN-d…
▽ More
Jet feedback from active galactic nuclei (AGN) harboured by brightest cluster galaxies is expected to play a fundamental role in regulating cooling in the intracluster medium (ICM). While observations and theory suggest energy within jet lobes balances ICM radiative losses, the modus operandi of energy communication with the ICM remains unclear. We present simulations of very high resolution AGN-driven jets launching in a live, cosmological galaxy cluster, within the moving mesh code AREPO. As the jet propagates through the ICM the majority of its energy, which is initially in the kinetic form, thermalizes quickly through internal shocks and inflates lobes of very hot gas. The jets effectively heat the cluster core, with $PdV$ work and weather-aided mixing being the main channels of energy transfer from the lobes to the ICM, while strong shocks and turbulence are subdominant. We additionally present detailed mock X-ray maps at different stages of evolution, revealing clear cavities surrounded by X-ray bright rims, with lobes being detectable for up to $\sim 10^8$ yr even when magnetic draping is ineffective. We find bulk motions in the cluster can significantly affect lobe propagation, offsetting them from the jet direction and imparting bulk velocities that can dominate over the buoyantly rising motion.
△ Less
Submitted 16 October, 2019; v1 submitted 30 January, 2019;
originally announced January 2019.
-
Fast and energetic AGN-driven outflows in simulated dwarf galaxies
Authors:
Sophie Koudmani,
Debora Sijacki,
Martin A. Bourne,
Matthew C. Smith
Abstract:
The systematic analysis of optical large-scale surveys has revealed a population of dwarf galaxies hosting AGN, which have been confirmed by X-ray follow-up observations. Recently, the MaNGA survey identified six dwarf galaxies that appear to have an AGN that is preventing on-going star formation. It is therefore timely to study the physical properties of dwarf galaxies, in particular whether the…
▽ More
The systematic analysis of optical large-scale surveys has revealed a population of dwarf galaxies hosting AGN, which have been confirmed by X-ray follow-up observations. Recently, the MaNGA survey identified six dwarf galaxies that appear to have an AGN that is preventing on-going star formation. It is therefore timely to study the physical properties of dwarf galaxies, in particular whether the presence of an AGN can affect their evolution. Using the moving mesh code AREPO, we have investigated different models of AGN activity, ranging from simple energy-driven spherical winds to collimated, mass-loaded, bipolar outflows in high resolution simulations of isolated dwarf galaxies hosting an active black hole. Our simulations also include a novel implementation of star formation and mechanical supernova (SN) feedback. We find that AGN outflows have a small but systematic effect on the central star formation rates (SFRs) for all set-ups explored, while substantial effects on the global SFR are only obtained with strong SNe and a sustained high-luminosity AGN with an isotropic wind. This suggests that AGN feedback in dwarf galaxies is unlikely to directly regulate their global SFRs. There is, however, a significant effect on outflow properties, which are notably enhanced by the AGN to much higher outflow temperatures and velocities, in agreement with kinematic signatures from the MaNGA survey. This indicates that AGN may play an indirect role in regulating the baryon cycle in dwarf galaxies by hindering cosmic gas inflows.
△ Less
Submitted 8 January, 2019; v1 submitted 11 December, 2018;
originally announced December 2018.
-
Cold Molecular Outflows in the Local Universe
Authors:
A. Fluetsch,
R. Maiolino,
S. Carniani,
A. Marconi,
C. Cicone,
M. A. Bourne,
T. Costa,
A. C. Fabian,
W. Ishibashi,
G. Venturi
Abstract:
We study molecular outflows in a sample of 45 local galaxies, both star forming and AGN, primarily by using CO data from the ALMA archive and from the literature. For a subsample we also compare the molecular outflow with the ionized and neutral atomic phases. We infer an empirical analytical function relating the outflow rate simultaneously to the SFR, $L_{\rm AGN}$, and galaxy stellar mass; this…
▽ More
We study molecular outflows in a sample of 45 local galaxies, both star forming and AGN, primarily by using CO data from the ALMA archive and from the literature. For a subsample we also compare the molecular outflow with the ionized and neutral atomic phases. We infer an empirical analytical function relating the outflow rate simultaneously to the SFR, $L_{\rm AGN}$, and galaxy stellar mass; this relation is much tighter than the relations with the individual quantities. The outflow kinetic power shows a larger scatter than in previous, more biased studies, spanning from 0.1 to 5 per cent of $L_{\rm AGN}$, while the momentum rate ranges from 1 to 30 times $L_{\rm AGN}/c$, indicating that these outflows can be both energy-driven, but with a broad range of coupling efficiencies with the ISM, and radiation pressure-driven. For about 10 per cent of the objects the outflow energetics significantly exceed the maximum theoretical values; we interpret these as 'fossil outflows' resulting from activity of a past strong AGN, which has now faded. We estimate that, in the stellar mass range probed here ($>$ 10$^{10}~\rm M_{\odot}$), less than 5 per cent of the outflowing gas escapes the galaxy. The molecular gas depletion time associated with the outflow can be as short as a few million years in powerful AGN, however, the total gas (H$_2$+HI) depletion times are much longer. Altogether, our findings suggest that even AGN-driven outflows might be relatively ineffective in clearing galaxies of their entire gas content, although they are likely capable of clearing and quenching the central region.
△ Less
Submitted 5 February, 2019; v1 submitted 14 May, 2018;
originally announced May 2018.
-
AGN jet feedback on a moving mesh: cocoon inflation, gas flows and turbulence
Authors:
Martin A. Bourne,
Debora Sijacki
Abstract:
In many observed galaxy clusters, jets launched by the accretion process on to supermassive black holes, inflate large-scale cavities filled with energetic, relativistic plasma. This process is thought to be responsible for regulating cooling losses, thus moderating the inflow of gas on to the central galaxy, quenching further star formation and maintaining the galaxy in a red and dead state. In t…
▽ More
In many observed galaxy clusters, jets launched by the accretion process on to supermassive black holes, inflate large-scale cavities filled with energetic, relativistic plasma. This process is thought to be responsible for regulating cooling losses, thus moderating the inflow of gas on to the central galaxy, quenching further star formation and maintaining the galaxy in a red and dead state. In this paper, we implement a new jet feedback scheme into the moving mesh-code arepo, contrast different jet injection techniques and demonstrate the validity of our implementation by comparing against simple analytical models. We find that jets can significantly affect the intracluster medium (ICM), offset the overcooling through a number of heating mechanisms, as well as drive turbulence, albeit within the jet lobes only. Jet-driven turbulence is, however, a largely ineffective heating source and is unlikely to dominate the ICM heating budget even if the jet lobes efficiently fill the cooling region, as it contains at most only a few per cent of the total injected energy. We instead show that the ICM gas motions, generated by orbiting substructures, while inefficient at heating the ICM, drive large-scale turbulence and when combined with jet feedback, result in line-of-sight velocities and velocity dispersions consistent with the Hitomi observations of the Perseus cluster.
△ Less
Submitted 1 November, 2017; v1 submitted 22 May, 2017;
originally announced May 2017.
-
Do AGN outflows quench or enhance star formation?
Authors:
Kastytis Zubovas,
Martin A. Bourne
Abstract:
AGN outflows can remove large quantities of gas from their host galaxy spheroids, potentially shutting off star formation. On the other hand, they can compress this gas, potentially enhancing or triggering star formation, at least for short periods. We present a set of idealised simulations of AGN outflows affecting turbulent gas spheres, and investigate the effect of the outflow and the AGN radia…
▽ More
AGN outflows can remove large quantities of gas from their host galaxy spheroids, potentially shutting off star formation. On the other hand, they can compress this gas, potentially enhancing or triggering star formation, at least for short periods. We present a set of idealised simulations of AGN outflows affecting turbulent gas spheres, and investigate the effect of the outflow and the AGN radiation field upon gas fragmentation. We show that AGN outflows of sufficient luminosity shut off fragmentation while the nucleus is active, but gas compression results in a burst of fragmentation after the AGN switches off. Self-shielding of gas against the AGN radiation field allows some fragmentation to occur during outbursts, but too much shielding results in a lower overall fragmentation rate. For our idealised simulation setup, there is a critical AGN luminosity which results in the highest fragmentation rate, with outflows being too efficient at removing gas when $L > L_{\rm crit}$ and not efficient enough to compress the gas to high densities otherwise. These results, although preliminary, suggest that the interaction between AGN and star formation in their host galaxies is particularly complex and requires careful study in order to interpret observations correctly.
△ Less
Submitted 26 April, 2017; v1 submitted 31 March, 2017;
originally announced March 2017.
-
Too small to succeed: the difficulty of sustaining star formation in low-mass haloes
Authors:
Claire Cashmore,
Mark Wilkinson,
Chris Power,
Martin Bourne
Abstract:
We present high resolution simulations of an isolated dwarf spheroidal (dSph) galaxy between redshifts $z\sim10$ and $z\sim 4$, the epoch when several Milky Way dSph satellites experienced extended star formation, in order to understand in detail the physical processes which affect a low-mass halo's ability to retain gas. It is well-established that supernova feedback is very effective at expellin…
▽ More
We present high resolution simulations of an isolated dwarf spheroidal (dSph) galaxy between redshifts $z\sim10$ and $z\sim 4$, the epoch when several Milky Way dSph satellites experienced extended star formation, in order to understand in detail the physical processes which affect a low-mass halo's ability to retain gas. It is well-established that supernova feedback is very effective at expelling gas from a $3\times 10^7$M$_\odot$ halo, the mass of a typical redshift 10 progenitor of a redshift 0 halo with mass $\sim10^9$M$_\odot$. We investigate the conditions under which such a halo is able to retain sufficient high-density gas to support extended star formation. In particular, we explore the effects of: an increased relative concentration of the gas compared to the dark matter; a higher concentration dark matter halo; significantly lower supernova rates; enhanced metal cooling due to enrichment from earlier supernovae. We show that disk-like gas distributions retain more gas than spherical ones, primarily due to the shorter gas cooling times in the disk. However, a significant reduction in the number of supernovae compared to that expected for a standard initial mass function is still needed to allow the retention of high density gas. We conclude that the progenitors of the observed dSphs would only have retained the gas required to sustain star formation if their mass, concentration and gas morphology were already unusual for those of a dSph-mass halo progenitor by a redshift of 10.
△ Less
Submitted 6 February, 2017;
originally announced February 2017.
-
A low upper mass limit for the central black hole in the late-type galaxy NGC 4414
Authors:
Sabine Thater,
Davor Krajnović,
Martin A. Bourne,
Michele Cappellari,
Tim de Zeeuw,
Eric Emsellem,
John Magorrian,
Richard M. McDermid,
Marc Sarzi,
Glenn van de Ven
Abstract:
We present our mass estimate of the central black hole in the isolated spiral galaxy NGC 4414. Using natural guide star adaptive optics assisted observations with the Gemini Near-Infrared Integral Field Spectrometer (NIFS) and the natural seeing Gemini Multi-Object Spectrographs-North (GMOS), we derived two-dimensional stellar kinematic maps of NGC 4414 covering the central 1.5 arcsec and 10 arcse…
▽ More
We present our mass estimate of the central black hole in the isolated spiral galaxy NGC 4414. Using natural guide star adaptive optics assisted observations with the Gemini Near-Infrared Integral Field Spectrometer (NIFS) and the natural seeing Gemini Multi-Object Spectrographs-North (GMOS), we derived two-dimensional stellar kinematic maps of NGC 4414 covering the central 1.5 arcsec and 10 arcsec, respectively, at a NIFS spatial resolution of 0.13 arcsec. The kinematic maps reveal a regular rotation pattern and a central velocity dispersion dip down to around 105 km/s. We constructed dynamical methods using two different methods: Jeans anisotropic dynamical modeling and axisymmetric Schwarzschild modeling. Both modeling methods give consistent results, but we cannot constrain the lower mass limit and only measure an upper limit for the black hole mass of Mbh= 1.56 x 10^6 Msun(at 3 sigma level) which is at least 1 sigma below the recent Mbh-sigma_e relations. Further tests with dark matter, mass-to-light ratio variation and different light models confirm that our results are not dominated by uncertainties. The derived upper limit mass is not only below the Mbh-sigma_e relation, but is also five times lower than the lower limit black hole mass anticipated from the resolution limit of the sphere of influence. This proves that via high quality integral field data we are now able to push black hole measurements down to at least five times less than the resolution limit.
△ Less
Submitted 24 October, 2016; v1 submitted 13 October, 2016;
originally announced October 2016.
-
A simple way to improve AGN feedback prescription in SPH simulations
Authors:
Kastytis Zubovas,
Martin A. Bourne,
Sergei Nayakshin
Abstract:
AGN feedback is an important ingredient in galaxy evolution, however its treatment in numerical simulations is necessarily approximate, requiring subgrid prescriptions due to the dynamical range involved in the calculations. We present a suite of SPH simulations designed to showcase the importance of the choice of a particular subgrid prescription for AGN feedback. We concentrate on two approaches…
▽ More
AGN feedback is an important ingredient in galaxy evolution, however its treatment in numerical simulations is necessarily approximate, requiring subgrid prescriptions due to the dynamical range involved in the calculations. We present a suite of SPH simulations designed to showcase the importance of the choice of a particular subgrid prescription for AGN feedback. We concentrate on two approaches to treating wide-angle AGN outflows: thermal feedback, where thermal and kinetic energy is injected into the gas surrounding the SMBH particle, and virtual particle feedback, where energy is carried by tracer particles radially away from the AGN. We show that the latter model produces a far more complex structure around the SMBH, which we argue is a more physically correct outcome. We suggest a simple improvement to the thermal feedback model - injecting the energy into a cone, rather than spherically symmetrically - and show that this markedly improves the agreement between the two prescriptions, without requiring any noticeable increase in the computational cost of the simulation.
△ Less
Submitted 17 December, 2015;
originally announced December 2015.
-
Simulating feedback from nuclear clusters: the impact of multiple sources
Authors:
Martin A. Bourne,
Chris Power
Abstract:
Nuclear star clusters (NCs) are found to exist in the centres of many galaxies and appear to follow scaling relations similar to those of super-massive black holes. Previous analytical work has suggested that such relations are a consequence of feedback regulated growth. We explore this idea using high resolution hydrodynamical simulations, focusing on the validity of the simplifying assumptions m…
▽ More
Nuclear star clusters (NCs) are found to exist in the centres of many galaxies and appear to follow scaling relations similar to those of super-massive black holes. Previous analytical work has suggested that such relations are a consequence of feedback regulated growth. We explore this idea using high resolution hydrodynamical simulations, focusing on the validity of the simplifying assumptions made in analytical models. In particular, we investigate feedback emanating from multiple stellar sources rather than from a single source, as is usually assumed, and show that collisions betweens shells of gas swept up by feedback leads to momentum cancellation and the formation of high density clumps and filaments. This high density material is resistant both to expulsion from the galaxy potential and to disruption by feedback; if it falls back onto the NC, we expect the gas to be available for further star formation or for feeding a central black hole. We also note our results may have implications for the evolution of globular clusters and stellar clusters in high redshift dark matter halos.
△ Less
Submitted 19 October, 2015;
originally announced October 2015.
-
The resolution bias: low resolution feedback simulations are better at destroying galaxies
Authors:
Martin A. Bourne,
Kastytis Zubovas,
Sergei Nayakshin
Abstract:
Feedback from super-massive black holes (SMBHs) is thought to play a key role in regulating the growth of host galaxies. Cosmological and galaxy formation simulations using smoothed particle hydrodynamics (SPH), which usually use a fixed mass for SPH particles, often employ the same sub-grid Active galactic nuclei (AGN) feedback prescription across a range of resolutions. It is thus important to a…
▽ More
Feedback from super-massive black holes (SMBHs) is thought to play a key role in regulating the growth of host galaxies. Cosmological and galaxy formation simulations using smoothed particle hydrodynamics (SPH), which usually use a fixed mass for SPH particles, often employ the same sub-grid Active galactic nuclei (AGN) feedback prescription across a range of resolutions. It is thus important to ask how the impact of the simulated AGN feedback on a galaxy changes when only the numerical resolution (the SPH particle mass) changes. We present a suite of simulations modelling the interaction of an AGN outflow with the ambient turbulent and clumpy interstellar medium (ISM) in the inner part of the host galaxy at a range of mass resolutions. We find that, with other things being equal, degrading the resolution leads to feedback becoming more efficient at clearing out all gas in its path. For the simulations presented here, the difference in the mass of the gas ejected by AGN feedback varies by more than a factor of ten between our highest and lowest resolution simulations. This happens because feedback-resistant high density clumps are washed out at low effective resolutions. We also find that changes in numerical resolution lead to undesirable artifacts in how the AGN feedback affects the AGN immediate environment.
△ Less
Submitted 29 July, 2015;
originally announced July 2015.
-
Black hole feedback in a multiphase interstellar medium
Authors:
Martin A. Bourne,
Sergei Nayakshin,
Alexander Hobbs
Abstract:
Ultrafast outflows (UFOs) from supermassive black holes (SMBHs) are thought to regulate the growth of SMBHs and host galaxies, resulting in a number of observational correlations. We present high-resolution numerical simulations of the impact of a thermalized UFO on the ambient gas in the inner part of the host galaxy. Our results depend strongly on whether the gas is homogeneous or clumpy. In the…
▽ More
Ultrafast outflows (UFOs) from supermassive black holes (SMBHs) are thought to regulate the growth of SMBHs and host galaxies, resulting in a number of observational correlations. We present high-resolution numerical simulations of the impact of a thermalized UFO on the ambient gas in the inner part of the host galaxy. Our results depend strongly on whether the gas is homogeneous or clumpy. In the former case all of the ambient gas is driven outward rapidly as expected based on commonly used energy budget arguments, while in the latter the flows of mass and energy decouple. Carrying most of the energy, the shocked UFO escapes from the bulge via paths of least resistance, taking with it only the low-density phase of the host. Most of the mass is however in the high-density phase, and is affected by the UFO much less strongly, and may even continue to flow inwards. We suggest that the UFO energy leakage through the pores in the multiphase interstellar medium (ISM) may explain why observed SMBHs are so massive despite their overwhelmingly large energy production rates. The multiphase ISM effects reported here are probably under-resolved in cosmological simulations but may be included in prescriptions for active galactic nuclei feedback in future simulations and in semi-analytical models.
△ Less
Submitted 22 May, 2014;
originally announced May 2014.
-
Inverse Compton X-ray signature of AGN feedback
Authors:
Martin A. Bourne,
Sergei Nayakshin
Abstract:
Bright AGN frequently show ultra-fast outflows (UFOs) with outflow velocities vout ~0.1c. These outflows may be the source of AGN feedback on their host galaxies sought by galaxy formation modellers. The exact effect of the outflows on the ambient galaxy gas strongly depends on whether the shocked UFOs cool rapidly or not. This in turn depends on whether the shocked electrons share the same temper…
▽ More
Bright AGN frequently show ultra-fast outflows (UFOs) with outflow velocities vout ~0.1c. These outflows may be the source of AGN feedback on their host galaxies sought by galaxy formation modellers. The exact effect of the outflows on the ambient galaxy gas strongly depends on whether the shocked UFOs cool rapidly or not. This in turn depends on whether the shocked electrons share the same temperature as ions (one temperature regime; 1T) or decouple (2T), as has been recently suggested. Here we calculate the Inverse Compton spectrum emitted by such shocks, finding a broad feature potentially detectable either in mid-to-high energy X-rays (1T case) or only in the soft X-rays (2T). We argue that current observations of AGN do not seem to show evidence for the 1T component. The limits on the 2T emission are far weaker, and in fact it is possible that the observed soft X-ray excess of AGN is partially or fully due to the 2T shock emission. This suggests that UFOs are in the energy-driven regime outside the central few pc, and must pump considerable amounts of not only momentum but also energy into the ambient gas. We encourage X-ray observers to look for the Inverse Compton components calculated here in order to constrain AGN feedback models further.
△ Less
Submitted 13 September, 2013; v1 submitted 11 June, 2013;
originally announced June 2013.