[abridged] Energy relaxation around a massive black hole (MBH) is key to establishing the dynamical state of galactic nuclei, and the nature of close stellar interactions with the MBH. The standard description of relaxation as diffusion provides a perturbative 2nd-order solution in the weak two-body interaction limit. We run N-body simulations and find that this solution fails to describe the non-Gaussian relaxation on short timescale, which is strongly influenced by extreme events even in the weak limit, and is thus difficult to characterize and measure. We derive a non-perturbative solution for relaxation as an anomalous diffusion process, and develop a robust estimation technique to measure it in simulations. These enable us to analyze and model our numerical results, and validate in detail, for the first time, this model of energy relaxation around an MBH on all timescales. We derive the relation between the energy diffusion time, t_E, and the time for a small perturbation to return to steady state, t_r, in a relaxed, single mass cusp around a MBH. We constrain the contribution of strong encounters, measure that of the weakest encounters, determine the value of the Coulomb logarithm, and provide a robust analytical estimate for t_E in a finite nuclear stellar cusp. We find that t_r ~ 10t_E ~(5/32)Q^2P_h/N_h log Q, where Q=M_bh/M_* is the MBH to star mass ratio, the orbital period P_h and number of stars N_h are evaluated at the energy scale corresponding to the MBH's sphere of influence, E_h=sigma_inf^2, where sigma_inf is the velocity dispersion far from the MBH. We conclude, using the observed cosmic M_bh/sigma correlation, that cusps around lower-mass MBHs (M_bh<10^7 Mo), which evolved passively over a Hubble time, should be relaxed. We consider the effects of anomalous energy diffusion on orbital perturbations of stars observed near the Galactic MBH.

In the braneworld scenario the four dimensional effective Einstein equation has extra source terms, which arise from the embedding of the 3-brane in the bulk. These non-local effects, generated by the free gravitational field of the bulk, may provide an explanation for the dynamics of the neutral hydrogen clouds at large distances from the galactic center, which is usually explained by postulating the existence of the dark matter. In the present paper we consider the asymptotic behavior of the galactic rotation curves in the brane world models, and we compare the theoretical results with observations of both High Surface Brightness and Low Surface Brightness galaxies. For the chosen sample of galaxies we determine first the baryonic parameters by fitting the photometric data to the adopted galaxy model; then we test the hypothesis of the Weyl fluid acting as dark matter on the chosen sample of spiral galaxies by fitting the tangential velocity equation of the combined baryonic-Weyl model to the rotation curves. We give an analytical expression for the rotational velocity of a test particle on a stable circular orbit in the exterior region to a galaxy, with Weyl fluid contributions included. The model parameter ranges for which the $\chi^2$ test provides agreement (within 1$\sigma$ confidence level) with observations on the velocity fields of the chosen galaxy sample are then determined. There is a good agreement between the theoretical predictions and observations, showing that extra-dimensional models can be effectively used as a viable alternative to the standard dark matter paradigm.

We use N-body simulations to study the evolution of the orbital eccentricities of stars deposited near (<0.05 pc) the Milky Way massive black hole (MBH), starting from initial conditions motivated by two competing models for their origin: formation in a disk followed by inward migration; and exchange interactions involving a binary star. The first model predicts modest eccentricities, lower than those observed in the S-star cluster, while the second model predicts higher eccentricities than observed. The N-body simulations include a dense cluster of 10 M_sun stellar black holes (SBHs), expected to accumulate near the MBH by mass segregation. Perturbations from the SBHs tend to randomize the stellar orbits, partially erasing the dynamical signatures of their origin. The eccentricities of the initially highly eccentric stars evolve, in 20 Myr (the S-star lifespan), to a distribution that is consistent at the ~95 % level with the observed eccentricity distribution. In contrast, the eccentricities of the initially more circular orbits fail to evolve to the observed values in 20 Myr, arguing against the disk migration scenario. We find that 20 % - 30 % of the S-stars are tidally disrupted by the MBH over their lifetimes, and that the S-stars are not likely to be ejected as hypervelocity stars outside the central 0.05 pc by close encounters with stellar black holes.

Resonant relaxation (RR) is a rapid relaxation process that operates in the nearly-Keplerian potential near a massive black hole (MBH). RR dominates the dynamics of compact remnants that inspiral into a MBH and emit gravitational waves (extreme mass ratio inspiral events, EMRIs). RR can either increase the EMRI rate, or strongly suppress it, depending on its still poorly-determined efficiency. We use small-scale Newtonian N-body simulations to measure the RR efficiency and to explore its possible dependence on the stellar number density profile around the MBH, and the mass-ratio between the MBH and a star (a single-mass stellar population is assumed). We develop an efficient and robust procedure for detecting and measuring RR in N-body simulations. We present a suite of simulations with a range of stellar density profiles and mass-ratios, and measure the mean RR efficiency in the near-Keplerian limit. We do not find a statistically significant dependence on the density profile or the mass-ratio. Our numerical determination of the RR efficiency in the Newtonian, single-mass population approximations, suggests that RR will likely enhance the EMRI rate by a factor of a few over the rates predicted assuming only slow stochastic two-body relaxation.

We give a gravitational upper limit for the mass of static degenerate fermionic dark matter objects. The treatment we use includes fully relativistic equations for describing the static solutions of these objects. We study the influence of the annihilation of the particles on this mass limit. We give the change of its value over the age of the Universe with annihilation cross sections relevant for such fermions constituting the dark matter. Our calculations take into account the possibility of Dirac as well Majorana spinors.

Recent observations of the Galactic center revealed a nuclear disk of young OB stars, in addition to many similar outlying stars with higher eccentricities and/or high inclinations relative to the disk (some of them possibly belonging to a second disk). Binaries in such nuclear disks, if they exist in non-negligible fractions, could have a major role in the evolution of the disks through binary heating of this stellar system. We suggest that interactions with/in binaries may explain some (or all) of the observed outlying young stars in the Galactic center. Such stars could have been formed in a disk, and later on kicked out from it through binary related interactions, similar to ejection of high velocity runaway OB stars in young clusters throughout the galaxy.

Due to the importance of accurate Fourier parameters, we devise a method that is more appropriate for deriving these parameters on low-quality data than the traditional Fourier fitting. Based on the accurate light curves of 248 fundamental mode RR Lyrae stars, we test the power of a full-fetched implementation of the template method in the computation of the Fourier decomposition. The applicability of the method is demonstrated also on datasets of filter passbands different from that of the template set. We examine in more detail the question of the estimation of Fourier- based iron abundance [Fe/H] and average brightness. We get, for example, for light curves sampled randomly in 30 data points with sigma=0.03 mag observational noise that optimized direct Fourier fits yield sigma_[Fe/H]=0.33, whereas the template fits result in sigma_[Fe/H]=0.18. Tests made on the RR Lyrae database of the Large Magellanic Cloud (LMC) of the Optical Gravitational Lensing Experiment (OGLE) support the applicability of the method on real photometric time series. These tests also show that the dominant part of error in estimating the average brightness comes from other sources, most probably from crowding effects, even for under-sampled light curves.

Compact object clusters are likely to exist in the centre of some galaxies because of mass segregation. The high densities and velocities reached in them deserves a better understanding. The formation of binaries and their subsequent merging by gravitational radiation emission is important to the evolution of such clusters. We address the evolution of such a system in a relativistic regime. The recurrent mergers at high velocities create an object with a mass much larger than the average. For this aim we modified the direct {\sc Nbody6}++ code to include post-Newtonian effects to the force during two-body encounters. We adjusted the equations of motion to include for the first time the effects of both periastron shift and energy loss by emission of gravitational waves and so to study the eventual decay and merger of radiating binaries. The method employed allows us to give here an accurate post-Newtonian description of the formation of a run-away compact object by successive mergers with surrounding particles, as well as the distribution of characteristic eccentricities in the events. This study should be envisaged as a first step towards a detailed, accurate study of possible gravitational waves sources thanks to the combination of the direct {\sc Nbody} numerical tool with the implementation of post-Newtonian terms on it.