Resonant relaxation near a massive black hole: the stellar distribution and gravitational wave sources

Resonant relaxation (RR) of orbital angular momenta occurs near massive black holes (MBHs) where the stellar orbits are nearly Keplerian and so do not precess significantly. The resulting coherent torques efficiently change the magnitude of the angular momenta and rotate the orbital inclination in all directions. As a result, many of the tightly bound stars very near the MBH are rapidly destroyed by falling into the MBH on low-angular momentum orbits, while the orbits of the remaining stars are efficiently randomized. We solve numerically the Fokker-Planck equation in energy for the steady state distribution of a single mass population with a RR sink term. We find that the steady state current of stars, which sustains the accelerated drainage close to the MBH, can be up to ~10 times larger than that due to non-coherent 2-body relaxation alone. RR mostly affects tightly bound stars, and so it increases only moderately the total tidal disruption rate, which is dominated by stars originating from less bound orbits farther away. We show that the event rate of gravitational wave (GW) emission from inspiraling stars, originating much closer to the MBH, is dominated by RR dynamics. The GW event rate depends on the uncertain efficiency of RR. The efficiency indicated by the few available simulations implies rates ~10 times higher than those predicted by 2-body relaxation, which would improve the prospects of detecting such events by future GW detectors, such as LISA. However, a higher, but still plausible RR efficiency can lead to the drainage of all tightly bound stars and strong suppression of GW events from inspiraling stars. We apply our results to the Galactic MBH, and show that the observed dynamical properties of stars there are consistent with RR.Comment: Accepted to ApJ; Minor revision

Galactic Centre stellar winds and Sgr A* accretion

(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.Comment: 16 pages, 18 colour figures. Accepted by MNRAS. Full resolution paper and movies available at http://www.mpa-garching.mpg.de/~jcuadra/Winds/ . (v2: minor changes

The Milky Way's Fermi Bubbles: Echoes of the Last Quasar Outburst?

{\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.Comment: 6 pages, 0 figures; accepted for publication in MNRA

On the origin of kinematic distribution of the sub-parsec young stars in the Galactic center

Within a half-parsec from the Galactic center (GC), there is a population of coeval young stars which appear to reside in a coherent disk. Surrounding this dynamically-cool stellar system, there is a population of stars with a similar age and much larger eccentricities and inclinations relative to the disk. We propose a hypothesis for the origin of this dynamical dichotomy. Without specifying any specific mechanism, we consider the possibility that both stellar populations were formed within a disk some 6 Myr ago. But this orderly structure was dynamically perturbed outside-in by an intruding object with a mass ~10^4 Msun, which may be an intermediate-mass black hole (IMBH) or a dark stellar cluster hosting an IMBH. We suggest that the perturber migrated inward to ~0.15-0.3pc from the GC under the action of dynamical friction. Along the way, it captured many stars in the outer disk region into its mean-motion resonance, forced them to migrate with it, closely encountered with them, and induced the growth of their eccentricity and inclination. But stars in the inner regions of the disk retain their initial coplanar structure. We predict that some of the inclined and eccentric stars surrounding the disk may have similar Galactocentric semimajor axis. Future precision determination of their kinematic distribution of these stars will not only provide a test for this hypothesis but also evidences for the presence of an IMBH or a dark cluster at the immediate proximity of the massive black hole at the GC. (abridged)Comment: 14 pages, including 13 figures, typo corrected, reference added, ApJ in pres

Evolution of the X-ray spectrum in the flare model of Active Galactic Nuclei

Nayakshin & Kazanas (2002) have considered the time-dependent illumination of an accretion disc in Active Galactic Nuclei, in the lamppost model. We extend their study to the flare model, which postulates the release of a large X-ray flux above a small region of the accretion disc. A fundamental difference with the lamppost model is that the region of the disc below the flare is not illuminated before the onset of the flare. A few test models show that the spectrum which follows immediately the increase in continuum flux should display the characteristics of a highly illuminated but dense gas, i.e. very intense X-ray emission lines and ionization edges in the soft X-ray range. The behaviour of the iron line is different in the case of a "moderate" and a ``strong'' flare: for a moderate flare, the spectrum displays a neutral component of the Fe K$\alpha$ line at 6.4 keV, gradually leading to more highly ionized lines. For a strong flare, the lines are already emitted by FeXXV (around 6.7 keV) after the onset, and have an equivalent width of several hundreds of eV. We find that the observed correlations between $R$, $\Gamma$, and the X-ray flux, are well accounted by a combination of flares having not achieved pressure equilibrium, strongly suggesting that the observed spectrum is dominated by regions in non-pressure equilibrium, typical of the onset of the flares. Finally a flare being confined to a small region of the disc, the spectral lines should be narrow (except for a weak Compton broadening), Doppler shifted, and moving.Comment: 14 pages, 13 figures, accepted in A & A, english corrected versio