Showing 1–50 of 77 results for author: Ruszkowski, M

On the Origin of the Ancient, Large-Scale Cold Front in the Perseus Cluster of Galaxies

Authors: Elena Bellomi, John ZuHone, Rainer Weinberger, Stephen Walker, Irina Zhuravleva, Mateusz Ruszkowski, Maxim Markevitch

Abstract: The intracluster medium of the Perseus Cluster exhibits spiral-shaped X-ray surface brightness discontinuities known as ``cold fronts'', which simulations indicate are caused by the sloshing motion of the gas after the passage of a subcluster. Recent observations of Perseus have shown that these fronts extend to large radii. In this work, we present simulations of the formation of sloshing cold fr… ▽ More The intracluster medium of the Perseus Cluster exhibits spiral-shaped X-ray surface brightness discontinuities known as ``cold fronts'', which simulations indicate are caused by the sloshing motion of the gas after the passage of a subcluster. Recent observations of Perseus have shown that these fronts extend to large radii. In this work, we present simulations of the formation of sloshing cold fronts in Perseus using the AREPO magnetohydrodynamics code, to produce a plausible scenario for the formation of the large front at a radius of 700 kpc. Our simulations explore a range of subcluster masses and impact parameters. We find that low-mass subclusters cannot generate a cold front that can propagate to such a large radius, and that small impact parameters create too much turbulence which leads to the disruption of the cold front before it reaches such a large distance. Subclusters which make only one core passage produce a stable initial front that expands to large radii, but without a second core passage of the subcluster other fronts are not created at a later time in the core region. We find a small range of simulations with subclusters with mass ratios of $R \sim 1:5$ and initial impact parameter of $θ~\sim~20-25^\circ$ which not only produce the large cold front but a second set in the core region at later times. These simulations indicate that the ``ancient' cold front is $\sim$6-8.5 Gyr old. For the simulations providing the closest match with observations, the subcluster has completely merged into the main cluster. △ Less

Submitted 13 October, 2023; originally announced October 2023.

Comments: 24 pages, 16 figures, submitted to ApJ

Shock-Driven Periodic Variability in a Low-Mass-Ratio Supermassive Black Hole Binary

Authors: K. Whitley, A. Kuznetsova, K. Gültekin, M. Ruszkowski

Abstract: We investigate the time-varying electromagnetic emission of a low-mass-ratio supermassive black hole binary (SMBHB) embedded in a circumprimary disk, with a particular interest in variability of shocks driven by the binary. We perform a 2D, locally isothermal hydrodynamics simulation of a SMBHB with mass ratio $q=0.01$ and separation $a=100\;R_g$, using a physically self-consistent steady disk mod… ▽ More We investigate the time-varying electromagnetic emission of a low-mass-ratio supermassive black hole binary (SMBHB) embedded in a circumprimary disk, with a particular interest in variability of shocks driven by the binary. We perform a 2D, locally isothermal hydrodynamics simulation of a SMBHB with mass ratio $q=0.01$ and separation $a=100\;R_g$, using a physically self-consistent steady disk model. We estimate the electromagnetic variability from the system by monitoring accretion onto the secondary and using an artificial viscosity scheme to capture shocks and monitor the energy dissipated. The SMBHB produces a wide, eccentric gap in the disk, previously only observed for larger mass ratios, which we attribute to our disk model being much thinner ($H/R\approx0.01$ near the secondary) than is typical of previous works. The eccentric gap drives periodic accretion onto the secondary SMBH on a timescale matching the orbital period of the binary, $t_{\rm{bin}}\approx0.1\;\rm{yr}$, implying that the variable accretion regime of the SMBHB parameter space extends to lower mass ratios than previously established. Shocks driven by the binary are periodic, with a period matching the orbital period, and the shocks are correlated with the accretion rate, with peaks in the shock luminosity lagging peaks in the accretion rate by $0.43\;t_{\rm{bin}}$. We propose that the correlation of these quantities represents a useful identifier of SMBHB candidates, via observations of correlated variability in X-ray and UV monitoring of AGN, rather than single-waveband periodicity alone. △ Less

Submitted 18 October, 2023; v1 submitted 10 October, 2023; originally announced October 2023.

Comments: 12 pages, 8 figures, accepted by MNRAS

Cosmic ray feedback in galaxies and galaxy clusters -- A pedagogical introduction and a topical review of the acceleration, transport, observables, and dynamical impact of cosmic rays

Authors: Mateusz Ruszkowski, Christoph Pfrommer

Abstract: Understanding the physical mechanisms that control galaxy formation is a fundamental challenge in contemporary astrophysics. Recent advances in the field of astrophysical feedback strongly suggest that cosmic rays (CRs) may be crucially important for our understanding of cosmological galaxy formation and evolution. The appealing features of CRs are their relatively long cooling times and relativel… ▽ More Understanding the physical mechanisms that control galaxy formation is a fundamental challenge in contemporary astrophysics. Recent advances in the field of astrophysical feedback strongly suggest that cosmic rays (CRs) may be crucially important for our understanding of cosmological galaxy formation and evolution. The appealing features of CRs are their relatively long cooling times and relatively strong dynamical coupling to the gas. In galaxies, CRs can be close to equipartition with the thermal, magnetic, and turbulent energy density in the interstellar medium, and can be dynamically very important in driving large-scale galactic winds. Similarly, CRs may provide a significant contribution to the pressure in the circumgalactic medium. In galaxy clusters, CRs may play a key role in addressing the classic cooling flow problem by facilitating efficient heating of the intracluster medium and preventing excessive star formation. Overall, the underlying physics of CR interactions with plasmas exhibit broad parallels across the entire range of scales characteristic of the interstellar, circumgalactic, and intracluster media. Here we present a review of the state-of-the-art of this field and provide a pedagogical introduction to cosmic ray plasma physics, including the physics of wave-particle interactions, acceleration processes, CR spatial and spectral transport, and important cooling processes. The field is ripe for discovery and will remain the subject of intense theoretical, computational, and observational research over the next decade with profound implications for the interpretation of the observations of stellar and supermassive black hole feedback spanning the entire width of the electromagnetic spectrum and multi-messenger data. △ Less

Submitted 6 October, 2023; v1 submitted 5 June, 2023; originally announced June 2023.

Comments: invited A&ARv review; revised version; accepted for publication; 238 pages

Line Emission Mapper (LEM): Probing the physics of cosmic ecosystems

Authors: Ralph Kraft, Maxim Markevitch, Caroline Kilbourne, Joseph S. Adams, Hiroki Akamatsu, Mohammadreza Ayromlou, Simon R. Bandler, Marco Barbera, Douglas A. Bennett, Anil Bhardwaj, Veronica Biffi, Dennis Bodewits, Akos Bogdan, Massimiliano Bonamente, Stefano Borgani, Graziella Branduardi-Raymont, Joel N. Bregman, Joseph N. Burchett, Jenna Cann, Jenny Carter, Priyanka Chakraborty, Eugene Churazov, Robert A. Crain, Renata Cumbee, Romeel Dave , et al. (85 additional authors not shown)

Abstract: The Line Emission Mapper (LEM) is an X-ray Probe for the 2030s that will answer the outstanding questions of the Universe's structure formation. It will also provide transformative new observing capabilities for every area of astrophysics, and to heliophysics and planetary physics as well. LEM's main goal is a comprehensive look at the physics of galaxy formation, including stellar and black-hole… ▽ More The Line Emission Mapper (LEM) is an X-ray Probe for the 2030s that will answer the outstanding questions of the Universe's structure formation. It will also provide transformative new observing capabilities for every area of astrophysics, and to heliophysics and planetary physics as well. LEM's main goal is a comprehensive look at the physics of galaxy formation, including stellar and black-hole feedback and flows of baryonic matter into and out of galaxies. These processes are best studied in X-rays, and emission-line mapping is the pressing need in this area. LEM will use a large microcalorimeter array/IFU, covering a 30x30' field with 10" angular resolution, to map the soft X-ray line emission from objects that constitute galactic ecosystems. These include supernova remnants, star-forming regions, superbubbles, galactic outflows (such as the Fermi/eROSITA bubbles in the Milky Way and their analogs in other galaxies), the Circumgalactic Medium in the Milky Way and other galaxies, and the Intergalactic Medium at the outskirts and beyond the confines of galaxies and clusters. LEM's 1-2 eV spectral resolution in the 0.2-2 keV band will make it possible to disentangle the faintest emission lines in those objects from the bright Milky Way foreground, providing groundbreaking measurements of the physics of these plasmas, from temperatures, densities, chemical composition to gas dynamics. While LEM's main focus is on galaxy formation, it will provide transformative capability for all classes of astrophysical objects, from the Earth's magnetosphere, planets and comets to the interstellar medium and X-ray binaries in nearby galaxies, AGN, and cooling gas in galaxy clusters. In addition to pointed observations, LEM will perform a shallow all-sky survey that will dramatically expand the discovery space. △ Less

Submitted 12 April, 2023; v1 submitted 17 November, 2022; originally announced November 2022.

Comments: 18 pages. White paper for a mission concept to be submitted for the 2023 NASA Astrophysics Probes opportunity. v2: All-sky survey figure expanded, references fixed. v3: Added energy resolution measurements for prototype detector array. v4: Author list and reference fixes

Turbulent Heating in a Stratified Medium

Authors: Chaoran Wang, S. Peng Oh, M. Ruszkowski

Abstract: There is considerable evidence for widespread subsonic turbulence in galaxy clusters, most notably from {\it Hitomi}. Turbulence is often invoked to offset radiative losses in cluster cores, both by direct dissipation and by enabling turbulent heat diffusion. However, in a stratified medium, buoyancy forces oppose radial motions, making turbulence anisotropic. This can be quantified via the Froude… ▽ More There is considerable evidence for widespread subsonic turbulence in galaxy clusters, most notably from {\it Hitomi}. Turbulence is often invoked to offset radiative losses in cluster cores, both by direct dissipation and by enabling turbulent heat diffusion. However, in a stratified medium, buoyancy forces oppose radial motions, making turbulence anisotropic. This can be quantified via the Froude number ${\rm Fr}$, which decreases inward in clusters as stratification increases. We exploit analogies with MHD turbulence to show that wave-turbulence interactions increase cascade times and reduces dissipation rates $ε\propto {\rm Fr}$. Equivalently, for a given energy injection/dissipation rate $ε$, turbulent velocities $u$ must be higher compared to Kolmogorov scalings. High resolution hydrodynamic simulations show excellent agreement with the $ε\propto {\rm Fr}$ scaling, which sets in for ${\rm Fr} < 0.1$. We also compare previously predicted scalings for the turbulent diffusion coefficient $D \propto {\rm Fr}^2$ and find excellent agreement, for ${\rm Fr} < 1$. However, we find a different normalization, corresponding to stronger diffusive suppression by more than an order of magnitude. Our results imply that turbulent diffusion is more heavily suppressed by stratification, over a much wider radial range, than turbulent dissipation. Thus, the latter potentially dominates. Furthermore, this shift implies significantly higher turbulent velocities required to offset cooling, compared to previous models. These results are potentially relevant to turbulent metal diffusion (which is likewise suppressed), and to planetary atmospheres. △ Less

Submitted 3 May, 2022; originally announced May 2022.

Fermi and eRosita bubbles as relics of the past activity of the Galactic black hole

Authors: H. -Y. Karen Yang, M. Ruszkowski, E. Zweibel

Abstract: The newly launched X-ray satellite, eRosita, has recently revealed two gigantic bubbles extending to ~80 degrees above and below the Galactic center. The morphology of these "eRosita bubbles" bears a remarkable resemblance to the Fermi bubbles previously discovered by the Fermi Gamma-ray Space Telescope and its counterpart, the microwave haze. The physical origin of these striking structures has b… ▽ More The newly launched X-ray satellite, eRosita, has recently revealed two gigantic bubbles extending to ~80 degrees above and below the Galactic center. The morphology of these "eRosita bubbles" bears a remarkable resemblance to the Fermi bubbles previously discovered by the Fermi Gamma-ray Space Telescope and its counterpart, the microwave haze. The physical origin of these striking structures has been intensely debated; however, because of their symmetry about the Galactic center, they likely originate from some energetic outbursts from the Galactic center in the past. Here we propose a theoretical model in which the eRosita bubbles, Fermi bubbles, and the microwave haze could be simultaneously explained by a single event of jet activity from the central supermassive black hole a few million years ago. Using numerical simulations, we show that this model could successfully reproduce the morphology and multi-wavelength spectra of the observed bubbles and haze, which allows us to derive critical constraints on the energetics and timescales of the outburst. This study serves as an important step forward in our understanding of the past Galactic center activity of our Milky Way Galaxy, and may bring valuable insights into the broader picture of supermassive black hole-galaxy co-evolution in the context of galaxy formation. △ Less

Submitted 4 March, 2022; originally announced March 2022.

Comments: 23 pages, 5 figures, matched to the published version in Nature Astronomy

Stress-Testing Cosmic Ray Physics: The Impact of Cosmic Rays on the Surviving Disk of Ram Pressure Stripped Galaxies

Authors: Ryan Jeffrey Farber, Mateusz Ruszkowski, Stephanie Tonnesen, Paco Holguin

Abstract: Cluster spiral galaxies suffer catastrophic losses of the cool, neutral gas component of their interstellar medium due to ram pressure stripping, contributing to the observed quenching of star formation in the disk compared to galaxies in lower density environments. However, the short term effects of ram pressure on the star formation rate and AGN activity of galaxies undergoing stripping remain u… ▽ More Cluster spiral galaxies suffer catastrophic losses of the cool, neutral gas component of their interstellar medium due to ram pressure stripping, contributing to the observed quenching of star formation in the disk compared to galaxies in lower density environments. However, the short term effects of ram pressure on the star formation rate and AGN activity of galaxies undergoing stripping remain unclear. Numerical studies have recently demonstrated cosmic rays can dramatically influence galaxy evolution for isolated galaxies, yet their influence on ram pressure stripping remains poorly constrained. We perform the first cosmic-ray magneto-hydrodynamic simulations of an $L_{*}$ galaxy undergoing ram pressure stripping, including radiative cooling, self-gravity of the gas, star formation, and stellar feedback. We find the microscopic transport of cosmic rays plays a key role in modulating the star formation enhancement experienced by spirals at the outskirts of clusters compared to isolated spirals. Moreover, we find that galaxies undergoing ram pressure stripping exhibit enhanced gas accretion onto their centers, which may explain the prevalence of AGN in these objects. In agreement with observations, we find cosmic rays significantly boost the global radio emission of cluster spirals. Although the gas removal rate is relatively insensitive to cosmic ray physics, we find that cosmic rays significantly modify the phase distribution of the remaining gas disk. These results suggest observations of galaxies undergoing ram pressure stripping may place novel constraints on cosmic-ray calorimetry and transport. △ Less

Submitted 11 January, 2022; originally announced January 2022.

Comments: 15 pages, 10 figures, submitted to MNRAS

Non-Kolmogorov turbulence in multiphase intracluster medium driven by cold gas precipitation and AGN jets

Authors: C. Wang, M. Ruszkowski, C. Pfrommer, S. Peng Oh, H. -Y. K. Yang

Abstract: AGN feedback is responsible for maintaining plasma in global thermal balance in extended halos of elliptical galaxies and galaxy clusters. Local thermal instability in the hot gas leads to the formation of precipitating cold gas clouds that feed the central supermassive black holes, thus heating the hot gas and maintaining global thermal equilibrium. We perform three dimensional MHD simulations of… ▽ More AGN feedback is responsible for maintaining plasma in global thermal balance in extended halos of elliptical galaxies and galaxy clusters. Local thermal instability in the hot gas leads to the formation of precipitating cold gas clouds that feed the central supermassive black holes, thus heating the hot gas and maintaining global thermal equilibrium. We perform three dimensional MHD simulations of self-regulated AGN feedback in a Perseus-like galaxy cluster with the aim of understanding the impact of the feedback physics on the turbulence properties of the hot and cold phases of the ICM. We find that, in general, the cold phase velocity structure function (VSF) is steeper than the prediction from Kolmogorov's theory. We attribute the physical origin of the steeper slope of the cold phase VSF to the driving of turbulent motions primarily by the gravitational acceleration acting on the ballistic clouds. We demonstrate that, in the pure hydrodynamical case, the precipitating cold filaments may be the dominant agent driving turbulence in the hot ICM. The arguments in favor of this hypothesis are that: (i) the cold phase mass dominates over hot gas mass in the inner cool core; (ii) hot and cold gas velocities are spatially correlated; (iii) both the cold and hot phase velocity distributions are radially biased. We show that, in the MHD case, the turbulence in the ambient hot medium (excluding the jet cone regions) can also be driven by the AGN jets. The driving is then facilitated by enhanced coupling due to magnetic fields of the ambient gas and the AGN jets. In the MHD case, turbulence may thus be driven by a combination of AGN jet stirring and filament motions. We conclude that future observations, including those from high spatial and spectral resolution X-ray missions, may help to constrain self-regulated AGN feedback by quantifying the multi-temperature VSF in the ICM. △ Less

Submitted 20 December, 2020; originally announced December 2020.

Comments: 12pages, 15figures, submitted to MNRAS

Whistler-regulated MHD: Transport equations for electron thermal conduction in the high $β$ intracluster medium of galaxy clusters

Authors: J. F. Drake, C. Pfrommer, C. S. Reynolds, M. Ruszkowski, M. Swisdak, A. Einarsson, T. Thomas, A. B. Hassam, G. T. Roberg-Clark

Abstract: Transport equations for electron thermal energy in the high $β_e$ intracluster medium (ICM) are developed that include scattering from both classical collisions and self-generated whistler waves. The calculation employs an expansion of the kinetic electron equation along the ambient magnetic field in the limit of strong scattering and assumes whistler waves with low phase speeds… ▽ More Transport equations for electron thermal energy in the high $β_e$ intracluster medium (ICM) are developed that include scattering from both classical collisions and self-generated whistler waves. The calculation employs an expansion of the kinetic electron equation along the ambient magnetic field in the limit of strong scattering and assumes whistler waves with low phase speeds $V_w\sim{v}_{te}/β_e\ll{v}_{te}$ dominate the turbulent spectrum, with $v_{te}$ the electron thermal speed and $β_e\gg1$ the ratio of electron thermal to magnetic pressure. We find: (1) temperature-gradient-driven whistlers dominate classical scattering when $L_c>L/β_e$, with $L_c$ the classical electron mean-free-path and $L$ the electron temperature scale length, and (2) in the whistler dominated regime the electron thermal flux is controlled by both advection at $V_w$ and a comparable diffusive term. The findings suggest whistlers limit electron heat flux over large regions of the ICM, including locations unstable to isobaric condensation. Consequences include: (1) the Field length decreases, extending the domain of thermal instability to smaller length-scales, (2) the heat flux temperature dependence changes from $T_e^{7/2}/L$ to $V_wnT_e\sim{T}_e^{1/2}$, (3) the magneto-thermal and heat-flux driven buoyancy instabilities are impaired or completely inhibited, and (4) sound waves in the ICM propagate greater distances, as inferred from observations. This description of thermal transport can be used in macroscale ICM models. △ Less

Submitted 27 August, 2021; v1 submitted 15 July, 2020; originally announced July 2020.

Comments: Accepted for publication in The Astrophysical Journal

Chaotic cold accretion in giant elliptical galaxies heated by AGN cosmic rays

Authors: Chaoran Wang, Mateusz Ruszkowski, H. -Y. Karen Yang

Abstract: Black hole feedback plays a central role in shaping the circumgalactic medium (CGM) of elliptical galaxies. We systematically study the impact of plasma physics on the evolution of ellipticals by performing three-dimensional non-ideal magneto-hydrodynamic simulations of the interactions of active galactic nucleus (AGN) jets with the CGM including magnetic fields, and cosmic rays (CRs) and their tr… ▽ More Black hole feedback plays a central role in shaping the circumgalactic medium (CGM) of elliptical galaxies. We systematically study the impact of plasma physics on the evolution of ellipticals by performing three-dimensional non-ideal magneto-hydrodynamic simulations of the interactions of active galactic nucleus (AGN) jets with the CGM including magnetic fields, and cosmic rays (CRs) and their transport processes. We find that the physics of feedback operating on large galactic scales depends very sensitively on plasma physics operating on small scales. Specifically, we demonstrate that: (i) in the purely hydrodynamical case, the AGN jets initially maintain the atmospheres in global thermal balance. However, local thermal instability generically leads to the formation of massive cold disks in the vicinity of the central black hole in disagreement with observations; (ii) including weak magnetic fields prevents the formation of the disks because local B-field amplification in the precipitating cold gas leads to strong magnetic breaking, which quickly extracts angular momentum from the accreting clouds. The magnetic fields transform the cold clouds into narrow filaments that do not fall ballistically; (iii) when plasma composition in the AGN jets is dominated by CRs, and CR transport is neglected, the atmospheres exhibit cooling catastrophes due to inefficient heat transfer from the AGN to CGM despite Coulomb/hadronic CR losses being present; (iv) including CR streaming and heating restores agreement with the observations, i.e., cooling catastrophes are prevented and massive cold central disks do not form. The AGN power is reduced as its energy is utilized efficiently. △ Less

Submitted 8 October, 2019; originally announced October 2019.

Comments: submitted to MNRAS

Supermassive Black Hole Feedback

Authors: Mateusz Ruszkowski, Daisuke Nagai, Irina Zhuravleva, Corey Brummel-Smith, Yuan Li, Edmund Hodges-Kluck, Hsiang-Yi Karen Yang, Kaustuv Basu, Jens Chluba, Eugene Churazov, Megan Donahue, Andrew Fabian, Claude-André Faucher-Giguère, Massimo Gaspari, Julie Hlavacek-Larrondo, Michael McDonald, Brian McNamara, Paul Nulsen, Tony Mroczkowski, Richard Mushotzky, Christopher Reynolds, Alexey Vikhlinin, Mark Voit, Norbert Werner, John ZuHone , et al. (1 additional authors not shown)

Abstract: Understanding the processes that drive galaxy formation and shape the observed properties of galaxies is one of the most interesting and challenging frontier problems of modern astrophysics. We now know that the evolution of galaxies is critically shaped by the energy injection from accreting supermassive black holes (SMBHs). However, it is unclear how exactly the physics of this feedback process… ▽ More Understanding the processes that drive galaxy formation and shape the observed properties of galaxies is one of the most interesting and challenging frontier problems of modern astrophysics. We now know that the evolution of galaxies is critically shaped by the energy injection from accreting supermassive black holes (SMBHs). However, it is unclear how exactly the physics of this feedback process affects galaxy formation and evolution. In particular, a major challenge is unraveling how the energy released near the SMBHs is distributed over nine orders of magnitude in distance throughout galaxies and their immediate environments. The best place to study the impact of SMBH feedback is in the hot atmospheres of massive galaxies, groups, and galaxy clusters, which host the most massive black holes in the Universe, and where we can directly image the impact of black holes on their surroundings. We identify critical questions and potential measurements that will likely transform our understanding of the physics of SMBH feedback and how it shapes galaxies, through detailed measurements of (i) the thermodynamic and velocity fluctuations in the intracluster medium (ICM) as well as (ii) the composition of the bubbles inflated by SMBHs in the centers of galaxy clusters, and their influence on the cluster gas and galaxy growth, using the next generation of high spectral and spatial resolution X-ray and microwave telescopes. △ Less

Submitted 22 March, 2019; originally announced March 2019.

Comments: 10 pages, submitted to the Astro2020 decadal

Report number: https://ui.adsabs.harvard.edu/link_gateway/2019BAAS...51c.326R/PUB_PDF

Journal ref: 2019BAAS...51c.326R

Physics of cosmic plasmas from high angular resolution X-ray imaging of galaxy clusters

Authors: Maxim Markevitch, Esra Bulbul, Eugene Churazov, Simona Giacintucci, Ralph Kraft, Matthew Kunz, Daisuke Nagai, Elke Roediger, Mateusz Ruszkowski, Alex Schekochihin, Reinout van Weeren, Alexey Vikhlinin, Stephen A. Walker, Qian H. S. Wang, Norbert Werner, Daniel Wik, Irina Zhuravleva, John ZuHone

Abstract: Galaxy clusters are massive dark matter-dominated systems filled with X-ray emitting, optically thin plasma. Their large size and relative simplicity (at least as astrophysical objects go) make them a unique laboratory to measure some of the interesting plasma properties that are inaccessible by other means but fundamentally important for understanding and modeling many astrophysical phenomena --… ▽ More Galaxy clusters are massive dark matter-dominated systems filled with X-ray emitting, optically thin plasma. Their large size and relative simplicity (at least as astrophysical objects go) make them a unique laboratory to measure some of the interesting plasma properties that are inaccessible by other means but fundamentally important for understanding and modeling many astrophysical phenomena -- from solar flares to black hole accretion to galaxy formation and the emergence of the cosmological Large Scale Structure. While every cluster astrophysicist is eagerly anticipating the direct gas velocity measurements from the forthcoming microcalorimeters onboard XRISM, Athena and future missions such as Lynx, a number of those plasma properties can best be probed by high-resolution X-ray imaging of galaxy clusters. Chandra has obtained some trailblazing results, but only grazed the surface of such studies. In this white paper, we discuss why we need arcsecond-resolution, high collecting area, low relative background X-ray imagers (with modest spectral resolution), such as the proposed AXIS and the imaging detector of Lynx. △ Less

Submitted 15 March, 2019; originally announced March 2019.

Comments: Science white paper submitted for Astro2020 Decadal survey. 9 pages, 4 figures

Galaxy Winds in the Age of Hyperdimensional Astrophysics

Authors: Grant R. Tremblay, Evan E. Schneider, Alexey Vikhlinin, Lars Hernquist, Mateusz Ruszkowski, Benjamin D. Oppenheimer, Ralph P. Kraft, John ZuHone, Michael A. McDonald, Massimo Gaspari, Megan Donahue, G. Mark Voit

Abstract: The past decade began with the first light of ALMA and will end at the start of the new era of hyperdimensional astrophysics. Our community-wide movement toward highly multiwavelength and multidimensional datasets has enabled immense progress in each science frontier identified by the 2010 Decadal Survey, particularly with regard to black hole feedback and the cycle of baryons in galaxies. Facilit… ▽ More The past decade began with the first light of ALMA and will end at the start of the new era of hyperdimensional astrophysics. Our community-wide movement toward highly multiwavelength and multidimensional datasets has enabled immense progress in each science frontier identified by the 2010 Decadal Survey, particularly with regard to black hole feedback and the cycle of baryons in galaxies. Facilities like ALMA and the next generation of integral field unit (IFU) spectrographs together enable mapping the physical conditions and kinematics of warm ionized and cold molecular gas in galaxies in unprecedented detail (Fig. 1). JWST's launch at the start of the coming decade will push this capability to the rest-frame UV at redshifts z > 6, mapping the birth of stars in the first galaxies at cosmic dawn. Understanding of their subsequent evolution, however, now awaits an ability to map the processes that transform galaxies directly, rather than the consequences of those processes in isolation. In this paper, we argue that doing so requires an equivalent revolution in spatially resolved spectroscopy for the hot plasma that pervades galaxies, the atmospheres in which they reside, and the winds that are the engines of their evolution. △ Less

Submitted 13 March, 2019; originally announced March 2019.

Comments: A Science White Paper submitted to the Astro2020 Decadal Survey (7 papers, 4 figures)

"SZ spectroscopy" in the coming decade: Galaxy cluster cosmology and astrophysics in the submillimeter

Authors: Kaustuv Basu, Jens Erler, Jens Chluba, Jacques Delabrouille, J. Colin Hill, Tony Mroczkowski, Michael D. Niemack, Mathieu Remazeilles, Jack Sayers, Douglas Scott, Eve M. Vavagiakis, Michael Zemcov, Manuel Aravena, James G. Bartlett, Nicholas Battaglia, Frank Bertoldi, Maude Charmetant, Sunil Golwala, Terry L. Herter, Pamela Klaassen, Eiichiro Komatsu, Benjamin Magnelli, Adam B. Mantz, P. Daniel Meerburg, Jean-Baptiste Melin , et al. (8 additional authors not shown)

Abstract: Sunyaev-Zeldovich (SZ) effects were first proposed in the 1970s as tools to identify the X-ray emitting hot gas inside massive clusters of galaxies and obtain their velocities relative to the cosmic microwave background (CMB). Yet it is only within the last decade that they have begun to significantly impact astronomical research. Thanks to the rapid developments in CMB instrumentation, measuremen… ▽ More Sunyaev-Zeldovich (SZ) effects were first proposed in the 1970s as tools to identify the X-ray emitting hot gas inside massive clusters of galaxies and obtain their velocities relative to the cosmic microwave background (CMB). Yet it is only within the last decade that they have begun to significantly impact astronomical research. Thanks to the rapid developments in CMB instrumentation, measurement of the dominant thermal signature of the SZ effects has become a routine tool to find and characterize large samples of galaxy clusters and to seek deeper understanding of several important astrophysical processes via high-resolution imaging studies of many targets. With the notable exception of the Planck satellite and a few combinations of ground-based observatories, much of this "SZ revolution" has happened in the photometric mode, where observations are made at one or two frequencies in the millimeter regime to maximize the cluster detection significance and minimize the foregrounds. Still, there is much more to learn from detailed and systematic analyses of the SZ spectra across multiple wavelengths, specifically in the submillimeter (>300 GHz) domain. The goal of this Science White Paper is to highlight this particular aspect of SZ research, point out what new and potentially groundbreaking insights can be obtained from these studies, and emphasize why the coming decade can be a golden era for SZ spectral measurements. △ Less

Submitted 12 March, 2019; originally announced March 2019.

Comments: Science white paper submitted for the Astro2020 decadal review, 5 pages + references

Cosmic Ray Acceleration of Cool Clouds in the Circumgalactic Medium

Authors: Joshua Wiener, Ellen G. Zweibel, Mateusz Ruszkowski

Abstract: We investigate a mechanism for accelerating cool (10$^4$ K) clouds in the circumgalactic medium (CGM) with cosmic rays (CRs), possibly explaining some characteristics of observed high velocity clouds (HVCs). Enforcing CRs to stream down their pressure gradient into a region of slow streaming speed results in significant buildup of CR pressure which can accelerate the CGM. We present the results of… ▽ More We investigate a mechanism for accelerating cool (10$^4$ K) clouds in the circumgalactic medium (CGM) with cosmic rays (CRs), possibly explaining some characteristics of observed high velocity clouds (HVCs). Enforcing CRs to stream down their pressure gradient into a region of slow streaming speed results in significant buildup of CR pressure which can accelerate the CGM. We present the results of the first two-dimensional magnetohydrodynamic (MHD) simulations of such `CR bottlenecks,' expanding on simpler simulations in 1D from \cite{wiener17a}. Although much more investigation is required, we find two main results. First, radiative cooling in the interfaces of these clouds is sufficient to keep the cloud intact to CR wave heating. Second, cloud acceleration depends almost linearly with the injected CR flux at low values (comparable to that expected from a Milky Way-like star formation rate), but scales sublinearly at higher CR fluxes in 1D simulations. 2D simulations show hints of sublinear dependence at high CR fluxes but are consistent with pure linear dependence up to the CR fluxes tested. It is therefore plausible to accelerate cool clouds in the CGM to speeds of hundreds of km s$^{-1}$. △ Less

Submitted 4 March, 2019; originally announced March 2019.

Role of Cosmic Ray Streaming and Turbulent Damping in Driving Galactic Winds

Authors: F. Holguin, M. Ruszkowski, A. Lazarian, R. Farber, H. -Y. K. Yang

Abstract: Large-scale galactic winds driven by stellar feedback are one phenomenon that influences the dynamical and chemical evolution of a galaxy, redistributing material throughout the circumgalatic medium. Non-thermal feedback from galactic cosmic rays (CRs) -high-energy charged particles accelerated in supernovae and young stars - can impact the efficiency of wind driving. The streaming instability lim… ▽ More Large-scale galactic winds driven by stellar feedback are one phenomenon that influences the dynamical and chemical evolution of a galaxy, redistributing material throughout the circumgalatic medium. Non-thermal feedback from galactic cosmic rays (CRs) -high-energy charged particles accelerated in supernovae and young stars - can impact the efficiency of wind driving. The streaming instability limits the speed at which they can escape. However, in the presence of turbulence, the streaming instability is subject to suppression that depends on the magnetization of turbulence given by its Alfvén Mach number. While previous simulations that relied on a simplified model of CR transport have shown that super-Alfvénic streaming of CRs enhances galactic winds, in the present paper we take into account a realistic model of streaming suppression. We perform three-dimensional magnetohydrodynamic simulations of a section of a galactic disk and find that turbulent damping dependent on local magnetization of turbulent interstellar medium (ISM) leads to more spatially extended gas and CR distributions compared to the earlier streaming calculations, and that scale-heights of these distributions increase for stronger turbulence. Our results indicate that the star formation rate increases with the level of turbulence in the ISM. We also find that the instantaneous wind mass loading is sensitive to local streaming physics with the mass loading dropping significantly as the strength of turbulence increases. △ Less

Submitted 12 September, 2019; v1 submitted 15 July, 2018; originally announced July 2018.

Comments: 11 pages, 5 figures, accepted for publication in MNRAS

AGN Feedback and Multi-phase Gas in Giant Elliptical Galaxies

Authors: Chaoran Wang, Yuan Li, Mateusz Ruszkowski

Abstract: Recent observations have found extended multi-phase gas in a significant fraction of massive elliptical galaxies. We perform high-resolution three-dimensional hydrodynamical simulations of two idealized elliptical galaxies -- one representing a typical galaxy characterized by initial conditions conducive to the development of thermal instability and the other one less likely to develop thermal ins… ▽ More Recent observations have found extended multi-phase gas in a significant fraction of massive elliptical galaxies. We perform high-resolution three-dimensional hydrodynamical simulations of two idealized elliptical galaxies -- one representing a typical galaxy characterized by initial conditions conducive to the development of thermal instability and the other one less likely to develop thermal instability -- in order to study the development of thermal instability and the formation of multi-phase structures. We analyze the interplay between radiative cooling, momentum-driven AGN feedback, star formation, and stellar feedback from both young and old stars. We find that in one class of elliptical galaxies, the entropy of the hot halo gas rises sharply as a function of radius, and the hot halo is thermally stable and run-away cooling can only happen in the very center of the galaxy. In other class of ellipticals, the hot halo gas has a cooling to free-fall time ratio close to 10, and the non-linear perturbation driven by AGN feedback can cause the hot gas to frequently precipitate into extended multi-phase filaments. Both multi- and single-phase elliptical galaxies experience cooling-driven AGN feedback cycles. Interestingly, AGN feedback maintains the multi- or single-phase nature of the halo but does not turn multi-phase galaxies into single-phase ones or vice versa. Some of the extended cold gas in the multi-phase galaxy also forms young stars. The level of star formation and its spatial distribution are in excellent agreement with {\it Hubble} observations of nearby elliptical galaxies. △ Less

Submitted 8 May, 2018; originally announced May 2018.

Comments: submitted to MNRAS

Acoustic Disturbances in Galaxy Clusters

Authors: Ellen G. Zweibel, Vladimir V. Mirnov, Mateusz Ruszkowski, Christopher S. Reynolds, H. -Y. Karen Yang, Andrew C. Fabian

Abstract: Galaxy cluster cores are pervaded by hot gas which radiates at far too high a rate to maintain any semblance of a steady state; this is referred to as the cooling flow problem. Of the many heating mechanisms that have been proposed to balance radiative cooling, one of the most attractive is dissipation of acoustic waves generated by Active Galactic Nuclei (AGN). Fabian (2005) showed that if the wa… ▽ More Galaxy cluster cores are pervaded by hot gas which radiates at far too high a rate to maintain any semblance of a steady state; this is referred to as the cooling flow problem. Of the many heating mechanisms that have been proposed to balance radiative cooling, one of the most attractive is dissipation of acoustic waves generated by Active Galactic Nuclei (AGN). Fabian (2005) showed that if the waves are nearly adiabatic, wave damping due to heat conduction and viscosity must be well below standard Coulomb rates in order to allow the waves to propagate throughout the core. Because of the importance of this result, we have revisited wave dissipation under galaxy cluster conditions in a way that accounts for the self limiting nature of dissipation by electron thermal conduction, allows the electron and ion temperature perturbations in the waves to evolve separately, and estimates kinetic effects by comparing to a semi-collisionless theory. While these effects considerably enlarge the toolkit for analyzing observations of wavelike structures and developing a quantitative theory for wave heating, the drastic reduction of transport coefficients proposed in Fabian (2005) remains the most viable path to acoustic wave heating of galaxy cluster cores. △ Less

Submitted 27 March, 2018; v1 submitted 13 February, 2018; originally announced February 2018.

Comments: Accepted for publication in Astrophys. J

Powering of H$α$ filaments by cosmic rays

Authors: Mateusz Ruszkowski, H. -Y. Karen Yang, Christopher S. Reynolds

Abstract: Cluster cool cores possess networks of line-emitting filaments. These filaments are thought to originate via uplift of cold gas from cluster centers by buoyant active galactic nuclei (AGN) bubbles, or via local thermal instability in the hot intracluster medium (ICM). Therefore, the filaments are either the signatures of AGN feedback or feeding of supermassive black holes. Despite being characteri… ▽ More Cluster cool cores possess networks of line-emitting filaments. These filaments are thought to originate via uplift of cold gas from cluster centers by buoyant active galactic nuclei (AGN) bubbles, or via local thermal instability in the hot intracluster medium (ICM). Therefore, the filaments are either the signatures of AGN feedback or feeding of supermassive black holes. Despite being characterized by very short cooling times, the filaments are significant H$α$ emitters, which suggests that some process continuously powers these structures. Many cool cores host diffuse radio mini halos and AGN injecting radio plasma, suggesting that cosmic rays (CRs) and magnetic fields are present in the ICM. We argue that the excitation of Alfvén waves by CR streaming, and the replenishment of CR energy via accretion onto the filaments of high plasma-$β$ ICM characterized by low CR pressure support, can provide the adequate amount of heating to power and sustain the emission from these filaments. This mechanism does not require the CRs to penetrate the filaments even if the filaments are magnetically isolated from the ambient ICM and it may operate irrespectively of whether the filaments are dredged up from the center or form in situ in the ICM. This picture is qualitatively consistent with non-thermal line ratios seen in the cold filaments. Future X-ray observations of the iron line complex with XARM, Lynx, or Athena could help to test this model by providing constraints on the amount of CRs in the hot plasma that is cooling and accreting onto the filaments. △ Less

Submitted 27 March, 2018; v1 submitted 12 February, 2018; originally announced February 2018.

Comments: ApJ, accepted

Unveiling the Origin of the Fermi Bubbles

Authors: H. -Y. K. Yang, M. Ruszkowski, E. G. Zweibel

Abstract: The Fermi bubbles, two giant structures above and below the Galactic center (GC), are among the most important discoveries of the Fermi Gamma-ray Space Telescope. Studying their physical origin has been providing valuable insights into cosmic-ray transport, the Galactic magnetic field, and past activity at the GC in the Milky Way galaxy. Despite their importance, the formation mechanism of the bub… ▽ More The Fermi bubbles, two giant structures above and below the Galactic center (GC), are among the most important discoveries of the Fermi Gamma-ray Space Telescope. Studying their physical origin has been providing valuable insights into cosmic-ray transport, the Galactic magnetic field, and past activity at the GC in the Milky Way galaxy. Despite their importance, the formation mechanism of the bubbles is still elusive. Over the past few years there have been numerous efforts, both observational and theoretical, to uncover the nature of the bubbles. In this article, we present an overview of the current status of our understanding of the bubbles' origin, and discuss possible future directions that will help to distinguish different scenarios of bubble formation. △ Less

Submitted 28 February, 2018; v1 submitted 11 February, 2018; originally announced February 2018.

Comments: Matched to the published version in the special issue of Galaxies

Journal ref: Galaxies 6 (2018) 29

Impact of Cosmic Ray Transport on Galactic Winds

Authors: R. Farber, M. Ruszkowski, H. -Y. K. Yang, E. G. Zweibel

Abstract: The role of cosmic rays generated by supernovae and young stars has very recently begun to receive significant attention in studies of galaxy formation and evolution due to the realization that cosmic rays can efficiently accelerate galactic winds. Microscopic cosmic ray transport processes are fundamental for determining the efficiency of cosmic ray wind driving. Previous studies focused on model… ▽ More The role of cosmic rays generated by supernovae and young stars has very recently begun to receive significant attention in studies of galaxy formation and evolution due to the realization that cosmic rays can efficiently accelerate galactic winds. Microscopic cosmic ray transport processes are fundamental for determining the efficiency of cosmic ray wind driving. Previous studies focused on modeling of cosmic ray transport either via constant diffusion coefficient or via streaming proportional to the Alfven speed. However, in predominantly cold, neutral gas, cosmic rays can propagate faster than in the ionized medium and the effective transport can be substantially larger; i.e., cosmic rays can decouple from the gas. We perform three-dimensional magnetohydrodynamical simulations of patches of galactic disks including the effects of cosmic rays. Our simulations include the decoupling of cosmic rays in the cold, neutral interstellar medium. We find that, compared to the ordinary diffusive cosmic ray transport case, accounting for the decoupling leads to significantly different wind properties such as the gas density and temperature, significantly broader spatial distribution of cosmic rays, and larger wind speed. These results have implications for X-ray, gamma-ray and radio emission, and for the magnetization and pollution of the circumgalactic medium by cosmic rays. △ Less

Submitted 23 February, 2018; v1 submitted 14 July, 2017; originally announced July 2017.

Comments: 14 pages, 6 figures, Accepted for publication in the Astrophysical Journal

The Spatially Uniform Spectrum of the Fermi Bubbles: the Leptonic AGN Jet Scenario

Authors: H. -Y. K. Yang, M. Ruszkowski

Abstract: The Fermi bubbles, two giant gamma-ray bubbles above and below the Galactic center (GC), are among the most important findings of the Fermi Gamma-ray Space Telescope. Because of the proximity, spatially resolved, multi-wavelength observations offer excellent opportunities to learn about the physical origin of the bubbles. One of the unique features of the observed bubbles is that their gamma-ray s… ▽ More The Fermi bubbles, two giant gamma-ray bubbles above and below the Galactic center (GC), are among the most important findings of the Fermi Gamma-ray Space Telescope. Because of the proximity, spatially resolved, multi-wavelength observations offer excellent opportunities to learn about the physical origin of the bubbles. One of the unique features of the observed bubbles is that their gamma-ray spectrum, including a high-energy cutoff at ~110 GeV and the overall shape of the spectrum, is nearly spatially uniform. The high-energy spectral cutoff is suggestive of a leptonic origin as it could be a signature of synchrotron and inverse-Compton (IC) cooling of cosmic-ray (CR) electrons; however, even for a leptonic model, it is not obvious why the spectrum should be spatially uniform. In this work, we investigate the formation of the Fermi bubble in the leptonic jet scenario using a newly implemented CRSPEC module in FLASH that allows us to track the evolution of CR spectrum on-the-fly during the simulations. We show that the high-energy cutoff is caused by fast synchrotron and IC cooling near the GC when the jets were first launched. After the jets propagate away from the GC, the dynamical timescale of the jets become the shortest among all relevant timescales, and therefore the spectrum is essentially advected with only mild cooling losses. This could explain why the bubble spectrum is nearly spatially uniform: the CRs from different parts of the bubbles as seen today all share the same origin near the GC at early stage of the bubble expansion. We find that the CR spatial and spectral distribution predicted by the leptonic jet model can simultaneously match the normalization, spectral shape, and high-energy cutoff of the observed gamma-ray spectrum and their spatial uniformity, suggesting that past AGN jet activity is a likely mechanism for the formation of the Fermi bubbles. △ Less

Submitted 30 November, 2017; v1 submitted 15 June, 2017; originally announced June 2017.

Comments: 14 pages, 7 figures, matched to the ApJ published version, added Fig. 4 of CR spectral evolution

Journal ref: 2017, ApJ, 850, 2

The Effects of Ram Pressure on the Cold Clouds in the Centers of Galaxy Clusters

Authors: Yuan Li, Mateusz Ruszkowski, Grant Tremblay

Abstract: We discuss the effect of ram pressure on the cold clouds in the centers of cool-core galaxy clusters, and in particular, how it reduces cloud velocity and sometimes causes an offset between the cold gas and young stars. The velocities of the molecular gas in both observations and our simulations fall in the range of $100-400$ km/s, much lower than expected if they fall from a few tens of kpc balli… ▽ More We discuss the effect of ram pressure on the cold clouds in the centers of cool-core galaxy clusters, and in particular, how it reduces cloud velocity and sometimes causes an offset between the cold gas and young stars. The velocities of the molecular gas in both observations and our simulations fall in the range of $100-400$ km/s, much lower than expected if they fall from a few tens of kpc ballistically. If the intra-cluster medium (ICM) is at rest, the ram pressure of the ICM only slightly reduces the velocity of the clouds. When we assume that the clouds are actually "fluffier" because they are co-moving with a warm-hot layer, the velocity becomes smaller. If we also consider the AGN wind in the cluster center by adding a wind profile measured from the simulation, the clouds are further slowed down at small radii, and the resulting velocities are in general agreement with the observations and simulations. Because ram pressure only affects gas but not stars, it can cause a separation between a filament and young stars that formed in the filament as they move through the ICM together. This separation has been observed in Perseus and also exists in our simulations. We show that the star-filament offset combined with line-of-sight velocity measurements can help determine the true motion of the cold gas, and thus distinguish between inflows and outflows. △ Less

Submitted 20 March, 2017; originally announced March 2017.

Comments: 9 pages, 10 figures. Submitted to ApJ

Cosmic ray feedback heating of the intracluster medium

Authors: M. Ruszkowski, H. -Y. K. Yang, C. S. Reynolds

Abstract: Self-regulating AGN feedback in the cool cores of galaxy clusters plays central role in solving the decades-old cooling flow problem, but one major problem remains unsolved - how is the AGN energy thermalized in the ICM and what are the effective black hole feeding rates in realistic systems? We perform a suite of 3D MHD AMR simulations of AGN feedback in a cool core cluster including cosmic ray (… ▽ More Self-regulating AGN feedback in the cool cores of galaxy clusters plays central role in solving the decades-old cooling flow problem, but one major problem remains unsolved - how is the AGN energy thermalized in the ICM and what are the effective black hole feeding rates in realistic systems? We perform a suite of 3D MHD AMR simulations of AGN feedback in a cool core cluster including cosmic ray (CR) physics. CRs are supplied to the ICM via collimated AGN jets and subsequently disperse in the magnetized ICM via streaming, and interact with the ICM via hadronic, Coulomb, and streaming instability heating. We find that CR transport is an essential model ingredient needed for AGN feedback to self-regulate, at least within the context of the physical model considered here. When CR streaming is neglected, the suppression of CR mixing with the ICM by magnetic fields significantly reduces ICM heating, which leads to cooling catastrophes. In the opposite case, CRs come into contact with the ambient ICM and efficiently heat it, which results in globally stable atmospheres. Moreover, the dynamical state and intermittency of the central AGN are dramatically altered when CR streaming is present. We find that CR streaming heating dominates over the heating due to Coulomb and hadronic processes. Importantly, in simulations that include CR streaming, CR pressure support in the central 100 kpc is very low and does not demonstrably violate observational constraints. On the contrary, when CR streaming is neglected, CR energy is not spent on the ICM heating and CR pressure builds up to the level that is in disagreement with the data. Overall, our models demonstrate that CR heating is a viable channel for the thermalization of AGN energy in clusters, and likely also in elliptical galaxies, and that CRs play an important role in determining AGN intermittency and the dynamical state of cool core atmospheres. △ Less

Submitted 25 January, 2017; originally announced January 2017.

Comments: submitted to ApJ

AGN Heating in Simulated Cool-Core Clusters

Authors: Yuan Li, Mateusz Ruszkowski, Greg L. Bryan

Abstract: We analyze heating and cooling processes in an idealized simulation of a cool-core cluster, where momentum-driven AGN feedback balances radiative cooling in a time-averaged sense. We find that, on average, energy dissipation via shock waves is almost an order of magnitude higher than via turbulence. Most of the shock waves in the simulation are very weak shocks with Mach numbers smaller than 1.5,… ▽ More We analyze heating and cooling processes in an idealized simulation of a cool-core cluster, where momentum-driven AGN feedback balances radiative cooling in a time-averaged sense. We find that, on average, energy dissipation via shock waves is almost an order of magnitude higher than via turbulence. Most of the shock waves in the simulation are very weak shocks with Mach numbers smaller than 1.5, but the stronger shocks, although rare, dissipate energy more effectively. We find that shock dissipation is a steep function of radius, with most of the energy dissipated within 30 kpc, while radiative cooling loses area less concentrated. However, adiabatic processes and mixing (of post-shock materials and the surrounding gas) are able to redistribute the heat throughout the core. A considerable fraction of the AGN energy also escapes the core region. The cluster goes through cycles of AGN outbursts accompanied by periods of enhanced precipitation and star formation, over Gyr timescales. The cluster core is under-heated at the end of each cycle, but over-heated at the peak of the AGN outburst. During the heating-dominant phase, turbulent dissipation alone is often able to balance radiative cooling at every radius but, when this is occurs, shock waves inevitably dissipate even more energy. Our simulation explains why some clusters, such as Abell 2029, are cooling dominated, while in some other clusters, such as Perseus, various heating mechanisms including shock heating, turbulent dissipation and bubble mixing can all individually balance cooling, and together, overheat the core. △ Less

Submitted 16 November, 2016; originally announced November 2016.

Comments: 12 pages, 11 figures, submitted to ApJ

The Efficiency of Magnetic Field Amplification at Shocks by Turbulence

Authors: Suoqing Ji, S. Peng Oh, Mateusz Ruszkowski, Maxim Markevitch

Abstract: Turbulent dynamo field amplification has often been invoked to explain the strong field strengths in thin rims in supernova shocks ($\sim 100 \, μ$G) and in radio relics in galaxy clusters ($\sim μ$G). We present high resolution MHD simulations of the interaction between pre-shock turbulence, clumping and shocks, to quantify the conditions under which turbulent dynamo amplification can be signific… ▽ More Turbulent dynamo field amplification has often been invoked to explain the strong field strengths in thin rims in supernova shocks ($\sim 100 \, μ$G) and in radio relics in galaxy clusters ($\sim μ$G). We present high resolution MHD simulations of the interaction between pre-shock turbulence, clumping and shocks, to quantify the conditions under which turbulent dynamo amplification can be significant. We demonstrate numerically converged field amplification which scales with Alfvén Mach number, $B/B_0 \propto {\mathcal M}_{\rm A}$, up to ${\mathcal M}_{\rm A} \sim 150$. This implies that the post-shock field strength is relatively independent of the seed field. Amplification is dominated by compression at low ${\mathcal M}_{\rm A}$, and stretching (turbulent amplification) at high ${\mathcal M}_{\rm A}$. For high $\mathcal{M}_{\rm A}$, the $B$-field grows exponentially and saturates at equipartition with turbulence, while the vorticity jumps sharply at the shock and subsequently decays; the resulting field is orientated predominately along the shock normal (an effect only apparent in 3D and not 2D). This agrees with the radial field bias seen in supernova remnants. By contrast, for low $\mathcal{M}_{\rm A}$, field amplification is mostly compressional, relatively modest, and results in a predominantly perpendicular field. The latter is consistent with the polarization seen in radio relics. Our results are relatively robust to the assumed level of gas clumping. Our results imply that the turbulent dynamo may be important for supernovae, but is only consistent with the field strength, and not geometry, for cluster radio relics. For the latter, this implies strong pre-existing $B$-fields in the ambient cluster outskirts. △ Less

Submitted 11 October, 2016; v1 submitted 28 March, 2016; originally announced March 2016.

Comments: 15 pages, 11 figures, published version on MNRAS

Global simulations of galactic winds including cosmic ray streaming

Authors: Mateusz Ruszkowski, H. -Y. Karen Yang, Ellen Zweibel

Abstract: Galactic outflows play an important role in galactic evolution. Despite their importance, a detailed understanding of the physical mechanisms responsible for the driving of these winds is lacking. In an effort to gain more insight into the nature of these flows, we perform global three-dimensional magneto-hydrodynamical simulations of an isolated Milky Way-size starburst galaxy. We focus on the dy… ▽ More Galactic outflows play an important role in galactic evolution. Despite their importance, a detailed understanding of the physical mechanisms responsible for the driving of these winds is lacking. In an effort to gain more insight into the nature of these flows, we perform global three-dimensional magneto-hydrodynamical simulations of an isolated Milky Way-size starburst galaxy. We focus on the dynamical role of cosmic rays injected by supernovae, and specifically on the impact of the streaming and anisotropic diffusion of cosmic rays along the magnetic fields. We find that these microphysical effects can have a significant effect on the wind launching and mass loading factors depending on the details of the plasma physics. Due to the cosmic ray streaming instability, cosmic rays propagating in the interstellar medium scatter on self-excited Alfven waves and couple to the gas. When the wave growth due to the streaming instability is inhibited by some damping process, such as the turbulent damping, the cosmic ray coupling to the gas is weaker and their effective propagation speed faster than the Alfven speed. Alternatively, cosmic rays could scatter from "extrinsic turbulence" that is driven by another mechanism. We demonstrate that the presence of moderately super-Alfvenic cosmic ray streaming enhances the efficiency of galactic wind driving. Cosmic rays stream away from denser regions near the galactic disk along partially ordered magnetic fields and, in the process, accelerate more tenuous gas away from the galaxy. For cosmic ray acceleration efficiencies broadly consistent with the observational constraints, cosmic rays reduce the galactic star formation rates and significantly aid in launching galactic winds. △ Less

Submitted 15 February, 2016; originally announced February 2016.

Comments: submitted to ApJ

Deep Chandra observation and numerical studies of the nearest cluster cold front in the sky

Authors: N. Werner, J. A. ZuHone, I. Zhuravleva, Y. Ichinohe, A. Simionescu, S. W. Allen, M. Markevitch, A. C. Fabian, U. Keshet, E. Roediger, M. Ruszkowski, J. S. Sanders

Abstract: We present the results of a very deep (500 ks) Chandra observation, along with tailored numerical simulations, of the nearest, best resolved cluster cold front in the sky, which lies 90 kpc (19 arcmin) to the north-west of M 87. The northern part of the front appears the sharpest, with a width smaller than 2.5 kpc (1.5 Coulomb mean free paths; at 99 per cent confidence). Everywhere along the front… ▽ More We present the results of a very deep (500 ks) Chandra observation, along with tailored numerical simulations, of the nearest, best resolved cluster cold front in the sky, which lies 90 kpc (19 arcmin) to the north-west of M 87. The northern part of the front appears the sharpest, with a width smaller than 2.5 kpc (1.5 Coulomb mean free paths; at 99 per cent confidence). Everywhere along the front, the temperature discontinuity is narrower than 4-8 kpc and the metallicity gradient is narrower than 6 kpc, indicating that diffusion, conduction and mixing are suppressed across the interface. Such transport processes can be naturally suppressed by magnetic fields aligned with the cold front. Interestingly, comparison to magnetohydrodynamic simulations indicates that in order to maintain the observed sharp density and temperature discontinuities, conduction must also be suppressed along the magnetic field lines. However, the northwestern part of the cold front is observed to have a non-zero width. While other explanations are possible, the broadening is consistent with the presence of Kelvin-Helmholtz instabilities (KHI) on length-scales of a few kpc. Based on comparison with simulations, the presence of KHI would imply that the effective viscosity of the intracluster medium is suppressed by more than an order of magnitude with respect to the isotropic Spitzer-like temperature dependent viscosity. Underneath the cold front, we observe quasi-linear features that are ~10 per cent brighter than the surrounding gas and are separated by ~15 kpc from each other in projection. Comparison to tailored numerical simulations suggests that the observed phenomena may be due to the amplification of magnetic fields by gas sloshing in wide layers below the cold front, where the magnetic pressure reaches ~5-10 per cent of the thermal pressure, reducing the gas density between the bright features. △ Less

Submitted 16 November, 2015; v1 submitted 21 June, 2015; originally announced June 2015.

Comments: Published in MNRAS 455, 846

Journal ref: 2016 MNRAS 455, 846

Cooling, AGN Feedback and Star Formation in Simulated Cool-Core Galaxy Clusters

Authors: Yuan Li, Greg L. Bryan, Mateusz Ruszkowski, G. Mark Voit, Brian W. O'Shea, Megan Donahue

Abstract: Numerical simulations of active galactic nuclei (AGN) feedback in cool-core galaxy clusters have successfully avoided classical cooling flows, but often produce too much cold gas. We perform adaptive mesh simulations that include momentum-driven AGN feedback, self-gravity, star formation and stellar feedback, focusing on the interplay between cooling, AGN heating and star formation in an isolated… ▽ More Numerical simulations of active galactic nuclei (AGN) feedback in cool-core galaxy clusters have successfully avoided classical cooling flows, but often produce too much cold gas. We perform adaptive mesh simulations that include momentum-driven AGN feedback, self-gravity, star formation and stellar feedback, focusing on the interplay between cooling, AGN heating and star formation in an isolated cool-core cluster. Cold clumps triggered by AGN jets and turbulence form filamentary structures tens of kpc long. This cold gas feeds both star formation and the supermassive black hole (SMBH), triggering an AGN outburst that increases the entropy of the ICM and reduces its cooling rate. Within 1-2 Gyr, star formation completely consumes the cold gas, leading to a brief shutoff of the AGN. The ICM quickly cools and redevelops multiphase gas, followed by another cycle of star formation/AGN outburst. Within 6.5 Gyr, we observe three such cycles. There is good agreement between our simulated cluster and the observations of cool-core clusters. ICM cooling is dynamically balanced by AGN heating, and a cool-core appearance is preserved. The minimum cooling time to free-fall time ratio typically varies between a few and $\gtrsim 20$. The star formation rate (SFR) covers a wide range, from 0 to a few hundred $\rm M_{\odot}\, yr^{-1}$, with an average of $\sim 40 \,\rm M_{\odot}\, yr^{-1}$. The instantaneous SMBH accretion rate shows large variations on short timescales, but the average value correlates well with the SFR. Simulations without stellar feedback or self-gravity produce qualitatively similar results, but a lower SMBH feedback efficiency (0.1% compared to 1%) results in too many stars. △ Less

Submitted 9 March, 2015; originally announced March 2015.

Comments: 16 pages, 13 figures, submitted to ApJ

Ram pressure stripping in elliptical galaxies: II. magnetic field effects

Authors: Min-Su Shin, Mateusz Ruszkowski

Abstract: We investigate the effects of magnetic fields and turbulence on ram pressure stripping in elliptical galaxies using ideal magnetohydrodynamics simulations. We consider weakly-magnetised interstellar medium (ISM) characterised by subsonic turbulence, and two orientations of the magnetic fields in the intracluster medium (ICM) - parallel and perpendicular to the direction of the galaxy motion throug… ▽ More We investigate the effects of magnetic fields and turbulence on ram pressure stripping in elliptical galaxies using ideal magnetohydrodynamics simulations. We consider weakly-magnetised interstellar medium (ISM) characterised by subsonic turbulence, and two orientations of the magnetic fields in the intracluster medium (ICM) - parallel and perpendicular to the direction of the galaxy motion through the ICM. While the stronger turbulence enhances the ram pressure stripping mass loss, the magnetic fields tend to suppress the stripping rates, and the suppression is stronger for parallel fields. However, the effect of magnetic fields on the mass stripping rate is mild. Nevertheless, the morphology of the stripping tails depends significantly on the direction of the ICM magnetic field. The effect of the magnetic field geometry on the tail morphology is much stronger than that of the level of the ISM turbulence. The tail has a highly collimated shape for parallel fields, while it has a sheet-like morphology in the plane of the ICM magnetic field for perpendicular fields. The magnetic field in the tail is amplified irrespectively of the orientation of the ICM field. More strongly magnetised regions in the ram pressure stripping tails are expected to have systematically higher metallicity due to the strong concentration of the stripped ISM than the less magnetised regions. Strong dependence of the morphology of the stripped ISM on the magnetic field could potentially be used to constrain the relative orientation of the ram pressure direction and the dominant component of the ICM magnetic field. △ Less

Submitted 11 September, 2014; originally announced September 2014.

Comments: 19 pages, 21 figures, accepted for publication in MNRAS

Fermi bubble simulations: black hole feedback in the Milky Way

Authors: M. Ruszkowski, H. -Y. K. Yang, E. Zweibel

Abstract: The $Fermi$ gamma-ray telescope discovered a pair of bubbles at the Galactic center. These structures are spatially-correlated with the microwave emission detected by the WMAP and Planck satellites. These bubbles were likely inflated by a jet launched from the vicinity of a supermassive black hole in the Galactic center. Using MHD simulations, which self-consistently include interactions between c… ▽ More The $Fermi$ gamma-ray telescope discovered a pair of bubbles at the Galactic center. These structures are spatially-correlated with the microwave emission detected by the WMAP and Planck satellites. These bubbles were likely inflated by a jet launched from the vicinity of a supermassive black hole in the Galactic center. Using MHD simulations, which self-consistently include interactions between cosmic rays and magnetic fields, we build models of the supersonic jet propagation, cosmic ray transport, and the magnetic field amplification within the $Fermi$ bubbles. Our key findings are that: (1) the synthetic $Fermi$ gamma-ray and WMAP microwave spectra based on our simulations are consistent with the observations, suggesting that a single population of cosmic ray leptons may simultaneously explain the emission across a range of photon energies; (2) the model fits the observed centrally-peaked microwave emission if a second, more recent, pair of jets embedded in the $Fermi$ bubbles is included in the model. This is consistent with the observationally-based suggestion made by Su & Finkbeiner (2012); (3) the radio emission from the bubbles is expected to be strongly polarized due to the relatively high level of field ordering caused by elongated turbulent vortices. This effect is caused by the interaction of the shocks driven by the jets with the preexisting interstellar medium turbulence; (4) a layer of enhanced rotation measure in the shock-compressed region could exist in the bubble vicinity but the level of this enhancement depends on the details of the magnetic topology. △ Less

Submitted 24 November, 2013; originally announced November 2013.

Comments: IAU Symposium No. 303, The Galactic Center: Feeding and Feedback in a Normal Galactic Nucleus

The origin of cold gas in giant elliptical galaxies and its role in fueling radio-mode AGN feedback

Authors: N. Werner, J. B. R. Oonk, M. Sun, P. E. J. Nulsen, S. W. Allen, R. E. A. Canning, A. Simionescu, A. Hoffer, T. Connor, M. Donahue, A. C. Edge, A. C. Fabian, A. von der Linden, C. S. Reynolds, M. Ruszkowski

Abstract: We present a multi-wavelength study of the interstellar medium in eight nearby, X-ray and optically bright, giant elliptical galaxies. Using Herschel PACS, we map the cold gas in the lines of [CII], [OI], and [OIb]. Additionally, we present Ha+[NII] imaging of warm ionized gas with the SOAR telescope, and a study of the hot X-ray emitting plasma with Chandra. All systems with extended Ha emission… ▽ More We present a multi-wavelength study of the interstellar medium in eight nearby, X-ray and optically bright, giant elliptical galaxies. Using Herschel PACS, we map the cold gas in the lines of [CII], [OI], and [OIb]. Additionally, we present Ha+[NII] imaging of warm ionized gas with the SOAR telescope, and a study of the hot X-ray emitting plasma with Chandra. All systems with extended Ha emission in our sample (6/8 galaxies) display significant [CII] line emission indicating the presence of cold gas. This emission is co-spatial with the Ha+[NII] emitting nebulae and the lowest entropy X-ray emitting plasma. The entropy profiles of the hot galactic atmospheres show a clear dichotomy, with the systems displaying extended emission line nebulae having lower entropies beyond r~1 kpc than the cold-gas-poor systems. We show that while the hot atmospheres of the cold-gas-poor galaxies are thermally stable outside of their innermost cores, the atmospheres of the cold-gas-rich systems are prone to cooling instabilities. This result indicates that the cold gas is produced chiefly by thermally unstable cooling from the hot phase. We show that cooling instabilities may develop more easily in rotating systems and discuss an alternative condition for thermal instability for this case. The hot atmospheres of cold-gas-rich galaxies display disturbed morphologies indicating that the accretion of clumpy multiphase gas in these systems may result in variable power output of the AGN jets, potentially triggering sporadic, larger outbursts. In the two cold-gas-poor, X-ray morphologically relaxed galaxies of our sample, NGC 1399 and NGC 4472, powerful AGN outbursts may have destroyed or removed most of the cold gas from the cores, allowing the jets to propagate and deposit most of their energy further out, increasing the entropy of the hot galactic atmospheres and leaving their cores relatively undisturbed. △ Less

Submitted 13 January, 2014; v1 submitted 21 October, 2013; originally announced October 2013.

Comments: MNRAS in press

The Fermi Bubbles: Gamma-ray, Microwave, and Polarization Signatures of Leptonic AGN Jets

Authors: H. -Y. K. Yang, M. Ruszkowski, E. Zweibel

Abstract: The origin of the Fermi bubbles and the microwave haze is yet to be determined. To disentangle different models requires detailed comparisons between theoretical predictions and multi-wavelength observations. Our previous simulations have demonstrated that the primary features of the Fermi bubbles could be successfully reproduced by recent jet activity from the central active galactic nucleus (AGN… ▽ More The origin of the Fermi bubbles and the microwave haze is yet to be determined. To disentangle different models requires detailed comparisons between theoretical predictions and multi-wavelength observations. Our previous simulations have demonstrated that the primary features of the Fermi bubbles could be successfully reproduced by recent jet activity from the central active galactic nucleus (AGN). In this work, we generate gamma-ray and microwave maps and spectra based on the simulated properties of cosmic rays (CRs) and magnetic fields in order to examine whether the observed bubble and haze emission could be explained by leptons contained in the AGN jets. We also investigate the model predictions of the polarization properties of the Fermi bubbles. We find that: (1) The same population of leptons can simultaneously explain the bubble and haze emission given that the magnetic fields within the bubbles are very close to the exponentially distributed ambient field, which can be explained by mixing in of the ambient field followed by turbulent field amplification; (2) The centrally peaked microwave profile suggests CR replenishment, which is consistent with the presence of a more recent second jet event; (3) The bubble interior exhibits a high degree of polarization because of ordered radial magnetic field lines stretched by elongated vortices behind the shocks; highly-polarized signals could also be observed inside the draping layer; (4) Enhancement of rotation measures could exist within the shock-compressed layer because of increased gas density and more amplified and ordered magnetic fields. We discuss the possibility that the deficient haze emission at b<-35 degrees is due to the suppression of magnetic fields, which is consistent with the existence of lower-energy CRs causing the polarized emission at 2.3 GHz. Possible AGN jet composition in the leptonic scenario is also discussed. △ Less

Submitted 29 October, 2013; v1 submitted 12 July, 2013; originally announced July 2013.

Comments: 15 pages, 9 figures, matched with MNRAS published version

Powering of cool filaments in cluster cores by buoyant bubbles. I. Qualitative model

Authors: E. Churazov, M. Ruszkowski, A. Schekochihin

Abstract: Cool-core clusters (e.g., Perseus or M87) often possess a network of bright gaseous filaments, observed in radio, infrared, optical and X-ray bands. We propose that these filaments are powered by the reconnection of the magnetic field in the wakes of buoyant bubbles. AGN-inflated bubbles of relativistic plasma rise buoyantly in the cluster atmosphere, stretching and amplifying the field in the wak… ▽ More Cool-core clusters (e.g., Perseus or M87) often possess a network of bright gaseous filaments, observed in radio, infrared, optical and X-ray bands. We propose that these filaments are powered by the reconnection of the magnetic field in the wakes of buoyant bubbles. AGN-inflated bubbles of relativistic plasma rise buoyantly in the cluster atmosphere, stretching and amplifying the field in the wake to values of $\displaystyle β=8πP_{gas}/B^2\sim 1$. The field lines in the wake have opposite directions and are forced together as the bubble motion stretches the filament. This setup bears strong similarity to the coronal loops on the Sun or to the Earth's magneto-tail. The reconnection process naturally explains both the required level of local dissipation rate in filaments and the overall luminosity of filaments. The original source of power for the filaments is the potential energy of buoyant bubbles, inflated by the central AGN. △ Less

Submitted 11 September, 2013; v1 submitted 10 April, 2013; originally announced April 2013.

Comments: 6 pages, 1 figure, MNRAS, accepted

Chaotic cold accretion onto black holes

Authors: M. Gaspari, M. Ruszkowski, S. Peng Oh

Abstract: Using 3D AMR simulations, linking the 50 kpc to the sub-pc scales over the course of 40 Myr, we systematically relax the classic Bondi assumptions in a typical galaxy hosting a SMBH. In the realistic scenario, where the hot gas is cooling, while heated and stirred on large scales, the accretion rate is boosted up to two orders of magnitude compared with the Bondi prediction. The cause is the nonli… ▽ More Using 3D AMR simulations, linking the 50 kpc to the sub-pc scales over the course of 40 Myr, we systematically relax the classic Bondi assumptions in a typical galaxy hosting a SMBH. In the realistic scenario, where the hot gas is cooling, while heated and stirred on large scales, the accretion rate is boosted up to two orders of magnitude compared with the Bondi prediction. The cause is the nonlinear growth of thermal instabilities, leading to the condensation of cold clouds and filaments when t_cool/t_ff < 10. Subsonic turbulence of just over 100 km/s (M > 0.2) induces the formation of thermal instabilities, even in the absence of heating, while in the transonic regime turbulent dissipation inhibits their growth (t_turb/t_cool < 1). When heating restores global thermodynamic balance, the formation of the multiphase medium is violent, and the mode of accretion is fully cold and chaotic. The recurrent collisions and tidal forces between clouds, filaments and the central clumpy torus promote angular momentum cancellation, hence boosting accretion. On sub-pc scales the clouds are channelled to the very centre via a funnel. A good approximation to the accretion rate is the cooling rate, which can be used as subgrid model, physically reproducing the boost factor of 100 required by cosmological simulations, while accounting for fluctuations. Chaotic cold accretion may be common in many systems, such as hot galactic halos, groups, and clusters, generating high-velocity clouds and strong variations of the AGN luminosity and jet orientation. In this mode, the black hole can quickly react to the state of the entire host galaxy, leading to efficient self-regulated AGN feedback and the symbiotic Magorrian relation. During phases of overheating, the hot mode becomes the single channel of accretion (with a different cuspy temperature profile), though strongly suppressed by turbulence. △ Less

Submitted 22 April, 2013; v1 submitted 14 January, 2013; originally announced January 2013.

Comments: Accepted by MNRAS: added comments and references. Your feedback is welcome

The nature of filamentary cold gas in the core of the Virgo Cluster

Authors: N. Werner, J. B. R. Oonk, R. E. A. Canning, S. W. Allen, A. Simionescu, J. Kos, R. J. van Weeren, A. C. Edge, A. C. Fabian, A. von der Linden, P. E. J. Nulsen, C. S. Reynolds, M. Ruszkowski

Abstract: We present a multi-wavelength study of the emission-line nebulae located southeast of the nucleus of M87, the central dominant galaxy of the Virgo Cluster. We report the detection of far-infrared (FIR) [CII] line emission from the nebulae using observations made with Herschel PACS. The infrared line emission is extended and cospatial with optical Hα+[NII], far-ultraviolet CIV lines, and soft X-ray… ▽ More We present a multi-wavelength study of the emission-line nebulae located southeast of the nucleus of M87, the central dominant galaxy of the Virgo Cluster. We report the detection of far-infrared (FIR) [CII] line emission from the nebulae using observations made with Herschel PACS. The infrared line emission is extended and cospatial with optical Hα+[NII], far-ultraviolet CIV lines, and soft X-ray emission. The filamentary nebulae evidently contain multi-phase material spanning a temperature range of at least 5 orders of magnitude, from ~100 K to ~10^7 K. This material has most likely been uplifted by the AGN from the center of M87. The thermal pressure of the 10^4 K phase appears to be significantly lower than that of the surrounding hot intra-cluster medium (ICM) indicating the presence of additional turbulent and magnetic pressure in the filaments. If the turbulence in the filaments is subsonic then the magnetic field strength required to balance the pressure of the surrounding ICM is B~30-70 μG. The spectral properties of the soft X-ray emission from the filaments indicate that it is due to thermal plasma with kT~0.5-1 keV, which is cooling by mixing with the cold gas and/or radiatively. Charge exchange can be ruled out as a significant source of soft X-rays. Both cooling and mixing scenarios predict gas with a range of temperatures. This is at first glance inconsistent with the apparent lack of X-ray emitting gas with kT<0.5 keV. However, we show that the missing very soft X-ray emission could be absorbed by the cold gas in the filaments with an integrated hydrogen column density of ~1.6x10^21 cm^-2, providing a natural explanation for the apparent temperature floor to the X-ray emission at kT~0.5 keV. The FIR through ultra-violet line emission is most likely primarily powered by the ICM particles penetrating the cold gas following a shearing induced mixing process. △ Less

Submitted 6 March, 2013; v1 submitted 28 November, 2012; originally announced November 2012.

Comments: accepted for publication in the ApJ (scheduled for the April 20, 2013 issue)

Simulations of bent-double radio sources in galaxy groups

Authors: Brian J. Morsony, Jacob J. Miller, Sebastian Heinz, Emily Freeland, Eric Wilcots, Marcus Brueggen, Mateusz Ruszkowski

Abstract: Bent-double radio sources have been used as a probe to measure the density of intergalactic gas in galaxy groups. We carry out a series of high-resolution, 3D simulations of AGN jets moving through an external medium with a constant density in order to develop a general formula for the radius of curvature of the jets, and to determine how accurately the density of the intra-group medium (IGM) can… ▽ More Bent-double radio sources have been used as a probe to measure the density of intergalactic gas in galaxy groups. We carry out a series of high-resolution, 3D simulations of AGN jets moving through an external medium with a constant density in order to develop a general formula for the radius of curvature of the jets, and to determine how accurately the density of the intra-group medium (IGM) can be measured. Our simulations produce curved jets ending in bright radio lobes with an extended trail of low surface brightness radio emission. The radius of curvature of the jets varies with time by only about 25%. The radio trail seen in our simulations is typically not detected in known sources, but may be detectable in lower resolution radio observations. The length of this tail can be used to determine the age of the AGN. We also use our simulation data to derive a formula for the kinetic luminosity of observed jets in terms of the radius of curvature and jet pressure. In characterizing how well observations can measure the IGM density, we find that the limited resolution of typical radio observations leads to a systematic under-estimate of the density of about 50%. The unknown angles between the observer and the direction of jet propagation and direction of AGN motion through the IGM leads to an uncertainty of about 50% in estimates of the IGM density. Previous conclusions drawn using these sources, indicating that galaxy groups contain significant reservoirs of baryons in their IGM, are still valid when considering this level of uncertainty. In addition, we model the X-ray emission expected from bent-double radio sources. We find that known sources in reasonably dense environments should be detectable in ~100 ks Chandra observations. X-ray observations of these sources would place constraints on the IGM density and AGN velocity that are complementary to radio observations. △ Less

Submitted 4 February, 2013; v1 submitted 4 October, 2012; originally announced October 2012.

Comments: 14 pages, 8 figures, accepted by MNRAS

Solving the Cooling Flow Problem through Mechanical AGN Feedback

Authors: M. Gaspari, F. Brighenti, M. Ruszkowski

Abstract: Unopposed radiative cooling of plasma would lead to the cooling catastrophe, a massive inflow of condensing gas, manifest in the core of galaxies, groups and clusters. The last generation X-ray telescopes, Chandra and XMM, have radically changed our view on baryons, indicating AGN heating as the balancing counterpart of cooling. This work reviews our extensive investigation on self-regulated heati… ▽ More Unopposed radiative cooling of plasma would lead to the cooling catastrophe, a massive inflow of condensing gas, manifest in the core of galaxies, groups and clusters. The last generation X-ray telescopes, Chandra and XMM, have radically changed our view on baryons, indicating AGN heating as the balancing counterpart of cooling. This work reviews our extensive investigation on self-regulated heating. We argue that the mechanical feedback, based on massive subrelativistic outflows, is the key to solving the cooling flow problem, i.e. dramatically quenching the cooling rates for several Gyr without destroying the cool-core structure. Using a modified version of the 3D hydrocode FLASH, we show that bipolar AGN outflows can further reproduce fundamental observed features, such as buoyant bubbles, weak shocks, metals dredge- up, and turbulence. The latter is an essential ingredient to drive nonlinear thermal instabilities, which cause the formation of extended cold gas, a residual of the quenched cooling flow and, later, fuel for the feedback engine. Compared to clusters, groups and galaxies require a gentler mechanical feedback, in order to avoid catastrophic overheating. We highlight the essential characteristics for a realistic AGN feedback, with emphasis on observational consistency. △ Less

Submitted 14 September, 2012; originally announced September 2012.

Comments: Accepted by AN; 4 pages, 2 figures

The Fermi Bubbles: Supersonic AGN Jets with Anisotropic Cosmic Ray Diffusion

Authors: H. -Y. K. Yang, M. Ruszkowski, P. M. Ricker, E. Zweibel, D. Lee

Abstract: The Fermi Gamma-ray Space Telescope reveals two large bubbles in the Galaxy, which extend nearly symmetrically ~50 degrees above and below the Galactic center (GC). Using three-dimensional (3D) magnetohydrodynamic (MHD) simulations that self-consistently include the dynamical interaction between cosmic rays (CR) and thermal gas, and anisotropic CR diffusion along the magnetic field lines, we show… ▽ More The Fermi Gamma-ray Space Telescope reveals two large bubbles in the Galaxy, which extend nearly symmetrically ~50 degrees above and below the Galactic center (GC). Using three-dimensional (3D) magnetohydrodynamic (MHD) simulations that self-consistently include the dynamical interaction between cosmic rays (CR) and thermal gas, and anisotropic CR diffusion along the magnetic field lines, we show that the key characteristics of the observed gamma-ray bubbles and the spatially-correlated X-ray features in ROSAT 1.5 keV map can be successfully reproduced by a recent jet activity from the central active galactic nucleus (AGN). We find that after taking into account the projection of the 3D bubbles onto the sky, the physical heights of the bubbles can be much smaller than previously thought, greatly reducing the formation time of the bubbles to about a Myr. This relatively small bubble age is needed to reconcile the simulations with the upper limit of bubbles ages estimated from the cooling time of high-energy electrons. No additional physical mechanisms are required to suppress large-scale hydrodynamic instabilities because the evolution time is too short for them to develop. The simulated CR bubbles are edge-brightened, which is consistent with the observed projected flat surface brightness distribution. Furthermore, we demonstrate that the sharp edges of the observed bubbles can be due to anisotropic CR diffusion along magnetic field lines that drape around the bubbles during their supersonic expansion, with suppressed perpendicular diffusion across the bubble surface. Possible causes of the slight bends of the Fermi bubbles to the west are also discussed. △ Less

Submitted 26 October, 2012; v1 submitted 17 July, 2012; originally announced July 2012.

Comments: 21 pages, 11 figures, version accepted by ApJ, added discussion on cosmic-ray cooling and heating processes

Ram pressure stripping in elliptical galaxies: I. the impact of the interstellar medium turbulence

Authors: Min-Su Shin, Mateusz Ruszkowski

Abstract: Elliptical galaxies contain X-ray emitting gas that is subject to continuous ram pressure stripping over timescales comparable to cluster ages. The gas in these galaxies is not in perfect hydrostatic equilibrium. Supernova feedback, stellar winds, or active galactic nuclei (AGN) feedback can significantly perturb the interstellar medium (ISM). Using hydrodynamical simulations, we investigate the e… ▽ More Elliptical galaxies contain X-ray emitting gas that is subject to continuous ram pressure stripping over timescales comparable to cluster ages. The gas in these galaxies is not in perfect hydrostatic equilibrium. Supernova feedback, stellar winds, or active galactic nuclei (AGN) feedback can significantly perturb the interstellar medium (ISM). Using hydrodynamical simulations, we investigate the effect of subsonic turbulence in the hot ISM on the ram pressure stripping process in early-type galaxies. We find that galaxies with more turbulent ISM produce longer, wider, and more smoothly distributed tails of the stripped ISM than those characterised by weaker ISM turbulence. Our main conclusion is that even very weak internal turbulence, at the level of <15% of the average ISM sound speed, can significantly accelerate the gas removal from galaxies via ram pressure stripping. The magnitude of this effect increases sharply with the strength of turbulence. As most of the gas stripping takes place near the boundary between the ISM and the intraclustermedium (ICM), the boost in the ISM stripping rate is due to the "random walk" of the ISM from the central regions of the galactic potential well to larger distances, where the ram pressure is able to permanently remove the gas from galaxies. The ICM can be temporarily trapped inside the galactic potential well due to the mixing of the turbulent ISM with the ICM. The galaxies with more turbulent ISM, yet still characterised by very weak turbulence, can hold larger amounts of the ICM. [Abridged] △ Less

Submitted 29 September, 2012; v1 submitted 29 June, 2012; originally announced June 2012.

Comments: 12 pages, 11 figures, 2 tables; accepted for publication in MNRAS

Cold Fronts and Gas Sloshing in Galaxy Clusters with Anisotropic Thermal Conduction

Authors: J. A. ZuHone, M. Markevitch, M. Ruszkowski, D. Lee

Abstract: (Abridged) Cold fronts in cluster cool cores should be erased on short timescales by thermal conduction, unless protected by magnetic fields that are "draped" parallel to the front surfaces, suppressing conduction perpendicular to the fronts. We present MHD simulations of cold front formation in the core of a galaxy cluster with anisotropic thermal conduction, exploring a parameter space of conduc… ▽ More (Abridged) Cold fronts in cluster cool cores should be erased on short timescales by thermal conduction, unless protected by magnetic fields that are "draped" parallel to the front surfaces, suppressing conduction perpendicular to the fronts. We present MHD simulations of cold front formation in the core of a galaxy cluster with anisotropic thermal conduction, exploring a parameter space of conduction strengths parallel and perpendicular to the field lines. Including conduction has a strong effect on the temperature of the core and the cold fronts. Though magnetic field lines are draping parallel to the front surfaces, the temperature jumps across the fronts are nevertheless reduced. The field geometry is such that the cold gas below the front surfaces can be connected to hotter regions outside via field lines along directions perpendicular to the plane of the sloshing motions and along sections of the front which are not perfectly draped. This results in the heating of this gas below the front on a timescale of a Gyr, but the sharpness of the density and temperature jumps may still be preserved. By modifying the density distribution below the front, conduction may indirectly aid in suppressing Kelvin-Helmholtz instabilities. If conduction along the field lines is unsuppressed, we find that the characteristic sharp jumps in X-ray emission seen in observations of clusters do not form. This suggests that the presence of sharp cold fronts in hot clusters could be used to place upper limits on conduction in the {\it bulk} of the ICM. Finally, the combination of sloshing and anisotropic thermal conduction can result in a larger flux of heat to the core than either process in isolation. While still not sufficient to prevent a cooling catastrophe in the very central ($r \sim$ 5 kpc) regions of the cool core, it reduces significantly the mass of cool gas that accumulates outside those radii. △ Less

Submitted 11 November, 2012; v1 submitted 26 April, 2012; originally announced April 2012.

Comments: 19 pages, 14 figures, "emulateapj" format. Updated version to match referee's comments and suggestions. Accepted by the Astrophysical Journal

Impact of magnetic fields on ram pressure stripping in disk galaxies

Authors: M. Ruszkowski, M. Bruggen, D. Lee, M. -S. Shin

Abstract: (abridged) Ram pressure can remove significant amounts of gas from galaxies in clusters, and thus has a large impact on the evolution of cluster galaxies. Recent observations have shown that key properties of ram pressure stripped tails of galaxies are in conflict with predictions by simulations. To increase the realism of existing simulations, we simulated for the first time a disk galaxy exposed… ▽ More (abridged) Ram pressure can remove significant amounts of gas from galaxies in clusters, and thus has a large impact on the evolution of cluster galaxies. Recent observations have shown that key properties of ram pressure stripped tails of galaxies are in conflict with predictions by simulations. To increase the realism of existing simulations, we simulated for the first time a disk galaxy exposed to a uniformly magnetized wind including radiative cooling and self-gravity of the gas. We find that B-fields have a strong effect on the morphology of the gas in the tail of the galaxy. While in the pure hydro case the tail is very clumpy, the MHD case shows very filamentary structures in the tail. The filaments can be strongly supported by magnetic pressure and, when this is the case, the B-field vectors tend to be aligned with the filaments. The ram pressure stripping may lead to the formation of magnetized density tails that appear as bifurcated in the plane of the sky and resemble the double tails observed in ESO 137-001 and ESO 137-002. Such tails can be formed under a variety of situations, both for the disks oriented face-on with respect to the ICM wind and for the tilted ones. While this bifurcation is due to the generic tendency for the B-fields to produce very filamentary tail morphology, the tail properties are further shaped by the combination of the B-field orientation and the sliding of the field past the disk surface exposed to the wind. Magnetic draping does not strongly change the rate of gas stripping. For a face-on galaxy, the field tends to reduce the amount of stripping compared to the pure hydro case, and is associated with the formation of a magnetic draping layer on the side of the galaxy exposed to the ICM wind. For significantly tilted disks, the stripping rate may be enhanced by the ``scraping'' of the disk surface by the B-fields sliding past the ISM/ICM interface. △ Less

Submitted 19 February, 2014; v1 submitted 6 March, 2012; originally announced March 2012.

Comments: ApJ in press. arXiv admin note: text overlap with arXiv:0909.3097 by other authors

Shallow Dark Matter Cusps in Galaxy Clusters

Authors: Chervin F. P. Laporte, Simon D. M. White, Thorsten Naab, Mateusz Ruszkowski, Volker Springel

Abstract: We study the evolution of the stellar and dark matter components in a galaxy cluster of $10^{15} \, \rm{M_{\odot}}$ from $z=3$ to the present epoch using the high-resolution collisionless simulations of Ruszkowski & Springel (2009). At $z=3$ the dominant progenitor halos were populated with spherical model galaxies with and without accounting for adiabatic contraction. We apply a weighting scheme… ▽ More We study the evolution of the stellar and dark matter components in a galaxy cluster of $10^{15} \, \rm{M_{\odot}}$ from $z=3$ to the present epoch using the high-resolution collisionless simulations of Ruszkowski & Springel (2009). At $z=3$ the dominant progenitor halos were populated with spherical model galaxies with and without accounting for adiabatic contraction. We apply a weighting scheme which allows us to change the relative amount of dark and stellar material assigned to each simulation particle in order to produce luminous properties which agree better with abundance matching arguments and observed bulge sizes at $z=3$. This permits the study of the effect of initial compactness on the evolution of the mass-size relation. We find that for more compact initial stellar distributions the size of the final Brightest Cluster Galaxy grows with mass according to $r\propto M^{2}$, whereas for more extended initial distributions, $r\propto M$. Our results show that collisionless mergers in a cosmological context can reduce the strength of inner dark matter cusps with changes in logarithmic slope of 0.3 to 0.5 at fixed radius. Shallow cusps such as those found recently in several strong lensing clusters thus do not necessarily conflict with CDM, but may rather reflect on the initial structure of the progenitor galaxies, which was shaped at high redshift by their formation process. △ Less

Submitted 10 February, 2012; originally announced February 2012.

Comments: 8 pages, 4 figures, submitted to MNRAS

Cause and Effect of Feedback: Multiphase Gas in Cluster Cores Heated by AGN Jets

Authors: M. Gaspari, M. Ruszkowski, P. Sharma

Abstract: Multiwavelength data indicate that the X-ray emitting plasma in the cores of galaxy clusters is not cooling catastrophically. To large extent, cooling is offset by heating due to active galactic nuclei (AGN) via jets. The cool-core clusters, with cooler/denser plasmas, show multiphase gas and signs of some cooling in their cores. These observations suggest that the cool core is locally thermally u… ▽ More Multiwavelength data indicate that the X-ray emitting plasma in the cores of galaxy clusters is not cooling catastrophically. To large extent, cooling is offset by heating due to active galactic nuclei (AGN) via jets. The cool-core clusters, with cooler/denser plasmas, show multiphase gas and signs of some cooling in their cores. These observations suggest that the cool core is locally thermally unstable while maintaining global thermal equilibrium. Using high-resolution, three-dimensional simulations we study the formation of multiphase gas in cluster cores heated by highly-collimated bipolar AGN jets. Our key conclusion is that spatially extended multiphase filaments form only when the instantaneous ratio of the thermal instability and free-fall timescales (t_TI/t_ff) falls below a critical threshold of \approx 10. When this happens, dense cold gas decouples from the hot ICM phase and generates inhomogeneous and spatially extended Halpha filaments. These cold gas clumps and filaments `rain' down onto the central regions of the core, forming a cold rotating torus and in part feeding the supermassive black hole. Consequently, the self-regulated feedback enhances AGN heating and the core returns to a higher entropy level with t_TI/t_ff > 10. Eventually the core reaches quasi-stable global thermal equilibrium, and cold filaments condense out of the hot ICM whenever t_TI/t_ff \lesssim 10. This occurs despite the fact that the energy from AGN jets is supplied to the core in a highly anisotropic fashion. The effective spatial redistribution of heat is enabled in part by the turbulent motions in the wake of freely-falling cold filaments. Increased AGN activity can locally reverse the cold gas flow, launching cold filamentary gas away from the cluster center. Our criterion for the condensation of spatially extended cold gas is in agreement with observations and previous idealized simulations. △ Less

Submitted 27 October, 2011; originally announced October 2011.

Comments: Submitted to ApJ; 16 pages, 11 figures

Cosmological MHD simulations of cluster formation with anisotropic thermal conduction

Authors: M. Ruszkowski, D. Lee, M. Bruggen, I. Parrish, S. Peng Oh

Abstract: (abridged) The ICM has been suggested to be buoyantly unstable in the presence of magnetic field and anisotropic thermal conduction. We perform first cosmological simulations of galaxy cluster formation that simultaneously include magnetic fields, radiative cooling and anisotropic thermal conduction. In isolated and idealized cluster models, the magnetothermal instability (MTI) tends to reorient t… ▽ More (abridged) The ICM has been suggested to be buoyantly unstable in the presence of magnetic field and anisotropic thermal conduction. We perform first cosmological simulations of galaxy cluster formation that simultaneously include magnetic fields, radiative cooling and anisotropic thermal conduction. In isolated and idealized cluster models, the magnetothermal instability (MTI) tends to reorient the magnetic fields radially. Using cosmological simulations of the Santa Barbara cluster we detect radial bias in the velocity and magnetic fields. Such radial bias is consistent with either the inhomogeneous radial gas flows due to substructures or residual MTI-driven field rearangements that are expected even in the presence of turbulence. Although disentangling the two scenarios is challenging, we do not detect excess bias in the runs that include anisotropic thermal conduction. The anisotropy effect is potentially detectable via radio polarization measurements with LOFAR and SKA and future X-ray spectroscopic studies with the IXO. We demonstrate that radiative cooling boosts the amplification of the magnetic field by about two orders of magnitude beyond what is expected in the non-radiative cases. At z=0 the field is amplified by a factor of about 10^6 compared to the uniform magnetic field evolved due to the universal expansion alone. Interestingly, the runs that include both radiative cooling and anisotropic thermal conduction exhibit stronger magnetic field amplification than purely radiative runs at the off-center locations. In these runs, shallow temperature gradients away from the cluster center make the ICM neutrally buoyant. The ICM is more easily mixed in these regions and the winding up of the frozen-in magnetic field is more efficient resulting in stronger magnetic field amplification. △ Less

Submitted 11 October, 2010; originally announced October 2010.

Comments: submitted to ApJ, higher resolution figures available at: http://www.astro.lsa.umich.edu/~mateuszr/

Galaxy Motions, Turbulence and Conduction in Clusters of Galaxies

Authors: M. Ruszkowski, S. Peng Oh

Abstract: Unopposed radiative cooling in clusters of galaxies results in excessive mass deposition rates. However, the cool cores of galaxy clusters are continuously heated by thermal conduction and turbulent heat diffusion due to minor mergers or the galaxies orbiting the cluster center. These processes can either reduce the energy requirements for AGN heating of cool cores, or they can prevent overcooling… ▽ More Unopposed radiative cooling in clusters of galaxies results in excessive mass deposition rates. However, the cool cores of galaxy clusters are continuously heated by thermal conduction and turbulent heat diffusion due to minor mergers or the galaxies orbiting the cluster center. These processes can either reduce the energy requirements for AGN heating of cool cores, or they can prevent overcooling altogether. We perform 3D MHD simulations including field-aligned thermal conduction and self-gravitating particles to model this in detail. Turbulence is not confined to the wakes of galaxies but is instead volume-filling, due to the excitation of large-scale g-modes. We systematically probe the parameter space of galaxy masses and numbers. For a wide range of observationally motivated galaxy parameters, the magnetic field is randomized by stirring motions, restoring the conductive heat flow to the core. The cooling catastrophe either does not occur or it is sufficiently delayed to allow the cluster to experience a major merger that could reset conditions in the intracluster medium. Whilst dissipation of turbulent motions is negligible as a heat source, turbulent heat diffusion is extremely important; it predominates in the cluster center. However, thermal conduction becomes important at larger radii, and simulations without thermal conduction suffer a cooling catastrophe. Conduction is important both as a heat source and to reduce stabilizing buoyancy forces, enabling more efficient diffusion. Turbulence enables conduction, and conduction enables turbulence. In these simulations, the gas vorticity---which is a good indicator of trapped g-modes--increases with time. The vorticity growth is approximately mirrored by the growth of the magnetic field, which is amplified by turbulence. △ Less

Submitted 30 August, 2010; originally announced August 2010.

Comments: Submitted to MNRAS

Swimming against the current: Simulations of central AGN evolution in dynamic galaxy clusters

Authors: Brian J. Morsony, Sebastian Heinz, Marcus Bruggen, Mateusz Ruszkowski

Abstract: We present a series of three-dimensional hydrodynamical simulations of central AGN driven jets in a dynamic, cosmologically evolved galaxy cluster. Extending previous work, we study jet powers ranging from L_jet = 10^44 erg/s to L_jet = 10^46 erg/s and in duration from 30 Myr to 200 Myr. We find that large-scale motions of cluster gas disrupt the AGN jets, causing energy to be distributed througho… ▽ More We present a series of three-dimensional hydrodynamical simulations of central AGN driven jets in a dynamic, cosmologically evolved galaxy cluster. Extending previous work, we study jet powers ranging from L_jet = 10^44 erg/s to L_jet = 10^46 erg/s and in duration from 30 Myr to 200 Myr. We find that large-scale motions of cluster gas disrupt the AGN jets, causing energy to be distributed throughout the centre of the cluster, rather than confined to a narrow angle around the jet axis. Disruption of the jet also leads to the appearance of multiple disconnected X-ray bubbles from a long-duration AGN with a constant luminosity. This implies that observations of multiple bubbles in a cluster are not necessarily an expression of the AGN duty cycle. We find that the "sphere of influence" of the AGN, the radial scale within which the cluster is strongly affected by the jet, scales as R ~ L_jet^(1/3). Increasing the duration of AGN activity does not increase the radius affected by the AGN significantly, but does change the magnitude of the AGN's effects. How an AGN delivers energy to a cluster will determine where that energy is deposited: a high luminosity is needed to heat material outside the core of the cluster, while a low-luminosity, long-duration AGN is more efficient at heating the inner few tens of kpc. △ Less

Submitted 16 October, 2010; v1 submitted 15 March, 2010; originally announced March 2010.

Comments: 23 pages, 14 figures, published in MNRAS

AGN pairs: chance superpositions or black hole binaries?

Authors: M. Dotti, M. Ruszkowski

Abstract: Several active galactic nuclei (AGN) with multiple sets of emission lines separated by over 2000 km/s have been observed recently. These have been interpreted as being due to massive black hole (MBH) recoil following a black hole merger, MBH binaries, or chance superpositions of AGN in galaxy clusters. Moreover, a number of double-peaked AGN with velocity offsets of ~ a few 100 km/s have also be… ▽ More Several active galactic nuclei (AGN) with multiple sets of emission lines separated by over 2000 km/s have been observed recently. These have been interpreted as being due to massive black hole (MBH) recoil following a black hole merger, MBH binaries, or chance superpositions of AGN in galaxy clusters. Moreover, a number of double-peaked AGN with velocity offsets of ~ a few 100 km/s have also been detected and interpreted as being due to the internal kinematics of the narrow line regions or MBH binary systems. Here we reexamine the superposition model. Using the Millennium Run we estimate the total number of detectable AGN pairs as a function of the emission line offset. We show that AGN pairs with high velocity line separations up to ~2000 km/s are very likely to be chance superpositions of two AGN in clusters of galaxies for reasonable assumptions about the relative fraction of AGN. No superimposed AGN pairs are predicted for velocity offsets in excess of ~3000 km/s as the required AGN fractions would violate observational constraints. The high velocity AGN pair numbers predicted here are competitive with those predicted from the models relying on MBH recoil or MBH binaries. However, the model fails to account for the largest emission line velocity offsets that require the presence of MBH binaries. △ Less

Submitted 29 December, 2009; originally announced December 2009.

Comments: 4 pages, 2 figures, submitted to ApJ Letters

AGN Feedback and Bimodality in Cluster Core Entropy

Authors: Fulai Guo, S. Peng Oh, M. Ruszkowski

Abstract: We investigate a series of steady-state models of galaxy clusters, in which the hot intracluster gas is efficiently heated by active galactic nucleus (AGN) feedback and thermal conduction, and in which the mass accretion rates are highly reduced compared to those predicted by the standard cooling flow models. We perform a global Lagrangian stability analysis. We show for the first time that the… ▽ More We investigate a series of steady-state models of galaxy clusters, in which the hot intracluster gas is efficiently heated by active galactic nucleus (AGN) feedback and thermal conduction, and in which the mass accretion rates are highly reduced compared to those predicted by the standard cooling flow models. We perform a global Lagrangian stability analysis. We show for the first time that the global radial instability in cool core clusters can be suppressed by the AGN feedback mechanism, provided that the feedback efficiency exceeds a critical lower limit. Furthermore, our analysis naturally shows that the clusters can exist in two distinct forms. Globally stable clusters are expected to have either: 1) cool cores stabilized by both AGN feedback and conduction, or 2) non-cool cores stabilized primarily by conduction. Intermediate central temperatures typically lead to globally unstable solutions. This bimodality is consistent with the recently observed anticorrelation between the flatness of the temperature profiles and the AGN activity (Dunn & Fabian 2008) and the observation by Rafferty et al. (2008) that the shorter central cooling times tend to correspond to significantly younger AGN X-ray cavities. △ Less

Submitted 30 November, 2009; originally announced November 2009.

Comments: 4 pages, to appear in the proceedings of "The Monster's Fiery Breath: Feedback in Galaxies, Groups, and Clusters", Eds. Sebastian Heinz, Eric Wilcots (AIP conference series)

Are black holes big enough to quench cooling in cluster cool cores?

Authors: M. Ruszkowski

Abstract: Total energy arguments (e.g., Fabian et al. 2002) suggest that black holes need to have masses significantly in excess of the prediction from the classic black hole mass - velocity dispersion relation (M-sigma) in order to offset the cooling losses in massive cool core clusters. This suggests that the black holes may be too small to power such clusters. However, Lauer et al. (2007) argue that th… ▽ More Total energy arguments (e.g., Fabian et al. 2002) suggest that black holes need to have masses significantly in excess of the prediction from the classic black hole mass - velocity dispersion relation (M-sigma) in order to offset the cooling losses in massive cool core clusters. This suggests that the black holes may be too small to power such clusters. However, Lauer et al. (2007) argue that the black hole mass - bulge luminosity relationship is a better predictor of black hole masses in high luminosity galaxies and that this relationship predicts significantly higher masses in BCGs. They find slow increase in the velocity dispersion with luminosity and a more rapid increase in effective radii with luminosity seen in BCGs as opposed to less luminous galaxies. Motivated by these results and the theoretical work of Boylan-Kolchin et al. (2006) on isolated mergers, we perform high-resolution cosmological simulations of dry mergers in a massive galaxy cluster identified in the Millennium Run including both the dark matter halos and stellar bulges of merging galaxies. We demonstrate that the BCG clearly evolves away from the size-luminosity relation as defined by the smaller galaxies (i.e., the relation bends) and we also see a bending in the luminosity-sigma relation. As black hole mass is expected to be proportional to the mass and luminosity of the stellar bulge of the BCGs (if they were formed in predominantly dissipationless mergers), our findings are consistent with those of Lauer et al. (2007) on a qualitative level and suggest that the black holes in BCGs may indeed be more massive than predicted from the standard M-sigma relation. △ Less

Submitted 27 November, 2009; originally announced November 2009.

Comments: submitted to the conference proceedings of "The Monster's Fiery Breath"