Dark rings with bright rims are the indirect signposts of planets embedded in protoplanetary discs. In a recent first, an azimuthally elongated AU-scale blob, possibly a planet, was resolved with ALMA in TW Hya. The blob is at the edge of a cliff-like rollover in the dust disc rather than inside a dark ring. Here we build time-dependent models of TW Hya disc. We find that the classical paradigm cannot account for the morphology of the disc and the blob. We propose that ALMA-discovered blob hides a Neptune mass planet losing gas and dust. We show that radial drift of mm-sized dust particles naturally explains why the blob is located on the edge of the dust disc. Dust particles leaving the planet perform a characteristic U-turn relative to it, producing an azimuthally elongated blob-like emission feature. This scenario also explains why a 10 Myr old disc is so bright in dust continuum. Two scenarios for the dust-losing planet are presented. In the first, a dusty pre-runaway gas envelope of about 40 Earth mass Core Accretion planet is disrupted, e.g., as a result of a catastrophic encounter. In the second, a massive dusty pre-collapse gas giant planet formed by Gravitational Instability is disrupted by the energy released in its massive core. Future modelling may discriminate between these scenarios and allow us to study planet formation in an entirely new way -- by analysing the flows of dust and gas recently belonging to planets, informing us about the structure of pre-disruption planetary envelopes.

Recent ALMA observations indicate that the majority of bright protoplanetary discs show signatures of young moderately massive planets. I show that this result is paradoxical. The planets should evolve away from their observed states by radial migration and gas accretion in about 1\% of the system age. These systems should then hatch tens of giant planets in their lifetime, and there should exist a very large population of bright planet-less discs; none of this is observationally supported. An alternative scenario, in which the population of bright ALMA discs is dominated by secondary discs recently rejuvenated by deposition of new gas, is proposed. The data are well explained if the gaseous mass of the discs is comparable to a Jovian planet mass, and they last a small fraction of a Million years. Self-disruptions of dusty gas giant protoplanets, previously predicted in the context of the Tidal Downsizing theory of planet formation, provide a suitable mechanism for such injections of new fuel, and yield disc and planet properties commensurate with ALMA observations. If this scenario is correct, then the secondary discs have gas-to-dust ratios considerably smaller than 100, and long look ALMA and NIR/optical observations of dimmer targets should uncover dusty, not yet disrupted, gas clumps with sizes of order an AU. Alternatively, secondary discs could originate from late external deposition of gas into the system, in which case we expect widespread signatures of warped outer discs that have not yet come into alignment with the planets.

Recent ALMA observations may indicate a surprising abundance of sub-Jovian planets on very wide orbits in protoplanetary discs that are only a few million years old. These planets are too young and distant to have been formed via the Core Accretion (CA) scenario, and are much less massive than the gas clumps born in the classical Gravitational Instability (GI) theory. It was recently suggested that such planets may form by the partial destruction of GI protoplanets: energy output due to the growth of a massive core may unbind all or most of the surrounding pre-collapse protoplanet. Here we present the first 3D global disc simulations that simultaneously resolve grain dynamics in the disc and within the protoplanet. We confirm that massive GI protoplanets may self-destruct at arbitrarily large separations from the host star provided that solid cores of mass around 10-20 Earth masses are able to grow inside them during their pre-collapse phase. In addition, we find that the heating force recently analysed by Masset and Velasco Romero (2017) perturbs these cores away from the centre of their gaseous protoplanets. This leads to very complicated dust dynamics in the protoplanet centre, potentially resulting in the formation of multiple cores, planetary satellites, and other debris such as planetesimals within the same protoplanet. A unique prediction of this planet formation scenario is the presence of sub-Jovian planets at wide orbits in Class 0/I protoplanetary discs.

Direct imaging (DI) surveys suggest that gas giants beyond 20 AU are rare around FGK stars. However, it is not clear what this means for the formation frequency of Gravitational Instability (GI) protoplanets due to uncertainties in gap opening and migration efficiency. Here we combine state-of-the-art calculations of homogeneous planet contraction with a population synthesis code. We find DI constraints to be satisfied if protoplanet formation by GI occurs in tens of percent of systems if protoplanets `super migrate' to small separations. In contrast, GI may occur in only a few percent of systems if protoplanets remain stranded at wide orbits because their migration is `quenched' by efficient gap opening. We then use the frequency of massive giants in radial velocity surveys inside 5 AU to break this degeneracy - observations recently showed that this population does not correlate with the host star metallicity and is therefore suspected to have formed via GI followed by inward migration. We find that only the super-migration scenario can sufficiently explain this population whilst simultaneously satisfying the DI constraints and producing the right mass spectrum of planets inside 5 AU. If massive gas-giants inside 5 AU formed via GI, then our models imply that migration must be efficient and that the formation of GI protoplanets occurs in at least a tens of percent of systems.

Thanks to recent high resolution ALMA observations, there is an accumulating evidence for presence of giant planets with masses from $\sim 0.01$ Jupiter mass to a few Jupiter mass with separations up to $ 100$~AU in the annular structures observed in young protoplanetary discs. We point out that these observations set unique "live" constraints on the process of gas accretion onto sub-Jovian planets that were not previously available. Accordingly, we use a population synthesis approach in a new way: we build time-resolved models and compare the properties of the synthetic planets with the ALMA data at the same age. Applying the widely used gas accretion formulae leads to a deficit of sub-Jovian planets and an over-abundance of a few Jupiter mass planets compared to observations. We find that gas accretion rate onto planets needs to be suppressed by about an order of magnitude to match the observed planet mass function. This slower gas giant growth predicts that the planet mass should correlate positively with the age of the protoplanetary disc, albeit with a large scatter. This effect is not clearly present in the ALMA data but may be confirmed in the near future with more observations.

Gas clumps formed within massive gravitationally unstable circumstellar discs are potential seeds of gas giant planets, brown dwarfs and companion stars. Simulations show that competition between three processes -- migration, gas accretion and tidal disruption -- establishes what grows from a given seed. Here we investigate the robustness of numerical modelling of clump migration and accretion with the codes PHANTOM, GADGET, SPHINX, SEREN, GIZMO-MFM, SPHNG and FARGO. The test problem comprises a clump embedded in a massive disc at an initial separation of 120 AU. There is a general qualitative agreement between the codes, but the quantitative agreement in the planet migration rate ranges from $\sim 10$% to $\sim 50$%, depending on the numerical setup. We find that the artificial viscosity treatment and the sink particle prescription may account for much of the differences between the codes. In order to understand the wider implications of our work, we also attempt to reproduce the planet evolution tracks from our hydrodynamical simulations with prescriptions from three previous population synthesis studies. We find that the disagreement amongst the population synthesis models is far greater than that between our hydrodynamical simulations. The results of our code comparison project are therefore encouraging in that uncertainties in the given problem are probably dominated by the physics not yet included in the codes rather than by how hydrodynamics is modelled in them.

Numerical simulations of pebble dynamics inside gas clumps formed by gravitational instability of protoplanetary discs are presented. We find that dust-mediated Rayleigh-Taylor instabilities transport pebbles inward rapidly via dense metal-rich "fingers". This speeds up sedimentation of small pebbles by up to two orders of magnitude and yet does not impede grain growth because grains of all sizes sediment at the same collective speed as long as Stokes number is less than unity. In simulations with a fixed pebble size, solid planetary cores form if pebble size exceeds a few cm. Pebble growth leads to core formation in some hundreds of years even when pebbles injected into clumps are of mm or smaller sizes. Properties of the gas clump dictate what kind of cores can be made. Low central temperature clumps allow formation of solid cores out of refractory materials, whereas in the highest temperature clumps pebbles of any composition are vaporised and make fuzzy cores only. These results confirm that gravitational instability of protoplanetary discs is a robust mechanism of hatching cores from sub-Earth to Neptune mass, as well as gas giants with massive cores, solid or fuzzy. This mode of planet formation is especially promising for environs too young and distant (such as the ALMA-observed HL Tau disc) or too violent (such as circum-binary planets) to form via the Core Accretion scenarios.

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.

We run numerical simulations to study the accretion of gas and dust grains onto gas giant planets embedded into massive protoplanetary discs. The outcome is found to depend on the disc cooling rate, planet mass, grain size and irradiative feedback from the planet. If radiative cooling is efficient, planets accrete both gas and pebbles rapidly, open a gap and usually become massive brown dwarfs. In the inefficient cooling case, gas is too hot to accrete onto the planet but pebble accretion continues and the planets migrate inward rapidly. Radiative feedback from the planet tends to suppress gas accretion. Our simulations predict that metal enrichment of planets by dust grain accretion inversely correlates with the final planet mass, in accordance with the observed trend in the inferred bulk composition of Solar System and exosolar giant planets. To account for observations, however, as much as ~30-50% of the dust mass should be in the form of large grains.

It has recently been suggested that in the presence of driven turbulence discs may be much less stable against gravitational collapse than their non turbulent analogs, due to stochastic density fluctuations in turbulent flows. This mode of fragmentation would be especially important for gas giant planet formation. Here we argue, however, that stochastic density fluctuations due to turbulence do not enhance gravitational instability and disc fragmentation in the long cooling time limit appropriate for planet forming discs. These fluctuations evolve adiabatically and dissipate away by decompression faster than they could collapse. We investigate these issues numerically in 2D via shearing box simulations with driven turbulence and also in 3D with a model of instantaneously applied turbulent velocity kicks. In the former setting turbulent driving leads to additional disc heating that tends to make discs more, rather than less, stable to gravitational instability. In the latter setting, the formation of high density regions due to convergent velocity kicks is found to be quickly followed by decompression, as expected. We therefore conclude that driven turbulence does not promote disc fragmentation in protoplanetary discs and instead tends to make the discs more stable. We also argue that sustaining supersonic turbulence is very difficult in discs that cool slowly.

Direct imaging observations constrain the fraction of stars orbited by gas giant planets with separations greater than 10 au to about 0.01 only. This is widely believed to indicate that massive protoplanetary discs rarely fragment on planetary mass objects. I use numerical simulations of gas clumps embedded in massive gas discs to show that these observations are consistent with $\sim 0.2 - 10$ planetary mass clumps per star being born in young gravitationally unstable discs. A trio of processes -- rapid clump migration, tidal disruption and runaway gas accretion -- destroys or transforms all of the simulated clumps into other objects, resulting in a desert of gas giants beyond separation of approximately 10 au. The cooling rate of the disc controls which of the three processes is dominant. For cooling rates faster than a few local dynamical times, clumps always grow rapidly and become massive brown dwarfs or low mass stars. For longer cooling times, post-collapse (high density) planets migrate inward to $\sim 10-20$ au where they open a gap in the disc and then continue to migrate inward much less rapidly. Pre-collapse (low density) planets are tidally disrupted and may leave massive solid cores behind. Gas giant planets observed inside the desert, such as those in HR 8799, must have followed an unusual evolutionary path, e.g., their host disc being dispersed in a catastrophic fashion.

Tidal Downsizing is the modern version of the Kuiper (1951) scenario of planet formation. Detailed simulations of self-gravitating discs, gas fragments, dust grain dynamics, and planet evolutionary calculations are summarised here and used to build a predictive planet formation model and population synthesis. A new interpretation of exoplanetary and debris disc data, the Solar System's origins, and the links between planets and brown dwarfs is offered. This interpretation is contrasted with the current observations and the predictions of the Core Accretion theory. Observations that can distinguish the two scenarios are pointed out. In particular, Tidal Downsizing predicts that presence of debris discs, sub-Neptune mass planets, planets more massive than $\sim 5$~Jupiter masses and brown dwarfs should not correlate strongly with the metallicity of the host. For gas giants of $\sim$ Saturn to a few Jupiter mass, a strong host star metallicity correlation is predicted only inwards of a few AU from the host. Composition of massive cores is predicted to be dominated by rock rather than ices. Debris discs made by Tidal Downsizing are distinct from those made by Core Accretion at birth: they have an innermost edge always larger than about 1 au, have smaller total masses and are usually in a dynamically excited state. It is argued that planet formation in surprisingly young or very dynamic systems such as HL Tau and Kepler-444 may be a signature of Tidal Downsizing. Open questions and potential weaknesses of the hypothesis are pointed out.

Recent observations of debris discs, believed to be made up of remnant planetesimals, brought a number of surprises. Debris disc presence does not correlate with the host star's metallicity, and may anti-correlate with the presence of gas giant planets. These observations contradict both assumptions and predictions of the highly successful Core Accretion model of planet formation. Here we explore predictions of the alternative Tidal Downsizing (TD) scenario of planet formation. In TD, small planets and planetesimal debris is made only when gas fragments, predecessors of giant planets, are tidally disrupted. We show that these disruptions are rare in discs around high metallicity stars but release more debris per disruption than their low [M/H] analogs. This predicts no simple relation between debris disc presence and host star's [M/H], as observed. A detected gas giant planet implies in TD that its predecessor fragment was not disputed, potentially explaining why DDs are less likely to be found around stars with gas giants. Less massive planets should correlate with DD presence, and sub-Saturn planets ($M_{p} \sim 50 \M_{\oplus}$) should correlate with DD presence stronger than sub-Neptunes ($M_{p} \leq 15 \M_{\oplus}$). These predicted planet-DD correlations will be diluted and weakened in observations by planetary systems' long term evolution and multi-fragment effects neglected here. Finally, although presently difficult to observe, DDs around M dwarf stars should be more prevalent than around Solar type stars.

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.

I argue that feedback is as important to formation of planets as it is to formation of stars and galaxies. Energy released by massive solid cores puffs up pre-collapse gas giant planets, making them vulnerable to tidal disruptions by their host stars. I find that feedback is the ultimate reason for some of the most robust properties of the observed exoplanet populations: the rarity of gas giants at all separations from $\sim 0.1$ to $\sim 100$~AU, the abundance of $\sim 10 M_\oplus$ cores but dearth of planets more massive than $\sim 20 M_\oplus$. Feedback effects can also explain (i) rapid assembly of massive cores at large separations as needed for Uranus, Neptune and the suspected HL Tau planets; (ii) the small core in Jupiter yet large cores in Uranus and Neptune; (iii) the existence of rare "metal monster" planets such as CoRoT-20b, a gas giant made of heavy elements by up to $\sim 50$\%.

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.

We present population synthesis calculations of the Tidal Downsizing (TD) hypothesis for planet formation. Our models address the following observations: (i) most abundant planets being Super Earths; (ii) cores more massive than $\sim 5-15 M_\oplus$ are enveloped by massive atmospheres; (iii) the frequency of occurrence of close-in gas giant planets correlates strongly with metallicity of the host star; (iv) no such correlation is found for sub-Neptune planets; (v) presence of massive cores in giant planets; (vi) gas giant planets are over-abundant in metals compared to their host stars; (vii) this over-abundance decreases with planet's mass; (viii) a deep valley in the planet mass function between masses of $\sim 10-20 M_\oplus$ and $\sim 100 M_\oplus$. A number of observational predictions distinguish the model from Core Accretion: (a) composition of the massive cores is always dominated by rocks not ices; (b) the core mass function is smooth with no minimum at $\sim 3 M_\oplus$ and has no ice-dominated cores; (c) gas giants beyond 10 AU are insensitive to the host star metallicity; (d) objects more massive than $\sim 10 M_{\rm Jup}$ do not correlate or even anti-correlate with metallicity. The latter prediction is consistent with observations of low mass stellar companions. TD can also explain formation of planets in close binary systems. TD model is a viable alternative to the Core Accretion scenario in explaining many features of the observed population of exoplanets.

We present numerical models of the gas dynamics in the inner parsec of the Galactic centre. We follow the gas from its origin as stellar winds of several observed young massive stars, until it is either captured by the central black hole, or leaves the system. Unlike our previous models, we include an outflow from the inner accretion flow. Two different kinds of outflows are modelled: (i) an instantaneous-response feedback mode, in which the outflow rate is directly proportional to the current black hole gas capture rate; and (ii) an outburst mode, which is stronger but lasts for a limited time. The latter situation may be particularly relevant to Sgr A*, since there is evidence that Sgr A* was much brighter in the recent past. We find that both types of outflow perturb the gas dynamics near the Bondi radius and the black hole capture rate significantly. The effects persist longer than the outflow itself. We also compare the effects of spherically symmetric and collimated outflows, and find that the latter are far less efficient in transferring its energy to the surrounding gas near the capture radius. Our results imply that accretion feedback is important for non-radiative accretion flows not only within but also outside the capture radius. Steady-state Bondi accretion rate estimates that do not account for feedback outflows over-predict not only the accretion rate onto the black hole but also the capture rate at the Bondi radius itself. Finally, the steady-state assumption under which non-radiative flows have been routinely studied in the literature may have to be abandoned if accretion feedback is bursty in nature.

Core Accretion (CA), the de-facto accepted theory of planet formation, requires formation of massive solid cores as a prerequisite for assembly of gas giant planets. The observed metallicity correlations of exoplanets are puzzling in the context of CA. While gas giant planets are found preferentially around metal-rich host stars, planets smaller than Neptune orbit hosts with a wide range of metallicities. We propose an alternative interpretation of these observations in the framework of a recently developed planet formation hypothesis called Tidal Downsizing (TD). We perform population synthesis calculations based on TD, and find that the connection between the populations of the gas giant and the smaller solid-core dominated planets is non linear and not even monotonic. While gas giant planets formed in the simulations in the inner few AU region follow a strong positive correlation with the host star metallicity, the smaller planets do not. The simulated population of these smaller planets shows a shallow peak in their formation efficiency at around the Solar metallicity. This result is driven by the fact that at low metallicities the solid core's growth is damped by the scarcity of metals, whereas at high metallicities the fragments within which the cores grow contract too quickly, cutting the core's growth time window short. Finally, simulated giant gas planets do not show a strong host star metallicity preference at large separations, which may explain why one of the best known directly imaged gas giant planet systems, HR 8799, is metal poor.

Frequency of detected giant planets is observed to increase rapidly with metallicity of the host star. This is usually interpreted as evidence in support of the Core Accretion (CA) theory, which assembles giant planets as a result of formation of a massive solid core. A strong positive planet-metallicity correlation for giant planets formed in the framework of Gravitational disc Instability (GI) model is found here. The key novelty of this work is "pebble accretion" onto GI fragments which has been recently demonstrated to accelerate contraction of GI fragments. Driven closer to the star by the inward migration, only the fragments that accrete metals rapidly enough collapse and survive the otherwise imminent tidal disruption. The survival fraction of simulated planets correlates strongly with the metallicity of the host star, as observed.

One of many challenges in forming giant gas planets via Gravitational disc Instability model (GI) is an inefficient radiative cooling of the pre-collapse fragments. Since fragment contraction times are as long at $10^5 -10^7$ years, the fragments may be tidally destroyed sooner than they contract into gas giant planets. Here we explore the role of "pebble accretion" onto the pre-collapse giant planets and find an unexpected result. Despite larger dust opacity at higher metallicities, addition of metals actually accelerates -- rather than slows down -- collapse of high opacity, relatively low mass giant gas planets ($M_p$ below a few Jupiter masses). A simple analytical theory that explains this result exactly in idealised simplified cases is presented. The theory shows that planets with the central temperature in the range between $\sim$ 1000 to 2000K are especially sensitive to pebble accretion: addition of just $\sim 5$ to 10 % of metals by weight is sufficient to cause their collapse. These results show that dust grain physics and dynamics is essential for an accurate modelling of self-gravitating disc fragments and their near environments in the outer massive and cold protoplanetary discs.

Tidal Downsizing (TD) is a recently developed planet formation theory that supplements the classical Gravitational disc Instability (GI) model with planet migration inward and tidal disruptions of GI fragments in the inner regions of the disc. Numerical methods for a detailed population synthesis of TD planets are presented here. As an example application, the conditions under which GI fragments collapse faster than they migrate into the inner $a\sim$ few AU disc are considered. It is found that most gas fragments are tidally or thermally disrupted unless (a) their opacity is $\sim 3$ orders of magnitude less than the interstellar dust opacity at metallicities typical of the observed giant planets, or (b) the opacity is high but the fragments accrete large dust grains (pebbles) from the disc. Case (a) models produce very low mass solid cores ($M_{\rm core} < 0.1$ Earth masses) and follow a negative correlation of giant planet frequency with host star metallicity. In contrast, case (b) models produce massive solid cores, correlate positively with host metallicity and explain naturally while giant gas planets are over-abundant in metals.

Ultrafast outflows (UFOs) from supermassive black holes (SMBHs) are thought to regulate the growth of SMBHs and host galaxies, resulting in a number of observational correlations. We present high-resolution numerical simulations of the impact of a thermalized UFO on the ambient gas in the inner part of the host galaxy. Our results depend strongly on whether the gas is homogeneous or clumpy. In the former case all of the ambient gas is driven outward rapidly as expected based on commonly used energy budget arguments, while in the latter the flows of mass and energy decouple. Carrying most of the energy, the shocked UFO escapes from the bulge via paths of least resistance, taking with it only the low-density phase of the host. Most of the mass is however in the high-density phase, and is affected by the UFO much less strongly, and may even continue to flow inwards. We suggest that the UFO energy leakage through the pores in the multiphase interstellar medium (ISM) may explain why observed SMBHs are so massive despite their overwhelmingly large energy production rates. The multiphase ISM effects reported here are probably under-resolved in cosmological simulations but may be included in prescriptions for active galactic nuclei feedback in future simulations and in semi-analytical models.

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.

Giant planet formation in the core accretion (CA) paradigm is predicated by the formation of a core, assembled by the coagulation of grains and later by planetesimals within a protoplanetary disc. In contrast, in the disc instability paradigm, giant planet formation is believed to be independent of core formation: massive self-gravitating gas fragments cool radiatively and collapse as a whole. We show that giant planet formation in the disc instability model may be also enhanced by core formation for reasons physically very similar to the CA paradigm. In the model explored here, efficient grain sedimentation within an initial fragment (rather than the disc) leads to the formation of a core composed of heavy elements. We find that massive atmospheres form around cores and undergo collapse as a critical core mass is exceeded, analogous to CA theory. The critical mass of the core to initiate such a collapse depends on the fragment mass and metallicity, as well as core luminosity, but ranges from less than 1 to as much as $\sim80$ Earth masses. We therefore suggest that there are two channels for the collapse of a gaseous fragment to planetary scales within the disc instability model: (i) H$_2$ dissociative collapse of the entire gaseous clump, and (ii) core-assisted gas capture, as presented here. We suggest that the first of these two is favoured in metal-poor environments and for fragments more massive than $\sim 5-10$ Jupiter masses, whereas the second is favored in metal-rich environments and fragments of lower mass. [Abridged]

Over the last decade, X-ray observations of Sgr A* have revealed a black hole in a deep sleep, punctuated roughly once per day by brief flares. The extreme X-ray faintness of this supermassive black hole has been a long-standing puzzle in black hole accretion. To study the accretion processes in the Galactic Center, Chandra (in concert with numerous ground- and space-based observatories) undertook a 3 Ms campaign on Sgr A* in 2012. With its excellent observing cadence, sensitivity, and spectral resolution, this Chandra X-ray Visionary Project (XVP) provides an unprecedented opportunity to study the behavior of the closest supermassive black hole. We present a progress report from our ongoing study of X-ray flares, including the brightest flare ever seen from Sgr A*. Focusing on the statistics of the flares and the quiescent emission, we discuss the physical implications of X-ray variability in the Galactic Center.

We study effects of AGN feedback outflows on multi-phase inter stellar medium (ISM) of the host galaxy. We argue that SMBH growth is dominated by accretion of dense cold clumps and filaments. AGN feedback outflows overtake the cold medium, compress it, and trigger a powerful starburst -- a positive AGN feedback. This predicts a statistical correlation between AGN luminosity and star formation rate at high luminosities. Most of the outflow's kinetic energy escapes from the bulge via low density voids. The cold phase is pushed outward only by the ram pressure (momentum) of the outflow. The combination of the negative and positive forms of AGN feedback leads to an $M-\sigma$ relation similar to the result of King (2003). Due to porosity of cold ISM in the bulge, SMBH influence on the low density medium of the host galaxy is significant even for SMBH well below the $M-\sigma$ mass. The role of SMBH feedback in our model evolves in space and time with the ISM structure. In the early gas rich phase, SMBH accelerates star formation in the bulge. During later gas poor (red-and-dead) phases, SMBH feedback is mostly negative everywhere due to scarcity of the cold ISM.

The large number of detected giant exoplanets offers the opportunity to improve our understanding of the formation mechanism, evolution, and interior structure of gas giant planets. The two main models for giant planet formation are core accretion and disk instability. There are substantial differences between these formation models, including formation timescale, favorable formation location, ideal disk properties for planetary formation, early evolution, planetary composition, etc. First, we summarize the two models including their substantial differences, advantages, and disadvantages, and suggest how theoretical models should be connected to available (and future) data. We next summarize current knowledge of the internal structures of solar- and extrasolar- giant planets. Finally, we suggest the next steps to be taken in giant planet exploration.

Bright AGN frequently show ultra-fast outflows (UFOs) with outflow velocities vout ~0.1c. These outflows may be the source of AGN feedback on their host galaxies sought by galaxy formation modellers. The exact effect of the outflows on the ambient galaxy gas strongly depends on whether the shocked UFOs cool rapidly or not. This in turn depends on whether the shocked electrons share the same temperature as ions (one temperature regime; 1T) or decouple (2T), as has been recently suggested. Here we calculate the Inverse Compton spectrum emitted by such shocks, finding a broad feature potentially detectable either in mid-to-high energy X-rays (1T case) or only in the soft X-rays (2T). We argue that current observations of AGN do not seem to show evidence for the 1T component. The limits on the 2T emission are far weaker, and in fact it is possible that the observed soft X-ray excess of AGN is partially or fully due to the 2T shock emission. This suggests that UFOs are in the energy-driven regime outside the central few pc, and must pump considerable amounts of not only momentum but also energy into the ambient gas. We encourage X-ray observers to look for the Inverse Compton components calculated here in order to constrain AGN feedback models further.

Most supermassive black holes (SMBHs) are accreting at very low levels and are difficult to distinguish from the galaxy centers where they reside. Our own Galaxy's SMBH provides a uniquely instructive exception, and we present a close-up view of its quiescent X-ray emission based on 3 mega-second of Chandra observations. Although the X-ray emission is elongated and aligns well with a surrounding disk of massive stars, we can rule out a concentration of low-mass coronally active stars as the origin of the emission based on the lack of predicted Fe Kalpha emission. The extremely weak H-like Fe Kalpha line further suggests the presence of an outflow from the accretion flow onto the SMBH. These results provide important constraints for models of the prevalent radiatively inefficient accretion state.

We present the first systematic analysis of the X-ray variability of Sgr A* during the Chandra X-ray Observatory's 2012 Sgr A* X-ray Visionary Project (XVP). With 38 High Energy Transmission Grating Spectrometer (HETGS) observations spaced an average of 7 days apart, this unprecedented campaign enables detailed study of the X-ray emission from this supermassive black hole at high spatial, spectral and timing resolution. In 3 Ms of observations, we detect 39 X-ray flares from Sgr A*, lasting from a few hundred seconds to approximately 8 ks, and ranging in 2-10 keV luminosity from ~1e34 erg/s to 2e35 erg/s. Despite tentative evidence for a gap in the distribution of flare peak count rates, there is no evidence for X-ray color differences between faint and bright flares. Our preliminary X-ray flare luminosity distribution dN/dL is consistent with a power law with index -1.9 (+0.3 -0.4); this is similar to some estimates of Sgr A*'s NIR flux distribution. The observed flares contribute one-third of the total X-ray output of Sgr A* during the campaign, and as much as 10% of the quiescent X-ray emission could be comprised of weak, undetected flares, which may also contribute high-frequency variability. We argue that flares may be the only source of X-ray emission from the inner accretion flow.

It is well known that massive protoplanetary disc are gravitationally unstable beyond tens of AU from their parent star. The eventual fate of the self-gravitating gas clumps born in the disc is currently not understood, although the range of uncertainty is well known. If clumps migrate inward rapidly, they are tidally disrupted, which may leave behind giant or terrestrial like planets. On the other hand, if clumps migrate less rapidly, they tend to accrete gas, becoming proto brown dwarfs or low mass companions to the parent star. Here we argue that radiative feedback of contracting clumps (protoplanets) on their discs is an important effect that has been overlooked in previous calculations. We show analytically that temperature in clump's vicinity may be high enough to support a quasi-static atmosphere if the clump mass is below a critical value, $M_{\rm cr} \sim 6$ Jupiter masses ($M_J$). This may arrest further gas accretion onto the clump and thus promote formation of planets rather than low mass companions. We use numerical simulations to evaluate these analytical conclusions by studying migration and accretion of gas clumps as a function of their initial mass, $M_i$. Simulations neglecting the radiative preheating effect show that gas clumps with mass less than $\sim 2 M_J$ migrate inward rapidly; more massive clumps result in low mass companions. In contrast, simulations that include radiative preheating from the clump show that clumps as massive as $8 M_J$ migrate inward rapidly and end up tidally disrupted in the inner disc. We conclude that, with other parameters being equal, previous simulations neglecting radiative feedback from self-gravitating clumps over-estimated the population of brown dwarfs and low mass stellar companions and under-estimated the population of planets.

We suggest a direct link between the two "missing" baryon problems of contemporary galaxy formation theory: (1) that large galaxies are known to contain too little gas and stars and (2) that too few dwarf satellite galaxies are observed around large galaxies compared with cosmological simulations. The former can be explained by invoking some energetic process -- most likely AGN or star formation feedback -- which expels to infinity a significant fraction of the gas initially present in the proto-galaxy, while the latter problem is usually explained by star formation feedback inside the dwarf or tidal and ram pressure stripping of the gas from the satellite galaxy by its parent. Here we point out that the host galaxy "missing" baryons, if indeed ejected at velocities of hundreds to a thousand km s$^{-1}$, must also affect smaller satellite galaxies by stripping or shocking the gas there. We estimate the fraction of gas removed from the satellites as a function of the satellite galaxy's properties. Applying these results to a Milky Way like dark matter halo, we find that this singular shock ram pressure stripping event may be quite efficient in removing the gas from the satellites provided that they are close enough. We also use the orbital and mass modelling data for eight Galactic dwarf spheroidal (dSph) satellites, and find that it is likely that many of them have been affected by the Galactic outflow. We point out that galactic outflows of the host may also trigger a starburst in the satellite galaxies by over-pressuring their gas discs. This process may be responsible for the formation of the globular clusters observed in some of the Milky Way's dSphs (e.g. the Fornax and Sagittarius dSphs).

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.

We point out that protoplanets created in the framework of the Tidal Downsizing (TD) theory for planet formation play a very important role for the evolution of accretion discs hosting them. Since all TD protoplanets are initially as massive as $\sim 10$ Jupiter masses, they are able to open very deep gaps in their discs, and even completely isolate the inner disc flows from the outer ones. Furthermore, in contrast to other planet formation theories, TD protoplanets are mass donors for their protostars. One potentially observable signature of planets being devoured by their protostars are FU Ori like outbursts, and episodic protostar accretion more generally, as discussed by a number of authors recently. Here we explore another observational implication of TD hypothesis: dust poor inner accretion flows, which we believe may be relevant to some of the observed mm-bright transitional discs around protostars. In our model, a massive protoplanet interrupts the flow of the outer dust-rich disc on its protostar, and at the same time looses a part of its dust-poor envelope into the inner disc. This then powers the observed gas-but-no-dust accretion onto the star. Upon a more detailed investigation, we find that this scenario is quite natural for young massive discs but is less so for older discs, e.g., those whose self-gravitating phase has terminated a fraction of a Million year or more ago. This stems from the fact that TD protoplanets of such an age should have contracted significantly, and so are unlikely to loose much mass. Therefore, we conclude that either (i) the population of "transition discs" with large holes and dust-poor accretion is much younger than generally believed; or (ii) there is a poorly understood stage for late removal of dust-poor envelopes from TD planets; (iii) another explanation for the observations is correct.

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.

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.

We argue that current observations of $M - \sigma$ relations for galaxies can be used to constrain theories of super-massive black holes (SMBH) feeding. In particular, assuming that SMBH mass is limited only by the feedback on the gas that feeds it, we show that SMBHs fed via a planar galaxy scale gas flow, such as a disc or a bar, should be much more massive than their counterparts fed by quasi-spherical inflows. This follows from the relative inefficiency of AGN feedback on a flattened inflow. We find that even under the most optimistic conditions for SMBH feedback on flattened inflows, the mass at which the SMBH expels the gas disc and terminates its own growth is a factor of several higher than the one established for quasi-spherical inflows. Any beaming of feedback away from the disc and any disc self-shadowing strengthens this result further. Contrary to this theoretical expectation, recent observations have shown that SMBH in pseudobulge galaxies (which are associated with barred galaxies) are typically under- rather than over-massive when compared with their classical bulge counterparts at a fixed value of $\sigma$. We conclude from this that SMBHs are not fed by large (100 pc to many kpc) scale gas discs or bars, most likely because such planar flows are turned into stars too efficiently to allow any SMBH growth. Based on this and other related observational evidence, we argue that most SMBHs grow by chaotic accretion of gas clouds with a small and nearly randomly distributed direction of angular momentum.

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.

Core Accretion, the most widely accepted scenario for planet formation, postulates existence of km-sized solid bodies, called planetesimals, arranged in a razor-thin disc in the earliest phases of planet formation. In the Tidal Downsizing hypothesis, an alternative scenario for formation of planets, grain growth, sedimentation and formation of planetary cores occur inside dense and massive gas clumps formed in the outer cold disc by gravitational instability. As a clump migrates inward, tidal forces of the star remove all or most of the gas from the clump, downsizing it to a planetary mass body. Here we argue that such a clump may form not only the planetary core but also numerous smaller bodies. As an example, we consider the simplest case of bodies on circular orbits around the planetary core in the centre of the gas clump. Bodies smaller than 1 km suffer a strong enough aerodynamic drag, spiral in and accrete onto the solid core rapidly; bodies in the planetesimal size range lose their centrifugal support very slowly. We find that planetesimals orbiting the protoplanetary core closely remain gravitationally bound to it; these may be relevant to formation of satellites of giant planets. Planetesimals on more distant orbits within the host clump are unbound from the protoplanet and are set on mildly eccentric heliocentric orbits, generically forming wide rings. These may correspond to debris discs around main sequence stars and the Kuiper belt in the Solar System. For the latter in particular, our hypothesis naturally explains the observed sharp outer edge and the "mass deficit" of the Kuiper belt.

Accretion onto supermassive black holes (SMBHs) in galaxy formation simulations is frequently modelled by the Bondi-Hoyle formalism. Here we examine the validity of this approach analytically and numerically. We argue that the character of the flow where one evaluates the gas properties is unlikely to satisfy the simple Bondi-Hoyle model. Only in the specific case of hot virialised gas with zero angular momentum and negligible radiative cooling is the Bondi-Hoyle solution relevant. In the opposite extreme, where the gas is in a state of free-fall at the evaluation radius due to efficient cooling and the dominant gravity of the surrounding halo, the Bondi-Hoyle formalism can be erroneous by orders of magnitude in either direction. This may impose artificial trends with halo mass in cosmological simulations by being wrong by different factors for different halo masses. We propose an expression for the sub-grid accretion rate which interpolates between the free-fall regime and the Bondi-Hoyle regime, therefore taking account of the contribution of the halo to the gas dynamics.

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

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 has been recently proposed that giant protoplanets migrating inward through the disc more rapidly than they contract could be tidally disrupted when they fill their Roche lobes ~ 0.1 AU away from their parent protostars. Here we consider the process of mass and angular momentum exchange between the tidally disrupted planet and the surrounding disc in detail. We find that the planet's adiabatic mass-radius relation and its ability to open a deep gap in the disc determine whether the disruption proceeds as a sudden runaway or a balanced quasi-static process. In the latter case the planet feeds the inner disc through its Lagrangian L1 point like a secondary star in a stellar binary system. As the planet loses mass it gains specific angular momentum and normally migrates in the outward direction until the gap closes. Numerical experiments show that planet disruption outbursts are preceded by long "quiescent" periods during which the disc inward of the planet is empty. The hole in the disc is created when the planet opens a deep gap, letting the inner disc to drain onto the star while keeping the outer one stalled behind the planet. We find that the mass-losing planet embedded in a realistic protoplanetary disc spawns an extremely rich set of variability patterns. In a subset of parameter space, there is a limit cycle behaviour caused by non-linear interaction between the planet mass loss and the disc hydrogen ionisation instability. We suggest that tidal disruptions of young massive planets near their stars may be responsible for the observed variability of young accreting protostars such as FU Ori, EXor and T Tauri stars in general.

The last decade has seen a dramatic confirmation that an in situ star formation is possible inside the inner parsec of the Milky Way. Here we suggest that giant planets, solid terrestrial-like planets, comets and asteroids may also form in these environments, and that this may have observational implications for Active Galactic Nuclei (AGN). Like in debris discs around main sequence stars, collisions of large solid objects should initiate strong fragmentation cascades. The smallest particles in such a cascade - the microscopic dust - may provide a significant opacity. We put a number of observational and physical constraints on AGN obscuring torii resulting from such fragmentation cascades. We find that torii fed by fragmenting asteroids disappear at both low and high AGN luminosities. At high luminosities, $L \sim L_{\rm Edd}$, where $L_{\rm Edd}$ is the Eddington limit, the AGN radiation pressure blows out the microscopic dust too rapidly. At low luminosities, on the other hand, the AGN discs may avoid gravitational fragmentation into stars and solids. We also note that these fragmentation cascades may be responsible for astrophysically "large" dust particles of approximately micrometer sizes that were postulated by some authors to explain unusual absorption properties of the AGN torii.

We study the dynamical evolution of stars and gas close to the centre of the Milky Way. Any plausible means of forming the young stars observed at the Galactic Centre leaves behind a residual gas disc at ~0.01pc radii. We show that the combined effects of viscous accretion and gravitational interactions with stars do not remove the residual gas efficiently, and that a substantial gas disc, interior to the stellar disc, persists for >10Myr after the stars form. Since no such disc is currently seen at the Galactic Centre we argue that it has been accreted by the super-massive black hole. This scenario offers an attractive connection between nuclear star formation and black hole feeding, and we suggest that the "missing" gas may have been used to power Sgr A*.

Recent Kepler observations revealed an unexpected abundance of "hot" Earth-size to Neptune-size planets in the inner ~0.02-0.2 AU from their parent stars. Another "oddity" of recent observations is the discovery of a hot jupiter planet around a very metal poor star of extragalactic origin. We point out that these observations have natural explanations in the framework of the Tidal Downsizing hypothesis, a new and an alternative model for planet formation. We show that disassociation of molecular hydrogen leads to the "second collapse" of the giant planet embryos. The collapse forms much denser fluid configurations, closely resembling isolated proto-giant planets on the Hayashi track before they enter the degenerate "main sequence" phase. We find that, in the early gas-rich phase of proto-planetary discs, these configurations migrate inward more rapidly than they contract. They are therefore tidally disrupted at characteristic distance of ~ 0.05 -0.1 AU from their parent stars. The disruption is likely to leave solid cores behind, which we identify with the hot Super-Earths observed by Kepler. In addition, we show that single hot jupiters (also found by the Kepler telescope) are not surprising in the context of TD hypothesis. Finally, formation of coreless giant planets by the second collapse of the giant embryos occurs for either very massive embryos (more than ~ 20 Jupiter masses) or very metal poor ones. We suggest that this process explains the extra-galactic hot jupiter.

It has long been recognised that the main obstacle to accretion of gas onto supermassive black holes (SMBHs) is large specific angular momentum. It is feared that the gas settles in a large scale disc, and that accretion would then proceed too inefficiently to explain the masses of the observed SMBHs. Here we point out that, while the mean angular momentum in the bulge is very likely to be large, the deviations from the mean can also be significant. Indeed, cosmological simulations show that velocity and angular momentum fields of gas flows onto galaxies are very complex. Furthermore, inside bulges the gas velocity distribution can be further randomised by the velocity kicks due to feedback from star formation. We perform hydrodynamical simulations of gaseous rotating shells infalling onto an SMBH, attempting to quantify the importance of velocity dispersion in the gas at relatively large distances from the black hole. We implement this dispersion by means of a supersonic turbulent velocity spectrum. We find that, while in the purely rotating case the circularisation process leads to efficient mixing of gas with different angular momentum, resulting in a low accretion rate, the inclusion of turbulence increases this accretion rate by up to several orders of magnitude. We show that this can be understood based on the notion of "ballistic" accretion, whereby dense filaments, created by convergent turbulent flows, travel through the ambient gas largely unaffected by hydrodynamical drag. We derive a simple analytical formula that captures the numerical results to within a factor of a few. Rescaling our results to astrophysical bulges, we argue that this "ballistic" mode of accretion could provide the SMBHs with a sufficient supply of fuel without the need to channel the gas via large-scale discs or bars, and therefore that star formation in bulges can be a strong catalyst for SMBH accretion.

We use SPH simulations with an approximate radiative cooling prescription to model evolution of a massive and large ($\sim 100$ AU) very young protoplanetary disc. We also model dust growth and gas-grain dynamics with a second fluid approach. It is found that the disc fragments onto a large number of $\sim 10$ Jupiter mass clumps that cool and contract slowly. Some of the clumps evolve onto eccentric orbits delivering them into the inner tens of AU, where they are disrupted by tidal forces from the star. Dust grows and sediments inside the clumps, displaying a very strong segregation, with the largest particles forming dense cores in the centres. The density of the dust cores may exceed that of the gas and is limited only by the numerical constraints, indicating that these cores should collapse into rocky planetary cores. One particular giant planet embryo migrates inward close enough to be disrupted at about 10 AU, leaving a self-bound solid core of about 7.5 $\mearth$ mass on a low eccentricity orbit at a radius of $\sim$ 8 AU. These simulations support the recent suggestions that terrestrial and giant planets may be the remnants of tidally disrupted giant planet embryos.

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

While gas accretion onto some massive black holes (MBHs) at the centers of galaxies actively powers luminous emission, the vast majority of MBHs are considered dormant. Occasionally, a star passing too near a MBH is torn apart by gravitational forces, leading to a bright panchromatic tidal disruption flare (TDF). While the high-energy transient Swift J164449.3+573451 ("Sw 1644+57") initially displayed none of the theoretically anticipated (nor previously observed) TDF characteristics, we show that the observations (Levan et al. 2011) suggest a sudden accretion event onto a central MBH of mass ~10^6-10^7 solar masses. We find evidence for a mildly relativistic outflow, jet collimation, and a spectrum characterized by synchrotron and inverse Compton processes; this leads to a natural analogy of Sw 1644+57 with a smaller-scale blazar. The phenomenologically novel Sw 1644+57 thus connects the study of TDFs and active galaxies, opening a new vista on disk-jet interactions in BHs and magnetic field generation and transport in accretion systems.

Comets are believed to be born in the outer Solar System where the temperature is assumed to have never exceeded T ~ 100 K. Surprisingly, observations and samples of cometary dust particles returned to Earth showed that they are in fact made of a mix of ices, as expected, but also of materials forged at high-temperatures (T ~ 1500 K). We propose a radically new view regarding the origin of the high-temperature processed materials in comets, based on the recent "Tidal Downsizing" (TD) hypothesis for planet formation. In the latter, the outer proto-planetary disc is gravitationally unstable and forms massive giant planet embryos (GEs). These hot (T ~ hundreds to 2000 K) and dense regions, immersed in the background cold and low density disc, are eventually disrupted. We propose that both planets and the high-T materials in comets are synthesised inside the GEs. Disruption of GEs separates planets and small solids as the latter are "frozen-in" into gas and are peeled off together with it. These small solids are then mixed with the ambient cold disc containing ices before being incorporated into comets. Several predictions of this picture may be testable with future observations of the Solar System and exoplanets.

Recent {\em Kepler} observations revealed an unexpected abundance of "hot" Earth-size to Neptune-size planets in the inner $0.02-0.2$ AU from their parent stars. We propose that these smaller planets are the remnants of massive giant planets that migrated inward quicker than they could contract. We show that such disruptions naturally occur in the framework of the Tidal Downsizing hypothesis for planet formation. We find that the characteristic planet-star separation at which such "hot disruptions" occur is $R \approx 0.03-0.2$ AU. This result is independent of the planet's embryo mass but is dependent on the accretion rate in the disc. At high accretion rates, $\dot M \simgt 10^{-6}\msun$ yr$^{-1}$, the embryo is unable to contract quickly enough and is disrupted. At late times, when the accretion rate drops to $\dot M \simlt 10^{-8} \msun$ yr$^{-1}$, the embryos migrate sufficiently slow to not be disrupted. These "late arrivals" may explain the well known population of hot jupiters. If type I migration regime is inefficient, then our model predicts a pile-up of planets at $R\sim 0.1$ AU as the migration rate suddenly switches from the type II to type I in that region.

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.

The standard picture of planet formation posits that giant gas planets are over-grown rocky planets massive enough to attract enormous gas atmospheres. It has been shown recently that the opposite point of view is physically plausible: the rocky terrestrial planets are former giant planet embryos dried of their gas "to the bone" by the influences of the parent star. Here we provide a brief overview of this "Tidal Downsizing" hypothesis in the context of the Solar System structure.

Black holes grow by accreting matter from their surroundings. However, angular momentum provides an efficient natural barrier to accretion and so only the lowest angular momentum material will be available to feed the black holes. The standard sub-grid model for black hole accretion in galaxy formation simulations - based on the Bondi-Hoyle method - does not account for the angular momentum of accreting material, and so it is unclear how representative the black hole accretion rate estimated in this way is likely to be. In this paper we introduce a new sub-grid model for black hole accretion that naturally accounts for the angular momentum of accreting material. Both the black hole and its accretion disc are modelled as a composite accretion disc particle. Gas particles are captured by the accretion disc particle if and only if their orbits bring them within its accretion radius R_acc, at which point their mass is added to the accretion disc and feeds the black hole on a viscous timescale t_visc. The resulting black hole accretion rate (dM/dt)_BH powers the accretion luminosity L_acc ~ (dM/dt)_BH, which drives black hole feedback. Using a series of controlled numerical experiments, we demonstrate that our new accretion disc particle method is more physically self-consistent than the Bondi-Hoyle method. We also discuss the physical implications of the accretion disc particle method for systems with a high degree of rotational support, and we argue that the M_BH-sigma relation in these systems should be offset from the relation for classical bulges and ellipticals, as appears to be observed.

It has been proposed recently that the first step in the formation of both rocky and gas giant planets is dust sedimentation into a solid core inside a gas clump (giant planet embryo). The clumps are then assumed to migrate closer to the star where their metal poor envelopes are sheared away by the tidal forces or by an irradiation-driven mass loss. We consider the implications of this hypothesis for natal rotation rates of both terrestrial and gas giant planets. It is found that both types of planets may rotate near their break up angular frequencies at birth. The direction of the spin should coincide with that of the parent disc and the star, except in cases of embryos that had close interactions or mergers with other embryos in the past. Furthermore, the large repository of specific angular momentum at birth also allows formation of close binary rocky planets inside the same embryos. We compare these predictions with rotation rates of planets in the Solar System and also question whether the Earth-Moon pair could have been formed within the same giant planet embryo.

Giant planet embryos are believed to be spawned by gravitational instability in massive extended (R ~ 100 AU) protostellar discs. In a recent paper we have shown that dust can sediment inside the embryos, as argued earlier by Boss (1998) in a slightly different model. Here we study numerically the next stage of this process -- the formation of a solid core. If conditions are conducive to solid core formation, the centre of the gas cloud goes through the following sequence of phases: (i) becomes grain (and metal) rich; (ii) forms a terrestrial mass solid core via a rapid collapse driven by self-gravity of the grains; (iii) starts to accrete a gaseous atmosphere when the solid core reaches mass of a few to 10 Earth masses. This sequence of events may build either terrestrial planet cores or metal rich giant planets inside the larger gas reservoir of the giant planet embryo. In a companion Letter we argue that tidal and irradiation effects from the parent star should disrupt the outer metal poor layers of the embryo, releasing nearly "ready to use" planets. We propose this as an alternative way to build planets.

In the context of massive fragmenting protoplanetary discs, Boss (1998) suggested that grains can grow and sediment inside giant planet embryos formed at R ~ 5 AU away from the star. Several authors since then criticised the suggestion. Convection may prevent grain sedimentation, and the embryos cannot even form so close to the parent star as cooling is too inefficient at these distances. Here we reconsider the grain sedimentation process suggested by Boss (1998) but inside an embryo formed, as expected in the light of the cooling constraints, at R ~ 100 AU. Such embryos are much less dense and are also cooler. We make analytical estimates of the process and also perform simple spherically symmetric radiation hydrodynamics simulations to test these ideas. We find that convection in our models does not become important before a somewhat massive (~ an Earth mass, this is clarified in a followup paper) solid core is built. Turbulent mixing slows down dust sedimentation but is overwhelmed by grain sedimentation when the latter grow to a centimetres size. The minimum time required for dust sedimentation to occur is a few thousand years, and is a strong function of the embryo's mass, dust content and opacity. An approximate analytical criterion is given to delineate conditions in which a giant embryo contracts and heats up faster than dust can sediment. As Boss et al (2002), we argue that core formation through grain sedimentation inside the giant planet embryos may yield an unexplored route to form giant gas and giant ice planets. The present model also stands at the basis of paper III, where we study the possibility of forming terrestrial planet cores by tidal disruption and photoevaporation of the planetary envelope.

We hypothesise that planets are made by tidal downsizing of migrating giant planet embryos. The proposed scheme for planet formation consists of these steps: (i) a massive young protoplanetary disc fragments at R ~ several tens to hundreds of AU on gaseous clumps with masses of a few Jupiter masses; (ii) the clumps cool and contract, and simultaneously migrate closer in to the parent star; (iii) as earlier suggested by Boss (1998), dust sediments inside the gas clumps to form terrestrial mass solid cores; (iv) if the solid core becomes more massive than ~ 10 Earth masses, a massive gas atmosphere collapses onto the solid core; (v) when the gas clumps reach the inner few AU from the star, tidal shear and evaporation due to stellar irradiation peel off the outer metal-poor envelope of the clump. If tidal disruption occurs quickly, while the system is still in stage (iii), a terrestrial planet core is left. If it happens later, in stage (iv), a metal rich gas giant planet with a solid core emerges from the envelope.

The observed super-massive black hole (SMBH) mass -- galaxy velocity dispersion ($M_{\rm cmo} - \sigma$) correlation, and the similar correlation for nuclear star clusters, may be established when winds/outflows from the CMO ("central massive object") drive gas out of the potential wells of classical bulges. Timescales of growth for these objects may explain why smaller bulges appear to host preferentially NCs while larger ones contain SMBHs only. Despite much recent progress, feedback processes in bulge/galaxy formation are far from being understood. Our numerical simulations show that that understanding how the CMO feeds is as important a piece of the puzzle as understanding how its feedback affects its host galaxy.

Numerical simulations for the non-linear development of Kelvin-Helmholtz instability in two different density layers have been performed with the particle-based method (Godunov SPH) developed by Inutsuka (2002). The Godunov SPH can describe the Kelvin-Helmholtz instability even with a high density contrast, while the standard SPH shows the absence of the instability across a density gradient (Agertz et al. 2007). The interaction of a dense blob with a hot ambient medium has been performed also. The Godunov SPH describes the formation and evolution of the fingers due to the combinations of Rayleigh-Taylor, Richtmyer-Meshkov, and Kelvin-Helmholtz instabilities. The blob test result coincides well with the results of the grid-based codes. An inaccurate handling of a density gradient in the standard SPH has been pointed out as the direct reason of the absence of the instabilities. An unphysical force happens at the density gradient even in a pressure equilibrium, and repulses particles from the initial density discontinuity. Therefore, the initial perturbation damps, and a gap forms at the discontinuity. The unphysical force has been studied in terms of the consistency of a numerical scheme. Contrary to the standard SPH, the momentum equation of the Godunov SPH doesnt use the particle approximation, and has been derived from the kernel convolution or a new Lagrangian function. The new Lagrangian function used in the Godunov SPH is more analogous to the real Lagrangian function for continuum. The momentum equation of the Godunov SPH has much better linear consistency, so the unphysical force is greatly reduced compared to the standard SPH in a high density contrast.

The observed super-massive black hole (SMBH) mass -- galaxy velocity dispersion ($M_{\rm bh} - \sigma$) correlation may be established when winds/outflows from the SMBH drive gas out of the potential wells of classical bulges. Here we present numerical simulations of this process in a static isothermal potential. Simple spherically symmetric models of SMBH feedback at the Eddington luminosity can successfully explain the $M_{\rm bh} - \sigma$ and nuclear cluster mass $M_{\rm NC}-\sigma$ correlations, as well as why larger bulges host SMBHs while smaller ones host nuclear star clusters. However these models do not specify how SMBHs feed on infalling gas whilst simultaneously producing feedback that drives gas out of the galaxy. More complex models with rotation and/or anisotropic feedback allow SMBHs to feed via a disc or regions not exposed to SMBH winds, but in these more realistic cases it is not clear why a robust $M_{\rm bh} - \sigma$ relation should be established. In fact, some of the model predictions contradict observations. For example, an isotropic SMBH wind impacting on a disc (rather than a shell) of aspect ratio $H/R \ll 1$ requires the SMBH mass to be larger by a factor $\sim R/H$, which is opposite to what is observed. We conclude that understanding how a SMBH feeds is as important a piece of the puzzle as understanding how its feedback affects its host galaxy. Finally, we note that in aspherical cases the SMBH outflows induce differential motions in the bulge. This may pump turbulence that is known to hinder star formation in star forming regions. SMBH feedback thus may not only drive gas out of the bulge but also reduce the fraction of gas turned into stars.

We present new observations of the nuclear star cluster in the central parsec of the Galaxy with the adaptive optics assisted, integral field spectrograph SINFONI on the ESO/VLT. Our work allows the spectroscopic detection of early and late type stars to m_K >= 16, more than 2 magnitudes deeper than our previous data sets. Our observations result in a total sample of 177 bona fide early-type stars. We find that most of these Wolf Rayet (WR), O- and B- stars reside in two strongly warped disks between 0.8" and 12" from SgrA*, as well as a central compact concentration (the S-star cluster) centered on SgrA*. The later type B stars (m_K>15) in the radial interval between 0.8" and 12" seem to be in a more isotropic distribution outside the disks. The observed dearth of late type stars in the central few arcseconds is puzzling, even when allowing for stellar collisions. The stellar mass function of the disk stars is extremely top heavy with a best fit power law of dN/dm ~ m^(-0.45+/-0.3). Since at least the WR/O-stars were formed in situ in a single star formation event ~6 Myrs ago, this mass function probably reflects the initial mass function (IMF). The mass functions of the S-stars inside 0.8" and of the early-type stars at distances beyond 12" are compatible with a standard Salpeter/Kroupa IMF (best fit power law of dN/dm ~ m^(-2.15+/-0.3)).

It is now well established that many galaxies have nuclear star clusters (NCs) whose total masses correlate with the velocity dispersion (sigma) of the galaxy spheroid in a very similar way to the well--known supermassive black hole (SMBH) M - sigma relation. Previous theoretical work suggested that both correlations can be explained by a momentum feedback argument. Observations further show that most known NCs have masses < 10^8 Msun, while SMBHs frequently have masses > 10^8 Msun, which remained unexplained in previous work. We suggest here that this changeover reflects a competition between the SMBH and nuclear clusters in the feedback they produce. When one of the massive objects reaches its limiting M-sigma value, it drives the gas away and hence cuts off its own mass and also the mass of the ``competitor''. The latter is then underweight with respect to the expected M-sigma mass (abridged).

We compute the effect of an orbiting gas disc in promoting the coalescence of a central supermassive black hole binary. Unlike earlier studies, we consider a finite mass of gas with explicit time dependence: we do not assume that the gas necessarily adopts a steady state or a spatially constant accretion rate, i.e. that the merging black hole was somehow inserted into a pre--existing accretion disc. We consider the tidal torque of the binary on the disc, and the binary's gravitational radiation. We study the effects of star formation in the gas disc in a simple energy feedback framework. The disc spectrum differs in detail from that found before. In particular, tidal torques from the secondary black hole heat the edges of the gap, creating bright rims around the secondary. These rims do not in practice have uniform brightness either in azimuth or time, but can on average account for as much as 50 per cent of the integrated light from the disc. This may lead to detectable high--photon--energy variability on the relatively long orbital timescale of the secondary black hole, and thus offer a prospective signature of a coalescing black hole binary. We also find that the disc can drive the binary to merger on a reasonable timescale only if its mass is at least comparable with that of the secondary black hole, and if the initial binary separation is relatively small, i.e. $a_0 \lesssim 0.05$ pc. Star formation complicates the merger further by removing mass from the disc. In the feedback model we consider, this sets an effective limit to the disc mass. As a result, binary merging is unlikely unless the black hole mass ratio is $\la 0.001$. Gas discs thus appear not to be an effective solution to the `last parsec' problem for a significant class of mergers.

We present a new framework for radiation hydrodynamics simulations. Gas dynamics is modelled by the Smoothed Particle Hydrodynamics (SPH) method, whereas radiation transfer is simulated via a time-dependent Monte-Carlo approach that traces photon packets. As a first step in the development of the method, in this paper we consider the momentum transfer between radiation field and gas, which is important for systems where radiation pressure is high. There is no fundamental limitations on the number of radiation sources, geometry or the optical depth of the problems that can be studied with the method. However, as expected for any Monte-Carlo transfer scheme, stochastic noise presents a serious limitation. We present a number of tests that show that the errors of the method can be estimated accurately by considering Poisson noise fluctuations in the number of photon packets that SPH particles interact with per dynamical time. It is found that for a reasonable accuracy the momentum carried by photon packets must be much smaller than a typical momentum of SPH particles. We discuss numerical limitations of the code, and future steps that can be taken to improve performance and applicability of the method.

The central parsec around the super-massive black hole in the Galactic Center hosts more than 100 young and massive stars. Outside the central cusp (R~1") the majority of these O and Wolf-Rayet (WR) stars reside in a main clockwise system, plus a second, less prominent disk or streamer system at large angles with respect to the main system. Here we present the results from new observations of the Galactic Center with the AO-assisted near-infrared imager NACO and the integral field spectrograph SINFONI on the ESO/VLT. These include the detection of 27 new reliably measured WR/O stars in the central 12" and improved measurements of 63 previously detected stars, with proper motion uncertainties reduced by a factor of four compared to our earlier work. We develop a detailed statistical analysis of their orbital properties and orientations. Half of the WR/O stars are compatible with being members of a clockwise rotating system. The rotation axis of this system shows a strong transition as a function of the projected distance from SgrA*. The main clockwise system either is either a strongly warped single disk with a thickness of about 10 degrees, or consists of a series of streamers with significant radial variation in their orbital planes. 11 out of 61 clockwise moving stars have an angular separation of more than 30 degrees from the clockwise system. The mean eccentricity of the clockwise system is 0.36+/-0.06. The distribution of the counter-clockwise WR/O star is not isotropic at the 98% confidence level. It is compatible with a coherent structure such as stellar filaments, streams, small clusters or possibly a disk in a dissolving state. The observed disk warp and the steep surface density distribution favor in situ star formation in gaseous accretion disks as the origin of the young stars.

Young massive stars in the central parsec of our Galaxy are best explained by star formation within at least one, and possibly two, massive self-gravitating gaseous discs. With help of numerical simulations, we here consider whether the observed population of young stars could have originated from a large angle collision of two massive gaseous clouds at R approx. 1 parsec from Sgr A*. In all the simulations performed, the post-collision gas flow forms an inner, nearly circular gaseous disc and one or two eccentric outer filaments, consistent with the observations. Furthermore, the radial stellar mass distribution is always very steep, Sigma proportional to R^-2, again consistent with the observations. All of our simulations produce discs that are warped by between 30 to 60 degrees, in accordance with the most recent observations. The 3D velocity structure of the stellar distribution is sensitive to initial conditions (e.g., the impact parameter of the clouds) and gas cooling details. For example, the runs in which the inner disc is fed intermittently with material possessing fluctuating angular momentum result in multiple stellar discs with different orbital orientations, contradicting the observed data. In all the cases the amount of gas accreted by our inner boundary condition is large, enough to allow Sgr A* to radiate near its Eddington limit over approx. 10^5 years. This suggests that a refined model would have physically larger clouds (or a cloud and a disc such as the circumnuclear disc) colliding at a distance of a few parsecs rather than 1 parsec as in our simulations.

Young massive stars in the central parsec of our Galaxy are best explained by star formation within at least one, and possibly two, massive self-gravitating gaseous discs. With help of numerical simulations, we here consider whether the observed population of young stars could have originated from a large angle collision of two massive gaseous clouds at R ~ 1 pc from Sgr A*. In all the simulations performed, the post-collision gas flow forms an inner nearly circular gaseous disc and one or two eccentric outer filaments, consistent with the observations. Furthermore, the radial stellar mass distribution is always very steep, sigma proportional to R^-2, again consistent with the data. The 3D velocity structure of the stellar distribution is however sensitive to initial conditions (e.g., the impact parameter of the clouds) and gas cooling details. In all the cases the amount of gas accreted by our inner boundary condition is large, enough to allow Sgr A* to radiate near its Eddington limit during ~ 10^5 years. This suggests that a refined model would have physically larger clouds (or a cloud and a disc such as the CND) colliding at a distance of a few pc rather than 1 pc as in our simulations.

We present numerical simulations of stellar wind dynamics in the central parsec of the Galactic centre, studying in particular the accretion of gas on to Sgr A*, the super-massive black hole. Unlike our previous work, here we use state-of-the-art observational data on orbits and wind properties of individual wind-producing stars. Since wind velocities were revised upwards and non-zero eccentricities were considered, our new simulations show fewer clumps of cold gas and no conspicuous disc-like structure. The accretion rate is dominated by a few close `slow wind stars' v_w < 750 km/s, and is consistent with the Bondi estimate, but variable on time-scales of tens to hundreds of years. This variability is due to the stochastic in-fall of cold clumps of gas, as in earlier simulations, and to the eccentric orbits of stars. The present models fail to explain the high luminosity of Sgr A* a few hundred years ago implied by Integral observations, but we argue that the accretion of a cold clump with a small impact parameter could have caused it. Finally, we show the possibility of constraining the total mass-loss rate of the `slow wind stars' using near infra-red observations of gas in the central few arcseconds.

Steady state accretion discs larger than ~ 0.01-0.1 pc are known to be gravitationally unstable for the accretion rates needed to explain super-massive black hole (SMBH) activity. We propose that SMBH are fed by a succession of mass deposition events with randomly directed angular momenta. Because of incomplete angular momentum cancellation a warped accretion disc forms in the inner few parsec. The orientation of the disc performs a random walk. Deposition of new material promotes SMBH accretion at rates much faster than viscous. Observational implications of this picture include: (i) lighter accretion discs that can fuel AGN and quasars and yet avoid star formation at R >> 0.1 pc; (ii) star formation inside the disc is not a function of mass accretion rate only. It can take place at high or low accretion rates, e.g., when too few clouds arrive in the inner region. An example of this might be the central parsec of our Galaxy. (iii) The discs can form Compton-thick obscuring structures of ~ parsec size as required in AGN unification models; (iv) faster black hole growth resulting from misalignment of the disc and the black hole spin in the early Universe; (v) Isotropic deposition of SMBH energy and momentum feedback in the galaxy bulge. This may help explain the high efficiency with which it seems to be operating in the Universe. (vi) No correlation between SMBH activity and the presence of kiloparsec scale bars or gaseous discs in galactic bulges; (vii) Bodily collisions between gaseous components of merging galaxies facilitate production of gas streams feeding the centre of the combined galaxy. Mergers should thus be catalysts of SMBH growth. (viii) Conversely, galaxies experiencing fewer mergers are more likely to form massive nuclear star clusters than feed their SMBHs.

The soft excess seen in many AGN is most probably due to partially ionized material moving at relativistic speeds close to the black hole. There are currently two potential geometries for this material, one where it is out of the line of sight, seen via reflection, e.g. the accretion disc, the other where it is in the line of sight, seen in absorption e.g. a wind above the disc. Both models require apparent fine tuning of the ionization parameter of this material in order to produce the large jump in opacity at $\sim 0.7$ keV associated with OVII/VIII, as required to make the soft excess. However, Chevalier et al (2006) show that these states rather naturally dominate the absorption spectrum for soft X-ray spectra if the illuminated material is (at least approximately) in pressure balance. Here we explore whether hydrostatic pressure equilibrium in a disc can likewise naturally select the required ionization states in reflection. We find the opposite. The soft excess X-ray excess is much weaker in the hydrostatic models than it is in the constant density models. Since even the constant density models cannot fit the largest soft excesses seen without the intrinsic continuum being hidden from view, this means that reflection from a hydrostatic disc cannot realistically match the data. Even if the disc structure is instead more like a constant density atmosphere, the required fine--tuning of the ionization parameter still remains a problem for reflection models.

Precise mechanisms by which Active Galactic Nuclei (AGN) receive their gaseous fuel is still a mystery. Here I draw attention to the extra ordinary star formation event that took place in the central ~ 0.5 parsec of our Galaxy. The most reliable explanation of the event seems to be that two somewhat massive nearly co-eval gaseous disks failed to accrete on Sgr A*, the super-massive black hole (SMBH) in our Galaxy, and instead cooled down and gravitationally collapsed, forming the stars observed now. This emphasises that star formation must be an important part of AGN feeding puzzle. I also discuss a model in which stellar winds create the observed obscuration of AGN. These winds are cold, clumpy and dusty, as required by the observations, but they are Compton-thin unless wind outflow rate is highly super-Eddington. This argument is in fact a general one, independent of the wind driving mechanism. I thus suggest that winds may be important for optically thin absorbers, and that a better model for optically thick AGN obscuration is a warped accretion/star forming disk.

We numerically model fragmentation of a gravitationally unstable gaseous disc under conditions that may be appropriate for the formation of the young massive stars observed in the central parsec of our Galaxy. In this study, we adopt a simple prescription with a locally constant cooling time. We find that, for cooling times just short enough to induce disc fragmentation, stars form with a top-heavy Initial Mass Function (IMF), as observed in the Galactic Centre (GC). For shorter cooling times, the disc fragments much more vigorously, leading to lower average stellar masses. Thermal feedback associated with gas accretion onto protostars slows down disc fragmentation, as predicted by some analytical models. We also simulate the fragmentation of a gas stream on an eccentric orbit in a combined Sgr A* plus stellar cusp gravitational potential. The stream precesses, self-collides and forms stars with a top-heavy IMF. None of our models produces large enough co-moving groups of stars that could account for the observed ``mini star cluster'' IRS13E in the GC. In all of the gravitationally unstable disc models that we explored, star formation takes place too fast to allow any gas accretion onto the central super-massive black hole. While this can help to explain the quiescence of `failed AGN' such as Sgr A*, it poses a challenge for understanding the high gas accretion rates infered for many quasars.

We summarize some of the new results from contributions made to the Galactic Centre Workshop that took place in Bad Honnef, Germany, on April 18-22, 2006.

Due to dynamical friction stellar mass black holes and neutron stars are expected to form high density cusps in the inner parsec of our Galaxy. These compact remnants, expected to number around 20000, may be accreting cold dense gas present there, and give rise to potentially observable X-ray emission. Here we build a simple but detailed time-dependent model of such emission. We find that at least several X-ray sources of this nature should be detectable with Chandra at any one time. Turning this issue around, we also ask a question of what current observational constraints might be telling us about the total number of compact remnants. A cusp with ~ 40 thousand black holes over-predicts the number of discrete sources and the total X-ray luminosity of the inner parsec, and is hence ruled out. Future observations of the distribution and orbits of the cold ionised gas in the inner parsec of Sgr A* will put tighter constraints on the cusp of compact remnants.

Straight-forward models of X-ray reflection in the inner region of accretion discs predict that primary X-ray flux and the flux reflected off the surface of the disc should vary together, albeit a short light travel time delay. Most of the observations, however, show that the X-ray flux can vary while the reflected features remain constant. Here we propose a simple explanation to this. In all likelihood, the emission of a moderately optically thick magnetic flare atop an accretion disc is anisotropic. A constant energy release rate in a flare will appear to produce a variable X-ray flux as the flare rotates with the accretion disc anchoring the magnetic tube. The reflector, on the other hand, receives a constant X-ray flux from the flare. Since the reflected emission is azimuthally symmetric, the observer will see a roughly constant reflected flux (neglecting relativistic effects). The model does not produce quasi-periodic oscillations (QPO) if magnetic flux tubes are sheared out faster than they complete one orbit.

Compact remnants -- stellar mass black holes and neutron stars formed in the inner few parsec of galactic centres are predicted to sink into the central parsec due to dynamical friction on low mass stars, forming a high concentration cusp (Morris 1993). Same physical region may also contain very high density molecular clouds and accretion discs that are needed to fuel SMBH activity. Here we estimate gas capture rates onto the cusp of stellar remnants, and the resulting X-ray luminosity, as a function of the accretion disc mass. At low disc masses, most compact objects are too dim to be observable, whereas in the high disc case most of them are accreting at their Eddington rates. We find that for low accretion disc masses, compact remnant cusps may be more luminous than the central SMBHs. This ``diffuse'' emission may be of importance for local moderately bright AGN, especially Low Luminosity AGN. We also briefly discuss how this expected emission can be used to put constraints on the black hole cusp near our Galactic Centre.

Momentum feedback from super-Eddington accretion offers a simple explanation for the observed M-sigma and M-M_{spheroid} relations between supermassive black holes (SMBH) and the spheroids of their host galaxies. Recently Ferrarese et al. and Wehner & Harris observed analogous relations between the masses of central star clusters and their hosts. We show that stellar winds and supernovae from such nuclear clusters (NC) give similar feedback explanations for this case also, and we discuss the connection to the Faber-Jackson relation for the spheroids themselves.

Observational evidence accumulated over the past decade indicates that accretion discs in X-ray binaries are viscously stable unless they accrete very close to the Eddington limit. This is at odds with the most basic standard accretion disc theory, but could be explained by either having the discs to be much cooler whereby they are not radiation pressure dominated, or by a more sophisticated viscosity law. Here we argue that the latter is taking place in practice, on the basis of a stability analysis that assumes that the magneto-rotational-instability (MRI) responsible for generating the turbulent stresses inside the discs is also the source for a magnetically dominated corona. We show that observations of stable discs in the high/soft states of black hole binaries, on the one hand, and of the strongly variable microquasar GRS 1915+105 on the other, can all be explained if the magnetic turbulent stresses inside the disc scale proportionally to the geometric mean of gas and total pressure with a constant of proportionality (viscosity parameter) having a value of a few times 10^{-2}. Implications for bright AGN are also briefly discussed

The role of star-formation driven outflows in the obscuration of the central source in the Active Galactic Nuclei (AGN) is discussed. The outflow from a sub-parsec scale accretion disc is numerically modelled for parameters appropriate to the Galactic Centre. The resulting obscuration pattern is very patchy, with some lines of sight becoming optically thick to Thomson scattering. A fixed observer would see column depth changing by factors of many over time scales of order months to hundreds of years, depending on the physical size of the outflow region. Such winds may be relevant for obscuration of some AGN and especially "changing look AGN". However, averaged over the sky as seen from the central source, these winds are always optically thin unless wind outflow rates are super-Eddington. A simple scaling argument shows that this is true not only for stellar-driven winds but for any AGN winds. We therefore conclude that AGN winds are unable to account for the vast majority of optically thick obscured AGN (a significant fraction of all AGN). We suggest that the most likely source of optically thick obscuration in AGN is a warped parsec scale accretion disc.

We report the definite spectroscopic identification of 41 OB supergiants, giants and main sequence stars in the central parsec of the Galaxy. Detection of their absorption lines have become possible with the high spatial and spectral resolution and sensitivity of the adaptive optics integral field spectrometer SPIFFI/SINFONI on the ESO VLT. Several of these OB stars appear to be helium and nitrogen rich. Almost all of the ~80 massive stars now known in the central parsec (central arcsecond excluded) reside in one of two somewhat thick (<|h|/R>~0.14) rotating disks. These stellar disks have fairly sharp inner edges (R~1") and surface density profiles that scale as R^{-2}. We do not detect any OB stars outside the central 0.5 pc. The majority of the stars in the clockwise system appear to be on almost circular orbits, whereas most of those in the `counter-clockwise' disk appear to be on eccentric orbits. Based on its stellar surface density distribution and dynamics we propose that IRS 13E is an extremely dense cluster (core density > 3x10^8 sunmass/pc^3), which has formed in the counter-clockwise disk. The stellar contents of both systems are remarkably similar, indicating a common age of ~6+/-2 Myr. The K-band luminosity function of the massive stars suggests a top-heavy mass function and limits the total stellar mass contained in both disks to ~1.5x10^4 sunmass. Our data strongly favor in situ star formation from dense gas accretion disks for the two stellar disks. This conclusion is very clear for the clockwise disk and highly plausible for the counter-clockwise system.

We consider the structure of self-gravitating marginally stable accretion disks in galactic centers in which a small fraction of the disk mass has been converted into proto-stars. We find that proto-stars accrete gaseous disk matter at prodigious rates. Mainly due to the stellar accretion luminosity, the disk heats up and geometrically thickens, shutting off further disk fragmentation. The existing proto-stars however continue to gain mass by gas accretion. As a results, the initial mass function for disk-born stars at distances R ~ 0.03-3 parsec from the super-massive black hole should be top-heavy. The effect is most pronounced at around R ~ 0.1 parsec. We suggest that this result explains observations of rings of young massive stars in our Galaxy and in M31, and predict that more of such rings will be discovered.

It is believed that young massive stars orbiting Sgr A* in two stellar discs on scales of 0.1-0.2 parsecs were formed either farther out in the Galaxy and then quickly migrated inward, or in situ in a massive self-gravitating disc. Comparing N-body evolution of stellar orbits with observational constraints, we set upper limits on the masses of the two stellar systems. These masses turn out to be few times lower than the expected total stellar mass estimated from the observed young high-mass stellar population and the standard galactic IMF. If these stars were formed in situ, in a massive self-gravitating disc, our results suggest that the formation of low-mass stars was suppressed by a factor of at least a few, requiring a top-heavy initial mass function (IMF) for stars formed near sgr A*.

(ABRIDGED) We present in detail our new 3D numerical models for the accretion of stellar winds on to Sgr A*. In our most sophisticated models, we put stars on realistic orbits around Sgr A*, include `slow' winds (300 km/s), and account for radiative cooling. We first model only one phase `fast' stellar winds (1000 km/s). For wind sources fixed in space, the accretion rate is Mdot ~ 1e-5 Msun/yr, fluctuates by < 10%, and is in a good agreement with previous models. In contrast, Mdot decreases by an order of magnitude for stars following circular orbits, and fluctuates by ~ 50%. Then we allow a fraction of stars to produce slow winds. Much of these winds cool radiatively, forming cold clumps immersed into the X-ray emitting gas. We test two orbital configurations for the stars in this scenario, an isotropic distribution and two rotating discs with perpendicular orientation. The morphology of cold gas is quite sensitive to the orbits. In both cases, however, most of the accreted gas is hot, with an almost constant Mdot ~ 3e-6 Msun/yr, consistent with Chandra observations. The cold gas accretes in intermittent, short but powerful episodes which may give rise to large amplitude variability in the luminosity of Sgr A* on time scales of 10s to 100s of years. The circularisation radii for the flows are ~ 1e3 and 1e4 Rsch, for the one and two-phase wind simulations, respectively, never forming the quasi-spherical accretion flows suggested in some previous work. Our work suggests that, averaged over time scales of 100s to 1000s of years, the radiative and mechanical luminosity of Sgr A* may be substantially higher than it is in its current state. Further improvements of the wind accretion modelling of Sgr A* will rely on improved observational constraints for the wind properties and stellar orbits.

Few dozens of young high mass stars orbit Sgr A* at distances as short as 0.1 parsec, where star formation should be quenched by the strong tidal shear from Sgr A*. The puzzling young stellar population is believed to come into existence in one of the two ways: (i) "normal" star formation at several tens of parsec in a very massive star cluster that then spiraled in, or (ii) star formation in situ in a massive self-gravitating disk. We propose to constrain these two scenarios via the expected X-ray emission from young low mass stars that should have formed alongside the massive stars. To this end we compare the recent Chandra observations of X-ray emission from young stars in the Orion Nebula, and the Chandra observations of Sgr A* field. We show that the cluster inspiral model is ruled out irrespectively of the initial mass function (IMF) of the young stars. In addition, for the in situ model, we find that no more than few thousand low-mass stars could have formed alongside the massive stars. This is more than a factor of ten fewer than expected if these stars were formed with the standard IMF as elsewhere in the Galaxy. The young stars in the GC are thus the first solid observational evidence for star formation in AGN disks and also require the IMF of these stars to be top-heavy. We briefly consider implication of these results for AGN in general.

Sgr A* is currently being fed by winds from a cluster of gravitationally bound young mass-loosing stars. Using observational constraints on the orbits, mass loss rates and wind velocities of these stars, we numerically model the distribution of gas in the ~ 0.1--10'' region around Sgr A*. We find that radiative cooling of recently discovered slow winds leads to the formation of many cool filaments and blobs, and to a thin and rather light accretion disc of about an arcsecond scale. The disc however does not extend all the way to our inner boundary. Instead, hot X-ray emitting gas dominates the inner arcsecond. In our simulations, cool streams of gas frequently enter this region on low angular momentum orbits, and are then disrupted and heated up to the ambient hot gas temperature. The accreting gas around Sgr A* is thus two-phase, with a hot component, observable at X-ray wavelengths, and a cool component, which may be responsible for the majority of time variability of Sgr A* emission on hundred and thousand years time-scales. We obtain an accretion rate of a few times 10^-6 msun/year, consistent with Chandra estimates, but variable on time-scales even shorter than hundred years. These results strongly depend on the chosen stellar orbits and wind parameters. Further observational input is thus key to a better modelling of Sgr A* wind accretion.

Orientation of parsec-scale accretion disks in AGN is likely to be nearly random for different black hole feeding episodes. Since AGN accretion disks are unstable to self-gravity on parsec scales, star formation in these disks will create young stellar disks, similar to those recently discovered in our Galactic Center. The disks blend into the quasi-spherical star cluster enveloping the AGN on time scales much longer than a likely AGN lifetime. Therefore, the gravitational potential within the radius of the black hole influence is at best axi-symmetric rather than spherically symmetric. Here we show that as a result, a newly formed accretion disk will be warped. For the simplest case of a potential resulting from a thin stellar ring, we calculate the disk precession rates, and the time dependent shape. We find that, for a realistic parameter range, the disk becomes strongly warped in few hundred orbital times. We suggest that this, and possibly other mechanisms of accretion disk warping, have a direct relevance to the problem of AGN obscuration, masing warped accretion disks, narrow Fe K-alpha lines, etc.

The extremely hot and tenuous accretion flow in the immediate vicinity of Sgr A* is believed to be invisible (too dim) in the X-ray band, except for short X-ray flares. Here we point out that during pericenter passages, close brightest stars irradiate the inner region of the accretion flow, providing a plenty of optical/UV photons. These seed photons are Compton up-scattered by the hot electrons of the accretion flow to higher frequencies, some into the X-ray band, potentially making the innermost accretion flow much brighter in X-rays than usual. We propose to use coordinated near infra-red and X-ray observations of close star passages to put constraints onto Sgr A* accretion theories. The absence of a noticeable change in the steady emission of Sgr A* as observed by Chandra in the year 2002, when the star named S2 passed through a pericenter of its orbit, already rules out the hotter of the ``standard'' Advection-Dominated Accretion Flows. The less dense accretion flows, in particular the model of Yuan et al. (2003), passes the test and is constrained to accretion rates no larger than $\sim 10^{-7}$ Solar masses per year.

Sgr A* is extra-ordinarily dim in all wavelengths requiring a very low accretion rate at the present time. However, at a radial distance of a fraction of a parsec from Sgr A*, two rings populated by young massive stars suggest a recent burst of star formation in a rather hostile environment. Here we explore two ways of creating such young stellar rings with a gaseous accretion disk: by self-gravity in a massive disk, and by capturing "old" low mass stars and growing them via gas accretion in a disk. The minimum disk mass is above 10^4 Msun for the first mechanism and is few tens times larger for the second one. The observed relatively small velocity dispersion of the stars rules out disks more massive than around 10^5 Msun: heavier stellar or gas disks would warp each other by orbital precession in an axisymmetric potential too strongly. The capture of "old" stars by a disk is thus unlikely as the origin of the young stellar disks. The absence of a massive nuclear gas disk in Sgr A* now implies that the disk was either accreted by the SMBH, which would then imply almost a quasar-like luminosity for Sgr A*, or was consumed in the star formation episode. The latter possibility appears to be more likely on theoretical grounds. We also consider whether accretion disk plane changes, expected to occur due to fluctuations in the angular momentum of gas infalling into the central parsec of a galaxy, would dislodge the embedded stars from the disk midplane. We find that the stars leave the disk midplane only if the disk orientation changes on time scales much shorter than the disk viscous time.

Quasar accretion disks are believed to form stars by self-gravity. Low Luminosity Active Galactic Nuclei (LLAGN) are much dimmer galactic centers, and are often believed to be quasars that ran out of gaseous fuel. LLAGN accretion disks should thus co-exist with thousands to millions of stars or proto-stars left from the previous stronger accretion activity. In principle, these stars may produce several important effects: (i) contribute to the optical/UV spectra of some LLAGN; (ii) reprocessing of the stellar radiation in the dusty disks could dominate the LLAGN infra-red spectra; (iii) deplete the (accretion) gas disk much faster than it can accrete onto the supper-massive black hole (SMBH); (iv) stars, individually or in groups, may slow down and modulate the accretion flow significantly due to their inertia. In this way they may produce the LLAGN cut-off disks; (v) alternatively, frequent enough stellar collisions and resulting stellar disruptions could keep the inner disk empty. Here we explore these ideas. We find that, despite ``low'' luminosities of LLAGN, unrealistically high stellar densities are required to make a sizable radiative contribution to the (HST) optical/UV spectra of these galactic nuclei. Stellar contribution to the infrared spectrum is more likely. Further, if LLAGN are in a quasi steady-state for as long as 10^7 years or more, too high stellar densities would again be required to significantly affect the dynamics of accretion flow. However, if LLAGN are ``short''-lived phenomena, e.g. t < 10^5 years, quiescent states of quiescence-outburst cycles, then embedded stars may be much more important through the mass effects (iii) -- (v). With observations of LLAGN becoming progressively better, it will be more and more difficult to neglect the presence of close stars in and around nuclear accretion disks.

We study the accretion flow of a hot gas captured by the black hole gravity in the presence of a thin cold accretion disk. Such geometrical arrangement is expected in Active Galactic Nuclei (AGN) and in galactic X-ray binary systems because both hot and cold gases are present in the black hole vicinity. Previous astrophysical literature concentrated on the evaporation of the cold disk in the classical heat conduction limit. Here we consider the inverse process, i.e. condensation of the hot gas onto the cold disk. We find two distinct condensation regimes. (i) In the classical thermal conduction limit, the radiative cooling in the hot gas itself force condensation above a certain critical accretion rate. Most of the flow energy in this case is re-emitted as X-ray radiation. (ii) Below a certain minimum accretion rate, the hot electrons are collisionless and the classical heat flux description becomes invalid. We use the ``non-local'' heat flux approach borrowed from the terrestrial laser heated plasma experiments. Due to their very large mean free path, the hot particles penetrate deep into the cold disk where the radiative losses are significant enough to enable condensation. In this case the hot flow energy is inconspicuously re-radiated by the transition layer in many UV and especially optical recombination lines (e.g., Ly$\alpha$, $H\alpha$, H$\beta$) as well as via the optically thick disk emission. We derive an approximate analytical solution for the dynamics of the hot condensing flow. If the cold disk is inactive, i.e. accumulating mass for a future accretion outburst, then the two-phase flows appear radiatively inefficient. These condensing solutions may be relevant to \sgra, low luminosity AGN, and transient binary accreting systems in quiescence.

Cold inactive disks are believed to exist in Low Luminosity AGN (LLAGN). They may also exist in the nuclei of inactive galaxies and in the center of our own Galaxy. These disks would then be embedded in the observed dense nuclear stellar clusters. Making the simplest assumption of an optically thick disk, we explore several ways to detect the disk presence through its interaction with the cluster. First of these is the eclipses of close bright stars by the disk. The second is the increase in the infrared flux of the disk due to illumination of its surface by such stars during close passages. Finally the surface brightness of the star cluster should show an anisotropy that depends on the inclination angle of the disk. We apply the first two of the methods to Sgr A*, the super massive black hole in our Galactic Center. Using the orbital parameters of the close star S2, we strongly rule out a disk optically thick in the near infrared unless it has a relatively large inner hole. For disks with no inner holes, we estimate that the data permit a disk with infrared optical depth no larger than about 0.01. Such a disk could also be responsible for the detected 3.8 microns excess in the spectrum of S2. The constraints on the disk that we obtain here can be reconciled with the disk parameters needed to explain the observed X-ray flares if dust particles in the disk have sizes greater than ~ 30 microns. The destruction of small dust particles by strong UV heating and shocks from star passages through the disk, and grain growth during "quiescent" times, are mentioned as possible mechanisms of creating the unusual grain size distribution. (...)

Sgr A*, the putative black hole in our Galactic Center (GC), is extraordinary dim in all frequencies. Apparently the black hole is unable to accrete at the Bondi accretion rate for some reason. Another mystery of Sgr A* is the recently discovered large magnitude and short duration X-ray flares (Baganoff et al. 2001). Here we point out that X-ray flares should be expected from star passages through an inactive (i.e. formerly accreting) disk. There are thousands of stars in Sgr A* stellar cluster, and each star will pass through the disk twice per orbit. A shock hot enough to emit X-rays is formed around the star when the latter is inside the disk. We develop a simple yet somewhat detailed model for the X-ray emission from the shock. The duration of the flares, their X-ray spectra, frequency of events, and weakness of emission in the radio and near infra-red appear to be consistent with the available observations of X-ray flares. We therefore suggest that at least a fraction of the observed flares is due to the star-disk passages. Such star-disk flares are also of interest in the nuclei of nearby inactive galaxies, especially in connection with perspective Constellation-X and XEUS X-ray missions.

The black hole in our Galactic Center is extremely underluminous for the amount of hot gas available for accretion. Theoretical understanding of this fact rests on a likely but not entirely certain assumption that the electrons in the accreting gas are much cooler than the protons. In this case the hot gas as a whole is too hot to accrete, and is too tenous to radiate away its gravitational energy. Here we propose a drastically different picture of the accretion process in Sgr A* not based on the unchecked two-temperature assumption. Namely, we argue that there should exist a very cold inactive disk -- a remnant of a past stronger accretion activity in Sgr A*. Such a disk would be a very efficient cooling surface for the hot flow. We show that under certain conditions the cooling due to thermal conduction cannot be balanced by the viscous heating in the hot flow. Along with the heat, the hot flow loses its viscosity and thus ability to accrete. It settles (condences) onto the cold disk slighlty inside of the circularization radius. If the latter is very large, then the liberated energy, and the luminosity emitted, is orders of magnitude less than naively expected. We build a simple analytical model for this flow and calculate the expected spectra that appear to be in a very reasonble agreement with observations. Strong additional support for the presence of the inactive disk comes from the recent observations of X-ray flares in Sgr A*. The properties of these flares are very similar to those produced by stars passing through a cold disk.

A cold neutral and extremely dim accretion disk may be present as a remnant of a past vigorous activity around the black hole in our Galactic Center (GC). Here we discuss ways to detect such a disk through its interaction with numerous stars present in the central ~0.1 parsec of the Galaxy. The first major effect expected is X-ray and near infrared (NIR) flares arising when stars pass through the disk. The second is eclipses of the stars by the disk. We point out conditions under which the properties of the expected X-ray flares are similar to those recently discovered by Chandra. Since orbits of bright stars are now being precisely measured, the combination of the expected flares and eclipses offers an invaluable tool for constraining the disk density, size, plane and even direction of rotation. The winds of the O-type stars are optically thick to free-free absorption in radio frequencies. If present near Sgr A* core, such powerful stellar winds can modulate and even occult the radio source.

We fit the broad-band RXTE PCA and HEXTE spectrum from 3--200 keV with reflection models which calculate the vertical ionization structure of an X-ray illuminated disc. We consider two geometries corresponding to a truncated disc/inner hot flow and magnetic flares above an untruncated disc. Both models are able to fit the PCA 3--20 keV data, but with very different spectral components. In the magnetic flare models the 3--20 keV PCA spectrum contains a large amount of highly ionized reflection while in the truncated disc models the amount of reflection is rather small. The Compton downscattering rollover in reflected emission means that the magnetic flare models predict a break in the spectrum at the high energies covered by the HEXTE bandpass which is not seen. By contrast the weakly illuminated truncated disc models can easily fit the 3--200 keV spectra.

We model the reflected spectrum expected from localized magnetic flares above an ionised accretion disk. We concentrate on the case of very luminous magnetic flares above a standard accretion disk extending down to the last stable orbit, and use a simple parameterisation to allow for an X-ray driven wind. Full disk spectra including relativistic smearing are calculated. When fit with the constant density reflection models, these spectra give both a low reflected fraction and a small line width as seen in the hard spectra from Galactic Black Hole Binaries and Active Galactic Nuclei. We fit our calculated spectra to real data from the low/hard state of Nova Muscae and Cyg X-1 and show that these models give comparable $\chi^2$ to those obtained from the constant density reflection models which implied a truncated disk. This explicitally demonstrates that the data are consistent with either magnetic flares above an ionized disk extending down to the last stable orbit around a black hole, or with non-ionized, truncated disks.

We perform a first study of time-dependent X-ray reflection in photo-ionized accretion disks. We assume a step-functional change in the X-ray flux and use a simplified prescription to describe the time evolution of the illuminated gas density profile in response to changes in the flux. We find that the dynamical time for re-adjustment of the hydrostatic balance is an important relaxation time scale of the problem since it affects evolution of the ionization state of the reflector. Because of this the Fe K-alpha line emissivity depends on the shape and intensity of the illuminating flux in prior times, and hence it is not a function of the instantaneous illuminating spectrum. Moreover, during the transition, a prominent Helium-like component of the Fe K-alpha line may appear. As a result, the line flux may appear to be completely uncorrelated with X-ray continuum flux on time scales shorter than the dynamical time. In addition, the time-dependence of the illuminating flux may leave imprints even on the time-averaged line spectra, which may be used as an additional test of accretion disk geometry. Our findings appear to be important for the proposed Fe K-alpha line reverberation studies in lamppost-like geometries for accretion rates exceeding about $\sim 1%$ of the Eddington value. However, most AGN do not show Helium-like lines that are prominent in such models, probably indicating that these models are not applicable to real sources.

We present the spectral analysis of an early XMM-Newton observation of the luminous Seyfert 1 galaxy Markarian 509. We find the hard (2-10 keV) continuum slope, including reflection, to be somewhat flatter ($\Gamma=1.75$) than for a typical BLS1. The most obvious feature in the hard X-ray spectrum is a narrow emission line near 6.4 keV, with an equivalent width of 50 eV. The energy and strength of this line is consistent with fluorescence from `neutral' iron in the molecular torus, and we note the emerging ubiquity of this feature in XMM-Newton and Chandra observations of Seyfert 1 galaxies over a wide luminosity range. We also find evidence for a second emission line at 6.7-6.9 keV, which we attempt to model by reflection from a highly ionised disc. A `soft excess', evident as an upward curvature in the continuum flux below 1.5 keV, cannot be explained solely by enhanced reflection from the ionised disc. The RGS spectrum shows only weak discrete emission and absorption features in the soft X-ray spectrum, supporting our conclusion that the soft excess emission in Mkn 509 represents the high energy portion of optically thick, thermal emission from the inner accretion disc.

It is currently believed that the ``standard'' accretion disk theory under-predicts the observed X-ray luminosity from Soft X-ray Transients (SXT) in quiescence by as much as 4 to 6 orders of magnitude. This failure of the standard model is considered to be an important argument for the existence of the alternative mode of accretion -- Advection Dominated Accretion Flows (ADAF) in astrophysics, since these flows allow a much higher level of X-ray emission in quiescence, in agreement with the observations. Here we point out that, in stark contrast to steady-state standard disks, such disks in quiescence (being non-steady) produce most of the X-ray emission very far from the last stable orbit. Taking this into account, these disks can accommodate the observed X-ray luminosities of SXTs rather naturally. Our theory predicts that Fe K-alpha lines from standard accretion disks in quiescence should be narrow even though the cold disk goes all the way down to the last stable orbit.

We point out that luminous magnetic flares, thought to occur in standard AGN accretion disks, cannot be located much higher than few pressure scale heights above the disk. Using this fact, we estimate the fraction of the disk surface illuminated by a typical flare. This fraction turns out to be very small for geometrically thin disks, which implies that the instantaneous Fe K-alpha emission line from a specific magnetic flare is narrow. The line is red- or blue-shifted depending on the position of the observer relative to the flare and sweeps across the line band with time. We present several examples of theoretical time-resolved line profiles from such flares for a non-rotating black hole. The observations of such moving features with future X-ray telescopes will present a powerful test of the accretion disk geometry and may also test General Relativity in the strong field limit.

Reflection of X-rays from cool material around a black hole is one of the few observational diagnostics of the accretion flow geometry. Models of this reflected spectrum generally assume that the accretion disk can be characterized by material in a single ionization state. However, several authors have recently stressed the importance of the classic ionization instability for X-ray irradiated gas in hydrostatic balance. This instability leads to a discontinuous transition in the vertical structure of the disk, resulting in a hot ionized skin above much cooler material. If the Compton temperature of the skin is high then even iron is completely ionized, and the skin does not produce any spectral features. These new models, where the ionization structure of the disk is calculated self-consistently, require an excessive amount of computing power and so are difficult to use in directly fitting observed X-ray spectra. Instead, we invert the problem by simulating X-ray spectra produced by the new reflection models, and then fit these with the old, single zone reflection models, to assess the extent to which the derived accretion geometry depends on the reflection model used. We find that the single zone ionization models can severely underestimate the covering fraction of the ``cold'' material as seen from the X-ray source if the optical depth in the ionized skin is of order unity, and that this can produce an apparent correlation between the covering fraction and the X-ray spectral index similar in nature to that reported by Zdziarski, Lubinski and Smith (1999).

We test the "reprocessing paradigm" of the optical -- UV AGN variability by detailed modeling of the correlated X-ray -- optical (3590 and 5510 Angstrem) variability of the recent multiwavelength campaign of NGC 3516. We produce model optical light curves by convolving the observed X-ray flux with the response function of an infinite, thin accretion disk, illuminated by a point-like X-ray source at a given height above the compact object (the lamp-post model) and compare their properties to those observed. Special attention is given to the correct computation of the X-ray albedo of the disk. We further compute the X-ray reflection response at two energies (E = 1 and 20 keV) and argue for the possibility of hard lags in their cross spectra. We also compute the continuum Optical -- UV and the X-ray reflection spectra as well as the Fe K-alpha fluorescent line profiles which we also compare to observations. Despite the large (~50 percent) amplitude excursions of the X-ray flux, the model optical light curves exhibit variability amplitudes of 3 -- 4 percent, not unlike those observed. However, the model CCF between the X-ray and the model optical variations show clear lags of 0.1 and 0.25 days for black hole masses M = 10^7 and 10^8 Solar masses, respectively, not apparent in the data. The synchrony of X-ray - optical variations points toward the BH mass < 10^7 Solar masses, whereas the X-ray spectra prefer BH mass > 10^8. Our conclusion is that the combination of the observed optical/UV/X-ray spectral and timing observations are inconsistent with the lamp-post model geometry for NGC 3516.

X-ray illumination of accretion disks is an invaluable diagnostic of the structure of these disks because of the associated iron K$\alpha$ emission. Here we point out that the resulting reflected spectra depend very sensitively on the geometry of the X-ray source, and that this fact can be efficiently used to test these models observationally. In particular, we discuss three different accretion disk geometries: the ``lamppost model'', accretion disks with magnetic flares, and the model with a full corona overlying a cold thin disk. We show that in the case of the lamppost model, unless the X-ray luminosity of the central source is larger than that of the cold disk by a factor of 10 or more, a significant fraction of iron in the ionized skin of the disk is in the hydrogen and helium-like ions. Because these ions have large fluorescence yields, the resulting reflected spectra look strongly ionized, with Equivalent Width (EW) of the line {\em increasing} with X-ray luminosity $L_x$ up to the maximum of $\sim 500$ eV. This situation contrasts to the magnetic flare model, where the large X-ray flux near flares completely ionizes the skin of the disk and thus the resulting spectra appear to be that from a neutral material. The line EW in this model {\em anti-correlates} with X-ray luminosity, and becomes arbitrarily small when $L_x$ is a good fraction of the Eddington luminosity. Finally, in the full corona case, due to the additional pressure and weight of the corona, the gas pressure (and its density) below the corona is always large enough to make the gas very cool and effectively neutral. No highly ionized skin forms in such a model. If the corona is Thomson thin, then EW of the line does not depend on the accretion disk or corona luminosities for the full corona model.

We present full disk X-ray reflection spectra for two currently popular accretion flow geometries for AGN -- the lamppost model frequently used to discuss the iron line reverberation in AGN, and the model where the X-rays are produced in magnetic flares above a cold accretion disk (AD). The lamppost spectra contain several spectroscopic features characteristic of highly ionized material that are not seen in the X-ray spectra of most AGN. The magnetic flare model, on the other hand, produces reflected spectra that are roughly a super-position of a power-law and a {\em neutral-like} reflection and iron K$\alpha$ line, and are thus more in line with typical AGN X-ray spectra. Furthermore, because of the difference in the ionization structure of the illuminated material in the two models, the line equivalent width increases with the X-ray luminosity, $L_x$, for the lamppost, and decreases with $L_x$ for the flare model. In light of these theoretical insights, recent iron line reverberation studies of AGN, the X-ray Baldwin effect, and the general lack of X-ray reflection features in distant quasars all suggest that, for high accretion rates, the cold accretion disk is covered by a Thomson thick, completely ionized skin. Because the latter is only possible when the X-rays are concentrated to small emitting regions, we believe that this presents a strong evidence for the magnetic flare origin of X-rays in AGN.

We present a simple analytical formula for the Thomson depth of the X-ray heated skin of accretion disks valid at any radius and for a broad range of spectral indices of the incident X-rays, accretion rates and black hole masses. We expect that this formula may find useful applications in studies of geometry of the inner part of accretion flows around compact objects, and in several other astrophysically important problems, such as the recently observed X-ray ``Baldwin'' effect (i.e., monotonic decrease of Fe line's equivalent width with the X-ray luminosity of AGN), the problem of missing Lyman edge in AGN, and line and continuum variability studies in accretion disks around compact objects. We compute the reflected X-ray spectra for several representative cases and show that for hard X-ray spectra and large ionizing fluxes the skin represents a perfect mirror that does not produce any Fe lines or absorption features. At the same time, for soft X-ray spectra or small ionizing fluxes, the skin produces very strong ionized absorption edge and highly ionized Fe lines that should be observable in the reflected spectra.

We study the X-ray illumination of an accretion disk. We relax the simplifying assumption of constant gas density used in most previous studies; instead we determine the density from hydrostatic balance. It is found that the thermal ionization instability prevents the illuminated gas from attaining temperatures at which the gas is unstable. In particular, the uppermost layers of the X-ray illuminated gas are found to be almost completely ionized and at the local Compton temperature ($\sim 10^7 - 10^8$ K); at larger depths, the gas temperature drops abruptly to form a thin layer with $T\sim 10^6$ K, while at yet larger depths it decreases sharply to the disk effective temperature. We find that most of the Fe K$\alpha$ line emission and absorption edge are produced in the coolest, deepest layers, while the Fe atoms in the hottest, uppermost layers are generally almost fully ionized, hence making a negligible contribution to reprocessing features in $\sim 6.4-10$ keV energy range. We provide a summary of how X-ray reprocessing features depend on parameters of the problem. The results of our self-consistent calculations are both quantitatively and qualitatively different from those obtained using the constant density assumption. Therefore, we conclude that X-ray reflection calculations should always utilize hydrostatic balance in order to provide a reliable theoretical interpretation of observed X-ray spectra of AGN and GBHCs.

We investigate different accretion disk models and viscosity prescriptions in order to provide a basic explanation for the exotic temporal behavior in GRS 1915+105. Based on the fact that the overall cycle times are very much longer than the rise/fall time scales in GRS 1915, we rule out the geometry of ADAF or a hot quasi-spherical region plus a cold outer disk for this source. We thus concentrate on geometrically thin Shakura-Sunyaev type disks (Shakura & Sunyaev 1973; hereafter SS73). We have devised a modified viscosity law that has a quasi-stable upper branch. Via numerical simulations, we show that the model does account for several gross observational features of GRS 1915+105. On the other hand, the rise/fall time scales are not short enough, and no rapid oscillations on time scales $\simlt$ 10 s emerge naturally from the model. We then consider and numerically test a more elaborate model that includes the cold disk, a corona, and plasma ejections from the inner disk region and show that this model allows us to reproduce several additional observed features of GRS 1915+105. We conclude that the most likely structure of the accretion flow in this source is that of a cold disk with a modified viscosity prescription, plus a corona that accounts for much of the X-ray emission, and unsteady plasma ejections that occur when the luminosity of the source is high.

Spectrum of Seyfert 1 Galaxies (S1G hereafter) is very similar to that of several Galactic Black Hole Candidates (GBHCs) in their hard state, suggestive that both classes of objects have similar physical processes. However, recent work has shown that reprocessing features make the patchy corona-disk model (PCD model: the best current explanation of S1G spectrum) problematic for GBHCs. To address the similarities and differences in spectrum of Seyferts and GBHCs, we consider the structure of the ionized X-ray skin near an active magnetic flare. We show that the X-ray skin is subject to a thermal instability, similar in nature to the well known ionization instability of quasar emission line regions. Due to the much higher ionizing X-ray flux in GBHCs, the only stable solution for the upper layer of the accretion disk is that in which it is highly ionized and is at the Compton temperature ($\sim $ few keV). We show that this accounts for the difference in spectrum of GBHCs and S1G. In addition, same instability, applied to S1G, leads to the X-ray skin temperature $T\sim 1-3\times 10^5$ Kelvin, which then may explain the observed spectral shape of BBB.

Recent progress in observations and understanding of spectra of Seyfert Galaxies shows that X-rays are most likely produced by magnetic flares on the surface of the disk, similar to Solar X-ray emission. However, while the model reproduces the shape of the spectrum well, the question of the overall normalization of the X-ray component relative to the bolometric luminosity has not been previously considered. Here we show that, in gas-dominated accretion disks, the magnetic energy transport can indeed power the corona in a way consistent with observations, {\it if} magnetic field in the disk is mostly contained to strong magnetic flux tubes. However, in radiation dominated disks radiation diffusion makes the field weak/diffuse, and thus the magnetic energy transport is less efficient, i.e., in such disks the disk intrinsic emission should be the dominant component in the overall spectrum. We compare our findings to observations, and conclude that our theory can account for the often observed ``steeper when brighter'' behavior in AGN and the hard-soft spectral transitions in GBHCs.

The spectrum of Seyfert 1 Galaxies is very similar to that of several Galactic Black Hole Candidates (GBHCs) in their hard state, suggestive that both classes of objects have similar physical processes. While it appears that an accretion disk corona (ADC) model, where the corona sandwiches the cold accretion disk, is capable of explaining the observations of Seyfert galaxies, recent work has shown that this model is problematic for GBHCs. To address the differences in the spectra of Seyferts and GBHCs, we consider the structure of the atmosphere of the accretion disk in a region near an active magnetic flare (we refer to this region as transition layer). We show that the transition layer is subject to a thermal instability. Due to the much higher ionizing X-ray flux in GBHCs, the only stable solution for the upper layer of the accretion disk is that in which it is highly ionized and is at the Compton temperature ($kT \sim $ a few keV). Using numerical simulations for a slab geometry ADC, we show that the presence of a transition layer can significantly alter the spectrum of escaping radiation for modest optical depths. The presence of the transition layer also leads to a reduction in reprocessing features, i.e., the iron line and the X-ray reflection hump. We find that due to the thermal instability of the surface of the accretion disk, which leads to different endpoints for GBHCs and Seyfert galaxies, the X-ray spectra from these two types of objects can be explained by a single unifying ADC model.

Rapid progress in multi-wavelength observations of Seyfert Galaxies in recent years is providing evidence that X-ray emission in these objects may be produced by magnetic flares occurring above a cold accretion disk. Here we attempt to develop a physically consistent model of accretion disks producing radiation via magnetic flares as well as the optically thick intrinsic disk emission. We place a special emphasis on comparison of the theory and the observations of Active Galactic Nuclei (AGN) and Galactic Black Hole Candidates (GBHCs). In particular, theoretical issues considered in this work provide tentative explanations for: (1) the origin of the Big Blue Bump of Seyfert 1 Galaxies. (2) the recently discovered X-ray Baldwin effect and the controversy over the existence of BBBs in quasars more luminous than typical Seyferts. (3) The differences and similarities in X-ray spectra of Galactic Black Hole Candidates and Seyfert 1s. (4) The hard/soft state transitions of GBHCs and AGN. Our conclusion is that the agreement between the theory and observations is very encouraging and we thus suggest that the physics of magnetic flares is the physics that should be added to the standard accretion disk theory in order to produce a more realistic description of radiation from accretion flows with a large angular momentum.

Some recent work has shown that the spectrum of Seyfert 1 Galaxies is very similar to that of several Galactic Black Hole Candidates (GBHCs) in their hard state. However, the smallness of the observed reflection component in Cygnus X-1, in addition to other constraints, seem to rule out the two-phase model (otherwise successful in the case of Seyfert Galaxies) for GBHCs. Here, we show that the latter conclusion is based on a number of key assumptions that probably are not valid when the overlying corona is patchy, e.g., when it is comprised of localized magnetic flares above the disk. We show that in GBHCs the energy deposited by the X-rays cannot be re-radiated fast enough to maintain equilibrium, unless the X-ray skin heats up to the Compton temperature, at which point the gas is mostly ionized. This leads to a substantially reduced cooling rate for the active regions due to the correspondingly smaller number of re-injected low-energy photons. We model this effect by introducing a transition layer situated between the corona and the cold disk, and find that the resulting spectrum is harder than that obtained with the standard (and unrealistic) two-phase model. We apply this model to Cygnus X-1 and show that it can account for its observed spectrum. This analysis therefore seems to provide a consistentpicture for both Seyfert Galaxies and GBHCs within the same framework, with differences arising due to the changing physical conditionsin the two categories of sources, rather than due to an ad hoc variation of the model parameters.

We here examine the conditions of the two-phase disk model under which magnetic flares arise above the cold accretion disk due to magnetic buoyancy and produce X-rays via Comptonization of the disk's soft radiation. We find that the disk's ability to produce strong magnetic flares is substantially diminished in its radiation dominated regions due to the diffusion of radiation into the magnetic flux tubes. Using a simplified, yet physically self-consistent, model that takes this effect into account, we show that the hard X-ray spectrum of some GBHCs can be explained as the X-ray emission by magnetic flares only when the disk's bolometric luminosity is a relatively small fraction ($\sim$ 5%) of the Eddington value, $L_{Edd}$. Further, we compute the hard ($20-200$ keV) and soft ($1-20$ keV) X-ray power as a function of the disk's luminosity, and find an excellent agreement with the available data for GBHC transient and persistent sources. We conclude that the observed high-energy spectrum of stellar-sized accretion disk systems can be explained by Comptonization of the disk's soft radiation by the hot gas trapped inside the magnetic flares when the luminosity falls in the range $\sim 10^{-3}-10^{-1}\times L_{Edd}$. For higher luminosities, another emission mechanism must be at work. For lower luminosities, the X-ray emissivity may still be dominated by magnetic flares, but this process is more likely to be thermal or non-thermal bremstrahlung, so that the X-ray spectrum below $\sim 10^{-3}L_{Edd}$ may be quite distinct from the typical hard spectrum for higher luminosities.

We discuss several physical constraints on the nature of the Active Regions (AR) in Seyfert 1 Galaxies, and show that a plausible model consistent with these constraints is one in which the ARs are magnetically confined and ``fed''. The unique X-ray index of these sources points to a large compactness parameter ($l\gg 1$). This, together with the conditions required to account for the observed optical depth being close to unity, suggests that the magnetic energy density in the AR should be comparable to the equipartition value in the accretion disk, and that it should be released in a flare-like event above the surface of the cold accretion disk. We consider the various issues pertaining to magnetic flares and attempt to construct a coherent picture, including a reason for the optical depth in the AR being $\sim 1$, and an understanding of the characteristics of the X-ray reflection component and the power density spectra associated with this high-energy emission.

We here consider the pressure equilibrium during an intense magnetic flare above the surface of a cold accretion disk. Under the assumption that the heating source for the plasma trapped within the flaring region is an influx of energy transported inwards with a group velocity close to $c$, e.g., by magnetohydrodynamic waves, this pressure equilibrium can constrain the Thomson optical depth $\tau_T$ to be of order unity. We suggest that this may be the reason why $\tau_T\sim 1$ in Seyfert Galaxies. We also consider whether current data can distinguish between the spectrum produced by a single X-ray emitting region with $\tau_T\sim 1$ and that formed by many different flares spanning a range of $\tau_T$. We find that the current observations do not yet have the required energy resolution to permit such a differentiation. Thus, it is possible that the entire X-ray/$\gamma$-ray spectrum of Seyfert Galaxies is produced by many independent magnetic flares with an optical depth $0.5<\tau_T<2$.

An important feature of the EUV spectrum (known as the Big Blue Bump, hereafter BBB) in Seyfert Galaxies is the narrow range in its cutoff energy E0 from source to source, even though the luminosity changes by 4 orders of magnitude. Here we show that if the BBB is due to accretion disk emission, then in order to account for this ``universality'' in the value of E0, the emission mechanism is probably optically thin bremsstrahlung. In addition, we demonstrate that the two-phase model with active regions localized on the surface of the cold disk is consistent with this constraint if the active regions are very compact and are highly transient, i.e., they evolve faster than one dynamical time scale.

High-energy, multi-component plasmas in which pair creation and annihilation, lepton-lepton scattering, lepton-proton scattering, and Comptonization all contribute to establishing the particle and photon distributions, are present in a broad range of compact astrophysical objects. Earlier work has included much of the microphysics needed to account for electron-photon and electron-proton interactions, but little has been done to handle the redistribution of the particles as a result of their Coulomb interaction with themselves in an arbitrary case. Our goal here is to use a Fokker-Planck approach in order to develop a fully self-consistent theory for the interaction of arbitrarily distributed particles and radiation to arrive at an accurate representation of the high-energy plasma in these sources. We conduct several tests representative of two dominant segments of parameter space and discuss physical implications of the non-Maxwellian distribution function. Approximate analytical forms for the electron distribution in the case of strongly non-Maxwellian plasmas are given.

Many models of cosmological gamma-ray bursts involve the sudden release of $\sim 10^{51}$ erg which produce shells which expand at relativistic speeds (Lorentz $\Gamma$ factors of $10^{2-3}$). We investigate the kinematic limits on the source size due to the observed time structure in three types of bursts: short spikes, FREDs (Fast Rise, Exponentail decay), and long complex bursts. The emitting shell keeps up with the photons it produces reducing apparent durations by $\Gamma^2$ so that source sizes can be very large ($2c\Gamma^2 T_{dur}). However, the thickness of the emitting region is not effected by the shell motion so it must always be small, $c\Delta T$ where $\Delta T$ is the subpeak time scale. We argue that one can only view the bulk motion head-on so it is inappropriate to treat GRBs as viewing the sides of a jet. Although photons come from a region within an angle $\Gamma^{-1}$, we show that the curvature of the shell within that angle creates delays comparable to those associated with the duration of the event. As a result, most bursts should be like FREDs with sharp rises related to how long the shell emits and power law decays related to how long the shell expanded before becoming gamma-ray active. Few bursts have the long decay phases required for large shells resulting in unacceptable high densities for ISM objects to cause the observed subpeaks. To be consistent with the observations, perhaps very thick shells (which act as parallel slabs) are required to avoid the effects of the curvature, or the duration is dictated by a central engine.