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.

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\'en Mach number. While previous simulations that relied on a simplified model of CR transport have shown that super-Alfv\'enic 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.

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.

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.

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.

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.

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}$.

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.

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.

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$\alpha$ 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\'en waves by CR streaming, and the replenishment of CR energy via accretion onto the filaments of high plasma-$\beta$ 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.

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.

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.

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.

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.

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.

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.

Turbulent dynamo field amplification has often been invoked to explain the strong field strengths in thin rims in supernova shocks ($\sim 100 \, \mu$G) and in radio relics in galaxy clusters ($\sim \mu$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\'en 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.

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.

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.

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.

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.

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

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.

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.

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.

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 \beta =8\pi 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.

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.

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{\alpha}+[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 {\mu}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.

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.

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

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.

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]

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.

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.

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.

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.

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

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.

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.

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

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.

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.

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.

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.

The large scale structure of the present Universe is determined by the growth of dark matter density fluctuations and by the dynamical action of dark energy and dark matter. While much progress has been made in recent years in constraining the cosmological parameters, and in reconstructing the evolution in the large--scale structure of the dark matter distribution, we still lack an understanding of the evolution of the baryonic component of the Universe. Located at nodes of the cosmic web, clusters of galaxies are the largest collapsed structures in the Universe with total masses up to 10$^{15}$ M$_{\sun}$. Over 80% of their mass resides in the form of dark matter. The remaining mass is composed of baryons, most of which (about 85%) is a diffuse, hot plasma that radiates primarily in X-rays. X-ray observations of the evolving cluster population provide a unique opportunity to address such open and fundamental questions as: How do hot diffuse baryons dynamically evolve in dark matter potentials? How and when was the excess energy which we observe in the intergalactic medium generated? What is the cosmic history of heavy-element production and circulation? Our current knowledge comes primarily from detailed studies of clusters in the relatively nearby Universe (z$<$0.5). Major advances will come from high throughput, high spectral and spatial resolution X-ray observations that measure the thermodynamical properties and metal content of the first low mass clusters emerging at z $\sim$ 2 and directly trace their evolution into today's massive clusters.

Using a resimulation technique, we perform high-resolution cosmological simulations of dry mergers in a massive galaxy cluster identified in the Millennium Run. Our initial conditions include well resolved compound galaxy models consisting of dark matter halos and stellar bulges that are used to replace the most massive cluster progenitor halos at redshift z=3. By construction, our galaxy models obey the stellar mass-size relation initially. We demonstrate that the brightest cluster galaxy (BCG) evolves away from the Kormendy relation as defined by the smaller mass galaxies (i.e., the relation bends). As a result of the comparatively large number of mergers the BCG experiences, its total mass-to-light ratio becomes significantly higher than in typical elliptical galaxies. We also show that the mixing processes between dark matter and stars lead to a small but numerically robust tilt in the fundamental plane and that the BCG lies on the tilted plane. Our model is consistent with the observed steepening of the logarithmic mass-to-light gradient as a function of the stellar mass. Surprisingly, we find only tentative weak distortion in the Faber-Jackson relation that depends on the aperture size, unlike expected based on studies of isolated merger simulations. This may be due to differences in the distribution of galaxy orbits, which is given in our approach directly by the cosmological context while it has to be assumed in isolated merger simulations, and the fact that the BCG is located deep in the cluster potential well.

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.

Using deep Chandra observations of M84 we study the energetics of the interaction between the black hole and the interstellar medium of this early-type galaxy. We perform a detailed two dimensional reconstruction of the properties of the X-ray emitting gas using a constrained Voronoi tessellation method, identifying the mean trends and carrying out the fluctuation analysis of the thermodynamical properties of the hot ISM. In addition to the PV work associated with the bubble expansion, we identify and measure the wave energy associated with the mildly supersonic bubble expansion. We show that, depending on the age of the cavity and the associated wave, the waves can have a substantial contribution to the total energy release from the AGN. The energy dissipated in the waves tends to be concentrated near the center of M84 and in the direction perpendicular to the bubble outflow, possibly due to the interference of the waves generated by the expansion of northern and southern bubbles. We also find direct evidence for the escape of radio plasma from the ISM of the host galaxy into the intergalactic medium.

Feedback by active galactic nuclei (AGN) is frequently invoked to explain the cut-off of the galaxy luminosity function at the bright end and the absence of cooling flows in galaxy clusters. Meanwhile, there are recent observations of shock fronts around radio-loud AGN. Using realistic 3D simulations of jets in a galaxy cluster, we address the question what fraction of the energy of active galactic nuclei is dissipated in shocks. We find that weak shocks that encompass the AGN have Mach numbers of 1.1-1.2 and dissipate at least 2% of the mechanical luminosity of the AGN. In a realistic cluster medium, even a continuous jet can lead to multiple shock structures, which may lead to an overestimate of the AGN duty cycles inferred from the spatial distribution of waves.

Most cool core clusters of galaxies possess active galactic nuclei (AGN) in their centers. These AGN inflate buoyant bubbles containing non-thermal radio emitting particles. If such bubbles efficiently confine cosmic rays (CR) then this could explain ``radio ghosts'' seen far from cluster centers. We simulate the diffusion of cosmic rays from buoyant bubbles inflated by AGN. Our simulations include the effects of the anisotropic particle diffusion introduced by magnetic fields. Our models are consistent with the X-ray morphology of AGN bubbles, with disruption being suppressed by the magnetic draping effect. We conclude that for such magnetic field topologies, a substantial fraction of cosmic rays can be confined inside the bubbles on buoyant rise timescales even when the parallel diffusivity coefficient is very large. For isotropic diffusion at a comparable level, cosmic rays would leak out of the bubbles too rapidly to be consistent with radio observations. Thus, the long confinement times associated with the magnetic suppression of CR diffusion can explain the presence of radio ghosts. We show that the partial escape of cosmic rays is mostly confined to the wake of the rising bubbles, and speculate that this effect could: (1) account for the excitation of the H$\alpha$ filaments trailing behind the bubbles in the Perseus cluster, (2) inject entropy into the metal enriched material being lifted by the bubbles and, thus, help to displace it permanently from the cluster center and (3) produce observable $\gamma$-rays via the interaction of the diffusing cosmic rays with the thermal intracluster medium (ICM).

There is growing consensus that feedback from AGN is the main mechanism responsible for stopping cooling flows in clusters of galaxies. AGN are known to inflate buoyant bubbles that supply mechanical power to the intracluster gas (ICM). High Reynolds number hydrodynamical simulations show that such bubbles get entirely disrupted within 100 Myr, as they rise in cluster atmospheres, which is contrary to observations. This artificial mixing has consequences for models trying to quantify the amount of heating and star formation in cool core clusters of galaxies. It has been suggested that magnetic fields can stabilize bubbles against disruption. We perform MHD simulations of fossil bubbles in the presence of tangled magnetic fields using the high order PENCIL code. We focus on the physically-motivated case where thermal pressure dominates over magnetic pressure and consider randomly oriented fields with and without maximum helicity and a case where large scale external fields drape the bubble.We find that helicity has some stabilizing effect. However, unless the coherence length of magnetic fields exceeds the bubble size, the bubbles are quickly shredded. As observations of Hydra A suggest that lengthscale of magnetic fields may be smaller then typical bubble size, this may suggest that other mechanisms, such as viscosity, may be responsible for stabilizing the bubbles. However, since Faraday rotation observations of radio lobes do not constrain large scale ICM fields well if they are aligned with the bubble surface, the draping case may be a viable alternative solution to the problem. A generic feature found in our simulations is the formation of magnetic wakes where fields are ordered and amplified. We suggest that this effect could prevent evaporation by thermal conduction of cold Halpha filaments observed in the Perseus cluster.

The physics of the intracluster medium, in particular the values for the thermal conductivity and the viscosity are largely unknown and subject to an ongoing debate. Here, we study the effect of viscosity on the thermal state of the intracluster medium using three-dimensional cosmological simulations of structure formation. It is shown that viscosity, provided it is not too far off from the unmagnetised Spitzer value, has a significant effect on cluster profiles. In particular, it aids in heating the cool cores of clusters. The central cooling time of the most massive clusters in our simulation is increased by more than an order of magnitude. In large clusters, viscous heating may help to establish an entropy floor and to prevent a cooling catastrophe.

The intergalactic magnetic field within filaments should be less polluted by magnetised outflows from active galaxies than magnetic fields in clusters. Therefore, filaments may be a better laboratory to study magnetic field amplification by structure formation than galaxy clusters which typically host many more active galaxies. We present highly resolved cosmological AMR simulations of magnetic fields in the cosmos and make predictions about the evolution and structure of magnetic fields in filaments. Comparing our results to observational evidence for magnetic fields in filaments suggests that amplification of seed fields by gravitational collapse is not sufficient to produce IGM fields. Finally, implications for cosmic ray transport are discussed.

(abridged) We investigate in detail the role of active galactic nuclei on the physical state of the gas in galaxy groups and clusters, and the implications for anisotropy in the CMB from Sunyaev-Zeldovich effect. We include the effect of thermal conduction, and find that the resulting profiles of temperature and entropy are consistent with observations. Unlike previously proposed models, our model predicts that isentropic cores are not an inevitable consequence of preheating. The model also reproduces the observational trend for the density profiles to flatten in lower mass systems. We deduce the energy E_agn required to explain the entropy observations as a function of mass of groups and clusters M_cl and show that E_agn is proportional to M_cl^alpha with alpha~1.5. We demonstrate that the entropy measurements, in conjunction with our model, can be translated into constraints on the cluster--black hole mass relation. The inferred relation is nonlinear and has the form M_bh\propto M_cl^alpha. This scaling is an analog and extension of a similar relation between the black hole mass and the galactic halo mass that holds on smaller scales. We show that the central decrement of the CMB temperature is reduced due to the enhanced entropy of the ICM, and that the decrement predicted from the plausible range of energy input from the AGN is consistent with available data of SZ decrement. We show that AGN heating, combined with the observational constraints on entropy, leads to suppression of higher multipole moments in the angular power spectrum and we find that this effect is stronger than previously thought.

Recent X-ray observations reveal growing evidence for heating by active galactic nuclei (AGN) in clusters and groups of galaxies. AGN outflows play a crucial role in explaining the riddle of cooling flows and the entropy problem in clusters. Here we study the effect of AGN on the intra-cluster medium in a cosmological simulation using the adaptive mesh refinement FLASH code. We pay particular attention to the effects of conductivity and viscosity on the dissipation of weak shocks generated by the AGN activity in a realistic galaxy cluster. Our 3D simulations demonstrate that both viscous and conductive dissipation play an important role in distributing the mechanical energy injected by the AGN, offsetting radiative cooling and injecting entropy to the gas. These processes are important even when the transport coefficients are at a level of 10% of the Spitzer value. Provided that both conductivity and viscosity are suppressed by a comparable amount, conductive dissipation is likely to dominate over viscous dissipation. Nevertheless, viscous effects may still affect the dynamics of the gas and contribute a significant amount of dissipation compared to radiative cooling. We also present synthetic Chandra observations. We show that the simulated buoyant bubbles inflated by the AGN, and weak shocks associated with them, are detectable with the Chandra observatory.

The thermal broadening of emission lines of heavy ions is small enough such that Doppler shifts due to bulk motions may be detected with the next generation of X-ray observatories. This opens up the possibility to study gas velocities in the intra-cluster medium. Here we study the effect of bulk motions induced by a central active galactic nucleus (AGN) on the emission lines around the FeXXV complex. We have modelled the evolution of AGN-induced bubbles in a realistic cosmological framework and studied the resulting FeXXV line profiles. We found that in clusters with AGN feedback, motions induced by the inflation of bubbles and their buoyant rise lead to distinct features in the iron emission lines that are detectable with a spectral resolution of about 10 eV. These observations will help to determine the mechanical energy that resides in the bubbles and thereby the kinetic luminosity of the AGN.

We present three-dimensional simulations of viscous dissipation of AGN induced gas motions and waves in clusters of galaxies. These simulations are motivated by recent detections of ripples in the Perseus and Virgo clusters. Although the sound waves generated by buoyant bubbles decay with distance from the cluster center, we show that these waves can contribute substantially to offsetting the radiative cooling at distances significantly exceeding the bubble size. The energy flux of the waves declines more steeply with radius than the inverse-square law predicted by energy conservation, implying that dissipation plays an important role in tapping the wave energy. We show that such dispersing sound waves/weak shocks are detectable as ripples on unsharp-masked X-ray cluster maps, and point out that the interfaces between the intracluster medium and old bubbles are also clearly detectable in unsharp-masked X-ray maps. This opens up the possibility of detecting fossil bubbles that are difficult to detect in radio emission. This mode of heating is consistent with other observational constraints, such as the presence of cool rims around the bubbles and the absence of strong shocks. Thus, the mechanism offers a way of heating clusters in a spatially distributed and gentle fashion. We also discuss the energy transfer between the central AGN and the surrounding medium. In our numerical experiments, we find that roughly 65 per cent of the energy injected by the AGN is transferred to the intracluster medium and approximately 25 percent of the injected energy is dissipated by viscous effects and contributes to heating of the gas. The overall transfer of heat from the AGN to the gas is comparable to the radiative cooling losses. The simulations were performed with the FLASH adaptive mesh refinement code.

Recent X-ray observations of clusters of galaxies have shown that the entropy of the intracluster medium (ICM), even at radii as large as half the virial radius, is higher than that expected from gravitational processes alone. This is thought to be the result of nongravitational processes influencing the physical state of the ICM. In this paper, we investigate whether heating by a central AGN can explain the distribution of excess entropy as a function of radius. The AGN is assumed to inject buoyant bubbles into the ICM, which heat the ambient medium by doing pdV work as they rise and expand. Several authors have suggested that this "effervescent heating" mechanism could allow the central regions of clusters to avoid the ``cooling catastrophe''. Here we study the effect of effervescent heating at large radii. Our calculations show that such a heating mechanism is able to solve the entropy problem. The only free parameters of the model are the time-averaged luminosity and the AGN lifetime. The results are mainly sensitive to the total energy injected into the cluster. Our model predicts that the total energy injected by AGN should be roughly proportional to the cluster mass. The expected correlation is consistent with a linear relation between the mass of the central black hole(s) and the mass of the cluster, which is reminiscent of the Magorrian relation between the black hole and bulge mass.

We simulate the effects of viscous dissipation of waves that are generated by AGN activity in clusters of galaxies. We demonstrate that the amount of viscous heating associated with the dissipation of these waves can offset radiative cooling rates in cooling flow clusters of galaxies. This heating mechanism leads to spatially distributed and approximately symmetrical dissipation. The heating waves reach a given distance from the cluster center on a timescale shorter than the cooling time. This means that this heating mechanism has the potential of quenching cooling flows in a quasi-stable fashion. Moreover, the heating is gentle as no strong shocks are present in the simulations. We first investigated whether a single continuous episode of AGN activity can lead to adequate dissipation to balance cooling rates. These simulations demonstrated that, whereas secondary waves generated by the interaction of the rising bubble with the intracluster medium are clearly present, viscous heating associated with the dissipation of these waves is insufficient to balance radiative cooling. It is only when the central source is intermittent that the viscous dissipation of waves associated with subsequent episodes of activity can offset cooling. This suggests that the ripples observed in the Perseus cluster can be interpreted as being due to the AGN duty cycle, i.e., they trace AGN activity history. The simulations were performed using the PPM adaptive mesh refinement code FLASH in two dimensions.

Jets and winds are significant channels for energy loss from accreting black holes. These outflows mechanically heat their surroundings, through shocks as well as gentler forms of heating. We discuss recent efforts to understand the nature and distribution of mechanical heating by central AGNs in clusters of galaxies, using numerical simulations and analytic models. Specifically, we will discuss whether the relatively gentle `effervescent heating' mechanism can compensate for radiative losses in the central regions of clusters, and account for the excess entropy observed at larger radii.

Radiation dominated accretion disks are likely to be subject to the ``photon bubble'' instability, which may lead to strong density inhomogeneities on scales much shorter than the disk scale height. Such disks -- and magnetized, radiation-dominated atmospheres in general -- could radiate well above the Eddington limit without being disrupted. When density contrasts become large over distances of order the photon mean free path, radiative transfer cannot be described adequately using either the standard diffusion approximation or existing prescriptions for flux-limited diffusion. Using analytical and Monte Carlo techniques, we consider the effects of strong density gradients deep within radiation- and scattering-dominated atmospheres. We find that radiation viscosity -- i.e., the off-diagonal elements of the radiation stress tensor -- has an important effect on radiative transfer under such conditions. We compare analytical and numerical results in the specific case of a plane-parallel density wave structure and calculate Eddington enhancement factors due to the porosity of the atmosphere. Our results can be applied to the study of dynamical coupling between radiation forces and density inhomogeneities in radiation dominated accretion disks in two or three dimensions.

I discuss the transfer of polarized synchrotron radiation in relativistic jets. I argue that the main mechanism responsible for the circular polarization properties of compact synchrotron sources is likely to be Faraday conversion and that, contrary to common expectation, a significant rate of Faraday rotation does not necessarily imply strong depolarization. The long-term persistence of the sign of circular polarization, observed in some sources, is most likely due to a small net magnetic flux generated in the central engine, carried along the jet axis and superimposed on a highly turbulent magnetic field. I show that the mean levels of circular and linear polarizations depend on the number of field reversals along the line of sight and that the gradient in Faraday rotation across turbulent regions can lead to "correlation depolarization''. The model is potentially applicable to a wide range of synchrotron sources. In particular, I demonstrate how the model can naturally explain the excess of circular over linear polarization in the Galactic Center (Sgr A*) and the low-luminosity AGN M81*.

There is mounting observational evidence from Chandra for strong interaction between keV gas and AGN in cooling flows. It is now widely accepted that the temperatures of cluster cores are maintained at a level of 1 keV and that the mass deposition rates are lower than earlier ROSAT/Einstein values. Recent theoretical results suggest that thermal conduction can be very efficient even in magnetized plasmas. Motivated by these discoveries, we consider a ``double heating model'' which incorporates the effects of simultaneous heating by both the central AGN and thermal conduction from the hot outer layers of clusters. Using hydrodynamical simulations, we demonstrate that there exists a family of solutions that does not suffer from the cooling catastrophe. In these cases, clusters relax to a stable final state, which is characterized by minimum temperatures of order 1 keV and density and temperature profiles consistent with observations. Moreover, the accretion rates are much reduced, thereby reducing the need for excessive mass deposition rates required by the standard cooling flow models.

The transfer of polarized radiation in stochastic synchrotron sources is explored by means of analytic treatment and Monte Carlo simulations. We argue that the main mechanism responsible for the circular polarization properties of compact synchrotron sources is likely to be Faraday conversion and that, contrary to common expectation, a significant rate of Faraday rotation does not necessarily imply strong depolarization. The long-term persistence of the sign of circular polarization, observed in many sources, is most likely due to a small net magnetic flux generated in the central engine, carried along the jet axis and superimposed on a highly turbulent magnetic field. We show that the mean levels of circular and linear polarizations depend on the number of field reversals along the line of sight and that the gradient in Faraday rotation across turbulent regions can lead to ``correlation depolarization''. Our model is potentially applicable to a wide range of synchrotron sources. In particular, we demonstrate how our model can naturally explain the excess of circular over linear polarization in the Galactic Center and the nearby spiral galaxy M81 and discuss its application to the quasar 3C 279, the intraday variable blazar PKS 1519-273 and the X-ray binary SS 433.

The fluorescent iron K$\alpha$ emission line profile provides an excellent probe of the innermost regions of active galactic nuclei. Fe XXV and Fe XXVI in diffuse plasma above the accretion disc can affect the X-ray spectrum by iron K$\alpha$ resonant absorption. This in turn can influence the interpretation of the data and the estimation of the accretion disc and black hole parameters. We embark on a fully relativistic computation of this effect and calculate the iron line profile in the framework of a specific model in which rotating, highly ionized and resonantly-absorbing plasma occurs close to the black hole. This can explain the features seen in the iron K$\alpha$ line profile recently obtained by Nandra et al. (1999) for the Seyfert 1 galaxy NGC 3516. We show that the redshift of this feature can be mainly gravitational in origin and accounted for without the need to invoke fast accretion of matter onto the black hole. New X-ray satellites such as XMM, ASTRO-E and Chandra provide excellent opportunities to test the model against high quality observational data.

The broad X-ray iron line, detected in many active galactic nuclei, is likely to be produced by fluorescence from the X-ray illuminated central parts of an accretion disc close to a supermassive black hole. The time-averaged shape of the line can be explained most naturally by a combination of special and general relativistic effects. Such line profiles contain information about the black hole spin and the accretion disc as well as the geometry of the emitting region and may help to test general relativity in the strong gravity regime. In this paper we embark on the computation of the temporal response of the line to the illuminating flux. Previous studies concentrated on the calculation of reverberation signatures from static sources illuminating the disc. In this paper we focus on the more physically justified case of flares located above the accretion disc and corotating with it. We compute the time dependent iron line taking into account all general relativistic effects and show that its shape is of very complex nature, and also present light curves accompanying the iron line variability. We suggest that future X-ray satellites like XMM or Constellation-X may be capable of detecting features present in the computed reverberation maps.

We report results from an 18-day ROSAT HRI monitoring campaign on the ultrasoft Narrow-Line Seyfert 1 (NLS1) class quasar PHL 1092. This luminous, radio-quiet quasar showed strong X-ray variability in a short ROSAT PSPC observation, and ROSAT HRI monitoring of the similar object IRAS 13224-3809 revealed extreme variability on intermediate timescales. We wanted to determine whether remarkable X-ray variability persistently occurs in PHL 1092, and we also wanted to search for outstanding variability events that constrain emission processes. Given the large luminosity of PHL 1092 (about 5x10^{45} erg/s in the HRI band), we detect extremely rapid and large-amplitude X-ray variability throughout our monitoring. The maximum observed variability amplitude is a factor of about 14, and in the most rapid variability event the HRI count rate increases by a factor of about 3.8 in a rest-frame time interval of < 3580 s. The most rapid event has a rate change of luminosity of > 1.3x10^{42} erg/s^2, making it the most extreme such event we are aware of from a radio-quiet quasar. Standard `radiative efficiency limit' arguments imply a radiative efficiency larger than can be achieved by accretion onto a Kerr black hole rotating at the maximum plausible rate, although we point out that such arguments depend upon the geometry of initial radiation release. Relativistic motions of the X-ray source are probably causing the radiative efficiency limit to break down; such relativistic motions have also been inferred in the similar NLS1-class quasar PKS 0558--504.

Gamma-ray bursters emit a small fraction of their flux in X rays, and because X-ray detectors are often very sensitive they may probe the gamma-ray burst universe more deeply than the current best gamma-ray instruments. On the reasonable assumptions that spectra of bursts observed by BATSE may be used to predict the X-ray fluxes of gamma-ray bursts, and that any corona of bursts around M31 is similar to the one around the Milky Way, we predict the rate at which the wide field cameras on board BeppoSAX should detect bursts from the Milky Way and M31. These rates are such that a one-month observation of M31 would have to either detect bursts from M31 or exclude most galactic models of gamma-ray bursts. (It is shown how the remainder can be dealt with.) Therefore such an observation would settle the long-standing dispute over their location.