Massive globular clusters lose stars via internal and external processes. Internal processes include mainly two-body relaxation, while external processes include interactions with the Galactic tidal field. We perform a suite of N-body simulations of such massive clusters using three different direct-summation N-body codes, exploring different Galactic orbits and particle numbers. By inspecting the rate at which a star's energy changes as it becomes energetically unbound from the cluster, we can neatly identify two populations we call kicks and sweeps, that escape through two-body encounters internal to the cluster and the external tidal field, respectively. We find that for a typical halo globular cluster on a moderately eccentric orbit, sweeps are far more common than kicks but the total mass loss rate is so low that these clusters can survive for hundreds of Hubble times. The different N-body codes give largely consistent results, but we find that numerical artifacts may arise in relation to the time-step parameter of the Hermite integration scheme, namely that the value required for convergent results is sensitive to the number of particles.

Although their existence is not yet confirmed observationally, intermediate mass black holes (IMBHs) may play a key role in the dynamics of galactic nuclei. In this paper, we neglect the effect of the nuclear star cluster itself and investigate only how a small reservoir of IMBHs influences the secular dynamics of stellar-mass black hole binaries, using N-body simulations. We show that our simplifications are valid and that the IMBHs significantly enhance binary evaporation by pushing the binaries into the Hill-unstable region of parameter space, where they are separated by the SMBH's tidal field. For binaries in the S-cluster region of the Milky Way, IMBHs drive the binaries to merge in up to 1-6% of cases, assuming five IMBHs within 5 pc of mass 10,000 solar masses each. Observations of binaries in the Galactic center may strongly constrain the population of IMBHs therein.

We investigate the internal dynamics of anisotropic, rotating globular clusters with a multimass stellar population by performing new direct N-body simulations. In addition to the well-known radial mass segregation effect, where heavy stars and stellar remnants sink toward the center of the cluster, we find a mass segregation in the distribution of orbital inclinations as well. This newly discovered anisotropic mass segregation leads to the formation of a disk-like structure of massive objects near the equatorial plane of a rotating cluster. This result has important implications on the expected spatial distribution of black holes in globular clusters.

Resonant relaxation has been discussed as an efficient process that changes the angular momenta of stars orbiting around a central supermassive black hole due to the fluctuating gravitational field of the stellar cluster. Other spherical stellar systems, such as globular clusters, exhibit a restricted form of this effect where enhanced relaxation rate only occurs in the directions of the angular momentum vectors, but not in their magnitudes; this is called vector resonant relaxation (VRR). To explore this effect, we performed a large set of direct N-body simulations, with up to 512k particles and ~500 dynamical times. Contrasting our simulations with Spitzer-style Monte Carlo simulations, that by design only exhibit 2-body relaxation, we show that the temporal behavior of the angular momentum vectors in $N$-body simulations cannot be explained by 2-body relaxation alone. VRR operates efficiently in globular clusters with $N>10^4$. The fact that VRR operates in globular clusters may open way to use powerful tools in statistical physics for their description. In particular, since the distribution of orbital planes relaxes much more rapidly than the distribution of the magnitude of angular momentum and the radial action, the relaxation process reaches an internal statistical equilibrium in the corresponding part of phase space while the whole cluster is generally out of equilibrium, in a state of quenched disorder. We point out the need to include effects of VRR in Monte Carlo simulations of globular clusters.

We perform high resolution direct $N$-body simulations to study the effect of an accretion disc on stellar dynamics in an active galactic nucleus (AGN). We show that the interaction of the nuclear stellar cluster (NSC) with the gaseous disc (AD) leads to formation of a stellar disc in the central part of the NSC. The accretion of stars from the stellar disc onto the super-massive black hole is balanced by the capture of stars from the NSC into the stellar disc, yielding a stationary density profile. We derive the migration time through the AD to be 3\% of the half-mass relaxation time of the NSC. The mass and size of the stellar disc are 0.7\% of the mass and 5\% of the influence radius of the super-massive black hole. An AD lifetime shorter than the migration time would result in a less massive nuclear stellar disc. The detection of such a stellar disc could point to past activity of the hosting galactic nucleus.

Collisional relaxation describes the stochastic process with which a self-gravitating system near equilibrium evolves in phase space due to the fluctuating gravitational field of the system. The characteristic timescale of this process is called the relaxation time. In this paper, we highlight the difference between two measures of the relaxation time in globular clusters: (i) the diffusion time with which the isolating integrals of motion (i.e. energy E and angular momentum magnitude L) of individual stars change stochastically and (ii) the asymptotic timescale required for a family of orbits to mix in the cluster. More specifically, the former corresponds to the instantaneous rate of change of a star's E or L, while the latter corresponds to the timescale for the stars to statistically forget their initial conditions. We show that the diffusion timescales of E and L vary systematically around the commonly used half-mass relaxation time in different regions of the cluster by a factor of ~10 and ~100, respectively, for more than 20% of the stars. We define the mixedness of an orbital family at any given time as the correlation coefficient between its E or L probability distribution functions and those of the whole cluster. Using Monte Carlo simulations, we find that mixedness converges asymptotically exponentially with a decay timescale that is ~10 times the half-mass relaxation time.

The Galactic center of the Milky-Way harbors a massive black hole (BH) orbited by a diverse population of young and old stars. A significant fraction of the youngest stars ($\sim4-7$ Myr) reside in a thin stellar disk with puzzling properties; the disk appears to be warped, shows asymmetries, and contains one or more clumpy structures (e.g. IRS 13). Models explaining the clumping invoked the existence of an intermediate-mass BH of $10^{3}-10^{4}$ M$_{\odot}$, but no kinematic evidence for such a BH has been found. Here we use extended $N$-body simulations and hybrid self-consistent field method models to show that naturally formed residual temporal asphericity of the hosting nuclear star cluster gives rise to torques on the disk, which lead to changes in its orientation over time, and to recurrent formation and dissolution of single spiral arm ($m=1$ modes) structures. The changing orientation leads to a flapping-like behavior of the disk and to the formation of a warped disk structure. The spiral arms may explain the over-densities in the disk (clumping) and its observed asymmetry, without invoking the existence of an intermediate-mass BH. The spiral arms are also important for the overall disk evolution and can be used to constrain the structure and composition of the nuclear stellar cluster.

The dynamical formation of stellar-mass black hole-black hole binaries has long been a promising source of gravitational waves for the Laser Interferometer Gravitational-Wave Observatory (LIGO). Mass segregation, gravitational focusing, and multibody dynamical interactions naturally increase the interaction rate between the most massive black holes in dense stellar systems, eventually leading them to merge. We find that dynamical interactions, particularly three-body binary formation, enhance the merger rate of black hole binaries with total mass M_tot roughly as ~M_tot^beta, with beta >~ 4. We find that this relation holds mostly independently of the initial mass function, but the exact value depends on the degree of mass segregation. The detection rate of such massive black hole binaries is only further enhanced by LIGO's greater sensitivity to massive black hole binaries with M_tot <~ 80 solar masses. We find that for power-law BH mass functions dN/dM ~ M^-alpha with alpha <~ 2, LIGO is most likely to detect black hole binaries with a mass twice that of the maximum initial black hole mass and a mass ratio near one. Repeated mergers of black holes inside the cluster result in about ~5% of mergers being observed between two and three times the maximum initial black hole mass. Using these relations, one may be able to invert the observed distribution to the initial mass function with multiple detections of merging black hole binaries.

We analyse the scatter in the correlation between super-massive black hole (SMBH) mass and bulge stellar mass of the host galaxy, and infer that it cannot be accounted for by mergers alone. The merger-only scenario, where small galaxies merge to establish a proportionality relation between the SMBH and bulge masses, leads to a scatter around the linear proportionality line that increases with the square root of the SMBH (or bulge) mass. By examining a sample of 103 galaxies we find that the intrinsic scatter increases more rapidly than expected from the merger-only scenario. The correlation between SMBH masses and their host galaxy properties is therefore more likely to be determined by a negative feedback mechanism that is driven by an active galactic nucleus. We find a hint that some galaxies with missing stellar mass reside close to the centre of clusters and speculate that ram-pressure stripping of gas off the young galaxy as it moves near the cluster centre, might explain the missing stellar mass at later times.

The Laser Interferometer Gravitational Wave Observatory (LIGO) has recently discovered gravitational waves (GWs) emitted by merging black hole binaries. We examine whether future GW detections may identify triple companions of merging binaries. Such a triple companion causes variations in the GW signal due to (1) the varying path length along the line of sight during the orbit around the center of mass, (2) relativistic beaming, Doppler, and gravitational redshift, (3) the variation of the "light"-travel time in the gravitational field of the triple companion, and (4) secular variations of the orbital elements. We find that the prospects for detecting the triple companion are the highest for low-mass compact object binaries which spend the longest time in the LIGO frequency band. In particular, for merging neutron star binaries, LIGO may detect a white dwarf or M-dwarf perturber at signal to noise ratio of 8, if it is within 0.4 solar radius distance from the binary and the system is within a distance of 100 Mpc. Stellar mass (supermassive) black hole perturbers may be detected at a factor 5x (1000x) larger separations. Such pertubers in orbit around the merging binary emit GWs at frequencies above 1 mHz detectable by the Laser Interferometer Space Antenna (LISA) in coincidence.

We examine the effect of an accretion disc on the orbits of stars in the central star cluster surrounding a central massive black hole by performing a suite of 39 high-accuracy direct N-body simulations using state-of-the art software and accelerator hardware, with particle numbers up to 128k. The primary focus is on the accretion rate of stars by the black hole (equivalent to their tidal disruption rate for black holes in the small to medium mass range) and the eccentricity distribution of these stars. Our simulations vary not only the particle number, but disc model (two models examined), spatial resolution at the centre (characterised by the numerical accretion radius) and softening length. The large parameter range and physically realistic modelling allow us for the first time to confidently extrapolate these results to real galactic centres. While in a real galactic centre both particle number and accretion radius differ by a few orders of magnitude from our models, which are constrained by numerical capability, we find that the stellar accretion rate converges for models with N > 32k. The eccentricity distribution of accreted stars, however, does not converge. We find that there are two competing effects at work when improving the resolution: larger particle number leads to a smaller fraction of stars accreted on nearly-circular orbits, while higher spatial resolution increases this fraction. We scale our simulations to some nearby galaxies and find that the expected boost in stellar accretion (or tidal disruption, which could be observed as X-ray flares) in the presence of a gas disc is about a factor of 10. Even with this boost, the accretion of mass from stars is still a factor of ~ 100 slower than the accretion of gas from the disc. Thus, it seems accretion of stars is not a major contributor to black hole mass growth.

Astrophysical research in recent decades has made significant progress thanks to the availability of various $N$-body simulation techniques. With the rapid development of high-performance computing technologies, modern simulations have been able to take the computing power of massively parallel clusters with more than $10^5$ GPU cores. While unprecedented accuracy and dynamical scales have been achieved, the enormous amount of data being generated continuously poses great challenges for the subsequent procedures of data analysis and archiving. As an urgent response to these challenges, in this paper we propose an adaptive storage scheme for simulation data, inspired by the block time step integration scheme found in a number of direct $N$-body integrators available nowadays. The proposed scheme, namely the block time step storage scheme, works by minimizing the data redundancy with assignments of data with individual output frequencies as required by the researcher. As demonstrated by benchmarks, the proposed scheme is applicable to a wide variety of simulations. Despite the main focus of developing a solution for direct $N$-body simulation data, the methodology is transferable for grid-based or tree-based simulations where hierarchical time stepping is used.

We present GPU implementations of two fast force calculation methods, based on series expansions of the Poisson equation. One is the Self-Consistent Field (SCF) method, which is a Fourier-like expansion of the density field in some basis set; the other is the Multipole Expansion (MEX) method, which is a Taylor-like expansion of the Green's function. MEX, which has been advocated in the past, has not gained as much popularity as SCF. Both are particle-field method and optimized for collisionless galactic dynamics, but while SCF is a "pure" expansion, MEX is an expansion in just the angular part; it is thus capable of capturing radial structure easily, where SCF needs a large number of radial terms. We show that despite the expansion bias, these methods are more accurate than direct techniques for the same number of particles. The performance of our GPU code, which we call ETICS, is profiled and compared to a CPU implementation. On the tested GPU hardware, a full force calculation for one million particles took ~ 0.1 seconds (depending on expansion cutoff), making simulations with as many as $10^8$ particles fast on a comparatively small number of nodes.

Supermassive black holes are expected to pair as a result of galaxy mergers, and form a bound binary at parsec or sub-parsec scales. These scales are unresolved even in nearby galaxies, and thus detection of non-active black hole binaries must rely on stellar dynamics. Here we show that these systems could be indirectly detected through the trail that the black holes leave as they spiral inwards. We analyze two numerical simulations of inspiralling black holes (equal masses and 10:1 mass ratio) in the stellar environment of a galactic centre. We studied the effect of the binary on the structure of the stellar population, with particular emphasis on projected kinematics and directly measurable moments of the velocity distribution. We present those moments as high-resolution 2D maps. As shown in past scattering experiments, a torus of stars counter-rotating with respect to the black holes exists in scales ~ 5 to 10 times larger than the binary separation. While this is seen in the average velocity map in the unequal mass case, it is obscured by a more strongly co-rotating outer region in the equal mass case; however, the inner counter-rotation could still be detected by studying the higher moments of the velocity distribution. Additionally, the maps reveal a dip in velocity dispersion in the inner region, as well as more pronounced signatures in the higher distribution moments. These maps could serve as templates for integral field spectroscopy observations of nearby galactic centres. The discovery of such signatures may help census the population of supermassive black hole binaries and refine signal rate predictions for future space-based low frequency gravitational wave detectors.

We examine the intrinsic scatter in the correlation between black hole masses and their host bulge masses, and find that it cannot be accounted for by mergers alone. A simple merger scenario of small galaxies leads to a proportionality relation between the late-time black hole and bulge masses, with intrinsic scatter (in linear scale) increasing along the ridge line of the relation as the square root of the mass. By examining a sample of 86 galaxies with well measured black hole masses, we find that the intrinsic scatter increases with mass more rapidly than expected from the merger-only scenario. We discuss the possibility that the feedback mechanism that operated during galaxy formation involved the presence of a cooling flow.

We present a new approach to studying the evolution of massive black hole binaries in a stellar environment. By imposing conservation of total energy and angular momentum in scattering experiments, we find the dissipation forces that are exerted on the black holes by the stars, and thus obtain the decaying path of the binary from the classical dynamical friction regime down to subparsec scales. Our scheme lies between scattering experiments and N-body simulations. While still resolving collisions between stars and black holes, it is fast enough and allows to use a large enough number of particles to reach a smooth and convergent result. We studied both an equal mass and a 10:1 mass ratio binaries under various initial conditions. We show that while an equal mass binary stalls at a nearly circular orbit, a runaway growth of eccentricity occurs in the unequal mass case. This effect reduces the timescale for black hole coalescence through gravitational radiation to well below the Hubble time, even in spherical and gasless systems formed by dry mergers.

We use a recent sample of 49 galaxies to show that there is a proportionality relation between the black hole mass M_BH and the quantity \mu =M_G*\sigma /c, where M_G is mass of the spheroidal stellar component and \sigma is the stellar velocity dispersion. \mu is called the momentum parameter and the ratio is M_BH/\mu ~3.3. This result is applied to the penetrating-jet feedback model which argues that the correlation that holds is with a momentum-like parameter, although this feedback mechanism is based on energy balance.

The stalling radius of a merging massive binary black hole (BBH) is expected to be below 0".1 even in nearby galaxies (Yu 2002), and thus BBHs are not expected to be spatially resolved in the near future. However, as we show below, a BBH may be detectable through the significantly anisotropic stellar velocity distribution it produces on scales 5-10 times larger than the binary separation. We calculate the velocity distribution of stable orbits near a BBH by solving the restricted three body problem for a BBH embedded in a bulge potential. We present high resolution maps of the projected velocity distribution moments, based on snapshots of ~ 10^8 stable orbits. The kinematic signature of a BBH in the average velocity maps is a counter rotating torus of stars outside the BBH Hill spheres. The velocity dispersion maps reveal a dip in the inner region, and an excess of 20-40% further out, compared to a single BH of the same total mass. More pronounced signatures are seen in the third and fourth Gauss-Hermite velocity moments maps. The detection of these signatures may indicate the presence of a BBH currently, or at some earlier time, which depends on the rate of velocity phase space mixing following the BBH merger.

We conduct 2D numerical simulations of jets expanding into the slow wind of asymptotic giant branch stars. We show that the post-shock jets' material can explain the observed extended X-ray emission from some planetary nebulae (PNs). Such jets are thought to shape many PNs, and therefore it is expected that this process will contribute to the X-ray emission from some PNs. In other PNs (not simulated in this work) the source of the extended X-ray emission is the shocked spherical wind blown by the central star. In a small fraction of PNs both sources might contribute, and a two-temperatures gas will fit better the X-ray properties than a one-temperature gas. A spacial separation between these two components is expected.