In the universe's most massive galaxies, kinetic feedback from a central supermassive black hole appears to limit star formation. Abundant circumstantial evidence suggests that accumulation of cold gas near the central black hole strongly boosts the feedback output, keeping the ambient medium in a state marginally unstable to condensation and formation of cold gas clouds. However, the ability of that mechanism to self-regulate may depend on numerous environmental factors, including the depth of the potential well and the pressure of the surrounding circumgalactic medium (CGM). Here we present a suite of numerical simulations that explores the dependence of cold-fuelled bipolar kinetic feedback on those environmental factors. Halo mass in this simulation suite ranges from $2 \times 10^{12} \, M_\odot$ to $8 \times 10^{14} \, M_\odot$. We include the spatially extended mass and energy input from the massive galaxy's old stellar population, which is capable of sweeping gas out of the galaxy and away from the central black hole if the confining CGM pressure is sufficiently low. Our simulations show that this feedback mechanism is tightly self-regulating in a massive galaxy with a deep central potential and low CGM pressure, permitting only small amounts of multiphase gas to accumulate and allowing almost no star formation. In a massive galaxy of similar mass but a shallower central potential and greater CGM pressure the same feedback mechanism is more episodic, producing extended multiphase gas and occasionally allowing small rates of star formation ($\sim 0.1 \, M_\odot \, {\rm yr}^{-1}$). At the low-mass end of the explored range the mechanism becomes implausibly explosive, perhaps because the ambient gas initially has no angular momentum, which would have reduced the amount of condensed gas capable of fueling feedback.

Star formation in the universe's most massive galaxies proceeds furiously early in time but then nearly ceases. Plenty of hot gas remains available but does not cool and condense into star-forming clouds. Active galactic nuclei (AGN) release enough energy to inhibit cooling of the hot gas, but energetic arguments alone do not explain why quenching of star formation is most effective in high-mass galaxies. In fact, optical observations show that quenching is more closely related to a galaxy's central stellar velocity dispersion ($\sigma_v$) than to any other characteristic. Here, we show that high $\sigma_v$ is critical to quenching because a deep central potential well maximizes the efficacy of AGN feedback. In order to remain quenched, a galaxy must continually sweep out the gas ejected from its aging stars. Supernova heating can accomplish this task as long as the AGN sufficiently reduces the gas pressure of the surrounding circumgalactic medium (CGM). We find that CGM pressure acts as the control knob on a valve that regulates AGN feedback and suggest that feedback power self-adjusts so that it suffices to lift the CGM out of the galaxy's potential well. Supernova heating then drives a galactic outflow that remains homogeneous if $\sigma_v \gtrsim 240 \, {\rm km \, s^{-1}}$. AGN feedback can effectively quench galaxies with a comparable velocity dispersion, but feedback in galaxies with a much lower velocity dispersion tends to result in convective circulation and accumulation of multiphase gas within the galaxy.

We study ram pressure stripping in simulated Milky Way-like halos at z>=2 from the Figuring Out Gas & Galaxies In Enzo (FOGGIE) project. These simulations reach exquisite resolution in their circumgalactic medium (CGM) gas owing to FOGGIE's novel refinement scheme. The CGM of each halo spans a wide dynamic range in density and velocity over its volume---roughly 6 dex and 1000 km/s, respectively---translating into a 5 dex range in ram pressure imparted to interacting satellites. The ram pressure profiles of the simulated CGM are highly stochastic, owing to kpc-scale variations of the density and velocity fields of the CGM gas. As a result, the efficacy of ram pressure stripping depends strongly on the specific path a satellite takes through the CGM. The ram-pressure history of a single satellite is generally unpredictable and not well correlated with its approach vector with respect to the host galaxy. The cumulative impact of ram pressure on the simulated satellites is dominated by only a few short strong impulses---on average, 90% of the total surface momentum gained through ram pressure is imparted in 20% or less of the total orbital time. These results reveal an erratic mode of ram pressure stripping in Milky-Way like halos at high redshift---one that is not captured by a smooth spherically-averaged model of the circumgalactic medium.

Recent observations show that the star formation rate (SFR) in the {\it Phoenix} cluster's central galaxy is $\sim 500$ M$_\odot$ yr$^{-1}$. Even though {\it Phoenix} is a massive cluster ($M_{200} \approx 2.0\times 10^{15}$ M$_\odot$; $z\approx 0.6$) such a high central SFR is not expected in a scenario in which feedback from an active galactic nucleus (AGN) maintains the intracluster medium (ICM) in a state of rough thermal balance. It has been argued that either AGN feedback saturates in very massive clusters or the central supermassive black hole (SMBH) is too small to produce enough kinetic feedback and hence is unable to quench the catastrophic cooling. In this work, we present an alternate scenario wherein intense short-lived cooling and star formation phases followed by strong AGN outbursts are part of the AGN feedback loop. Using results from a 3D hydrodynamic simulation of a standard cool-core cluster ($M_{200}\sim 7\times10^{14}$ M$_\odot$; $z=0$), scaled to account for differences in mass and redshift, we argue that {\it Phoenix} is at the end of a cooling phase in which an AGN outburst has begun but has not yet arrested core cooling. This state of high cooling rate and star formation is expected to last for $\lesssim$ 100 Myr in {\it Phoenix}.

During the last decade, numerous and varied observations, along with increasingly sophisticated numerical simulations, have awakened astronomers to the central role the circumgalactic medium (CGM) plays in regulating galaxy evolution. It contains the majority of the baryonic matter associated with a galaxy, along with most of the metals, and must continually replenish the star forming gas in galaxies that continue to sustain star formation. And while the CGM is complex, containing gas ranging over orders of magnitude in temperature and density, a simple emergent property may be governing its structure and role. Observations increasingly suggest that the ambient CGM pressure cannot exceed the limit at which cold clouds start to condense out and precipitate toward the center of the potential well. If feedback fueled by those clouds then heats the CGM and causes it to expand, the pressure will drop and the "rain" will diminish. Such a feedback loop tends to suspend the CGM at the threshold pressure for precipitation. The coming decade will offer many opportunities to test this potentially fundamental principle of galaxy evolution.

Recent analysis shows that it is important to explicitly include the gravitational potential of the central brightest central galaxy (BCG) to infer the acceleration due to gravity ($g$) and the free-fall time ($t_{\rm ff} \equiv [2r/g]^{1/2}$) in cool cluster cores. Accurately measuring $t_{\rm ff}$ is crucial because according to numerical simulations cold gas condensation and strong feedback occur in cluster cores with min($t_{\rm cool}/t_{\rm ff}$) below a threshold value close to 10. Recent observations which include the BCG gravity show that the observed threshold in min($t_{\rm cool}/t_{\rm ff}$) lies at a somewhat higher value, close to 10-30; there are only a few clusters in which this ratio falls much below 10. In this paper we compare numerical simulations of feedback AGN (Active Galactic Nuclei) jets interacting with the intracluster medium (ICM), with and without a BCG potential. We find that, for a fixed feedback efficiency, the presence of a BCG does not significantly affect the temperature but increases (decreases) the core density (entropy) on average. Most importantly, min($t_{\rm cool}/t_{\rm ff}$) is only affected slightly by the inclusion of the BCG gravity. Also notable is that the lowest value of min($t_{\rm cool}/t_{\rm ff}$) in the NFW+BCG runs are about twice larger than in the NFW runs. We also look at the role of depletion of cold gas due to star formation and show that it only affects the rotationally dominant component (torus), while the radially dominant component (which regulates the feedback cycle) remains largely unaffected. The distribution of metals due to AGN jets in our simulations is predominantly along the jet direction and the radial spread of metals is less. We also show that the turbulence in cool core clusters is weak, consistent with recent Hitomi results on Perseus cluster.

An important ingredient in numerical modelling of high temperature magnetised astrophysical plasmas is the anisotropic transport of heat along magnetic field lines from higher to lower temperatures.Magnetohydrodynamics (MHD) typically involves solving the hyperbolic set of conservation equations along with the induction equation. Incorporating anisotropic thermal conduction requires to also treat parabolic terms arising from the diffusion operator. An explicit treatment of parabolic terms will considerably reduce the simulation time step due to its dependence on the square of the grid resolution ($\Delta x$) for stability. Although an implicit scheme relaxes the constraint on stability, it is difficult to distribute efficiently on a parallel architecture. Treating parabolic terms with accelerated super-time stepping (STS) methods has been discussed in literature but these methods suffer from poor accuracy (first order in time) and also have difficult-to-choose tuneable stability parameters. In this work we highlight a second order (in time) Runge Kutta Legendre (RKL) scheme (first described by Meyer et. al. 2012) that is robust, fast and accurate in treating parabolic terms alongside the hyperbolic conversation laws. We demonstrate its superiority over the first order super time stepping schemes with standard tests and astrophysical applications. We also show that explicit conduction is particularly robust in handling saturated thermal conduction. Parallel scaling of explicit conduction using RKL scheme is demonstrated up to more than $10^4$ processors.

Several arguments suggest that stochastic condensation of cold gas and its accretion onto the central supermassive black hole (SMBH) is essential for active galactic nuclei (AGN) feedback to work in the most massive galaxies that lie at the centres of galaxy clusters. Our 3-D hydrodynamic AGN jet-ICM (intracluster medium) simulations, looking at the detailed angular momentum distribution of cold gas and its time variability for the first time, show that the angular momentum of the cold gas crossing $\lesssim 1$ kpc is essentially isotropic. With almost equal mass in clockwise and counter-clockwise orientations, we expect a cancellation of angular momentum on roughly the dynamical time. This means that a compact accretion flow with a short viscous time ought to form, through which enough accretion power can be channeled into jet mechanical energy sufficiently quickly to prevent a cooling flow. The inherent stochasticity, expected in feedback cycles driven by cold gas condensation, gives rise to a large variation in the cold gas mass at the centres of galaxy clusters, for similar cluster and SMBH masses, in agreement with the observations. Such correlations are expected to be much tighter for the smoother hot/Bondi accretion. The weak correlation between cavity power and Bondi power obtained from our simulations also match observations.

Using high-resolution 3-D and 2-D (axisymmetric) hydrodynamic simulations in spherical geometry, we study the evolution of cool cluster cores heated by feedback-driven bipolar active galactic nuclei (AGN) jets. Condensation of cold gas, and the consequent enhanced accretion, is required for AGN feedback to balance radiative cooling with reasonable efficiencies, and to match the observed cool core properties. A feedback efficiency (mechanical luminosity $\approx \epsilon \dot{M}_{\rm acc} c^2$; where $\dot{M}_{\rm acc}$ is the mass accretion rate at 1 kpc) as small as $5 \times 10^{-5}$ is sufficient to reduce the cooling/accretion rate by $\sim 10$ compared to a pure cooling flow. This value is smaller compared to the ones considered earlier, and is consistent with the jet efficiency and the fact that only a small fraction of gas at 1 kpc is accreted on to the supermassive black hole (SMBH). We find hysteresis cycles in all our simulations with cold mode feedback: {\em condensation} of cold gas when the ratio of the cooling-time to the free-fall time ($t_{\rm cool}/t_{\rm ff}$) is $\lesssim 10$ leads to a sudden enhancement in the accretion rate; a large accretion rate causes strong jets and {\em overheating} of the hot ICM such that $t_{\rm cool}/t_{\rm ff} > 10$; further condensation of cold gas is suppressed and the accretion rate falls, leading to slow cooling of the core and condensation of cold gas, restarting the cycle. Therefore, there is a spread in core properties, such as the jet power, accretion rate, for the same value of core entropy or $t_{\rm cool}/t_{\rm ff}$. A fewer number of cycles are observed for higher efficiencies and for lower mass halos because the core is overheated to a longer cooling time. The 3-D simulations show the formation of a few-kpc scale, rotationally-supported, massive ($\sim 10^{11} M_\odot$) cold gas torus. (abstract abridged)

We report on the results of infrared spectroscopic mapping observations carried out in the nuclear region of Centaurus A (NGC5128). The 500 pc bipolar dust shell discovered by Quillen et al.(2006) is even more clearly seen in the 11.3 micron dust emission feature than previous broad band imaging. The pure rotational lines of molecular hydrogen other than the S(0) line are detected above the dusty disk and associated with the oval dust shell. The molecular hydrogen transitions indicate the presence of warm gas at temperatures 250--720K. The ratio of the dust emission features at 7.7 and 11.3 micron and the ratio of the [NeII](12.8) and 11.3 dust emission feature are lower in the 500 pc dust shell than in the star forming disk. The clearer shell morphology at 11.3 micron, warm molecular hydrogen emission in the shell, and variation in line ratios in the shell compared to those in the disk, confirm spectroscopically that this shell is a separate coherent entity and is unlikely to be a chance superposition of dust filaments. The physical conditions in the shell are most similar to Galactic supernova remnants where blast waves encounter molecular clouds. The lines requiring the highest level of ionization, [NeV] and [OIV], are detected 20--25 arcsec north-east and south-west of the nucleus and at position angles near the radio jet axis. Fine structure line ratios and limits from this region suggest that the medium is low density and illuminated by a hard radiation field at low ionization parameter. These higher S molecular hydrogen pure rotational transitions are also particularly bright in the same region as the [OIV] and [NeV] emission. This suggests that the gas associated with the dust shell has been excited near the jet axis and is part of an ionization cone.