Strong scaling relations between host galaxy properties (such as stellar mass, bulge mass, luminosity, effective radius etc) and their nuclear supermassive black hole's mass point towards a close co-evolution. In this work, we first review previous efforts supporting the fundamental importance of the relation between supermassive black hole mass and stellar velocity dispersion ($M_{\rm BH}$-$\sigma_{\rm e}$). We then present further original work supporting this claim via analysis of residuals and principal component analysis applied to some among the latest compilations of local galaxy samples with dynamically measured supermassive black hole masses. We conclude with a review of the main physical scenarios in favour of the existence of a $M_{\rm BH}$-$\sigma_{\rm e}$ relation, with a focus on momentum-driven outflows.

We present a semi-analytic model of isolated dwarf galaxy evolution and use it to study the build-up of observed correlations between dwarf galaxy properties. We analyse the evolution using models with averaged and individual halo mass assembly histories in order to determine the importance of stochasticity on the present-day properties of dwarf galaxies. The model has a few free parameters, but when these are calibrated using the halo mass - stellar mass and stellar mass-metallicity relations, the results agree with other observed dwarf galaxy properties remarkably well. Redshift evolution shows that even isolated galaxies change significantly over the Hubble time and that 'fossil dwarf galaxies' with properties equivalent to those of high-redshift analogues should be extremely rare, or non-existent, in the Local Universe. A break in most galaxy property correlations develops over time, at a stellar mass $M_* \simeq 10^7 M_\odot$. It is caused predominantly by the ionizing background radiation and can therefore in principle be used to constrain the properties of reionization.

There possibly was an AGN episode in the Galactic Centre about 6 Myr ago, powerful enough to produce the Fermi bubbles. We present numerical simulations of a possible scenario giving rise to an activity episode: a collision between a central gas ring surrounding the supermassive black hole (SMBH) and an infalling molecular cloud. We investigate different initial collision angles between the cloud and the ring. We follow the hydrodynamical evolution of the system following the collision using Gadget-3 hybrid N-body/SPH code and calculate the feeding rate of the SMBH accretion disc. This rate is then used as an input for a 1D thin $\alpha$-disc model in order to calculate the AGN luminosity. By varying the disc feeding radii we determine the limiting values for possible AGN accretion disc luminosity. Small angle collisions do not result in significant mass transport to the centre of the system, while models with highest collision angles transport close to $40\%$ of the initial matter to the accretion disc. Even with ring and cloud masses equal to $10^4 \, M_\odot$, which is the lower limit of present-day mass of the Circumnuclear ring in the Galactic Centre, the energy released over an interval of 1.5 Myr can produce $\sim 10\%$ of that required to inflate the Fermi bubbles. If the gas ring in the Galactic Centre 6 Myr ago had a mass of at least $10^5 \, M_\odot$, our proposed scenario can explain the formation of the Fermi bubbles. We estimate that such high-impact collisions might occur once every $\sim 10^8$ yr in our Galaxy.

Active galactic nuclei (AGN) probably control the growth of their host galaxies via feedback in the form of wide-angle wind-driven outflows. These establish the observed correlations between supermassive black hole (SMBH) masses and host galaxy properties, e.g. the spheroid velocity dispersion $\sigma$. In this paper we consider the growth of the SMBH once it starts driving a large-scale outflow through the galaxy. To clear the gas and ultimately terminate further growth of both the SMBH and the host galaxy, the black hole must continue to grow its mass significantly, by up to a factor of a few, after reaching this point. The mass increment $\Delta M_{\rm BH}$ depends sensitively on both galaxy size and SMBH spin. The galaxy size dependence leads to $\Delta M_{\rm BH} \propto \sigma^5$ and a steepening of the $M-\sigma$ relation beyond the analytically calculated $M \propto \sigma^4$, in agreement with observation. Slowly--spinning black holes are much less efficient in producing feedback, so at any given $\sigma$ the slowest--spinning black holes should be the most massive. Current observational constraints are consistent with this picture, but insufficient to test it properly; however, this should change with upcoming surveys.

Supermassive black holes (SMBHs) are found in the centres of most galaxies. Their masses, and hence their gravitational potentials, are negligible compared with those of the host galaxy. However, several strong correlations between SMBH masses and host galaxy properties have been observed, notably the $M-\sigma$ relation connecting the SMBH mass to the characteristic velocity of stars in the galaxy. The existence of these correlations implies that the SMBH influences the evolution of its host galaxy. In this review, we present the most promising physical model of this influence, known as the Active galactic nucleus (AGN) wind feedback model. Winds launched from the accretion disc around the SMBH can drive powerful outflows, provided that the SMBH is massive enough - this condition establishes the $M-\sigma$ relation. Outflows can have a profound influence on the evolution of the host galaxy, by compressing its gas and driving it out, affecting the star formation rate. We present the current status of the model and the observational evidence for it, as well as the directions of future research.

Warm absorbers are found in many AGN and consist of clouds moving at moderate radial velocities, showing complex ionization structures and having moderate to large column densities. Using 1D numerical calculations, we confirm earlier suggestions that the energy released by an AGN pushes the surrounding gas outward in a bubble until this reaches transparency. Typical AGN episode durations of $5\times 10^4$ yr supply enough energy for this, except in very gas-rich and/or very compact galaxies, such as those in the early Universe. In those galaxies, the AGN might remain hidden for many periods of activity, hiding the black hole growth. The typical radii of $0.1-1$ kpc, velocities of $100-1000$ km s$^{-1}$ and resulting optical depths are consistent with observations of warm absorbers. The resulting structure is a natural outcome of outflows driven by AGN buried in an optically thick gas envelope, and has a total mass comparable to the final $M -\sigma$ mass the central supermassive black hole will eventually reach.These results suggest that AGN can feed very efficiently by agitating this surrounding dense material. This may not be easy to observe, as this gas is Compton thick along many sightlines. The infall may produce episodic star formation in the centre, building up nuclear star clusters simultaneously with the growth of the central black hole.

In recent years, numerous active galactic nuclei have been discovered in ever smaller galaxies, questioning the paradigm that dwarf galaxies do not harbour central massive black holes. Even if such black holes exist, feeding them by gas streams is difficult, since star formation should be more efficient than AGN feeding in dwarf galaxies. In this paper, I investigate the possibility that tidal disruptions of stars are responsible for the observed AGN in dwarf galaxies. I show that the expected duty cycles of TDE-powered AGN, $f_{\rm AGN} \geq 0.5\%$, are consistent with observed AGN fractions assuming that the occupation fraction in dwarf galaxies is close to unity. Furthermore, I calculate the properties of outflows driven by TDE-powered AGN under idealised conditions and find that they might have noticeable effects on the host galaxies. Outflows themselves might not be detectable, except in gas-poor galaxies, where they can accelerate to $v_{\rm out} > 100$~km/s, but increased gas turbulence, more diffuse density profile and lower star formation efficiency can be discovered and used to constrain the black hole occupancy fraction and more nuanced effects on dwarf galaxy evolution. If massive black holes form from seeds that are much more massive than stellar black holes, then their outflows should be easily detectable; this result, aided by observations of high-redshift dwarf galaxies, provides a potential way of determining seed masses of black holes.

Galaxy evolution is affected by competing feedback processes. Stellar feedback dominates in low-mass galaxies, while AGN feedback predominantly affects massive ones. Recent observational results reveal the dependence of black hole accretion rate (BHAR) and star formation rate (SFR) on galaxy stellar mass, and give information on the galaxy mass at which the changeover between dominant feedback mechanisms occurs. I use this information to derive an empirical estimate of the coupling efficiency, $f_{\rm AGN}$, between AGN luminous energy output and AGN-driven galactic outflows, and the momentum loading factor $f_{\rm p,AGN}$ between the momentum of AGN radiation field and the outflow. The results are independent of any particular model of AGN feedback and show that AGN feedback must be very efficient and/or have very large momentum loading in order to explain current observations. I discuss possible ways of reaching the required efficiency and loading factor, and the selection effects that might result in only weak outflows being observed, while the most powerful ones may be generally obscured. There are significant uncertainties involved in the derivation of the result; I suggest ways of reducing them. In the near future, better estimates of coupling efficiency can help distinguish among AGN feedback models, investigate the redshift evolution and mass dependence of feedback efficiency.

Powerful outflows with a broad range of properties (such as velocity, ionization, radial scale and mass loss rate) represent a key feature of active galactic nuclei (AGN), even more so since they have been simultaneously revealed also in individual objects. Here we revisit in a simple analytical framework the recent remarkable cases of two ultraluminous infrared quasars, IRAS F11119+3257 and Mrk 231, which allow us to investigate the physical connection between multi-phase AGN outflows across the ladder of distance from the central supermassive black hole (SMBH). We argue that any major deviations from the standard outflow propagation models might encode unique information on the past SMBH accretion history, and briefly discuss how this could help address some controversial aspects of the current picture of AGN feedback.

Sgr A$^*$ is the super massive black hole residing in the centre of the Milky Way. There is plenty of observational evidence that a massive gas cloud fell into the central parsec of the Milky Way $\sim 6$ million years ago, triggering formation of a disc of young stars and activating Sgr A$^*$. In addition to the disc, there is an unexplained population of young stars on randomly oriented orbits. Here we hypothesize that these young stars were formed by fragmentation of a massive quasi-spherical gas shell driven out from Sgr A$^*$ potential well by an energetic outflow. To account for the properties of the observed stars, the shell must be more massive than $10^5$ Solar masses, be launched from inside $\sim 0.01$~pc, and the feedback outflow has to be highly super-Eddington albeit for a brief period of time, producing kinetic energy of at least $10^{55}$~erg. The young stars in the central parsec of the Galaxy may be a unique example of stars formed from atomic rather than molecular hydrogen, and forged by extreme pressure of black hole outflows.

Irregular dwarf galaxies of the Local Group have very varied properties and star formation histories. Some of them formed the majority of their stars very late compared to the others. Extreme examples are Leo A and Aquarius which reached the peak of star formation at $z<1$ ( > 6 Gyr after BB). This fact seemingly challenges the LCDM cosmology because the dark matter halos of these galaxies on average should assemble the majority of their masses before z~2 (<3 Gyr after BB). In this work we investigate whether the delayed star formation histories of some irregular dwarf galaxies could be explained purely by the stochasticity of their mass assembly histories coupled with the effect of cosmic reionization. We develop a semi-analytic model to follow the accretion of baryonic matter, star formation and stellar feedback in dark matter halos with present day virial masses 10^9 M_Sun < M < 10^11 M_Sun and with different stochastic growth histories obtained using the PINOCCHIO code based on Lagrangian perturbation theory. We obtain the distributions of observable parameters and the evolution histories for these galaxies. Accretion of baryonic matter is strongly suppressed after the epoch of reionization in some models but they continue to accrete dark matter and eventually reach enough mass for accretion of baryonic matter to begin again. These "reborn" model galaxies show very similar delayed star formation histories to those of Leo A and Aquarius. We find that the stochasticity caused by mass assembly histories is enhanced in systems with virial masses ~10^10 M_Sun because of their sensitivity to the photoionizing intergalactic radiation field after the epoch of reionization. This results in qualitatively different star formation histories in late- and early-forming galaxies and it might explain the peculiar star formation histories of irregular dwarf galaxies such as Leo A and Aquarius.

Massive large-scale AGN outflows are an important element of galaxy evolution, being a way through which the AGN can affect most of the host galaxy. However, outflows evolve on timescales much longer than typical AGN episode durations, therefore most AGN outflows are not observed simultaneously with the AGN episode that inflated them. It is therefore remarkable that rather tight correlations between outflow prop- erties and AGN luminosity exist. In this paper, I show that such correlations can be preserved during the fading phase of the AGN episode, provided that the AGN lumi- nosity evolves as a power law with exponent $\alpha_{\rm d} \sim 1$ at late times. I also show that subsequent AGN episodes that illuminate an ongoing outflow are unlikely to produce outflow momentum or energy rates rising above the observed correlations. However, there may be many difficult-to-detect outflows with momentum and energy rates lower than expected from the current AGN luminosity. Detailed observations of AGN out- flow properties might help constrain the activity histories of typical and/or individual AGN.

AGN outflows can remove large quantities of gas from their host galaxy spheroids, potentially shutting off star formation. On the other hand, they can compress this gas, potentially enhancing or triggering star formation, at least for short periods. We present a set of idealised simulations of AGN outflows affecting turbulent gas spheres, and investigate the effect of the outflow and the AGN radiation field upon gas fragmentation. We show that AGN outflows of sufficient luminosity shut off fragmentation while the nucleus is active, but gas compression results in a burst of fragmentation after the AGN switches off. Self-shielding of gas against the AGN radiation field allows some fragmentation to occur during outbursts, but too much shielding results in a lower overall fragmentation rate. For our idealised simulation setup, there is a critical AGN luminosity which results in the highest fragmentation rate, with outflows being too efficient at removing gas when $L > L_{\rm crit}$ and not efficient enough to compress the gas to high densities otherwise. These results, although preliminary, suggest that the interaction between AGN and star formation in their host galaxies is particularly complex and requires careful study in order to interpret observations correctly.

Observations of massive outflows with detectable central AGN typically find them within radii $\lesssim 10$ kpc. We show that this apparent size restriction is a natural result of AGN driving if this process injects total energy only of order the gas binding energy to the outflow, and the AGN varies over time (`flickers') as suggested in recent work. After the end of all AGN activity the outflow continues to expand to larger radii, powered by the thermal expansion of the remnant shocked AGN wind. We suggest that on average, outflows should be detected further from the nucleus in more massive galaxies. In massive gas--rich galaxies these could be several tens of kpc in radius. We also consider the effect that pressure of such outflows has on a galaxy disc. In moderately gas--rich discs, with gas-to-baryon fraction $< 0.2$, the outflow may induce star formation significant enough to be distinguished from quiescent by an apparently different normalisation of the Kennicutt-Schmidt law. The star formation enhancement is probably stronger in the outskirts of galaxy discs, so coasting outflows might be detected by their effects upon the disc even after the driving AGN has shut off. We compare our results to the recent inference of inside--out quenching of star formation in galaxy discs.

AGN feedback is an important ingredient in galaxy evolution, however its treatment in numerical simulations is necessarily approximate, requiring subgrid prescriptions due to the dynamical range involved in the calculations. We present a suite of SPH simulations designed to showcase the importance of the choice of a particular subgrid prescription for AGN feedback. We concentrate on two approaches to treating wide-angle AGN outflows: thermal feedback, where thermal and kinetic energy is injected into the gas surrounding the SMBH particle, and virtual particle feedback, where energy is carried by tracer particles radially away from the AGN. We show that the latter model produces a far more complex structure around the SMBH, which we argue is a more physically correct outcome. We suggest a simple improvement to the thermal feedback model - injecting the energy into a cone, rather than spherically symmetrically - and show that this markedly improves the agreement between the two prescriptions, without requiring any noticeable increase in the computational cost of the simulation.

Feedback from super-massive black holes (SMBHs) is thought to play a key role in regulating the growth of host galaxies. Cosmological and galaxy formation simulations using smoothed particle hydrodynamics (SPH), which usually use a fixed mass for SPH particles, often employ the same sub-grid Active galactic nuclei (AGN) feedback prescription across a range of resolutions. It is thus important to ask how the impact of the simulated AGN feedback on a galaxy changes when only the numerical resolution (the SPH particle mass) changes. We present a suite of simulations modelling the interaction of an AGN outflow with the ambient turbulent and clumpy interstellar medium (ISM) in the inner part of the host galaxy at a range of mass resolutions. We find that, with other things being equal, degrading the resolution leads to feedback becoming more efficient at clearing out all gas in its path. For the simulations presented here, the difference in the mass of the gas ejected by AGN feedback varies by more than a factor of ten between our highest and lowest resolution simulations. This happens because feedback-resistant high density clumps are washed out at low effective resolutions. We also find that changes in numerical resolution lead to undesirable artifacts in how the AGN feedback affects the AGN immediate environment.

The Central Molecular Zone (CMZ) of the Milky Way shows several peculiar properties: a large star formation rate, some of the most massive young star clusters and molecular clouds in the Galaxy, and a twisted ring morphology in molecular gas. In this paper, I use SPH simulations to show that most of these properties can be explained as due to a recent outburst of AGN activity in Sgr A*, the central supermassive black hole of the Milky Way. In particular, the narrow ring of dense gas, massive gas clouds, young star clusters and an elevated SFR can all be caused by the passage of an AGN outflow through the system, which compresses the gas and triggers fragmentation. Furthermore, I show that the asymmetric distribution of gas, as observed in the CMZ, can be produced by outflow-induced instabilities from an initially axisymmetric gas disc. Angular momentum mixing in the disc produces some low angular momentum material, which can subsequently feed Sgr A*. These processes can occur in any galaxy that experiences an AGN episode, leading to bursts of nuclear star formation much stronger than pure bar-driven mass inflows would predict.

Recent analytical and numerical models show that AGN outflows and jets create ISM pressure in the host galaxy that is several orders of magnitude larger than in quiescent systems. This pressure increase can confine and compress molecular gas, thus accelerating star formation. In this paper, we model the effects of increased ambient ISM pressure on spherically symmetric turbulent molecular clouds. We find that large external pressure confines the cloud and drives a shockwave into it, which, together with instabilities behind the shock front, significantly accelerates the fragmentation rate. The compressed clouds therefore convert a larger fraction of their mass into stars over the cloud lifetime, and produce clusters that are initially more compact. Neither cloud rotation nor shear against the ISM affect this result significantly, unless the shear velocity is higher than the sound speed in the confining ISM. We conclude that external pressure is an important element in the star formation process, provided that it dominates over the internal pressure of the cloud.

It is usually assumed that outflows from luminous AGN are either in the energy-conserving (non-radiative) or in the momentum-conserving (radiative) regime. We show that in a non-spherical geometry the effects of both regimes may manifest at the same time, and that it is the momentum of the outflow that sets the $M_{\rm BH}-\sigma$ relation. Considering an initially elliptical distribution of gas in the host galaxy, we show that a non-radiative outflow opens up a wide ``escape route'' over the paths of least resistance. Most of the outflow energy escapes in that direction. At the same time, in the directions of higher resistance, the ambient gas is affected mainly by the incident momentum from the outflow. Quenching SMBH growth requires quenching gas delivery along the paths of highest resistance, and therefore, it is the momentum of the outflow that limits the black hole growth. We present an analytical argument showing that such energy-conserving feedback bubbles driving leaky ambient shells will terminate SMBH growth once its mass reaches roughly the $M_\sigma$ mass derived earlier by King (2003) for momentum-conserving AGN outflows. Our simulations also have potentially important implications for observations of AGN jet feedback and starburst galaxy feedback. The collimation of the wide angle AGN outflow away from the symmetry plane, as found in our simulations, indicates that credit for work done by such outflows may sometimes be mistakenly given to AGN jets or star formation feedback since wide angle $v \sim 0.1 c$ outflows are harder to observe and the phase when they are present may be short.

Several active galaxies show strong evidence for fast ($v_{\rm out} \sim 1000~{\rm km\,s}^{-1}$) massive ($\dot{M} =$ several $\times 1000~\msun\,{\rm yr}^{-1}$) gas outflows. Such outflows are expected on theoretical grounds once the central supermassive black hole reaches the mass set by the $M - \sigma$ relation, and may be what makes galaxies become red and dead. Despite their high velocities, which imply temperatures far above those necessary for molecule dissociation, the outflows contain large amounts of molecular gas. To understand this surprising result, we investigate the gas cooling and show that the outflows cannot stably persist in high--temperature states. Instead the outflowing gas forms a two--phase medium, with cold dense molecular clumps mixed with hot tenuous gas, as observed. We also show that efficient cooling leads to star formation, providing an observable outflow signature. The central parts of the outflows can be intrinsically luminous gamma--ray sources, provided that the central black hole is still strongly accreting. We note also that these outflows can persist for $\sim 10^8$~yr after the central AGN has turned off, so that many observed outflows (particularly with high speeds) otherwise assumed to be driven by starbursts might also be of this type.

Recent well resolved numerical simulations of AGN feedback have shown that its effects on the host galaxy may be not only negative but also positive. In the late gas poor phase, AGN feedback blows the gas away and terminates star formation. However, in the gas-rich phase(s), AGN outflows trigger star formation by over-compressing cold dense gas and thus provide positive feedback on their hosts. In this paper we study this AGN-triggered starburst effect. We show that star formation rate in the burst increases until the star formation feedback counteracts locally the AGN outflow compression. Globally, this predicts a strong nearly linear statistical correlation between the AGN and starburst bolometric luminosities in disc galaxies, L_* \propto L_{AGN}^{5/6}. The correlation is statistical only because AGN activity may fluctuate on short time scales (as short as tens of years), and because AGN may turn off but its effects on the host may continue to last until the AGN-driven outflow leaves the host, which may be up to 10 times longer than the duration of the AGN activity. The coefficient in front of this relation depends on the clumpiness and morphology of the cold gas in the galaxy. A "maximum starburst" takes place in am azimuthally uniform gas disc, for which we derive an upper limit of L_* \sim 50 times larger than L_{AGN} for typical quasars. For more clumpy and/or compact cold gas distributions, the starburst luminosity decreases. We also suggest that similar AGN-triggerred starbursts are possible in hosts of all geometries, including during galaxy mergers, provided the AGN is activated. Finally, we note that due to the short duration of the AGN activity phase the accelerating influence of AGN on starbursts may be much more common than observations of simultaneous AGN and starbursts would suggest.

Several tens of hyper-velocity stars (HVSs) have been discovered escaping our Galaxy. These stars share a common origin in the Galactic centre and are distributed anisotropically in Galactic longitude and latitude. We examine the possibility that HVSs may be created as the result of supernovae occurring within binary systems in a disc of stars around Sgr A* over the last 100 Myr. Monte Carlo simulations show that the rate of binary disruption is ~10^-4 yr^-1, comparable to that of tidal disruption models. The supernova-induced HVS production rate (\Gamma_HVS) is significantly increased if the binaries are hardened via migration through a gaseous disc. Moderate hardening gives \Gamma_HVS ~ 2*10^-7 yr^-1 and an estimated population of ~20 HVSs in the last 100 Myr. Supernova-induced HVS production requires the internal and external orbital velocity vectors of the secondary binary component to be aligned when the binary is disrupted. This leaves an imprint of the disc geometry on the spatial distribution of the HVSs, producing a distinct anisotropy.

We suggest a common physical origin connecting the fast, highly ionized winds (UFOs) seen in nearby AGN, and the slower and less ionized winds of BAL QSOs. The primary difference is the mass loss rate in the wind, which is ultimately determined by the rate at which mass is fed towards the central supermassive black hole (SMBH) on large scales. This is below the Eddington accretion rate in most UFOs, and slightly super-Eddington in extreme UFOs such as PG1211+143, but ranges up to $\sim 10-50$ times this in BAL QSOs. For UFOs this implies black hole accretion rates and wind mass loss rates which are at most comparable to Eddington, giving fast, highly-ionized winds. In contrast BAL QSO black holes have mildly super-Eddington accretion rates, and drive winds whose mass loss rates are significantly super-Eddington, and so are slower and less ionized. This picture correctly predicts the velocities and ionization states of the observed winds, including the recently-discovered one in SDSS J1106+1939. We suggest that luminous AGN may evolve through a sequence from BAL QSO through LoBAL to UFO-producing Seyfert or quasar as their Eddington factors drop during the decay of a bright accretion event. LoBALs correspond to a short-lived stage in which the AGN radiation pressure largely evacuates the ionization cone, but before the large-scale accretion rate has dropped to the Eddington value. We show that sub-Eddington wind rates would produce an $M - \sigma$ relation lying above that observed. We conclude that significant SMBH mass growth must occur in super-Eddington phases, either as BAL QSOs, extreme UFOs, or obscured from direct observation.

Quasar feedback outflows are commonly invoked to drive gas out of galaxies in the early gas-rich epoch to terminate growth of galaxies. Here we present simulations that show that AGN feedback may drive not only gas but also stars out of their host galaxies under certain conditions. The mechanics of this process is as following: (1) AGN-driven outflows accelerate and compress gas filling the host galaxy; (2) the accelerated dense shells become gravitationally unstable and form stars on radial trajectories. For the spherically symmetric initial conditions explored here, the black hole needs to exceed the host's M_sigma mass by a factor of a few to accelerate the shells and the new stars to escape velocities. We discuss potential implications of these effects for the host galaxies: (i) radial mixing of bulge stars with the rest of the host; (ii) contribution of quasar outflows to galactic fountains as sources of high-velocity clouds; (iii) wholesale ejection of hyper velocity stars out of their hosts, giving rise to type II supernovae on galactic outskirts, and contributing to reionization and metal enrichment of the Universe; (iv) bulge erosion and even complete destruction in extreme cases resulting in overweight or bulgeless SMBHs.

Galaxies become red and dead when the central supermassive black hole (SMBH) becomes massive enough to drive an outflow beyond the virial radius of the halo. We show that this final SMBH mass is larger than the final SMBH mass in the bulge of a spiral galaxy by up to an order of magnitude. The M - \sigma relations in the two galaxy types are almost parallel (M \propto \sigma^{4+\beta}, with \beta < 1) but offset in normalization, with the extra SMBH mass supplied by the major merger transforming the galaxy into an elliptical, or by mass gain in a galaxy cluster. This agrees with recent findings that SMBH in two Brightest Cluster Galaxies are \sim 10\times the expected M-\sigma mass. We show that these results do not strongly depend on the assumed profile of the dark matter halo, so analytic estimates found for an isothermal potential are approximately valid in all realistic cases. Our results imply that there are in practice actually {\it three} M - \sigma relations, corresponding to spiral galaxies with evolved bulges, field elliptical galaxies and cluster centre elliptical galaxies. A fourth relation, corresponding to cluster spiral galaxies, is also possible, but such galaxies are expected to be rare. All these relations have the form M_{\rm BH} = C_n\sigma^4, with only slight difference in slope between field and cluster galaxies, but with slightly different coefficients C_n. Conflating data from galaxies of different types and fitting a single relation to them tends to produce a higher power of \sigma.

Growing Supermassive Black Holes (SMBH) are believed to influence their parent galaxies in a negative way, terminating their growth by ejecting gas out before it could turn into stars. Here we present some of the most sophisticated SMBH feedback simulations to date showing that quasar's effects on galaxies are not always negative. We find that when the ambient shocked gas cools rapidly, the shocked gas is compressed into thin cold dense shells, filaments and clumps. Driving these high density features out is much more difficult than analytical models predict since dense filaments are resilient to the feedback. However, in this regime quasars have another way of affecting the host -- by triggering a massive star formation burst in the cold gas by over-pressurising it. Under these conditions SMBHs actually accelerate star formation in the host, having a positive rather than negative effect on their host galaxies. The relationship between SMBH and galaxies is thus even more complex and symbiotic than currently believed. We also suggest that the instabilities found here may encourage the chaotic AGN feeding mode.

Deposition of a massive ($10^4$ to $10^5 \msun$) giant molecular cloud (GMC) into the inner parsec of the Galaxy is widely believed to explain the origin of over a hundred unusually massive young stars born there $\sim 6$ Myr ago. An unknown fraction of that gas could have been accreted by Sgr A*, the supermassive black hole (SMBH) of the Milky Way. It has been recently suggested that two observed $\gamma$-ray-emitting bubbles emanating from the very center of our Galaxy were inflated by this putative activity of Sgr A*. We run a suite of numerical simulations to test whether the observed morphology of the bubbles could be due to the collimation of a wide angle outflow from Sgr A* by the disc-like Central Molecular Zone (CMZ), a well known massive repository of molecular gas in the central $\sim 200$ pc. We find that an Eddington-limited outburst of Sgr A* lasting $\simeq 1$ Myr is required to reproduce the morphology of the {\it Fermi} bubbles, suggesting that the GMC mass was $\sim 10^5 \msun$ and it was mainly accreted by Sgr A* rather than used to make stars. We also find that the outflow from Sgr A* enforces strong angular momentum mixing in the CMZ disc, robustly sculpting it into a much narrower structure -- a ring -- perhaps synonymous with the recently reported "Herschel ring". In addition, we find that Sgr A* outflow is likely to have induced formation of massive star-forming GMCs in the CMZ. In this scenario, the Arches and Quintuplet clusters, the two observed young star clusters in the central tens of parsecs of the Galaxy, and also GMCs such as Sgr B2, owe their existence to the recent Sgr A* activity.

Fermi bubbles, the recently observed giant (~10 kpc high) gamma-ray emitting lobes on either side of our Galaxy (Su et al. 2010), appear morphologically connected to the Galactic center, and thus offer a chance to test several models of supermassive black hole (SMBH) evolution, feedback and relation with their host galaxies. We use a physical feedback model (King 2003, 2010) and novel numerical techniques (Nayakshin et al. 2009) to simulate a short burst of activity in Sgr A*, the central SMBH of the Milky Way, ~6 Myr ago, temporally coincident with a star formation event in the central parsec. We are able to reproduce the bubble morphology and energetics both analytically (Zubovas et al. 2011) and numerically (Zubovas & Nayakshin, in prep). These results provide strong support to the model, which was also used to simulate more extreme environments (Nayakshin & Power 2010).

During the last decade, wide-angle powerful outflows from AGN, both on parsec and kpc scales, have been detected in many galaxies. These outflows are widely suspected to be responsible for sweeping galaxies clear of their gas. We present the analytical model describing the propagation of such outflows and calculate their observable properties. Large-scale AGN-driven outflows should have kinetic luminosities \sim {\eta}L_Edd/2 \sim 0.05L_Edd and momentum rates \sim 20L_Edd/c, where L_Edd is the Eddington luminosity of the central black hole and {\eta} \sim 0.1 its radiative accretion efficiency. This creates an expanding two-phase medium in which molecular species coexist with hot gas, which can persist after the central AGN has switched off. This picture predicts outflow velocities \sim 1000 - 1500 km/s and mass outflow rates up to 4000 M_\odot/yr on kpc scales, fixed mainly by the host galaxy velocity dispersion (or equivalently black hole mass). We compare our prediction with recent observational data, finding excellent agreement, and suggest future observational tests of this picture.

It is widely suspected that AGN activity ultimately sweeps galaxies clear of their gas. We work out the observable properties required to achieve this. Large-scale AGN-driven outflows should have kinetic luminosities $\sim \eta\le/2 \sim 0.05\le$ and momentum rates $\sim 20\le/c$, where $\le$ is the Eddington luminosity of the central black hole and $\eta\sim 0.1$ its radiative accretion efficiency. This creates an expanding two-phase medium in which molecular species coexist with hot gas, which can persist after the central AGN has switched off. This picture predicts outflow velocities $\sim 1000 - 1500$ km\,s$^{-1}$ and mass outflow rates up to $4000 \msun\,{\rm yr}^{-1}$ on kpc scales, fixed mainly by the host galaxy velocity dispersion (or equivalently black hole mass). All these features agree with those of outflows observed in galaxies such as Mrk231. This strongly suggests that AGN activity is what sweeps galaxies clear of their gas on a dynamical timescale and makes them red and dead. We suggest future observational tests of this picture.

It is theoretically expected that a supermassive black hole (SMBH) in the centre of a typical nearby galaxy disrupts a Solar-type star every ~ 10^5 years, resulting in a bright flare lasting for months. Sgr A*, the resident SMBH of the Milky Way, produces (by comparison) tiny flares that last only hours but occur daily. Here we explore the possibility that these flares could be produced by disruption of smaller bodies - asteroids. We show that asteroids passing within an AU of Sgr A* could be split into smaller fragments which then vaporise by bodily friction with the tenuous quiescent gas accretion flow onto Sgr A*. The ensuing shocks and plasma instabilities may create a transient population of very hot electrons invoked in several currently popular models for Sgr A* flares, thus producing the required spectra. We estimate that asteroids larger than ~ 10 km in size are needed to power the observed flares, with the maximum possible luminosity of the order 10^39 erg s^-1. Assuming that the asteroid population per parent star in the central parsec of the Milky Way is not too dissimilar from that around stars in the Solar neighbourhood, we estimate the asteroid disruption rates, and the distribution of the expected luminosities, finding a reasonable agreement with the observations. We also note that planets may be tidally disrupted by Sgr A* as well, also very infrequently. We speculate that one such disruption may explain the putative increase in Sgr A* luminosity ~ 300 yr ago.

{\it Fermi}-LAT has recently detected two gamma ray bubbles disposed symmetrically with respect to the Galactic plane. The bubbles have been suggested to be in a quasi-steady state, inflated by ongoing star formation over the age of the Galaxy. Here we propose an alternative picture where the bubbles are the remnants of a large-scale wide-angle outflow from \sgra, the SMBH of our Galaxy. Such an outflow would be a natural consequence of a short but bright accretion event on to \sgra\ if it happened concurrently with the well known star formation event in the inner 0.5 pc of the Milky Way $\sim 6$ Myr ago. We find that the hypothesised near-spherical outflow is focussed into a pair of symmetrical lobes by the greater gas pressure along the Galactic plane. The outflow shocks against the interstellar gas in the Galaxy bulge. Gamma--ray emission could be powered by cosmic rays created by either \sgra\ directly or accelerated in the shocks with the external medium. The Galaxy disc remains unaffected, agreeing with recent observational evidence that supermassive black holes do not correlate with galaxy disc properties. We estimate that an accreted mass $\sim 2 \times 10^3\msun$ is needed for the accretion event to power the observed {\it Fermi}--LAT lobes. Within a factor of a few this agrees with the mass of the young stars born during the star formation event. This estimate suggests that roughly 50% of the gas was turned into stars, while the rest accreted onto \sgra. One interpretation of this is a reduced star formation efficiency inside the \sgra\ accretion disc due to stellar feedback, and the other a peculiar mass deposition geometry that resulted in a significant amount of gas falling directly inside the inner $\sim 0.03$ pc of the Galaxy.

We discuss massive outflows in galaxy bulges, particularly ones driven by accretion episodes where the central supermassive black hole reaches the Eddington limit. We show that the quasar radiation field Compton--cools the wind shock until this reaches distances $\sim 1$ kpc from the black hole, but becomes too dilute to do this at larger radii. Radiative processes cannot cool the shocked gas within the flow time at any radius. Outflows are therefore momentum-driven at small radii (as required to explain the $M - \sigma$ relation). At large radii they are energy-driven, contrary to recent claims. We solve analytically the motion of an energy--driven shell after the central source has turned off. This shows that the thermal energy in the shocked wind can drive further expansion for a time $\sim 10$ times longer than the active time of the central source. Outflows observed at large radii with no active central source probably result from an earlier short (few Myr) active phase of this source.

We show that star formation in galaxy bulges is self-regulating through momentum feedback, limiting the stellar bulge mass to M_b ~ sigma^4. Together with a black hole mass M_BH ~ sigma^4 set by AGN momentum feedback, this produces a linear M_BH - M_b relation. At low redshift this gives M_BH/M_b ~ 0.001, close to the observed ratio. We show that AGN feedback can remove any remaining gas from the bulge and terminate star formation once the central black hole reaches the M_BH - sigma value, contrary to earlier claims. We find a mild upward deviation from the sigma^4 law at higher redshift and at higher sigma.