In this paper I present an overview of the main observational properties of a special class of exotic objects (the so-called Blue Straggler Stars, BSSs) in Galactic Globular Clusters (GCs). The BSS specific frequency and their radial distribution are discussed in the framework of using this stellar population as probe of GC internal dynamics. In particular, the shape of the BSS radial distribution has been found to be a powerful tracer of the dynamical evolution of stellar systems, thus allowing the definition of an empirical “clock”able to measure the dynamical age of stellar aggregates from pure observational properties.

The study of intermediate-mass black holes (IMBHs) is a young and promising field of research. If IMBHs exist, they could explain the rapid growth of supermassive black holes by acting as seeds in the early stage of galaxy formation. Formed by runaway collisions of massive stars in young and dense stellar clusters, intermediate-mass black holes could still be present in the centers of globular clusters, today. Our group investigated the presence of intermediate-mass black holes for a sample of 10 Galactic globular clusters. We measured the inner kinematic profiles with integral-field spectroscopy and determined masses or upper limits of central black holes in each cluster. In combination with literature data we further studied the positions of our results on known black-hole scaling relations (such as M• − σ) and found a similar but flatter correlation for IMBHs. Applying cluster evolution codes, the change in the slope could be explained with the stellar mass loss occurring in clusters in a tidal field over its life time. Furthermore, we present results from several numerical simulations on the topic of IMBHs and integral field units (IFUs). We ran N-body simulations of globular clusters containing IMBHs in a tidal field and studied their effects on mass-loss rates and remnant fractions and showed that an IMBH in the center prevents core collapse and ejects massive objects more rapidly. These simulations were further used to simulate IFU data cubes. For the specific case of NGC 6388 we simulated two different IFU techniques and found that velocity dispersion measurements from individual velocities are strongly biased towards lower values due to blends of neighboring stars and background light. In addition, we use the Astrophysical Multipurpose Software Environment (AMUSE) to combine gravitational physics, stellar evolution and hydrodynamics to simulate the accretion of stellar winds onto a black hole.

I present an overview of our ongoing project aimed at building a new generation of velocity dispersion profiles ad rotation curves for a representative sample of Galactic globular clusters, from the the radial velocity of hundreds of individual stars distributed at different distances from the cluster center. The innermost portion of the profiles will be used to constrain the possible presence of intermediate-mass black holes. The adopted methodology consists of combining spectroscopic observations acquired with three different instruments at the ESO-VLT: the adaptive-optics assisted, integral field unit (IFU) spectrograph SINFONI for the innermost and highly crowded cluster cores, the multi-IFU spectrograph KMOS for the intermediate regions, and the multi-fiber instrument FLAMES/GIRAFFE-MEDUSA for the outskirts. The case of NGC 6388, representing the pilot project that motivated the entire program, is described in some details.

Finding an intermediate-mass black hole (IMBH) in a globular cluster (GC), or proving its absence, is a crucial ingredient in our understanding of galaxy formation and evolution. The challenge is to identify a unique signature of an IMBH that cannot be accounted for by other processes. Observational claims of IMBH detection are often based on analyses of the kinematics of stars, such as a rise in the velocity dispersion profile towards the centre. In this contribution we discuss the degeneracy between this IMBH signal and pressure anisotropy in the GC. We show that that by considering anisotropic models it is possible to partially explain the innermost shape of the projected velocity dispersion profile, even though models that do not account for an IMBH do not exhibit a cusp in the centre.

The recent studies show that the formation and evolution process of the nearby galaxies are still unclear. By using the Canada France Hawaii Telescope (CFHT) 3.6m telescope, the PanDAS shows complicated substructures (dwarf satellite galaxies, halo globular clusters, extended clusters, star streams, etc.) in the halo of M31 to ~150 kpc from the center of galaxy and M31-M33 interaction has been studied. In our work, we would like to investigate formation, evolution and interaction of M31 and M33, which are the nearest two spiral galaxies in Local Group. The star cluster systems of the two galaxies are good tracers to study the dynamics of the substructures and the interaction. Since 2010, the Xinglong 2.16m, Lijiang 2.4m and MMT 6.5m telescopes have been used for our spectroscopic observations. The radial velocities and Lick absorption-line indices can thus be measured with the spectroscopy and then ages, metallicities and masses of the star clusters can be fitted with the simple stellar population models. These parameters could be used as the input physical parameters for numerical simulations of M31-M33 interaction.

Dynamical interactions that take place between objects in dense stellar systems lead to frequent formation of exotic stellar objects, unusual binaries, and systems of higher multiplicity. They are most important for the formation of binaries with neutron stars and black holes, which are usually observationally revealed in mass-transferring binaries. Here we review the current understanding of compact object's retention, of the metallicity dependence on the formation of low-mass X-ray binaries with neutron stars, and how mass-transferring binaries with a black hole and a white dwarf can be formed. We discuss as well one old unsolved puzzle and two new puzzles posed by recent observations: what descendants do ultra-compact X-ray binaries produce, how are very compact triples formed, and how can black hole low-mass X-ray binaries acquire non-degenerate companions?

We will discuss the evolution of star clusters with a large initial binary fraction, up to 95%. The initial binary population is chosen to follow the invariant orbital-parameter distributions suggested by Kroupa (1995). The Monte Carlo MOCCA simulations of star cluster evolution are compared to the observations of Milone et al. (2012) for photometric binaries. It is demonstrated that the observed dependence on cluster mass of both the binary fraction and the ratio of the binary fractions inside and outside of the half mass radius are well recovered by the MOCCA simulations. This is due to a rapid decrease in the initial binary fraction due to the strong density-dependent destruction of wide binaries described by Marks, Kroupa & Oh (2011). We also discuss a new scenario for the formation of intermediate mass black holes in dense star clusters. In this scenario, intermediate mass black holes are formed as a result of dynamical interactions of hard binaries containing a stellar mass black hole, with other stars and binaries. We will discuss the necessary conditions to initiate the process of intermediate mass black hole formation and the dependence of its mass accretion rate on the global cluster properties.

The detection of intermediate mass black holes (IMBHs) in globular clusters has been hotly debated, with different observational methods delivering different outcomes for the same object. In order to understand these discrepancies, we construct detailed mock integral field spectroscopy (IFU) observations of globular clusters, starting from realistic Monte Carlo cluster simulations. The output is a data cube of spectra in a given field-of-view that can be analyzed in the same manner as real observations and compared to other (resolved) kinematic measurement methods. We show that the main discrepancies arise because the luminosity-weighted IFU observations can be strongly biased by the presence of a few bright stars that introduce a scatter in velocity dispersion measurements of several km s−1. We show that this intrinsic scatter can prevent a sound assessment of the central kinematics, and therefore should be fully taken into account to correctly interpret the signature of an IMBH.

We present ETICS, a collisionless N-body code based on two kinds of series expansions of the Poisson equation, implemented for graphics processing units (GPUs). The code is publicly available and can be used as a standalone program or as a library (an AMUSE plugin is included). One of the two expansion methods available 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 methods and optimized for collisionless galactic dynamics, but while SCF is a “pure” expansion, MEX is an expansion in just the angular part; thus, MEX is capable of capturing radial structure easily, while SCF needs a large number of radial terms.

The role of stability in the general three-body problem is investigated with regard to the tidal radius of a globular cluster (GC) in a galactic potential. This proceedings is a summary of two papers which outline the stability method (Kennedy 2014a) and compare the predicted stability boundary radius to observations of velocity dispersion profiles in Milky Way GCs (Kennedy 2014b).

In the solar neighborhood, where the typical relaxation timescale is larger than the cosmic age, at least 10% to 15% of Sun-like stars have planetary systems with Jupiter-mass planets. In contrast, dense star clusters, characterized by frequent close encounters, have been found to host very few planets. We carry out numerical simulations with different initial conditions to investigate the dynamical stability of planetary systems in star cluster environments.

We have performed a survey of N-body simulations to explore the dynamics of primordial binaries in multiple-population globular clusters. We show that, as a consequence of the initial differences between the spatial distribution of first-generation (FG) and second-generation (SG) stars, SG binaries are disrupted more efficiently than FG binaries. The effects of dynamical evolution on the surviving binaries produces a difference between the SG and the FG binary binding energy distribution with the SG population characterized by a larger fraction of high binding energy (more bound) binaries. We also explore the evolution of the radial variation of the SG-to-FG binary number ratio and find that although the global binary fraction decreases more rapidly for the SG population, the local binary fraction measured in the cluster inner regions may still be dominated by SG binaries.

We present evolutionary models of rotating self-gravitating systems (e.g. globular clusters, galaxy cores). These models are characterized by the presence of initial axisymmetry due to rotation. Central black hole seeds are alternatively included in our models, and black hole growth due to consumption of stellar matter is simulated until the central potential dominates the kinematics in the core. Goal is to study the long-term evolution (~ Gyr) of relaxed dense stellar systems, which deviate from spherical symmetry, their morphology and final kinematics. With this purpose, we developed a 2D Fokker-Planck analytical code, which results we confirm by detailed N-Body techniques, applying a high performance code, developed for GPU machines. We compare our models to available observations of galactic rotating globular clusters, and conclude that initial rotation modifies significantly the shape and lifetime of these systems, and can not be neglected in studying the evolution of globular clusters, and the galaxy itself.

The Arches cluster is a young, compact, and massive star cluster located in ~30pc away from the Galactic Center in projection. The cluster is located in the extreme environment of the Galactic Center, making it an excellent target for understanding the effects of star-forming environment on the mass function of star clusters. In this study, we estimate the initial condition (mass, concentration parameter, and galactocentric radius) of the Arches cluster by comparing Fokker-Planck calculations with observed velocity dispersion, surface density and mass function data.

We use photometry in the F220W, F250W, F330W, F435W filters from the High Resolution Channel of the Advanced Camera for Surveys and photometry in the F555W, F675W, and F814W filters from the Wide Field and Planetary Camera 2 aboard the Hubble Space Telescope to derive individual stellar reddenings and extinctions for member stars in the HD 97950 cluster in the giant H ii region NGC 3603. Within the standard deviation associated with E(λ−F555W)/E(F435W−F555W) in each filter, the cluster extinction curve at ultraviolet wavelengths tends to be greyer than the average Galactic extinction laws from Cardelli et al. (1989) and Fitzpatrick et al. (1999). It is closer to the extinction law derived by Calzetti et al. (2000) for starburst galaxies, where the 0.2175 μm bump is absent.

Globular clusters (GCs) is dominated by behaviors of high-mass components such as neutron stars (NSs) and black holes (BHs). In this work, we perform direct N-body simulations of a GC assuming two BH masses. We used the BH masses of 10 and 20 M⊙ with total math ratio between those two populations is assumed to be 2:1 and 5:1. Our results show that the heavier BHs (20 M⊙) are depleted in the early stage of cluster evolution. The existence of heavier BH components increase the retention rate of lighter BHs during the cluster evolution. About 30% of ejected BHs are in binaries.

We interpret the radio/X-ray correlation of LR ∝ LX~1.4 for LX/LEdd ≳ 10−3 with a detailed disk corona-jet model, in which the accretion flow and the jet are connected by a parameter, η, describing the fraction of the matter in the accretion flow ejected outward to form the jet. We calculate LR and LX at different Ṁ, adjusting η to fit the observed radio/X-ray correlation of the black hole X-ray transient H1743-322 for LX/LEdd > 10−3. It is found that the value of η for this radio/X-ray correlation for LX/LEdd > 10−3, is systematically less than that of the case for LX/LEdd < 10−3, which is consistent with the general idea that the jet is often relatively suppressed at the high luminosity phase in black hole X-ray binaries.

We explore X-ray spectral evolution and radio–X-ray correlation simultaneously for four X-ray binaries (XRBs). We find that hard X-ray photon indices, Γ, are anti- and positively correlated to X-ray fluxes when the X-ray flux, F3–9keV, is below and above a critical flux, FX,crit, which may be regulated by ADAF and disk-corona respectively. We find that the data points with anti-correlation of Γ-F3–9keV follow the universal radio–X-ray correlation of FR ∝ FXb (b ~ 0.5-0.7), while the data points with positive X-ray spectral evolution follow a steeper radio–X-ray correlation (b ~ 1.4, the so-called ‘outliers track’). The bright active galactic nuclei (AGNs) share similar X-ray spectral evolution and radio–X-ray correlation as XRBs in ‘outliers’ track, and we present a new fundamental plane of log LR=1.59+0.28−0.22 log LX−0.22+0.19−0.20 log MBH−28.97+0.45−0.45 for these radiatively efficient BH sources.

The direct-summation of N gravitational forces is a complex problem for which there is no analytical solution. Dense stellar systems such as galactic nuclei and stellar clusters are the loci of different interesting problems. In this work we present a new GPU, direct-summation N-body integrator written from scratch and based on the Hermite scheme. The first release of the code consists of the Hermite integrator for a system of N bodies with softening. We find an acceleration factor of about ≈ 90 of the GPU version in a single node as compared to the Serial-Single-CPU one. We additionally investigate the impact of using softening in the dynamics of a dense cluster. We study how it affects the two body relaxation, as compared with another code, NBODY6, which uses KS regularization, so as to understand the role of softening in the evolution of the system. This initial release is the first step towards more and more realistic scenarios, starting for a proper treatment for binary evolution, close encounters and the role of a massive black hole.

We are developing a software library to calculate gravitational interaction for the Hermite scheme on parallel computing systems supported by OpenCL API. Our library is partly compatible with a standard GRAPE-6A interface and is easily usable in existing N-body codes. Since our library is based on OpenCL standard API, our library is working on many parallel computing systems such as a multi-core CPU, a GPU, and a many-core architecture. We report the performance evaluation of our library on computing platforms from various vendors.