There is both theoretical expectation and some observational clues that intermediate mass black holes reside in nuclei of globular clusters. In order to find an independent indicator for their existence, we investigate in this paper how an IMBH manifests itself through its dynamical interaction with a binary rich globular cluster of moderate extension and mass. By means of direct $N$-body integration we follow the dynamical evolution of models of such a system over a time span of $\approx$ 0.8 Gyr and compare the cases with and without the primordial binaries as well as with and without the IMBH. In accord with previous results, we show that when present the IMBH develops a power-law density cusp of stars around it, regardless of the binary population in the cluster. If, however, binaries are present, their interaction with the IMBH leads to the production of high velocity escapers at a rate of the order of 0.1 Myr$^{-1}$. These stars may contribute to the population of high-velocity stars observed in the Galaxy. Clusters hosting the IMBH together with high number of binaries also form a denser halo of marginally unbound stars than clusters that lack either the IMBH or the rich binary population. Finally, we show that the binary population leads to an increased rate of direct interactions of stars with the IMBH, potentially observable as tidal disruption events.

Hypervelocity stars (HVSs) that are observed in the Galactic halo, are believed to be accelerated to large velocities by a process of tidal disruption of binary stars passing close to a supermassive black hole (SMBH) which resides in the center of the Galaxy. It is, however, still unclear, where these relatively young stars were born and which dynamical process pushed them to nearly radial orbits around the SMBH. In this paper we investigate the possibility that the young binaries originated from a thin eccentric disc, similar to the one observed in the Galactic center nowadays. By means of direct N-body simulations, we follow the dynamical evolution of an initially thin and eccentric disc of stars with a 100% binary fraction orbiting around the SMBH. Such a configuration leads to Kozai-Lidov oscillations of orbital elements, bringing considerable amount of binaries to close vicinity of the black hole. Subsequent tidal disruption of these binaries accelerates one of their component to velocities well above the escape velocity from the SMBH while the second component becomes tightly bound to the SMBH. We describe the main kinematic properties of the escaping and tightly bound stars within our model and compare them qualitatively to the properties of the observed HVSs and S-stars, respectively. The most prominent feature is a strong anisotropy in the directions of the escaping stars which is observed for the Galactic HVSs but not explained yet.

There is growing evidence of star formation in the vicinity of supermassive black holes (SMBH) in galactic nuclei. A viable scenario for this process assumes infall of a massive gas cloud towards the SMBH and subsequent formation of a dense accretion disc which gives birth to the young stars. Numerical hydrodynamical models indicate that this star formation process is rather fast and it precedes full circularization of the accretion flow, i.e. the new stars are born on elliptic orbits. By means of direct numerical N-body modeling, we show in this paper that the non-zero eccentricity of the stellar discs around the SMBH leads to an onset of various types of the Kozai-Lidov oscillations of a non-negligible subset of individual orbits in the disc, showing a remarkable robustness of this classical mechanism. Among others, we demonstrate that under certain circumstances, presence of an additional spherical cluster (which is generally known to damp Kozai-Lidov oscillations) may trigger such oscillations due to affecting the internal flow of the angular momentum through the disc. We conclude that the Kozai-Lidov oscillations are capable to substantially modify the initial structure of the disc (its thickness and distribution of eccentricities, in particular).

The centre of our Galaxy hosts almost two hundreds of very young stars, a subset of which is orbiting the central supermassive black hole (SMBH) in a relatively thin disc-like structure. First analyses indicated a power-law surface density profile of the disc, $\Sigma \propto R^\beta$ with $\beta = -2$. Recently, however, speculations about this profile arose. In particular, it now seems to be better described by a sort of broken power-law. By means of both analytical arguments and numerical N-body modelling, we show that such a broken power-law profile is a natural consequence of the two-body relaxation of the disc which is, due to small relative velocities of stars on nearby co-planar Keplerian orbits around the SMBH, effective enough to affect the evolution of the disc on time-scales comparable to its estimated age. In the inner, densest part of the disc, the profile becomes rather flat ($\beta ~ -1$) while the outer parts keep imprints of the initial state. Our numerical models show that the observed projected surface density profile of the young stellar disc can result from two-body relaxation driven evolution of a disc with initial single power-law profile with $-2 < \beta < -1.5$. In addition, we suggest that two-body relaxation may have caused a significant radial migration of the S-stars towards the central SMBH, thus playing an important role in their formation scenario.

Runaway stars are stars observed to have large peculiar velocities. Two mechanisms are thought to contribute to the ejection of runaway stars, both involve binarity (or higher multiplicity). In the binary supernova scenario a runaway star receives its velocity when its binary massive companion explodes as a supernova (SN). In the alternative dynamical ejection scenario, runaway stars are formed through gravitational interactions between stars and binaries in dense, compact clusters or cluster cores. Here we study the ejection scenario. We make use of extensive N-body simulations of massive clusters, as well as analytic arguments, in order to to characterize the expected ejection velocity distribution of runaways stars. We find the ejection velocity distribution of the fastest runaways (>~80 km s^-1) depends on the binary distribution in the cluster, consistent with our analytic toy model, whereas the distribution of lower velocity runaways appears independent of the binaries properties. For a realistic log constant distribution of binary separations, we find the velocity distribution to follow a simple power law; Gamma(v) goes like v^(-8/3) for the high velocity runaways and v^(-3/2) for the low velocity ones. We calculate the total expected ejection rates of runaway stars from our simulated massive clusters and explore their mass function and their binarity. The mass function of runaway stars is biased towards high masses, and depends strongly on their velocity. The binarity of runaways is a decreasing function of their ejection velocity, with no binaries expected to be ejected with v>150 km s^-1. We also find that hyper-runaways with velocities of hundreds of km s^-1 can be dynamically ejected from stellar clusters, but only at very low rates, which cannot account for a significant fraction of the observed population of hyper-velocity stars in the Galactic halo.

We investigate the dynamical evolution of the Orion Nebula Cluster (ONC) by means of direct N-body integrations. A large fraction of residual gas was probably expelled when the ONC formed, so we assume that the ONC was much more compact when it formed compared to its current size, in agreement with the embedded cluster radius-mass relation from Marks & Kroupa (2012). Hence, we assume that few-body relaxation played an important role during the initial phase of evolution of the ONC. In particular, three body interactions among OB stars likely led to their ejection from the cluster and, at the same time, to the formation of a massive object via runaway physical stellar collisions. The resulting depletion of the high mass end of the stellar mass function in the cluster is one of the important points where our models fit the observational data. We speculate that the runaway-mass star may have collapsed directly into a massive black hole (Mbh > 100Msun). Such a dark object could explain the large velocity dispersion of the four Trapezium stars observed in the ONC core. We further show that the putative massive black hole is likely to be a member of a binary system with appr. 70 per cent probability. In such a case, it could be detected either due to short periods of enhanced accretion of stellar winds from the secondary star during pericentre passages, or through a measurement of the motion of the secondary whose velocity would exceed 10 km/s along the whole orbit.

We investigate the orbital evolution of a system of N mutually interacting stars on initially circular orbits around the dominating central mass. We include perturbative influence of a distant axisymmetric source and an extended spherical potential. In particular, we focus on the case when the secular evolution of orbital eccentricities is suppressed by the spherical perturbation. By means of standard perturbation methods, we derive semi-analytic formulae for the evolution of normal vectors of the individual orbits. We find its two qualitatively different modes. Either the orbits interact strongly and, under such circumstances, they become dynamically coupled, precessing synchronously in the potential of the axisymmetric perturbation. Or, if their mutual interaction is weaker, the orbits precess independently, interchanging periodically their angular momentum, which leads to oscillations of inclinations. We argue that these processes may have been fundamental for the evolution of the disc of young stars orbiting the supermassive black hole in the centre of the Milky Way.

The Galactic centre hosts, according to observations, a number of early-type stars. About one half of those which are orbiting the central supermassive black hole on orbits with projected radii $\gtrsim$ 0.03 pc form a coherently rotating disc. Observations further reveal a massive gaseous torus and a significant population of late-type stars. In this paper, we investigate, by means of numerical N-body computations, the orbital evolution of the stellar disc, which we consider to be initially thin. We include the gravitational influence of both the torus and the late-type stars, as well as the self-gravity of the disc. Our results show that, for a significant set of system parameters, the evolution of the disc leads, within the lifetime of the early-type stars, to a configuration compatible with the observations. In particular, the disc naturally reaches a specific - perpendicular - orientation with respect to the torus, which is indeed the configuration observed in the Galactic centre. We, therefore, suggest that all the early-type stars may have been born within a single gaseous disc.

Models of the origin of young stars in the Galactic Centre are facing various problems. The most promissing scenario of the star formation in a thin self-gravitating disc naturally forms stars on coherently rotating orbits, but it fails to explain origin of several tens of stars that evidently do not belong to any of the disc-like structures in the GC. One possible solution lies in rather complicated initial conditions, assuming at least two infalling and interacting gas clouds. We present alternative solution showing that a single thin stellar disc may have given birth to all young stars in the GC. The outliers are explained as stars that have been stripped from the parent structure due to the gravitational interaction with the gaseous circum-nuclear disc.

Within the central parsec of the Galaxy, several tens of young stars orbiting a central supermassive black hole are observed. A subset of these stars forms a coherently rotating disc. Other observations reveal a massive molecular torus which lies at a radius ~1.5pc from the centre. In this paper we consider the gravitational influence of the molecular torus upon the stars of the stellar disc. We derive an analytical formula for the rate of precession of individual stellar orbits and we show that it is highly sensitive upon the orbital semi-major axis and inclination with respect to the plane of the torus as well as on the mass of the torus. Assuming that both the stellar disc and the molecular torus are stable on the time-scale >6Myr, we constrain the mass of the torus and its inclination with respect to the young stellar disc. We further suggest that all young stars observed in the Galactic Centre may have a common origin in a single coherently rotating structure with an opening angle <5deg, which was partially destroyed (warped) during its lifetime by the gravitational influence of the molecular torus.

We examine influence of the circum-nuclear disc (CND) upon the orbital evolution of young stars in the Galactic Centre. We show that gravity of the CND causes precession of the orbits which is highly sensitive upon the semi-major axis and inclination. We consider such a differential precession within the context of an ongoing discussion about the origin of the young stars and suggest a possibility that all of them have originated in a thin disc which was partially destroyed due to the influence of the CND during the period of ~6Myr.

Galactic nuclei are unique laboratories for the study of processes connected with the accretion of gas onto supermassive black holes. At the same time, they represent challenging environments from the point of view of stellar dynamics due to their extreme densities and masses involved. There is a growing evidence about the importance of the mutual interaction of stars with gas in galactic nuclei. Gas rich environment may lead to stellar formation which, on the other hand, may regulate accretion onto the central mass. Gas in the form of massive torus or accretion disc further influences stellar dynamics in the central parsec either via gravitational or hydrodynamical interaction. Eccentricity oscillations on one hand and energy dissipation on the other hand lead to increased rate of infall of stars into the supermassive black hole. Last, but not least, processes related to the stellar dynamics may be detectable with forthcoming gravitational waves detectors.

Mass segregation stands as one of the most robust features of the dynamical evolution of self-gravitating star clusters. In this paper we formulate parametrised models of mass segregated star clusters in virial equilibrium. To this purpose we introduce mean inter-particle potentials for statistically described unsegregated systems and suggest a single-parameter generalisation of its form which gives a mass-segregated state. We describe an algorithm for construction of appropriate star cluster models. Their stability over several crossing-times is verified by following the evolution by means of direct N-body integration.

We study the probability of close encounters between stars from a nuclear cluster and a massive black hole. The gravitational field of the system is dominated by the black hole in its sphere of influence. It is further modified by the cluster mean field (a spherical term) and a gaseous disc/torus (an axially symmetric term) causing a secular evolution of stellar orbits via Kozai oscillations. Intermittent phases of large eccentricity increase the chance that stars become damaged inside the tidal radius of the central hole. Such events can produce debris and lead to recurring episodes of enhanced accretion activity. We introduce an effective loss cone and associate it with tidal disruptions during the high-eccentricity phases of the Kozai cycle. By numerical integration of the trajectories forming the boundary of the loss cone we determine its shape and volume. We also include the effect of relativistic advance of pericentre. The potential of the disc has the efffect of enlarging the loss cone and, therefore, the predicted number of tidally disrupted stars should grow by factor of ~10^2. On the other hand, the effect of the cluster mean potential together with the relativistic pericentre advance act against the eccentricity oscillations. In the end we expect the tidal disruption events to be approximately ten times more frequent in comparison with the model in which the three effects -- the cluster mean field, the relativistic pericentre advance, and the Kozai mechanism -- are all ignored. The competition of different influences suppresses the predicted star disruption rate as the black hole mass increases. Hence, the process under consideration is more important for intermediate-mass black holes, M_bh~10^4M_s.

We discuss the structure of a central stellar cluster whose dynamics is influenced by gravitation of a supermassive black hole and by the dissipative interaction of orbiting stars with an accretion disc. We also take the effect of disc self-gravity into account. We show that the cluster properties are determined predominantly by the radial profile of the disc surface density. To this aim we develop a simple steady-state model of the central cluster and we estimate the rate at which stars migrate to the centre. This model is relevant for central regions of an active galactic nucleus (AGN) containing a rather dense accretion medium. In passing we also briefly mention a possibility that a fossil (inactive) disc could exist in the centre of our own Galaxy. Such a hypothetical disc could perturb the motion of stars and set them on highly elliptic trajectories with small pericentre distances. The required mass of the disc is less than one percent of the central black hole mass, i.e. below the upper limit permitted by present accuracy of the orbital parameters of S-stars in the Galactic Centre.

We propose that Kozai's phenomenon is responsible for the long-term evolution of stellar orbits near a supermassive black hole. We pursue the idea that this process may be driven by a fossil accretion disc in the centre of our Galaxy, causing the gradual orbital decay of stellar trajectories, while setting some stars on highly elliptic orbits. We evolve model orbits that undergo repetitive transitions across the disc over the period of ~10^7 years. We assume that the disc mass is small compared to the central black hole, and its gravitational field comparatively weak, yet non-zero, and we set the present values of orbital parameters of the model star consistent with those reported for the S2 star in Sagittarius A*. We show how a model trajectory decays and circularizes, but at some point the mean eccentricity is substantially increased by Kozai's resonance. In consequence the orbital decay of highly eccentric orbits is accelerated. A combination of an axially symmetric gravitational field and dissipative environment can provide a mechanism explaining the origin of stars on highly eccentric orbits tightly bound to the central black hole. In the context of other S-stars, we can conclude that an acceptable mass of the disc (i.e., M_d<=1 percent of the black hole mass) is compatible with their surprisingly young age and small pericentre distances, provided these stars were formed at r<=10^5 gravitational radii.

Structure of a quasi-stationary stellar cluster is modelled assuming that it is embedded in the gravitational field of a super-massive black hole. Gradual orbital decay of stellar trajectories is caused by the dissipative interaction with an accretion disc. Gravitational field of the disc is constructed and its effect on the cluster structure is taken into account as an axially symmetric perturbation. Attention is focused on a circumnuclear region (r<10^4 gravitational radii) where the effects of the central black hole and the disc dominate over the influence of an outer galaxy. It is shown how the stellar system becomes gradually flattened towards the disc plane. For certain combinations of the model parameters, a toroidal structure is formed by a fraction of stars. Growing anisotropy of stellar velocities as well as their segregation occur. The mass function of the inner cluster is modified and it progressively departs from the asymptotic form assumed in the outer cluster. A new stationary distribution can be characterized in terms of velocity dispersion of the stellar sample in the central region of the modified cluster.

Motion of stellar-mass satellites is studied around a massive compact body which is surrounded by a gaseous slab of a stationary accretion disc. The satellites suffer an orbital decay due to hydrodynamical interaction with the disc medium (transitions across the disc, gap opening in the disc, density waves) and gravitational radiation. Arbitrary orbital eccentricities and inclinations are considered, and it is observed how the competing effects depend on the parameters of the model, namely, the mass and compactness of the orbiters, the osculating elements of their trajectories, and surface density of the disc. These effects have a visible impact on the satellites long-term motion, and they can produce observational consequences with respect to galactic central clusters. It is shown that the satellite-disc collisions do not impose serious restrictions on the results of gravitational wave experiments if the disc medium is diluted and the orbiter is compact but they are important in the case of environments with relatively high density. We thus concentrate on application to accretion flows in which the density is not negligible. We discuss the expected quasi-stationary structure of the cluster that is established on sub-parsec scales within the sphere of gravitational influence of the central object. Relevant to this region, we give the power-law slopes defining the radial profile of modified clusters and we show that their values are determined by satellite interaction with the accretion flow rather than their initial distribution.

We have explored spectral line profiles due to spiral patterns in accretion discs around black holes. A parametrization was employed for the shape and emissivity of spiral waves, which can be produced by non-axisymmetric perturbations affecting the disc density and ionization structure. The effects of the light-travel time, energy shift, and gravitational focusing near to a rotating black hole were taken into account. A high-resolution ray-tracing code was used to follow the time variations of the synthetic line profile. A variety of expected spectral features were examined and the scheme applied to a broad iron line observed in MCG-6-30-15.

We further investigate the long-term evolution of a trajectory of a stellar-mass orbiter which is gravitationally bound to a massive central body. The orbiter undergoes repetitionary collisions with an accretion disc. We consider eccentric orbits intersecting the disc once or twice per each revolution, and we solve the equations for osculating elements. We find terminal radii of the orbits and time needed to bring the orbit into the disc plane as function of initial parameters, and we show that previous simplified estimates (derived for the case of low eccentricity) remain valid within factor of two. The discussion presented in this paper offers a toy model of the orbital evolution of a satellite which passes through the rarefied gaseous environment of a stationary disc. We demonstrate that the drag on the satellite should be taken into account in calculations of the stellar distribution near the super-massive black hole in galactic centers. On the other hand, it does not impose serious restrictions on the results of gravitational wave experiments, except for the case of discs with rather high surface density.