BZPEERPapers, people, evidence and research paths
Menu

Discover

SearchFind papers, people, topics, institutions and opportunities from one place.GraphSee how papers, authors, topics and evidence connect.

Research tools

HomeContinue from the papers, people and topics that matter now.CollectionsTurn saved papers into reading lists, evidence maps and shareable context.ResearchReconstruct the path that led you from search to paper to question.CollaborationMove from a paper, topic or expert request into a focused thread.PublishCreate full-text research articles with metadata, formulas and versions.

Network

ExpertsFind specialists whose work and reviews explain why they can help.CommunitiesJoin research groups organized around papers, topics and shared evidence.

Opportunities

ListingsBrowse roles, calls and collaborations connected to real research context.

Account

Log inReturn to your saved papers, graph views and active threads.Create accountStart saving papers, building a research trail and joining teams.
Log inCreate account

Source author record

M. Ruszkowski

M. Ruszkowski appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record
astro-ph.COastro-ph.HEastro-ph.GAphysics.plasm-ph

Catalog footprint

What is connected

22works
4topics
4close collaborators

Actions

Connect this record

Log in to claim
Open graphBrowse works

Research graph

See the researcher in context

Open full explorer

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

22 published item(s)

preprint2022arXiv

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

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.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

Turbulent Heating in a Stratified Medium

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.

0 citations0 reviewsNo code linkedOpen access
preprint2021arXiv

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

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.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

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

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.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

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

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.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Chaotic cold accretion onto black holes

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.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Fermi bubble simulations: black hole feedback in the Milky Way

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.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

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

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.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

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

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.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

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

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.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

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

(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.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Solving the Cooling Flow Problem through Mechanical AGN Feedback

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.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

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

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.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

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

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.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Cosmological MHD simulations of cluster formation with anisotropic thermal conduction

(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.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Galaxy Motions, Turbulence and Conduction in Clusters of Galaxies

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.

0 citations0 reviewsNo code linkedOpen access
preprint2009arXiv

AGN Feedback and Bimodality in Cluster Core Entropy

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.

0 citations0 reviewsNo code linkedOpen access
preprint2009arXiv

AGN pairs: chance superpositions or black hole binaries?

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.

0 citations0 reviewsNo code linkedOpen access
preprint2009arXiv

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

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.

0 citations0 reviewsNo code linkedOpen access
preprint2009arXiv

Dual black holes in merger remnants. I: linking accretion to dynamics

We study the orbital evolution and accretion history of massive black hole (MBH) pairs in rotationally supported circumnuclear discs up to the point where MBHs form binary systems. Our simulations have high resolution in mass and space which, for the first time, makes it feasible to follow the orbital decay of a MBH either counter- or co-rotating with respect to the circumnuclear disc. We show that a moving MBH on an initially counter-rotating orbit experiences an "orbital angular momentum flip" due to the gas-dynamical friction, i.e., it starts to corotate with the disc before a MBH binary forms. We stress that this effect can only be captured in very high resolution simulations. Given the extremely large number of gas particles used, the dynamical range is sufficiently large to resolve the Bondi-Hoyle-Lyttleton radii of individual MBHs. As a consequence, we are able to link the accretion processes to the orbital evolution of the MBH pairs. We predict that the accretion rate is significantly suppressed and extremely variable when the MBH is moving on a retrograde orbit. It is only after the orbital angular momentum flip has taken place that the secondary rapidly "lights up" at which point both MBHs can accrete near the Eddington rate for a few Myr. The separation of the double nucleus is expected to be around ~10 pc at this stage. We show that the accretion rate can be highly variable also when the MBH is co-rotating with the disc (albeit to a lesser extent) provided that its orbit is eccentric. Our results have significant consequences for the expected number of observable double AGNs at separations of <100 pc.

0 citations0 reviewsNo code linkedOpen access
preprint2009arXiv

Dual black holes in merger remnants. II: spin evolution and gravitational recoil

Using high resolution hydrodynamical simulations, we explore the spin evolution of massive dual black holes orbiting inside a circumnuclear disc, relic of a gas-rich galaxy merger. The black holes spiral inwards from initially eccentric co or counter-rotating coplanar orbits relative to the disc's rotation, and accrete gas that is carrying a net angular momentum. As the black hole mass grows, its spin changes in strength and direction due to its gravito-magnetic coupling with the small-scale accretion disc. We find that the black hole spins loose memory of their initial orientation, as accretion torques suffice to align the spins with the angular momentum of their orbit on a short timescale (<1-2 Myr). A residual off-set in the spin direction relative to the orbital angular momentum remains, at the level of <10 degrees for the case of a cold disc, and <30 degrees for a warmer disc. Alignment in a cooler disc is more effective due to the higher coherence of the accretion flow near each black hole that reflects the large-scale coherence of the disc's rotation. If the massive black holes coalesce preserving the spin directions set after formation of a Keplerian binary, the relic black hole resulting from their coalescence receives a relatively small gravitational recoil. The distribution of recoil velocities inferred from a simulated sample of massive black hole binaries has median <70 km/s much smaller than the median resulting from an isotropic distribution of spins.

0 citations0 reviewsNo code linkedOpen access
preprint2009arXiv

Shaken and stirred: conduction and turbulence in clusters of galaxies

(abridged) Uninhibited radiative cooling in clusters of galaxies would lead to excessive mass accretion rates contrary to observations. One of the key proposals to offset radiative energy losses is thermal conduction from outer, hotter layers of cool core clusters to their centers. However, conduction is sensitive to magnetic field topology. In cool-core clusters the heat buoyancy instability (HBI) leads to B-fields ordered preferentially in the direction perpendicular to that of gravity, which significantly reduces the level of conduction below the classical Spitzer-Braginskii value. However, the cluster cool cores are rarely in perfect hydrostatic equilibrium. Sloshing motions due to minor mergers, galaxy motions or AGN can significantly perturb the gas and affect the level of thermal conduction. We perform 3D AMR MHD simulations of the effect of turbulence on the properties of the anisotropic thermal conduction in cool core clusters. We show that very weak subsonic motions, well within observational constraints, can randomize the magnetic field and significantly boost effective thermal conduction beyond the saturated values expected in the pure unperturbed HBI case. We find that the turbulent motions can essentially restore the conductive heat flow to the cool core to level comparable to the theoretical maximum of 1/3 Spitzer for a highly tangled field. Runs with radiative cooling show that the cooling catastrophe can be averted and the cluster core stabilized. Above a critical Froude number, these same turbulent motions also eliminate the tangential bias in the velocity and magnetic field that is otherwise induced by the trapped g-modes. Our results can be tested with future radio polarization measurements, and have implications for efficient metal dispersal in clusters.

0 citations0 reviewsNo code linkedOpen access