-
Active galactic nucleus jet feedback in hydrostatic halos
Authors:
Rainer Weinberger,
Kung-Yi Su,
Kristian Ehlert,
Christoph Pfrommer,
Lars Hernquist,
Greg L. Bryan,
Volker Springel,
Yuan Li,
Blakesley Burkhart,
Ena Choi,
Claude-André Faucher-Giguère
Abstract:
Feedback driven by jets from active galactic nuclei is believed to be responsible for reducing cooling flows in cool-core galaxy clusters. We use simulations to model feedback from hydrodynamic jets in isolated halos. While the jet propagation converges only after the diameter of the jet is well resolved, reliable predictions about the effects these jets have on the cooling time distribution funct…
▽ More
Feedback driven by jets from active galactic nuclei is believed to be responsible for reducing cooling flows in cool-core galaxy clusters. We use simulations to model feedback from hydrodynamic jets in isolated halos. While the jet propagation converges only after the diameter of the jet is well resolved, reliable predictions about the effects these jets have on the cooling time distribution function only require resolutions sufficient to keep the jet-inflated cavities stable. Comparing different model variations, as well as an independent jet model using a different hydrodynamics code, we show that the dominant uncertainties are the choices of jet properties within a given model. Independent of implementation, we find that light, thermal jets with low momentum flux tend to delay the onset of a cooling flow more efficiently on a $50$ Myr timescale than heavy, kinetic jets. The delay of the cooling flow originates from a displacement and boost in entropy of the central gas. If the jet kinetic luminosity depends on accretion rate, collimated, light, hydrodynamic jets are able to reduce cooling flows in halos, without a need for jet precession or wide opening angles. Comparing the jet feedback with a `kinetic wind' implementation shows that equal amounts of star formation rate reduction can be achieved by different interactions with the halo gas: the jet has a larger effect on the hot halo gas while leaving the denser, star forming phase in place, while the wind acts more locally on the star forming phase, which manifests itself in different time-variability properties.
△ Less
Submitted 5 May, 2023; v1 submitted 21 November, 2022;
originally announced November 2022.
-
Self-regulated AGN feedback of light jets in cool-core galaxy clusters
Authors:
Kristian Ehlert,
Rainer Weinberger,
Christoph Pfrommer,
Rüdiger Pakmor,
Volker Springel
Abstract:
Heating from active galactic nuclei (AGN) is thought to stabilize cool-core clusters, limiting star formation and cooling flows. We employ radiative magneto-hydrodynamic (MHD) simulations to model light AGN jet feedback with different accretion modes (Bondi-Hoyle-Lyttleton and cold accretion) in an idealised Perseus-like cluster. Independent of the probed accretion model, accretion efficiency, jet…
▽ More
Heating from active galactic nuclei (AGN) is thought to stabilize cool-core clusters, limiting star formation and cooling flows. We employ radiative magneto-hydrodynamic (MHD) simulations to model light AGN jet feedback with different accretion modes (Bondi-Hoyle-Lyttleton and cold accretion) in an idealised Perseus-like cluster. Independent of the probed accretion model, accretion efficiency, jet density and resolution, the cluster self-regulates with central entropies and cooling times consistent with observed cool-core clusters in this non-cosmological setting. We find that increased jet efficiencies lead to more intermittent jet powers and enhanced star formation rates. Our fiducial low-density jets can easily be deflected by orbiting cold gaseous filaments, which redistributes angular momentum and leads to more extended cold gas distributions and isotropic bubble distributions. In comparison to our fiducial low momentum-density jets, high momentum-density jets heat less efficiently and enable the formation of a persistent cold-gas disc perpendicular to the jets that is centrally confined. Cavity luminosities measured from our simulations generally reflect the cooling luminosities of the intracluster medium (ICM) and correspond to averaged jet powers that are relatively insensitive to short periods of low-luminosity jet injection. Cold gas structures in our MHD simulations with low momentum-density jets generally show a variety of morphologies ranging from discy to very extended filamentary structures. In particular, magnetic fields are crucial to inhibit the formation of unrealistically massive cold gas discs by redistributing angular momentum between the hot and cold phases and by fostering the formation of elongated cold filaments that are supported by magnetic pressure.
△ Less
Submitted 25 March, 2023; v1 submitted 4 April, 2022;
originally announced April 2022.
-
Turning AGN bubbles into radio relics with sloshing: modeling CR transport with realistic physics
Authors:
John ZuHone,
Kristian Ehlert,
Rainer Weinberger,
Christoph Pfrommer
Abstract:
Radio relics are arc-like synchrotron sources at the periphery of galaxy clusters, produced by cosmic-ray electrons in a $μ$G magnetic field which are believed to have been (re-)accelerated by merger shock fronts. However, not all relics appear at the same location as shocks as seen in the X-ray. In a previous work, we suggested that the shape of some relics may result from the pre-existing spatia…
▽ More
Radio relics are arc-like synchrotron sources at the periphery of galaxy clusters, produced by cosmic-ray electrons in a $μ$G magnetic field which are believed to have been (re-)accelerated by merger shock fronts. However, not all relics appear at the same location as shocks as seen in the X-ray. In a previous work, we suggested that the shape of some relics may result from the pre-existing spatial distribution of cosmic-ray electrons, and tested this hypothesis using simulations by launching AGN jets into a cluster atmosphere with sloshing gas motions generated by a previous merger event. We showed that these motions could transport the cosmic ray-enriched material of the AGN bubbles to large radii and stretch it in a tangential direction, producing a filamentary shape resembling a radio relic. In this work, we improve our physical description for the cosmic rays by modeling them as a separate fluid which undergoes diffusion and Alfvén losses. We find that including this additional cosmic ray physics significantly diminishes the appearance of these filamentary features, showing that our original hypothesis is sensitive to the modeling of cosmic ray physics in the intracluster medium.
△ Less
Submitted 29 October, 2021; v1 submitted 8 October, 2021;
originally announced October 2021.
-
How merger-driven gas motions in galaxy clusters can turn AGN bubbles into radio relics
Authors:
John A. ZuHone,
Maxim Markevitch,
Rainer Weinberger,
Paul Nulsen,
Kristian Ehlert
Abstract:
Radio relics in galaxy clusters are extended synchrotron sources produced by cosmic-ray electrons in the $μ$G magnetic field. Many relics are found in the cluster periphery and have a cluster-centric, narrow arc-like shape, which suggests that the electrons are accelerated or re-accelerated by merger shock fronts propagating outward in the intracluster plasma. In the X-ray, some relics do exhibit…
▽ More
Radio relics in galaxy clusters are extended synchrotron sources produced by cosmic-ray electrons in the $μ$G magnetic field. Many relics are found in the cluster periphery and have a cluster-centric, narrow arc-like shape, which suggests that the electrons are accelerated or re-accelerated by merger shock fronts propagating outward in the intracluster plasma. In the X-ray, some relics do exhibit such shocks at the location of the relic, but many do not. We explore the possibility that radio relics trace not the shock fronts but the shape of the underlying distribution of seed relativistic electrons, lit up by a recent shock passage. We use magnetohydrodynamic simulations of cluster mergers and include bubbles of relativistic electrons injected by jets from the central AGN or from an off-center radio galaxy. We show that the merger-driven gas motions (a) can advect the bubble cosmic rays to very large radii, and (b) spread the relativistic seed electrons preferentially in tangential direction -- along the gravitational equipotential surfaces -- producing extended, filamentary or sheet-like regions of intracluster plasma enriched with aged cosmic rays, which resemble radio relics. Once a shock front passes across such a region, the sharp radio emission edges would trace the sharp boundaries of these enriched regions rather than the front. We also show that these elongated cosmic ray features are naturally associated with magnetic fields stretched tangentially along their long axis, which could help explain the high polarization of relics.
△ Less
Submitted 27 May, 2021; v1 submitted 3 December, 2020;
originally announced December 2020.
-
Connecting turbulent velocities and magnetic fields in galaxy cluster simulations with active galactic nuclei jets
Authors:
Kristian Ehlert,
Rainer Weinberger,
Christoph Pfrommer,
Volker Springel
Abstract:
The study of velocity fields of the hot gas in galaxy clusters can help to unravel details of microphysics on small-scales and to decipher the nature of feedback by active galactic nuclei (AGN). Likewise, magnetic fields as traced by Faraday rotation measurements (RMs) inform about their impact on gas dynamics as well as on cosmic ray production and transport. We investigate the inherent relations…
▽ More
The study of velocity fields of the hot gas in galaxy clusters can help to unravel details of microphysics on small-scales and to decipher the nature of feedback by active galactic nuclei (AGN). Likewise, magnetic fields as traced by Faraday rotation measurements (RMs) inform about their impact on gas dynamics as well as on cosmic ray production and transport. We investigate the inherent relationship between large-scale gas kinematics and magnetic fields through non-radiative magnetohydrodynamical simulations of the creation, evolution and disruption of AGN jet-inflated lobes in an isolated Perseus-like galaxy cluster, with and without pre-existing turbulence. In particular, we connect cluster velocity measurements with mock RM maps to highlight their underlying physical connection, which opens up the possibility of comparing turbulence levels in two different observables. For single jet outbursts, we find only a local impact on the velocity field, i.e. the associated increase in velocity dispersion is not volume-filling. Furthermore, in a setup with pre-existing turbulence, this increase in velocity dispersion is largely hidden. We use mock X-ray observations to show that at arcmin resolution, the velocity dispersion is therefore dominated by existing large-scale turbulence and is only minimally altered by the presence of a jet. For the velocity structure of central gas uplifted by buoyantly rising lobes, we find fast, coherent outflows with low velocity dispersion. Our results highlight that projected velocity distributions show complex structures which pose challenges for the interpretation of observations.
△ Less
Submitted 27 November, 2020;
originally announced November 2020.
-
Interaction of a cold cloud with a hot wind: the regimes of cloud growth and destruction and the impact of magnetic fields
Authors:
Martin Sparre,
Christoph Pfrommer,
Kristian Ehlert
Abstract:
Multiphase galaxy winds, the accretion of cold gas through galaxy haloes, and gas stripping from jellyfish galaxies are examples of interactions between cold and hot gaseous phases. There are two important regimes in such systems. A sufficiently small cold cloud is destroyed by the hot wind as a result of Kelvin-Helmholtz instabilities, which shatter the cloud into small pieces that eventually mix…
▽ More
Multiphase galaxy winds, the accretion of cold gas through galaxy haloes, and gas stripping from jellyfish galaxies are examples of interactions between cold and hot gaseous phases. There are two important regimes in such systems. A sufficiently small cold cloud is destroyed by the hot wind as a result of Kelvin-Helmholtz instabilities, which shatter the cloud into small pieces that eventually mix and dissolve in the hot wind. On the contrary, stripped cold gas from a large cloud mixes with the hot wind to intermediate temperatures, and then becomes thermally unstable and cools, causing a net accretion of hot gas to the cold tail. Using the magneto-hydrodynamical code AREPO, we perform cloud crushing simulations and test analytical criteria for the transition between the growth and destruction regimes to clarify a current debate in the literature. We find that the hot-wind cooling time sets the transition radius and not the cooling time of the mixed phase. Magnetic fields modify the wind-cloud interaction. Draping of wind magnetic field enhances the field upstream of the cloud and fluid instabilities are suppressed by a turbulently magnetised wind beyond what is seen for a wind with a uniform magnetic field. We furthermore predict jellyfish galaxies to have ordered magnetic fields aligned with their tails. We finally discuss how the results of idealised simulations can be used to provide input to subgrid models in cosmological (magneto-)hydrodynamical simulations, which cannot resolve the detailed small-scale structure of cold gas clouds in the circum-galactic medium.
△ Less
Submitted 12 October, 2020; v1 submitted 20 August, 2020;
originally announced August 2020.
-
The Sunyaev-Zel'dovich effect of simulated jet-inflated bubbles in clusters
Authors:
Kristian Ehlert,
Christoph Pfrommer,
Rainer Weinberger,
Rüdiger Pakmor,
Volker Springel
Abstract:
Feedback by active galactic nuclei (AGNs) is essential for regulating the fast radiative cooling of low-entropy gas at the centers of galaxy clusters and for reducing star formation rates of central ellipticals. The details of self-regulation depend critically on the unknown contents of AGN-inflated bubbles. Observations of the Sunyaev-Zeldovich (SZ) signal of AGN bubbles provide us with the abili…
▽ More
Feedback by active galactic nuclei (AGNs) is essential for regulating the fast radiative cooling of low-entropy gas at the centers of galaxy clusters and for reducing star formation rates of central ellipticals. The details of self-regulation depend critically on the unknown contents of AGN-inflated bubbles. Observations of the Sunyaev-Zeldovich (SZ) signal of AGN bubbles provide us with the ability to directly measure the lobe electron pressure given a bubble morphology. Here we compute the SZ signal of jet-inflated bubbles in three-dimensional magnetohydrodynamical simulations of the galaxy cluster MS0735.6+7421 with the Arepo code, and compare our synthetic SZ results to inferences obtained with popular modelling approaches. We find that cutting out ellipsoidal bubbles from a double-beta pressure profile only matches the inner bubble edges in the simulations and fails to account for the emission of the shock-enhanced pressure cocoon outside the bubbles. This additional contribution significantly worsens the accuracy of the cut-out method for jets with small inclinations with respect to the line of sight. Also, the kinetic SZ effect of the bubbles, a previously neglected contribution, becomes relevant at these smaller inclinations due to entrainment and mixing of the intracluster medium with low-density jet material. Fortunately, the different signs of the kinetic SZ signal in opposite lobes allow modelling this effect. We present an approximate method to determine the jet inclination, which combines jet power and lifetime estimates, the stand-off distance between jet head and bow shock, and the kinetic SZ effect, thereby helping to correctly infer the bubble contents.
△ Less
Submitted 17 December, 2018;
originally announced December 2018.
-
Simulations of the dynamics of magnetised jets and cosmic rays in galaxy clusters
Authors:
Kristian Ehlert,
Rainer Weinberger,
Christoph Pfrommer,
Rüdiger Pakmor,
Volker Springel
Abstract:
Feedback processes by active galactic nuclei in the centres of galaxy clusters appear to prevent large-scale cooling flows and impede star formation. However, the detailed heating mechanism remains uncertain. One promising heating scenario invokes the dissipation of Alfvén waves that are generated by streaming cosmic rays (CRs). In order to study this idea, we use three-dimensional magneto-hydrody…
▽ More
Feedback processes by active galactic nuclei in the centres of galaxy clusters appear to prevent large-scale cooling flows and impede star formation. However, the detailed heating mechanism remains uncertain. One promising heating scenario invokes the dissipation of Alfvén waves that are generated by streaming cosmic rays (CRs). In order to study this idea, we use three-dimensional magneto-hydrodynamical simulations with the AREPO code that follow the evolution of jet-inflated bubbles that are filled with CRs in a turbulent cluster atmosphere. We find that a single injection event produces the CR distribution and heating rate required for a successful CR heating model. As a bubble rises buoyantly, cluster magnetic fields drape around the leading interface and are amplified to strengths that balance the ram pressure. Together with helical magnetic fields in the bubble, this initially confines the CRs and suppresses the formation of interface instabilities. But as the bubble continues to rise, bubble-scale eddies significantly amplify radial magnetic filaments in its wake and enable CR transport from the bubble to the cooling intracluster medium. By varying the jet parameters, we obtain a rich and diverse set of jet and bubble morphologies ranging from Fanaroff-Riley type I-like (FRI) to FRII-like jets. We identify jet energy as the leading order parameter (keeping the ambient density profiles fixed), whereas jet luminosity is primarily responsible for setting the Mach numbers of shocks around FRII-like sources. Our simulations also produce FRI-like jets that inflate bubbles without detectable shocks and show morphologies consistent with cluster observations.
△ Less
Submitted 27 November, 2020; v1 submitted 14 June, 2018;
originally announced June 2018.
-
Constraining the coherence scale of the interstellar magnetic field using TeV gamma-ray observations of supernova remnants
Authors:
Matteo Pais,
Christoph Pfrommer,
Kristian Ehlert,
Maria Werhahn,
Georg Winner
Abstract:
Galactic cosmic rays (CRs) are believed to be accelerated at supernova remnant (SNR) shocks. In the hadronic scenario the TeV gamma-ray emission from SNRs originates from decaying pions that are produced in collisions of the interstellar gas and CRs. Using CR-magnetohydrodynamic simulations, we show that magnetic obliquity-dependent shock acceleration is able to reproduce the observed TeV gamma-ra…
▽ More
Galactic cosmic rays (CRs) are believed to be accelerated at supernova remnant (SNR) shocks. In the hadronic scenario the TeV gamma-ray emission from SNRs originates from decaying pions that are produced in collisions of the interstellar gas and CRs. Using CR-magnetohydrodynamic simulations, we show that magnetic obliquity-dependent shock acceleration is able to reproduce the observed TeV gamma-ray morphology of SNRs such as Vela Jr. and SN1006 solely by varying the magnetic morphology. This implies that gamma-ray bright regions result from quasi-parallel shocks (i.e., when the shock propagates at a narrow angle to the upstream magnetic field), which are known to efficiently accelerate CR protons, and that gamma-ray dark regions point to quasi-perpendicular shock configurations. Comparison of the simulated gamma-ray morphology to observations allows us to constrain the magnetic coherence scale $λ_B$ around Vela Jr. and SN1006 to $λ_B \simeq 13_{-4.3}^{+13}$ pc and $λ_B >200_{-10}^{+80}$ pc, respectively, where the ambient magnetic field of SN1006 is consistent with being largely homogeneous. We find consistent pure hadronic and mixed hadronic-leptonic models that both reproduce the multi-frequency spectra from the radio to TeV gamma rays and match the observed gamma-ray morphology. Finally, to capture the propagation of a SNR shock in a clumpy interstellar medium, we study the interaction of a shock with a dense cloud with numerical simulations and analytics. We construct an analytical gamma-ray model for a core collapse SNR propagating through a structured interstellar medium, and show that the gamma-ray luminosity is only biased by 30% for realistic parameters.
△ Less
Submitted 9 June, 2020; v1 submitted 8 May, 2018;
originally announced May 2018.
-
The effect of cosmic-ray acceleration on supernova blast wave dynamics
Authors:
M. Pais,
C. Pfrommer,
K. Ehlert,
R. Pakmor
Abstract:
Non-relativistic shocks accelerate ions to highly relativistic energies provided that the orientation of the magnetic field is closely aligned with the shock normal (quasi-parallel shock configuration). In contrast, quasi-perpendicular shocks do not efficiently accelerate ions. We model this obliquity-dependent acceleration process in a spherically expanding blast wave setup with the moving-mesh c…
▽ More
Non-relativistic shocks accelerate ions to highly relativistic energies provided that the orientation of the magnetic field is closely aligned with the shock normal (quasi-parallel shock configuration). In contrast, quasi-perpendicular shocks do not efficiently accelerate ions. We model this obliquity-dependent acceleration process in a spherically expanding blast wave setup with the moving-mesh code {\sc arepo} for different magnetic field morphologies, ranging from homogeneous to turbulent configurations. A Sedov-Taylor explosion in a homogeneous magnetic field generates an oblate ellipsoidal shock surface due to the slower propagating blast wave in the direction of the magnetic field. This is because of the efficient cosmic ray (CR) production in the quasi-parallel polar cap regions, which softens the equation of state and increases the compressibility of the post-shock gas. We find that the solution remains self-similar because the ellipticity of the propagating blast wave stays constant in time. This enables us to derive an effective ratio of specific heats for a composite of thermal gas and CRs as a function of the maximum acceleration efficiency. We finally discuss the behavior of supernova remnants expanding into a turbulent magnetic field with varying coherence lengths. For a maximum CR acceleration efficiency of about 15 per cent at quasi-parallel shocks (as suggested by kinetic plasma simulations), we find an average efficiency of about 5 per cent, independent of the assumed magnetic coherence length.
△ Less
Submitted 30 April, 2018;
originally announced May 2018.
-
Simulating the interaction of jets with the intra-cluster medium
Authors:
Rainer Weinberger,
Kristian Ehlert,
Christoph Pfrommer,
Rüdiger Pakmor,
Volker Springel
Abstract:
Jets from supermassive black holes in the centres of galaxy clusters are a potential candidate for moderating gas cooling and subsequent star formation through depositing energy in the intra-cluster gas. In this work, we simulate the jet-intra-cluster medium interaction using the moving-mesh magnetohydrodynamics code Arepo. Our model injects supersonic, low density, collimated and magnetised outfl…
▽ More
Jets from supermassive black holes in the centres of galaxy clusters are a potential candidate for moderating gas cooling and subsequent star formation through depositing energy in the intra-cluster gas. In this work, we simulate the jet-intra-cluster medium interaction using the moving-mesh magnetohydrodynamics code Arepo. Our model injects supersonic, low density, collimated and magnetised outflows in cluster centres, which are then stopped by the surrounding gas, thermalise and inflate low-density cavities filled with cosmic-rays. We perform high-resolution, non-radiative simulations of the lobe creation, expansion and disruption, and find that its dynamical evolution is in qualitative agreement with simulations of idealised low-density cavities that are dominated by a large-scale Rayleigh-Taylor instability. The buoyant rising of the lobe does not create energetically significant small-scale chaotic motion in a volume-filling fashion, but rather a systematic upward motion in the wake of the lobe and a corresponding back-flow perpendicular to it. We find that, overall, 50 per cent of the injected energy ends up in material which is not part of the lobe, and about 25 per cent remains in the inner 100 kpc. We conclude that jet-inflated, buoyantly rising cavities drive systematic gas motions which play an important role in heating the central regions, while mixing of lobe material is sub-dominant. Encouragingly, the main mechanisms responsible for this energy deposition can be modelled already at resolutions within reach in future, high-resolution cosmological simulations of galaxy clusters.
△ Less
Submitted 27 March, 2017;
originally announced March 2017.