by Anders, P. and Baumgardt, H. and Gaburov, E. and Zwart, S. Portegies
14 pages incl. 3 pages with figures and 4 pages of tables (analysis results), MNRAS in press

Most recent progress in understanding the dynamical evolution of star clusters relies on direct N-body simulations. Owing to the computational demands, and the desire to model more complex and more massive star clusters, hardware calculational accelerators, such as GRAPE special-purpose hardware or, more recently, GPUs (i.e. graphics cards), are generally utilised. In addition, simulations can be accelerated by adjusting parameters determining the calculation accuracy (i.e. changing the internal simulation time step used for each star).

We extend our previous thorough comparison (Anders et al. 2009) of basic quantities as derived from simulations performed either with STARLAB/KIRA or NBODY6. Here we focus on differences arising from using different hardware accelerations (including the increasingly popular graphic card accelerations/GPUs) and different calculation accuracy settings.

We use the large number of star cluster models (for a fixed stellar mass function, without stellar/binary evolution, primordial binaries, external tidal fields etc) already used in the previous paper, evolve them with STARLAB/KIRA (and NBODY6, where required), analyse them in a consistent way and compare the averaged results quantitatively. For this quantitative comparison, we apply the bootstrap algorithm for functional dependencies developed in our previous study.

In general we find very high comparability of the simulation results, independent of the used computer hardware (including the hardware accelerators) and the used N-body code. For the tested accuracy settings we find that for reduced accuracy (i.e. time step at least a factor 2.5 larger than the standard setting) most simulation results deviate significantly from the results using standard settings. The remaining deviations are comprehensible and explicable.

by Nissanke, Samaya and Vallisneri, Michele and Nelemans, Gijs and Prince, Thomas A.
17 pages, 3 figures, 5 tables, to be submitted to the Astrophysical Journal

Compact Galactic binaries where at least one member is a white dwarf (WD) or neutron star constitute the majority of individually detectable sources for future low-frequency space-based gravitational-wave (GW) observatories; in addition, they form an unresolved continuum, the dominant Galactic foreground at frequencies below a few mHz. A handful of ultra-compact binaries, observed at optical, ultraviolet and X-ray wavelengths, are known verification sources for space-based GW interferometers. Due to the paucity of electromagnetic observations, the majority of studies of Galactic-binary populations so far have been based on population-synthesis simulations. However, recent surveys have reported several new detections of compact binaries including double WDs, providing new constraints for population estimates. In this article, we evaluate the impact of revised local densities of interacting WD binaries on future low-frequency GW observations. Specifically: we consider five scenarios that explain these densities with different assumptions on the formation of interacting systems; we simulate corresponding populations of detached and interacting WD binaries; we estimate the number of individually detectable GW sources and the magnitude of the confusion-noise foreground, in the case of two GW interferometers with armlengths of 1 and 5 Mkm. We confirm earlier estimates of thousands of detached-binary detections, but project only a few ten to a few hundred detections of interacting systems. We also confirm estimates for the confusion-noise foreground (except in one scenario that explains smaller local densities of interacting systems with fewer progenitor detached systems). Last, we provide a general scaling argument that shows that the magnitude of the GW foreground can be derived robustly from the merger rate of Galactic WD binaries, and depends only weakly on the structure of the Galaxy.

by Li, Shuo and Liu, F. K. and Berczik, Peter and Chen, Xian and Spurzem, Rainer
38 pages, 10 figues; accepted for publication in ApJ

Supermassive black hole binaries (SMBHBs) are the products of frequent galaxy mergers. The coalescence of the SMBHBs is a distinct source of gravitational wave (GW) radiation. The detections of the strong GW radiation and their possible electromagnetic counterparts are essential. Numerical relativity suggests that the post-merger supermassive black hole (SMBH) gets a kick velocity up to 4000 km/s due to the anisotropic GW radiations. Here we investigate the dynamical co-evolution and interaction of the recoiling SMBHs and their galactic stellar environments with one million direct N-body simulations including the stellar tidal disruption by the recoiling SMBHs. Our results show that the accretion of disrupted stars does not significantly affect the SMBH dynamical evolution. We investigate the stellar tidal disruption rates as a function of the dynamical evolution of oscillating SMBHs in the galactic nuclei. Our simulations show that most of stellar tidal disruptions are contributed by the unbound stars and occur when the oscillating SMBHs pass through the galactic center. The averaged disruption rate is ~10^{-6} M_\odot yr^{-1}, which is about an order of magnitude lower than that by a stationary SMBH at similar galactic nuclei. Our results also show that a bound star cluster is around the oscillating SMBH of about ~ 0.7% the black hole mass. In addition, we discover a massive cloud of unbound stars following the oscillating SMBH. We also investigate the dependence of the results on the SMBH masses and density slopes of the galactic nuclei.

by Lousto, Carlos O. and Zlochower, Yosef and Dotti, Massimo and Volonteri, Marta
17 pages, 10 tables, 14 figures, revtex 4

We explore the newly discovered “hangup-kick” effect, which greatly amplifies the recoil for configuration with partial spin- orbital-angular momentum alignment, by studying a set of 48 new simulations of equal-mass, spinning black-hole binaries. We propose a phenomenological model for the recoil that takes this new effect into account and then use this model, in conjunction with statistical distributions for the spin magnitude and orientations, based on accretion simulations, to find the probabilities for observing recoils of several thousand km/s. In addition, we provide initial parameters, eccentricities, radiated linear and angular momentum, precession rates and remnant mass, spin, and recoils for all 48 configurations. Our results indicate that surveys exploring peculiar (redshifted or blueshifted) differential line-of-sight velocities should observe at least one case above 2000 km/s out of four thousand merged galaxies. The probability that a remnant BH receives a total recoil exceeding the ~2000 km/s escape velocity of large elliptical galaxies is ten times larger. Probabilities of recoils exceeding the escape velocity quickly rise to 5% for galaxies with escape velocities of 1000 km/s and nearly 20% for galaxies with escape velocities of 500 km/s. In addition the direction of these large recoils is strongly peaked toward the angular momentum axis, with very low probabilities of recoils exceeding 350 km/s for angles larger than 45 deg. with respect to the orbital angular momentum axis.

by Zubovas, Kastytis and King, Andrew R.
6 pages, 1 figure, 2 tables; accepted for publication in ApJ Letters

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

by Johannsen, Tim and Psaltis, Dimitrios and Gillessen, Stefan and Marrone, Daniel P. and Ozel, Feryal and Doeleman, Sheperd S. and Fish, Vincent L.
7 pages, 4 figures, 2 tables, submitted to ApJ

Dynamical mass measurements to date have allowed determinations of the mass M and the distance D of the galactic center black hole Sgr A* as well as those of other nearby supermassive black holes. In the case of Sgr A*, these measurements are limited by a degeneracy between the mass and distance scaling roughly as M ~ D^2. Future very-long baseline interferometric observations will image a bright and narrow ring surrounding the shadow of the supermassive black hole, if its accretion flow is optically thin. In this paper, we show that the combination of dynamical measurements and VLBI imaging of the ring of Sgr A* breaks the degeneracy between mass and distance. We estimate the signal to noise ratio of near-future VLBI arrays consisting of five to six stations and simulate measurements of the mass and distance of Sgr A* using the expected size of the ring image and existing data of stellar ephemerides. We demonstrate that VLBI observations at 1 mm can already improve the error on the mass by a factor of three compared to the results from the monitoring of stellar orbits alone; observations at 0.5 mm can reduce the error by as much as a factor of 7.5. In addition, we calculate the angular sizes of the bright rings of a number of other nearby supermassive black holes and identify the optimal targets besides Sgr A* that could be imaged by a ground-based VLBI array or a future space-VLBI mission allowing for refined mass measurements.

by Zocchi, A. and Bertin, G. and Varri, A. L.
18 pages, 7 figures. Accepted for publication in Astronomy & Astrophysics

We perform a systematic combined photometric and kinematic analysis of a sample of globular clusters under different relaxation conditions, based on their core relaxation time (as listed in available catalogs), by means of two well-known families of spherical stellar dynamical models. Systems characterized by shorter relaxation time scales are expected to be better described by isotropic King models, while less relaxed systems might be interpreted by means of non-truncated, radially-biased anisotropic f^(\nu) models, originally designed to represent stellar systems produced by a violent relaxation formation process and applied here for the first time to the study of globular clusters. The comparison between dynamical models and observations is performed by fitting simultaneously surface brightness and velocity dispersion profiles. For each globular cluster, the best-fit model in each family is identified, along with a full error analysis on the relevant parameters. Detailed structural properties and mass-to-light ratios are also explicitly derived. We find that King models usually offer a good representation of the observed photometric profiles, but often lead to less satisfactory fits to the kinematic profiles, independently of the relaxation condition of the systems. For some less relaxed clusters, f^(\nu) models provide a good description of both observed profiles. Some derived structural characteristics, such as the total mass or the half-mass radius, turn out to be significantly model-dependent. The analysis confirms that, to answer some important dynamical questions that bear on the formation and evolution of globular clusters, it would be highly desirable to acquire larger numbers of accurate kinematic data-points, well distributed over the cluster field.

by Burkert, Andreas and Schartmann, Mark and Alig, Christian and Gillessen, Stefan and Genzel, Reinhard and Fritz, Tobias and Eisenhauer, Frank
22 pages, 13 figures, submitted to ApJ

The origin, structure and evolution of the small gas cloud, G2, is investigated, that is on an orbit almost straight into the Galactic central supermassive black hole (SMBH). G2 is a sensitive probe of the hot accretion zone of Sgr A*, requiring gas temperatures and densities that agree well with models of captured shock-heated stellar winds. Its mass is equal to the critical mass below which cold clumps would be destroyed quickly by evaporation. Its mass is also constrained by the fact that at apocenter its sound crossing timescale was equal to its orbital timescale. Our numerical simulations show that the observed structure and evolution of G2 can be well reproduced if it formed in pressure equilibrium with the surrounding in 1995 at a distance from the SMBH of 7.6e16 cm. If the cloud would have formed at apocenter in the ‘clockwise’ stellar disk as expected from its orbit, it would be torn into a very elongated spaghetti-like filament by 2011 which is not observed. This problem can be solved if G2 is the head of a larger, shell-like structure that formed at apocenter. Our numerical simulations show that this scenario explains not only G2’s observed kinematical and geometrical properties but also the Br_gamma observations of a low surface brightness gas tail that trails the cloud. In 2013, while passing the SMBH G2 will break up into a string of droplets that within the next 30 years mix with the surrounding hot gas and trigger cycles of AGN activity.

by Vasiliev, E. and Athanassoula, E.
13 pages, 10 figures

We use both N-body simulations and integration in fixed potentials to explore the stability and the long-term secular evolution of self-consistent, equilibrium, non-rotating, triaxial spheroidal galactic models. More specifically, we consider Dehnen models built with the Schwarzschild method. We show that short-term stability depends on the degree of velocity anisotropy (radially anisotropic models are subject to rapid development of radial-orbit instability). Long-term stability, on the other hand, depends mainly on the properties of the potential, and in particular, on whether it admits a substantial fraction of strongly chaotic orbits. We show that in the case of a weak density cusp (gamma=1 Dehnen model) the N-body model is remarkably stable, while the strong-cusp (gamma=2) model exhibits substantial evolution of shape away from triaxiality, which we attribute to the effect of chaotic diffusion of orbits. The different behaviour of these two cases originates from the different phase space structure of the potential; in the weak-cusp case there exist numerous resonant orbit families that impede chaotic diffusion. We also find that it is hardly possible to affect the rate of this evolution by altering the fraction of chaotic orbits in the Schwarzschild model, which is explained by the fact that the chaotic properties of an orbit are not preserved by the N-body evolution. There are, however, parameters in Schwarzschild modelling that do affect the stability of an N-body model, so we discuss the recipes how to build a `good’ Schwarzschild model.

by Crocker, Roland M.
30 pages, 35 figures

We construct a simple model of the star-formation- (and resultant supernova-) driven mass and energy flows through the inner ~200 pc (in diameter) of the Galaxy. Our modelling is constrained, in particular, by the non-thermal radio continuum and {\gamma}-ray signals detected from the region. The modelling points to a current star-formation rate of 0.04 – 0.12 M\msun/year at 2{\sigma} confidence within the region with best-fit value in the range 0.08 – 0.12 M\msun/year which – if sustained over 10 Gyr – would fill out the ~ 10^9 M\msun stellar population of the nuclear bulge. Mass is being accreted on to the Galactic centre (GC) region at a rate ~0.3M\msun/year. The region’s star-formation activity drives an outflow of plasma, cosmic rays, and entrained, cooler gas. Neither the plasma nor the entrained gas reaches the gravitational escape speed, however, and all this material fountains back on to the inner Galaxy. The system we model can naturally account for the recently-observed ~> 10^6 ‘halo’ of molecular gas surrounding the Central Molecular Zone out to 100-200 pc heights. The injection of cooler, high-metallicity material into the Galactic halo above the GC may catalyse the subsequent cooling and condensation of hot plasma out of this region and explain the presence of relatively pristine, nuclear-unprocessed gas in the GC. The plasma outflow from the GC reaches a height of a few kpc and is compellingly related to the recently-discovered Fermi Bubbles. Our modelling demonstrates that ~ 10^9 M\msun of hot gas is processed through the GC over 10 Gyr. We speculate that the continual star-formation in the GC over the age of the Milky Way has kept the SMBH in a quiescent state thus preventing it from significantly heating the coronal gas, allowing for the continual accretion of gas on to the disk and the sustenance of star formation on much wider scales in the Galaxy [abridged].

by Giersz, Mirek and Heggie, Douglas C. and Hurley, Jarrod and Hypki, Arkadiusz
15 pages, 24 figures

We describe a major upgrade of a Monte Carlo code which has previously been used for many studies of dense star clusters. We outline the steps needed in order to calibrate the results of the new Monte Carlo code against N-body simulations for large $latex N$ systems, up to N=200000. The new version of the Monte Carlo code (called MOCCA), in addition to the old version, incorporates direct FewBody integrator for three- and four-body interactions, and new treatment of the escape process based on Fokushige and Heggie (2000). Now stars which fulfil the escape criterion are not removed immediately, but can stay in the system for a certain time which depends on the excess of the energy of a star above the critical energy. They are called potential escapers. FewBody integrator allows to follow all interaction channels, which are important for the rate of creation of various types of objects observed in star clusters, and assures that the energy generation by binaries is treated in a meaner similar to the N-body model.

There are at most three parameters which have to be adjusted against N-body simulations for large N. Two (or one, depends on the chosen approach) connected with the escape process and one responsible for determination of the interaction probabilities. The adopted free parameters are independent on N. They allow MOCCA code to reproduce N-body results, in a reasonably precision, not only for the rate of cluster evolution and the cluster mass distribution, but also for the detailed distributions of mass and binding energy of binaries.

The MOCCA code is at present the most advanced code for simulations of real star clusters. It can follow the cluster evolution in details comparable to N-body code, but orders of magnitude faster.

by Nakasato, Naohito
Journal of Computational Science, 2011; See our recent update at http://galaxy.u-aizu.ac.jp/trac/note/wiki/Octree_On_GPU

The kd-tree is a fundamental tool in computer science. Among other applications, the application of kd-tree search (by the tree method) to the fast evaluation of particle interactions and neighbor search is highly important, since the computational complexity of these problems is reduced from O(N^2) for a brute force method to O(N log N) for the tree method, where N is the number of particles. In this paper, we present a parallel implementation of the tree method running on a graphics processing unit (GPU). We present a detailed description of how we have implemented the tree method on a Cypress GPU. An optimization that we found important is localized particle ordering to effectively utilize cache memory. We present a number of test results and performance measurements. Our results show that the execution of the tree traversal in a force calculation on a GPU is practical and efficient.

by Gillessen, S. and Genzel, R. and Fritz, T. K. and Quataert, E. and Alig, C. and Burkert, A. and Cuadra, J. and Eisenhauer, F. and Pfuhl, O. and Dodds-Eden, K. and Gammie, C. F. and Ott, T.
in press at Nature

Measurements of stellar orbits provide compelling evidence that the compact radio source Sagittarius A* at the Galactic Centre is a black hole four million times the mass of the Sun. With the exception of modest X-ray and infrared flares, Sgr A* is surprisingly faint, suggesting that the accretion rate and radiation efficiency near the event horizon are currently very low. Here we report the presence of a dense gas cloud approximately three times the mass of Earth that is falling into the accretion zone of Sgr A*. Our observations tightly constrain the cloud’s orbit to be highly eccentric, with an innermost radius of approach of only ~3,100 times the event horizon that will be reached in 2013. Over the past three years the cloud has begun to disrupt, probably mainly through tidal shearing arising from the black hole’s gravitational force. The cloud’s dynamic evolution and radiation in the next few years will probe the properties of the accretion flow and the feeding processes of the super-massive black hole. The kilo-electronvolt X-ray emission of Sgr A* may brighten significantly when the cloud reaches pericentre. There may also be a giant radiation flare several years from now if the cloud breaks up and its fragments feed gas into the central accretion zone.

by Alexander, David M. and Hickox, Ryan C.
42 pages, 8 figures; Extensive review to appear in New Astronomy Reviews

Massive black holes (BHs) are at once exotic and yet ubiquitous, residing in the centers of massive galaxies in the local Universe. Recent years have seen remarkable advances in our understanding of how these BHs form and grow over cosmic time, during which they are revealed as active galactic nuclei (AGN). However, despite decades of research, we still lack a coherent picture of the physical drivers of BH growth, the connection between the growth of BHs and their host galaxies, the role of large-scale environment on the fueling of BHs, and the impact of BH-driven outflows on the growth of galaxies. In this paper we review our progress in addressing these key issues, motivated by the science presented at the “What Drives the Growth of Black Holes?” workshop held at Durham on 26th-29th July 2010, and discuss how these questions may be tackled with current and future facilities.

by Liu, K. and Wex, N. and Kramer, M. and Cordes, J. M. and Lazio, T. J. W.
12 pages, 10 Figures, accepted for publication in ApJ

The discovery of radio pulsars in compact orbits around Sgr A* would allow an unprecedented and detailed investigation of the spacetime of the supermassive black hole. This paper shows that pulsar timing, including that of a single pulsar, has the potential to provide novel tests of general relativity, in particular its cosmic censorship conjecture and no-hair theorem for rotating black holes. These experiments can be performed by timing observations with 100 micro-second precision, achievable with the Square Kilometre Array for a normal pulsar at frequency above 15 GHz. Based on the standard pulsar timing technique, we develop a method that allows the determination of the mass, spin, and quadrupole moment of Sgr A*, and provides a consistent covariance analysis of the measurement errors. Furthermore, we test this method in detailed mock data simulations. It seems likely that only for orbital periods below ~0.3 yr is there the possibility of having negligible external perturbations. For such orbits we expect a ~10^-3 test of the frame dragging and a ~10^-2 test of the no-hair theorem within 5 years, if Sgr A* is spinning rapidly. Our method is also capable of identifying perturbations caused by distributed mass around Sgr A*, thus providing high confidence in these gravity tests. Our analysis is not affected by uncertainties in our knowledge of the distance to the Galactic center, R0. A combination of pulsar timing with the astrometric results of stellar orbits would greatly improve the measurement precision of R0.

by Martín, S. and Martín-Pintado, J. and Montero-Castaño, M. and Ho, P. T. P. and Blundell, R.
12 pages, 22 figures. Accepted for publication in Astronomy and Astrophysics

The interstellar region within the few central parsecs around the super-massive black hole, Sgr A* at the very Galactic center is composed by a number of overlapping molecular structures which are subject to one of the most hostile physical environments in the Galaxy. We present high resolution (4″x3″~0.16×0.11 pc) interferometric observations of CN, 13CN, H2CO, SiO, c-C3H2 and HC3N emission at 1.3 mm towards the central ~4 pc of the Galactic center region. Strong differences are observed in the distribution of the different molecules. The UV resistant species CN, the only species tracing all previously identified circumnuclear disk (CND) structures, is mostly concentrated in optically thick clumps in the rotating filaments around Sgr A*. H2CO emission traces a shell-like structure that we interpret as the expansion of Sgr A East against the 50 km/s and 20 km/s giant molecular clouds (GMCs). We derive isotopic ratios 12C/13C~15-45 across most of the CND region. The densest molecular material, traced by SiO and HC3N, is located in the southern CND. The observed c-C3H2/HC3N ratio observed in the region is more than an order of magnitude lower than in Galactic PDRs. Toward the central region only CN was detected in absorption. Apart from the known narrow line-of-sight absorptions, a 90 km/s wide optically thick spectral feature is observed. We find evidences of an even wider (>100 km/s) absorption feature. Around 70-75% of the gas mass, concentrated in just the 27% densest molecular clumps, is associated with rotating structures and show evidences of association with each of the arcs of ionized gas in the mini-spiral structure. Chemical differentiation has been proven to be a powerful tool to disentangle the many overlapping molecular components in this crowded and heavily obscured region.

by Barausse, Enrico and Buonanno, Alessandra and Tiec, Alexandre Le
11 pages, 2 figures

Using the main result of a companion paper, in which the binding energy of a circular-orbit non-spinning compact binary system is computed at leading-order beyond the test-particle approximation, the exact expression of the effective-one-body (EOB) metric component $latex g^\text{eff}_{tt}$ is obtained through first order in the mass ratio. Combining these results with the recent gravitational self-force calculation of the periastron advance for circular orbits in the Schwarzschild geometry, the EOB metric component $latex g^\text{eff}_{rr}$ is also determined at linear order in the mass ratio. These results assume that the mapping between the real and effective Hamiltonians at the second and third post-Newtonian (PN) orders holds at all PN orders. Our findings also confirm the advantage of resumming the PN dynamics around the test-particle limit if the goal is to obtain a flexible model that can smoothly connect the test-mass and equal-mass limits.

by Jalali, B. and Baumgardt, H. and Kissler-Patig, M. and Gebhardt, K. and Noyola, E. and Lützgendorf, N. and de Zeeuw, P. T.
Accepted for publication in A&A

Supermassive black holes (SMBHs) are fundamental keys to understand the formation and evolution of their host galaxies. However, the formation and growth of SMBHs are not yet well understood. One of the proposed formation scenarios is the growth of SMBHs from seed intermediate-mass black holes (IMBHs, 10^2 to 10^5 M_{\odot}) formed in star clusters. In this context, and also with respect to the low mass end of the M-sigma relation for galaxies, globular clusters are in a mass range that make them ideal systems to look for IMBHs. Among Galactic star clusters, the massive cluster $latex \omega$ Centauri is a special target due to its central high velocity dispersion and also its multiple stellar populations. We study the central structure and dynamics of the star cluster $latex \omega$ Centauri to examine whether an IMBH is necessary to explain the observed velocity dispersion and surface brightness profiles. We perform direct N-body simulations to follow the dynamical evolution of $latex \omega$ Centauri. The simulations are compared to the most recent data-sets in order to explain the present-day conditions of the cluster and to constrain the initial conditions leading to the observed profiles. We find that starting from isotropic spherical multi-mass King models and within our canonical assumptions, a model with a central IMBH mass of 2% of the cluster stellar mass, i.e. a 5×10^4 M_{\odot} IMBH, provides a satisfactory fit to both the observed shallow cusp in surface brightness and the continuous rise towards the center of the radial velocity dispersion profile. In our isotropic spherical models, the predicted proper motion dispersion for the best-fit model is the same as the radial velocity dispersion one. (abridged)

by Barausse, Enrico and Buonanno, Alessandra and Hughes, Scott A. and Khanna, Gaurav and O’Sullivan, Stephen and Pan, Yi
19 pages, 14 figures, 6 tables

Using the effective-one-body (EOB) formalism and a time-domain Teukolsky code, we generate inspiral, merger, and ringdown waveforms in the small mass-ratio limit. We use EOB inspiral and plunge trajectories to build the Teukolsky equation source term, and compute full coalescence waveforms for a range of black hole spins. By comparing EOB waveforms that were recently developed for comparable mass binary black holes to these Teukolsky waveforms, we improve the EOB model for the (2,2), (2,1), (3,3), and (4,4) modes. Our results can be used to quickly and accurately extract useful information about merger waves for binaries with spin, and should be useful for improving analytic models of such binaries. Although in this analysis we only consider equatorial inspirals (orbital angular momentum parallel to spin), there is no issue of principle preventing us from considering inclined binaries. We will extend this analysis to examine misaligned spin-orbit configurations in future work.

by Banerjee, Sambaran and Kroupa, Pavel
14 pages, 4 figures. Published in The Astrophysical Journal Letters

Among the most explored directions in the study of dense stellar systems is the investigation of the effects of the retention of supernova remnants, especially that of the massive stellar remnant black holes (BHs), in star clusters. By virtue of their eventual high central concentration, these stellar mass BHs potentially invoke a wide variety of physical phenomena, the most important ones being emission of gravitational waves (GWs), formation of X-ray binaries, and modification of the dynamical evolution of the cluster. Here we propose, for the first time, that rapid removal of stars from the outer parts of a cluster by the strong tidal field in the inner region of our Galaxy can unveil its BH sub-cluster, which appears as a star cluster that is gravitationally bound by an invisible mass. We study the formation and properties of such systems through direct N-body computations and estimate that they can be present in significant numbers in the inner region of the Milky Way. We call such objects “dark star clusters” (DSCs) as they appear dimmer than normal star clusters of similar mass and they comprise a predicted, new class of entities. The finding of DSCs will robustly cross-check BH retention; they will not only constrain the uncertain natal kicks of BHs, thereby the widely debated theoretical models of BH formation, but will also pinpoint star clusters as potential sites for GW emission for forthcoming ground-based detectors such as the Advanced LIGO. Finally, we also discuss the relevance of DSCs for the nature of IRS 13E.

by Mancini, Luigi and Feoli, Antonio
11 pages, 8 figures, to appear in Astronomy & Astrophysics

Thanks to the angular resolution of modern telescopes and kinematic models, the existence of supermassive black holes (SMBHs) in the inner part of galaxies has been established on quite solid grounds. A possible correlation between the mass of SMBHs and the evolutionary state of their host galaxies is expected. Based on the recent 2D decomposition of mid-infrared Spiter/IRAC images of local galaxies with M_bh measurements, we investigated various scaling laws, studying what the best predictor of the mass of the central SMBHs is. We focused on the M_bh-M_G sigma^2 law, the relation between the mass of SMBHs and the kinetic energy of random motions of the corresponding host galaxies. In order to find the best fit for each of the scaling laws examined, we performed a least-squares regression of M_bh on x for the considered sample of galaxies, x being a whatever known parameter of the galaxy bulge. Our analysis shows that M_bh-M_G sigma^2 law fits the examined experimental data successfully as much as the other known scaling laws and shows a value of chi^2 better than the others, a result which is consistent with previous determinations. This means that a combination of sigma and M_G could be necessary to drive the correlations between M_bh and other bulge properties. This issue has been investigated by a careful analysis of the residuals of the various relations. In order to avoid rushed conclusions on galaxy activity and evolution, the indirect inferring of M_bh from the kinetic energy of random motions should be considered, especially when applied to higher redshift galaxies. This statement is suggested by a reanalysis of the SDSS data used to study the SMBH growth in the nearby Universe. Adopting the M_bh-M_G sigma^2 relation instead of the M_bh-sigma, a radio-quiet/radio-loud dichotomy appears in the SMBH mass distribution of the corresponding SDSS early-type AGN galaxies.

by Touma, Jihad R. and Sridhar, S.
41 pages, 5 figures, Submitted to MNRAS

We formulate the collisionless Boltzmann equation (CBE) for dense star clusters that lie within the radius of influence of a massive black hole in galactic nuclei. Our approach to these nearly Keplerian systems follows that of Sridhar and Touma (1999): Delaunay canonical variables are used to describe stellar orbits and we average over the fast Keplerian orbital phases. The stellar distribution function (DF) evolves on the longer time scale of precessional motions, whose dynamics is governed by a Hamiltonian, given by the orbit-averaged self-gravitational potential of the cluster. We specialize to razor-thin, planar discs and consider two counter-rotating (”$latex \pm$”) populations of stars. To describe discs of small eccentricities, we expand the $latex \pm$ Hamiltonian to fourth order in the eccentricities, with coefficients that depend self-consistently on the $latex \pm$ DFs. We construct approximate $latex \pm$ dynamical invariants and use Jeans’ theorem to construct time-dependent $latex \pm$ DFs, which are completely described by their centroid coordinates and shape matrices. When the centroid eccentricities are larger than the dispersion in eccentricities, the $latex \pm$ centroids obey a set of 4 autonomous ordinary differential equations. We show that these can be cast as a two-degree of freedom Hamiltonian system which is nonlinear, yet integrable. We study the linear instability of initially circular discs and derive a criterion for the counter-rotating instability. We then explore the rich nonlinear dynamics of counter-rotating discs, with focus on the variety of steadily precessing eccentric configurations that are allowed. The stability and properties of these configurations are studied as functions of parameters such as the disc mass ratios and angular momentum.

by Zubovas, Kastytis and Nayakshin, Sergei and Markoff, Sera
11 pages. MNRAS submitted

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.

by Antonini, Fabio and Capuzzo-Dolcetta, Roberto and Mastrobuono-Battisti, Alessandra and Merritt, David
15 pages, 14 figure. Submitted to ApJ

In one widely discussed model for the formation of nuclear star clusters (NSCs), massive globular clusters spiral into the center of a galaxy and merge to form the nucleus. It is now known that at least some NSCs coexist with supermassive black holes (SBHs); this is the case, for instance, in the Milky Way (MW). In this paper, we investigate how the presence of a SMBH at the center of the MW impacts the merger hypothesis for the formation of its NSC. Starting from a model consisting of a low-density nuclear stellar disk and the SMBH, we use N-body simulations to follow the successive inspiral and merger of (12) globular clusters. The clusters are started on circular orbits of radius 20 pc, and their initial masses and radii are set up in such a way as to be consistent with the galactic tidal field at that radius. The total accumulated mass is about 1.5×10^7 Solar masses. Each cluster is disrupted by the SMBH at a distance of roughly one parsec. The density profile that results after the final inspiral event is characterized by a core of roughly this radius, and an envelope with density that falls off rho \sim r^-2. These properties are similar to those of the MW NSC, with the exception of the core size, which in the MW is a little smaller. But by continuing the evolution of the model after the final inspiral event, we find that the core shrinks substantially via gravitational encounters in a time (when scaled to the MW) of 10 Gyr as the stellar distribution evolves toward a Bahcall-Wolf cusp. We also show that the luminosity function of the MW NSC is consistent with the hypothesis that a large fraction of the mass comes from (~10Gyr) old stars, brought in by globular clusters. We conclude that a model in which a large fraction of the mass of the MW NSC arose from infalling globular clusters is consistent with existing observational constraints.

by Li, Gongjie and Conroy, Charlie and Loeb, Abraham
8 pages, 3 figures, 2 tables

We investigate the evolution of the MBH-{\sigma} relation by examining the relationship between the intrinsic scatter in the MBH-{\sigma} relation and galaxy bolometric nuclear luminosity, the latter being a probe of the accretion rate of the black hole (BH). Our sample is composed of galaxies with classical bulges when possible, of which 38 have dynamically measured BHs masses, and 17 have BHs masses measured by reverberation mapping. In order to obtain the bolometric nuclear luminosity for galaxies with low nuclear luminosity, we convert the X-ray nuclear luminosity measured by Chandra to bolometric luminosity. We find that the scatter in the MBH-{\sigma} relation is uncorrelated with nuclear luminosity over seven orders of magnitude in luminosity, with the high luminosity end approaching the Eddington luminosity. This suggests that at the present epoch galaxies evolve along the MBH-{\sigma} relation. This conclusion is consistent with the standard paradigm that BHs grow contemporaneously with their host stellar spheroids.

by Melia, Fulvio and Falanga, Maurizio and Goldwurm, Andrea
Accepted for Publication in MNRAS, September 26, 2011

The Galaxy’s supermassive black hole, Sgr A*, produces an outburst of infrared radiation about once every 6 hours, sometimes accompanied by an even more energetic flurry of X-rays. The NIR photons are produced by nonthermal synchrotron processes, but we still don’t completely understand where or why these flares originate, nor exactly how the X-rays are emitted. The power-law electrons radiating the infrared light may be partially cooled, so the distribution may be a broken power law with a (”cooling break”) transition frequency. In addition, the emission region appears to be rather compact, possibly restricted to the inner edge of the accretion disk. In that case, the X-ray outburst may itself be due to synchrotron processes by the most energetic particles in this population. In this paper, we examine several key features of this proposal, producing relativistically correct polarimetric images of Sgr A*’s NIR and X-ray flare emission, in order to determine (1) whether the measured NIR polarization fraction is consistent with this geometry, and (2) whether the predicted X-ray to NIR peak fluxes are confirmed by the currently available multi-wavelength observations. We also calculate the X-ray polarization fraction and position angle (relative to that of the NIR photons) in anticipation of such measurements in the coming years. We show that whereas the polarization fraction and position angle of the X-rays are similar to those of the NIR component for synchrotron-cooled emission, these quantities are measurably different when the X-rays emerge from a scattering medium. It is clear, therefore, that the development of X-ray polarimetry will represent a major new tool for studying the spacetime near supermassive black holes.

by Stone, Nicholas and Loeb, Abraham
4 pages, 4 figures

When a star is tidally disrupted by a supermassive black hole (SMBH), the streams of liberated gas form an accretion disk after their return to pericenter. We demonstrate that Lense-Thirring precession in the spacetime around a rotating SMBH can produce significant time evolution of the disk angular momentum vector, due to both the periodic precession of the disk and the nonperiodic, differential precession of the bound debris streams. Jet precession and periodic modulation of disk luminosity are possible consequences. The persistence of the jetted X-ray emission in the Swift J164449.3+573451 flare suggests that the jet axis was aligned with the spin axis of the SMBH during this event.

by Vincent, F. H. and Paumard, T. and Perrin, G. and Gourgoulhon, E. and Eisenhauer, F. and Gillessen, S.
5 pages, 3 figures, to appear in the Proceedings of the French Society of Astronomy and Astrophysics (SF2A)

The ability of the near future second generation VLTI instrument GRAVITY to constrain the properties of the Galactic center black hole is investigated. The Galactic center infrared flares are used as probes of strong-field gravity, within the framework of the hot spot model according to which the flares are the signature of a blob of gas orbiting close to the black hole’s innermost stable circular orbit. Full general relativistic computations are performed, together with realistic observed data simulations, that lead to conclude that GRAVITY could be able to constrain the black hole’s inclination parameter.

by Alexander, R. D. and Smedley, S. L. and Nayakshin, S. and King, A. R.
7 pages, 3 figures. Accepted for publication in MNRAS

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

by Yusef-Zadeh, F. and Bushouse, H. and Wardle, M.
16 pages, 7 figures, ApJ, in press

We present HST/NICMOS data to study the surface brightness distribution of stellar light within the inner 10″ of Sgr A* at 1.4, 1.7 and 1.9 microns. We use these data to independently examine the surface brightness distribution that had been measured previously with NICMOS and to determine whether there is a drop in the surface density of stars very near Sgr A*. Our analysis confirms that a previously reported drop in the surface brightness within 0.8″ of Sgr A* is an artifact of bright and massive stars near that radius. We also show that the surface brightness profile within 5″ or ~0.2 pc of Sgr A* can be fitted with broken power laws. The power laws are consistent with previous measurements, in that the profile becomes shallower at small radii. For radii > 0.7″ the slope is beta=-0.34\pm0.04 where Sigma is proportional to r^beta and becomes flatter at smaller radii with beta=-0.13\pm0.04. Modeling of the surface brightness profile gives a stellar density that increases roughly as r^-1 within the inner 1″ of Sgr A*. This slope confirms earlier measurements in that it is not consistent with that expected from an old, dynamically-relaxed stellar cluster with a central supermassive black hole. Assuming that the diffuse emission is not contaminated by a faint population of young stars down to the 17.1 magnitude limit of our imaging data at 1.70$latex \mu$, the shallow cusp profile is not consistent with a decline in stellar density in the inner arcsecond. In addition, converting our measured diffuse light profile to a stellar mass profile, with the assumption that the light is dominated by K0 dwarfs, the enclosed stellar mass within radius r < 0.1 pc of Sgr A* is ~ 3.2×10^4 M_solar (r/0.1 {pc})^2.1.

by Gil-Merino, Rodrigo and Goicoechea, Luis J. and Shalyapin, Vyacheslav N. and Braga, Vittorio F.
21 pages, 7 text pages, 13 figures, 2 tables, accepted by The Astrophysical Journal

We present the results of our X-ray, UV and optical monitoring campaign of the first gravitationally lensed AGN from late 2009 to mid 2010. The trailing (B) image of the AGN 0957+561 shows the intrinsic continuum variations that were predicted in advance based on observations of the leading (A) image in the gr optical bands. This multiwavelength variability of the B image allows us to carry out a reverberation mapping analysis in the radio-loud AGN 0957+561 at redshift z = 1.41. We find that the U-band and r-band light curves are highly correlated with the g-band record, leading and trailing it by 3 +/- 1 days (U band) and 4 +/- 1 days (r band). These 1-sigma measurements are consistent with a scenario in which flares originated in the immediate vicinity of the supermassive black hole are thermally reprocessed in a standard accretion disk at about 10-20 Schwarzschild radii from the central dark object. We also report that the light curve for the X-ray emission with power-law spectrum is delayed with respect to those in the Ugr bands by about 32 days. Hence, the central driving source can not be a standard corona emitting the observed power-law X-rays. This result is also supported by X-ray reprocessing simulations and the absence of X-ray reflection features in the spectrum of 0957+561. We plausibly interpret the lack of reflection and the 32-day delay as evidence for a power-law X-ray source in the base of the jet at a typical height of about 200 Schwarzschild radii. A central EUV source would drive the variability of 0957+561.

by Gualandris, Alessia and Dotti, Massimo and Sesana, Alberto
5 pages, 5 figures, submitted to MNRAS

We study the evolution of the orientation of the orbital plane of massive black hole binaries (BHBs) in rotating stellar systems in which the total angular momentum of the stellar cusp is misaligned with respect to that of the binary. We compare results from direct summation N-body simulations with predictions from a simple theoretical model. We find that the same encounters between cusp stars and the BHB that are responsible for the hardening and eccentricity evolution of the binary, lead to a reorientation of the binary orbital plane. In particular, binaries whose angular momentum is initially misaligned with respect to that of the stellar cusp tend to realign their orbital planes with the angular momentum of the cusp on a timescale of a few hardening times. This is due to angular momentum exchange between stars and the BHB during close encounters, and may have important implications for the relative orientation of host galaxies and radio jets.

by Amaro-Seoane, Pau and Brem, Patrick and Cuadra, Jorge and Armitage, Philip J.
Submitted

Measurements of gravitational waves from the inspiral of a stellar-mass compact object into a massive black hole (MBH) are unique probes to test General Relativity (GR) and MBH properties, as well as the stellar distribution about these holes in galactic nuclei. Current data analysis techniques can provide us with parameter estimation with very narrow errors. However, an EMRI is not a two-body problem, since other stellar bodies orbiting nearby will influence the capture orbit. Any deviation from the isolated inspiral of the binary will induce a small, though observable deviation from the idealised waveform which could be misinterpreted as a failure of GR. Based on conservative analysis of mass segregation in a Milky Way like nucleus, we estimate that the possibility that a star has a semi-major axis comparable to that of the EMRI is non-negligible. This star introduces an observable perturbation in the orbit in the case in which we consider only loss of energy via gravitational radiation at periapsis. When considering the two first-order non-dissipative post-Newtonian contributions (the periapsis shift of the orbit) the evolution of the orbital elements of the EMRI turns out to be chaotic in nature. The implications of this study are twofold. From the one side, the application to testing GR and measuring MBHs parameters with the detection of EMRIs in galactic nuclei with a millihertz mission will be even more challenging than believed. From the other side, this behaviour could in principle be used as a signature of mass segregation in galactic nuclei.

by Konstantinidis, Symeon and Amaro-Seoane, Pau and Kokkotas, Kostas D.
Submitted

Contrary to supermassive and stellar-mass black holes (SBHs), the existence of intermediate-mass black holes (IMBHs) with masses ranging between 100 and 10,000 Msun has not yet been confirmed. The main problem in the detection is that the innermost stellar kinematics of globular clusters (GCs), the natural loci to IMBHs, are very difficult to resolve. However, if IMBHs reside in the center of GCs, a possibility is that they interact dynamically with their enviroment. A binary formed with the IMBH and a compact object of the GC would naturally lead to a prominent source of gravitational radiation, detectable with future observatories. We run for the first time direct-summation integrations of GCs with an IMBH including the dynamical evolution of the IMBH with the stellar system and relativistic effects, such as energy loss in gravitational waves (GWs) and periapsis shift, and gravitational recoil. We find in one of our models an intermediate-mass ratio inspiral (IMRI), which leads to a merger with a recoiling velocity higher than the escape velocity of the GC. The GWs emitted fall in the range of frequencies that a LISA-like observatory could detect, like the European eLISA or in mission options considered in the recent preliminary mission study conducted in China. The merger has an impact on the global dynamics of the cluster, as an important heating source is removed when the merged system leaves the GC. The detection of one IMRI would constitute a test of GR, as well as an irrefutable proof of the existence of IMBHs.

by Nakano, Hiroyuki and Zlochower, Yosef and Lousto, Carlos O. and Campanelli, Manuela
23 pages, 35 figures, revtex4

We revisit the scenario of small-mass-ratio (q) black-hole binaries; performing new, more accurate, simulations of mass ratios 10:1 and 100:1 for initially nonspinning black holes. We propose fitting functions for the trajectories of the two black holes as a function of time and mass ratio (in the range 1/100 < q < 1/10$) that combine aspects of post-Newtonian trajectories at smaller orbital frequencies and plunging geodesics at larger frequencies. We then use these trajectories to compute waveforms via black hole perturbation theory. Using the advanced LIGO noise curve, we see a match of ~99.5% for the leading (l,m)=(2,2) mode between the numerical relativity and perturbative waveforms. Nonleading modes have similarly high matches. We thus prove the feasibility of efficiently generating a bank of gravitational waveforms in the intermediate-mass-ratio regime using only a sparse set of full numerical simulations.

by Noyola, Eva and Baumgardt, Holger
Accepted for publication in ApJ

Surface photometry is a necessary tool to establish the dynamical state of stars clusters. We produce realistic HST-like images from N-body models of star clusters with and without central intermediate-mass black holes (IMBHs) in order to measure their surface brightness profiles. The models contain ~600,000 individual stars, black holes of various masses between 0% to 2% of the total mass, and are evolved for a Hubble time. We measure surface brightness and star count profiles for every constructed image in order to test the effect of intermediate mass black holes on the central logarithmic slope, the core radius, and the half-light radius. We use these quantities to test diagnostic tools for the presence of central black holes using photometry. We find that the the only models that show central shallow cusps with logarithmic slopes between -0.1 and -0.4 are those containing central black holes. Thus, the central logarithmic slope seems to be a good way to choose clusters suspect of containing intermediate-mass black holes. Clusters with steep central cusps can definitely be ruled out to host an IMBH. The measured r_c/r_h ratio has similar values for clusters that have not undergone core-collapse, and those containing a central black hole. We notice that observed Galactic globular clusters have a larger span of values for central slope and r_c/r_h than our modeled clusters, and suggest possible reasons that could account for this and contribute to improve future models.

by Fiestas, Jose and Porth, Oliver and Berczik, Peter and Spurzem, Rainer
15 pages, 7 figures,accepted by MNRAS

NBody realizations of axisymmetric collisional galaxy cores (e.g. M32, M33, NGC205, Milky Way) with embedded growing black holes are presented. Stars which approach the disruption sphere are disrupted and accreted to the black hole. We measure the zone of influence of the black hole and disruption rates in relaxation time scales. We show that secular gravitational instabilities dominate the initial core dynamics, while the black hole is small and growing due to consumption of stars. Later, the black hole potential dominates the core, and loss cone theory can be applied. Our simulations show that central rotation in galaxies can not be neglected for relaxed systems, and compare and discuss our results with the standard theory of spherically symmetric systems.

by Rice, W. K. M. and Armitage, P. J. and Mamatsashvili, G. R. and Lodato, G. and Clarke, C. J.
MNRAS, in press

Self-gravity becomes competitive as an angular momentum transport process in accretion discs at large radii, where the temperature is low enough that external irradiation likely contributes to the thermal balance. Irradiation is known to weaken the strength of disc self-gravity, and can suppress it entirely if the disc is maintained above the threshold for linear instability. However, its impact on the susceptibility of the disc to fragmentation is less clear. We use two-dimensional numerical simulations to investigate the evolution of self-gravitating discs as a function of the local cooling time and strength of irradiation. In the regime where the disc does not fragment, we show that local thermal equilibrium continues to determine the stress – which can be represented as an effective viscous alpha – out to very long cooling times (at least 240 dynamical times). In this regime, the power spectrum of the perturbations is uniquely set by the effective viscous alpha and not by the cooling rate. Fragmentation occurs for cooling times tau < beta_crit / Omega, where beta_crit is a weak function of the level of irradiation. We find that beta_crit declines by approximately a factor of two, as irradiation is increased from zero up to the level where instability is almost quenched. The numerical results imply that irradiation cannot generally avert fragmentation of self-gravitating discs at large radii; if other angular momentum transport sources are weak mass will build up until self-gravity sets in, and fragmentation will ensue.

by Antonini, Fabio and Merritt, David
28 pages, 23 figures

The density of stars in galactic bulges is often observed to be flat or slowly rising inside the influence radius of the supermassive black hole (SMBH). Chandrasekhar’s dynamical friction formula predicts little or no frictional force on a test body in such a core, regardless of its density, due to the absence of stars moving more slowly than the local circular velocity. We have tested this prediction using large-scale $latex N$-body experiments. The rate of orbital decay never drops precisely to zero, because stars moving faster than the test body also contribute to the frictional force. When the contribution from the fast-moving stars is included in the expression for the dynamical friction force, and the changes induced by the massive body on the stellar distribution are taken into account, Chandrasekhar’s theory is found to reproduce the rate of orbital decay remarkably well. However, this rate is still substantially smaller than the rate predicted by Chandrasekhar’s formula in its most widely-used forms, implying longer time scales for inspiral. Motivated by recent observations that suggest a parsec-scale core around the Galactic center SMBH, we investigate the evolution of a population of stellar-mass black holes (BHs) as they spiral in to the center of the Galaxy. After $latex \sim 10$ Gyr, we find that the density of BHs can remain substantially less than the density in stars at all radii; we conclude that it would be unjustified to assume that the spatial distribution of BHs at the Galactic center is well described by steady-state models. One consequence is that rates of capture of BHs by the SMBH at the Galactic center (EMRIs) may be much lower than in standard models. When capture occurs, inspiraling BHs often reach the gravitational-radiation-dominated regime while on orbits that are still highly eccentric.

by Brockamp, M. and Baumgardt, H. and Kroupa, P.
19 pages, 11 figures, accepted for publication in MNRAS

The disruption rate of stars by supermassive black holes (SMBHs) is calculated numerically with a modified version of Aarseth’s NBODY6 code. The initial stellar distribution around the SMBH follows a S\’{e}rsic n=4 profile representing bulges and early type galaxies. In order to infer relaxation driven effects and to increase the statistical significance, a very large set of N-body integrations with different particle numbers N, ranging from 10^{3} to 0.5 \cdot 10^{6} particles, is performed. Three different black hole capture radii are taken into account, enabling us to scale these results to a broad range of astrophysical systems with relaxation times shorter than one Hubble time, i.e. for SMBHs up to M_bh \approx 10^{7} M_sun. The computed number of disrupted stars are driven by diffusion in angular momentum space into the loss cone of the black hole and the rate scales with the total number of particles as dN/dt \propto N^{b}, where b is as large as 0.83. This is significantly steeper than the expected scaling dN/dt \propto ln(N) derived from simplest energy relaxation arguments. Only a relatively modest dependence of the tidal disruption rate on the mass of the SMBH is found and we discuss our results in the context of the M_bh/sigma relation. The number of disrupted stars contribute a significant part to the mass growth of black holes in the lower mass range as long as a significant part of the stellar mass becomes swallowed by the SMBH. This also bears direct consequences for the search and existence of IMBHs in globular clusters. For SMBHs similar to the galactic center black hole SgrA*, a tidal disruption rate of 55 \pm 27 events per Myr is deduced. Finally relaxation driven stellar feeding can not account for the masses of massive black holes M_bh \geq 10^{7} M_sun. (abridged)

by Wardle, Mark and Yusef-Zadeh, Farhad
Submitted to ApJ Letters

Several thin, Keplerian, sub-parsec megamaser disks have been discovered in the nuclei of active galaxies and used to precisely determine the mass of their host black holes. We show that there is an empirical linear correlation between the disk radius and the black hole mass. We demonstrate that such disks are naturally formed by the partial capture of molecular clouds passing through the galactic nucleus and temporarily engulfing the central supermassive black hole. Imperfect cancellation of the angular momenta of the cloud material colliding after passing on opposite sides of the hole leads to the formation of a compact disk. The radial extent of the disk is determined by the efficiency of this process and the Bondi-Hoyle capture radius of the black hole, and naturally produces the empirical linear correlation of the radial extent of the maser distribution with black hole mass. The disk has sufficient column density to allow X-ray irradiation from the central source to generate physical and chemical conditions conducive to the formation of 22 GHz H2O masers. For initial cloud column densities less than ~10^{23.5} cm^-2 the disk is non-self gravitating, consistent with the ordered kinematics of the edge-on megamaser disks; for higher cloud columns the disk would fragment and produce a compact stellar disk similar to that observed around Sgr A* at the galactic centre.

by Lousto, Carlos O. and Zlochower, Yosef
4 pages, 3 figures, revtex 4

We revisit the scenario of the gravitational radiation recoil acquired by the final remnant of a black-hole-binary merger by studying a set of configurations that have components of the spin both aligned with the orbital angular momentum and in the orbital plane. We perform a series of 24 new full numerical simulations for equal-mass and equal-spin-magnitude binaries, but with different spin orientations.

We extend previous recoil fitting formulas to include nonlinear terms in the spins and successfully include both the new and known results. For this new formula the predicted maximum velocity approaches 5000km/s. More importantly, from the astrophysical point of view, it reaches this maximum for spins partially aligned with the orbital angular momentum. The optimal configuration is near an equipartition of the hangup and superkick contributions. This newly discovered contribution to the recoil leads to an important increase of the probabilities of large recoils in generic astrophysical mergers. We measure these probabilities for the case of accretion-aligned spins and find non-negligible probabilities for supermassive black hole encounters leading to recoil velocities of several thousand km/s.

by Gultekin, Kayhan and Richstone, Douglas O. and Gebhardt, Karl and Faber, S. M. and Lauer, Tod R. and Bender, Ralf and Kormendy, John and Pinkney, Jason
10 pages, 8 figures, accepted by ApJ

We present HST STIS observations of the galaxy NGC 4382 (M85) and axisymmetric models of the galaxy to determine mass-to-light ration (M/L, V-band) and central black hole mass (M_BH). We find M/L = 3.74 +/- 0.1 (solar units) and M_BH = 1.3 (+5.2, -1.2) \times 10^7 M_sun at an assumed distance of 17.9 Mpc, consistent with no black hole. The upper limit, M_BH < 9.6 \times 10^7 M_sun (2{\sigma}) or M_BH < 1.4 \times 10^8 M_sun (3{\sigma}) is consistent with the current M-{\sigma} relation, which predicts M_BH = 8.8 \times 10^7 M_sun at {\sigma}_e = 182 km/s, but low for the current M-L relation, which predicts M_BH = 7.8 \times 10^8 M_sun at L_V = 8.9 \times 10^10 L_sun,V. HST images show the nucleus to be double, suggesting the presence of a nuclear eccentric stellar disk, in analogy to the Tremaine disk in M31. This conclusion is supported by the HST velocity dispersion profile. Despite the presence of this non-axisymmetric feature and evidence of a recent merger, we conclude that the reliability of our black hole mass determination is not hindered. The inferred low black hole mass may explain the lack of nuclear activity.

by Nieuwenhuizen, Theo M.
5 pages latex

It is assumed that a galaxy starts as a dark halo of a few million Jeans clusters (JCs), each of which consists of nearly a trillion micro brown dwarfs, MACHOs of earth mass. JCs in the galaxy center heat up their MACHOs by tidal forces, which makes them expand, so that coagulation and star formation occurs. Being continuously fed by matter from bypassing JCs, the star(s) may transform into a super massive black hole. It has a fast $latex t^3$ growth during the first mega years, and a slow $latex t^{1/3}$ growth at giga years. JCs disrupted by a close encounter can provide matter for the bulge. Those that survive can be so agitated that they form stars and become globular star clusters. Thus black holes mostly arise together with galactic bulges in their own environment and are about as old as the oldest globular clusters. The age 13.2 Gyr of the star HE 1523-0901 (Frebel et al. 2007) puts forward that the Galactic halo was fully assembled at that moment. In case of merging super massive black holes the JCs passing near the galactic center provide ideal assistance to overcome the last parsec.

by Sazonov, S. and Sunyaev, R. and Revnivtsev, M.
16 pages, 7 figures. Submitted to MNRAS

Chandra has detected optically thin, thermal X-ray emission with a size of ~1 arcsec and luminosity ~10^33 erg/s from the direction of the Galactic supermassive black hole (SMBH), Sgr A*. We suggest that a significant or even dominant fraction of this signal may be produced by several thousand late-type main-sequence stars that possibly hide in the central ~0.1 pc region of the Galaxy. As a result of tidal spin-ups caused by close encounters with other stars and stellar remnants, these stars should be rapidly rotating and hence have hot coronae, emitting copious amounts of X-ray emission with temperatures kT<~ a few keV. The Chandra data thus place an interesting upper limit on the space density of (currently unobservable) low-mass main-sequence stars near Sgr A*. This bound is close to and consistent with current constraints on the central stellar cusp provided by infrared observations. If coronally active stars do provide a significant fraction of the X-ray luminosity of Sgr A*, it should be variable on hourly and daily time scales due to giant flares occurring on different stars. Another consequence is that the quiescent X-ray luminosity and accretion rate of the SMBH are yet lower than believed before.

by Lützgendorf, N. and Kissler-Patig, M. and Noyola, E. and Jalali, B. and de Zeeuw, P. T. and Gebhardt, K. and Baumgardt, H.
12 pages, 12 figures, Accepted for publication in A&A

Intermediate-mass black holes (IMBHs) are of interest in a wide range of astrophysical fields. In particular, the possibility of finding them at the centers of globular clusters has recently drawn attention. IMBHs became detectable since the quality of observational data sets, particularly those obtained with HST and with high resolution ground based spectrographs, advanced to the point where it is possible to measure velocity dispersions at a spatial resolution comparable to the size of the gravitational sphere of influence for plausible IMBH masses. We present results from ground based VLT/FLAMES spectroscopy in combination with HST data for the globular cluster NGC 6388. The aim of this work is to probe whether this massive cluster hosts an intermediate-mass black hole at its center and to compare the results with the expected value predicted by the $latex M_{\bullet} – \sigma$ scaling relation. The spectroscopic data, containing integral field unit measurements, provide kinematic signatures in the center of the cluster while the photometric data give information of the stellar density. Together, these data sets are compared to dynamical models and present evidence of an additional compact dark mass at the center: a black hole. Using analytical Jeans models in combination with various Monte Carlo simulations to estimate the errors, we derive (with 68% confidence limits) a best fit black-hole mass of $latex (17 \pm 9) \times 10^3 M_{\odot}$ and a global mass-to-light ratio of $latex M/L_V = (1.6 \pm 0.3) \ M_{\odot}/L_{\odot}$.

by Gualandris, Alessia and Merritt, David
22 pages, 23 figures, submitted to ApJ

We simulate mergers between galaxies containing collisionally-relaxed nuclei around massive black holes (BHs). Our galaxies contain four mass groups, representative of old stellar populations; a primary goal is to understand the distribution of stellar-mass BHs after the merger. Mergers are followed using direct-summation N-body simulations, assuming a mass ratio of 1:3 and two different orbits. Evolution of the massive BH binary is followed until its separation has shrunk by a factor of 20 below the hard-binary separation. During the galaxy merger, large cores are carved out in the stellar distribution, with radii several times the influence radius of the massive BH. Much of the pre-existing mass segregation is erased during this phase. We follow the evolution of the merged galaxies for approximately three, central relaxation times after coalescence of the massive binary; both standard, and top-heavy, mass functions are considered. The cores that were formed in the stellar distribution persist, and the distribution of the stellar-mass black holes evolves against this essentially fixed background. Even after three central relaxation times, these models look very different from the relaxed, multi-mass models that are often assumed to describe the distribution of stars and stellar remnants near a massive BH; in particular, the density of stellar BHs is much smaller than in those models. We discuss the implications of our results for the EMRI problem and for the existence of Bahcall-Wolf cusps.

by Lang, Meagan and Holley-Bockelmann, Kelly and Bogdanovic, Tamara and Amaro-Seoane, Pau and Sesana, Alberto
9 pages, 1 figure, submitted to ApJ. Comments are welcome and may be incorporated into the paper with attribution

Observations of the Galactic Center (GC) have accumulated a multitude of “forensic” evidence indicating that several million years ago the center of the Milky Way galaxy was teaming with starforming and accretion-powered activity — this paints a rather different picture from the GC as we understand it today. We examine a possibility that this epoch of activity could have been triggered by the infall of a satellite galaxy into the Milky Way which began at the redshift of 10 and ended few million years ago with a merger of the Galactic supermassive black hole with an intermediate mass black hole brought in by the inspiralling satellite.

by Iorio, Lorenzo
LaTex2e, 21 pages, 4 tables, no figures

Empirically determining the averaged variations of the orbital parameters of the stars orbiting the Supermassive Black Hole (SBH) hosted by the Galactic Centre (GC) in Sgr A* is, in principle, a valuable tool to put on the test the General Theory of Relativity (GTR), in regimes far stronger than those tested so far, and certain key predictions of it like the no-hair theorems. We analytically work out the long-term variations of all the six osculating Keplerian orbital elements of a test particle orbiting a non-spherical, rotating body with quadrupole moment Q_2 and angular momentum S for a generic spatial orientation of its spin axis k. This choice is motivated by the fact that, basically, we do not know the position in the sky of the spin axis of the SBH in Sgr A* with sufficient accuracy. We apply our results to S2, which is the closest star discovered so far having an orbital period P_b = 15.98 yr, and to a hypothetical closer star X with P_b = 0.5 yr. Our calculations are quite general, not being related to any specific parameterization of k, and can be applied also to astrophysical binary systems, stellar planetary systems, and planetary satellite geodesy in which different reference frames, generally not aligned with the primary’s rotational axis, are routinely used.

by Volonteri, Marta and Stark, Daniel P.
MNRAS in press

Motivated by recent observational results that focus on high redshift black holes, we explore the effect of scatter and observational biases on the ability to recover the intrinsic properties of the black hole population at high redshift. We find that scatter and selection biases can hide the intrinsic correlations between black holes and their hosts, with ‘observable’ subsamples of the whole population suggesting, on average, positive evolution even when the underlying population is characterized by no- or negative evolution. We create theoretical mass functions of black holes convolving the mass function of dark matter halos with standard relationships linking black holes with their hosts. Under these assumptions, we find that the local MBH – sigma correlation is unable to fit the z = 6 black hole mass function proposed by Willott et al. (2010), overestimating the number density of all but the most massive black holes. Positive evolution or including scatter in the MBH – sigma correlation makes the discrepancy worse, as it further increases the number density of observable black holes. We notice that if the MBH – sigma correlation at z = 6 is steeper than today, then the mass function becomes shallower. This helps reproducing the mass function of z = 6 black holes proposed by Willott et al. (2010). Alternatively, it is possible that very few halos (of order 1/1000) host an active massive black hole at z = 6, or that most AGN are obscured, hindering their detection in optical surveys. Current measurements of the high redshift black hole mass function might be underestimating the density of low mass black holes if the active fraction or luminosity are a function of host or black hole mass. Finally, we discuss physical scenarios that can possibly lead to a steeper MBH – sigma relation at high redshift.