Source author record

Kastytis Zubovas

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

ResearcherUnclaimed source record

Catalog footprint

What is connected

23works
4topics
4close collaborators

Actions

Connect this record

Log in to claim

Research graph

See the researcher in context

Open full explorer

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

Building this map preview

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

Published work

23 published item(s)

preprint2022arXiv

Determining AGN luminosity histories using present-day outflow properties: a neural-network based approach

Large-scale outflows driven by active galactic nuclei (AGN) can have a profound influence on their host galaxies. The outflow properties themselves depend sensitively on the history of AGN energy injection during the lifetime of the outflow. Most observed outflows have dynamical timescales longer than the typical AGN episode duration, i.e. they have been inflated by multiple AGN episodes. Here, we present a neural-network based approach to inferring the most likely duty cycle and other properties of AGN based on the observable properties of their massive outflows. Our model recovers the AGN parameters of simulated outflows with typical errors $< 25\%$. We apply the method to a sample of 59 real molecular outflows and show that a large fraction of them have been inflated by AGN shining with a rather high duty cycle $δ_{\rm AGN} > 0.2$. This result suggests that nuclear activity in galaxies is clustered hierarchically in time, with long phases of more frequent activity composed of many short activity episodes. We predict that $\sim \! 19\%$ of galaxies should have AGN-driven outflows, but half of them are fossils - this is consistent with currently available data. We discuss the possibilities to investigate AGN luminosity histories during outflow lifetimes and suggest ways to use our software to test other physical models of AGN outflows. The source code of all of the software used here is made public.

preprint2022arXiv

Improving Black Hole Accretion Treatment in Hydrodynamical Simulations

The large galactic scales are connected to the many orders of magnitude smaller supermassive black hole (SMBH) scales by an episodic cycle of feeding and feedback. Active galactic nuclei (AGN) are powered by accretion onto SMBH and the majority of AGN energy, in near-Eddington regime, is produced in thin sub-pc accretion discs. Currently, it is very difficult to model processes that occur on vastly different scales, ranging from the circumnuclear gas reservoirs at tens to hundreds of parsecs, down to the accretion disc scales at <0.01 pc. While sub-grid prescriptions used in large-scale or cosmological simulations are able to reproduce large-scale feedback, we propose using a more realistic model in parsec-scale simulations, where it is important to get accurate timescales to understand how feedback affects gas dynamics and star formation in the vicinity of the AGN. To test our approach we use a sub-resolution thin accretion disc model, coupled to the SMBH, in a set of hydrodynamical simulations of a retrograde collision between a gas ring and a molecular cloud in an environment similar to the Galactic centre using the SPH code Gadget-3. The disc-mediated feeding of the SMBH is relatively smooth and delayed compared to an instantaneous feeding prescription. While the reduction of accretion due to feedback is present in both accretion disc and instantaneous feeding simulations, a clear central cavity appears only in accretion disc runs - hinting that a less volatile accretion phase could have a greater impact on the surrounding gas.

preprint2021arXiv

High-redshift SMBHs can grow from stellar-mass seeds via chaotic accretion

Extremely massive black holes, with masses $M_{\rm BH} > 10^9 M_\odot$, have been observed at ever higher redshifts. These results create ever tighter constraints on the formation and growth mechanisms of early black holes. Here we show that even the most extreme black hole known, Pōniuā'ena, can grow from a $10 M_\odot$ seed black hole via Eddington-limited luminous accretion, provided that accretion proceeds almost continuously, but is composed of a large number of episodes with individually-uncorrelated initial directions. This chaotic accretion scenario ensures that the growing black hole spins slowly, with the dimensionless spin parameter $a \lesssim 0.2$, so its radiative efficiency is also low, $ε\simeq 0.06$. If accretion is even partially aligned, with $20-40\%$ of accretion events happening in the same direction, the black hole spin and radiative efficiency are much higher, leading to significantly slower growth. We suggest that the chaotic accretion scenario can be completely falsified only if a $10^9 M_\odot$ black hole is discovered at $z \geq 9.1$, approximately $150$~Myr before Pōniuā'ena. The space density of extreme quasars suggests that only a very small fraction, roughly one in $4 \times 10^7$, of seed black holes need to encounter favourable growth conditions to produce the observed extreme quasars. Other seed black holes grow much less efficiently, mainly due to lower duty cycles, so are much more difficult to detect.

preprint2020arXiv

Evolution of dwarf galaxy observable parameters

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.

preprint2020arXiv

The case for the fundamental $M_{\rm BH}$-$σ$ relation

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}$-$σ_{\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}$-$σ_{\rm e}$ relation, with a focus on momentum-driven outflows.

preprint2016arXiv

The small observed scale of AGN--driven outflows, and inside--out disc quenching

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.

preprint2015arXiv

A simple way to improve AGN feedback prescription in SPH simulations

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.

preprint2015arXiv

AGN activity and nuclear starbursts: Sgr A* activity shapes the Central Molecular Zone

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.

preprint2015arXiv

The resolution bias: low resolution feedback simulations are better at destroying galaxies

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.

preprint2014arXiv

Collapse and fragmentation of molecular clouds under pressure

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.

preprint2014arXiv

Energy- and momentum-conserving AGN feedback outflows

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

preprint2014arXiv

Galaxy-wide outflows: cold gas and star formation at high speeds

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

preprint2013arXiv

AGN outflows trigger starbursts in gas-rich galaxies

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.

preprint2013arXiv

BAL QSOs and Extreme UFOs: the Eddington connection

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

preprint2013arXiv

Outflows of stars due to quasar feedback

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.

preprint2013arXiv

Supernovae in the Central Parsec: A Mechanism for Producing Spatially Anisotropic Hypervelocity Stars

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 (Γ_HVS) is significantly increased if the binaries are hardened via migration through a gaseous disc. Moderate hardening gives Γ_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.

preprint2012arXiv

AGN Winds and the Black-Hole - Galaxy Connection

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 η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 η \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.

preprint2012arXiv

Clearing Out a Galaxy

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 η\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 $η\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.

preprint2012arXiv

Fermi Bubbles in the Milky Way: the closest AGN feedback laboratory courtesy of Sgr A*?

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

preprint2012arXiv

Quasar feedback: accelerated star formation and chaotic accretion

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.

preprint2012arXiv

The M - σrelation in different environments

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 - σrelations in the two galaxy types are almost parallel (M \propto σ^{4+β}, with β< 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-σ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 - σ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σ^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 σ.

preprint2012arXiv

What's in a Fermi Bubble: a quasar episode in the Galactic centre

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

preprint2011arXiv

Sgr A* flares: tidal disruption of asteroids and planets?

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.