Nuclear stellar clusters (NSCs) are known to exist around massive black holes (MBHs) in galactic nuclei. Two formation scenarios were suggested for their origin: build-up of NSCs and Continuous in-situ star-formation. Here we study the effects of star formation on the build-up of NSCs and its implications for their long term evolution and their resulting structure. We show that continuous star-formation can lead to the build-up of an NSC with properties similar to those of the Milky-way NSC. We also find that the general structure of the old stellar population in the NSC with in-situ star-formation could be very similar to the steady-state Bahcall-Wolf cuspy structure. However, its younger stellar population does not yet achieve a steady state. In particular, formed/evolved NSCs with in-situ star-formation contain differential age-segregated stellar populations which are not yet fully mixed. Younger stellar populations formed in the outer regions of the NSC have a cuspy structure towards the NSC outskirts, while showing a core-like distribution inwards; with younger populations having larger core sizes.

We report on GalevNB (Galev for N-body simulations), an integrated software solution that provides N-body users direct access to the software package GALEV (GALaxy EVolutionary synthesis models). GalevNB is developed for the purpose of a direct comparison between N-body simulations and observations. It converts the fundamental stellar properties of N-body simulations, i.e., stellar mass, temperature, stellar luminosity and metallicity, into observational magnitudes for a variety of filters of widely used instruments/telescopes (HST, ESO, SDSS, 2MASS), and into spectra that span from far-UV (90 Å) to near-IR (160 μm).

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.

We present the estimates of the black hole spins of five quasars. The peaks of the spectra of the accretion discs surrounding massive black holes in quasars are in the far-UV or soft X-ray band, which are usually not observed. However, in the disc corona model, the soft photons from the disc are Comptonized to high energy in the hot corona, and the hard X-ray spectra (luminosity and spectral shape) contain the information of the incident spectra from the disc. The values of black hole spin parameter a are inferred from the spectral fitting, which spread over a large range, ~ −0.94 to 0.998.

Discovery of ultra-compact dwarfs (UCDs) in the past 15 years blurs the once thought clear division between classic globular clusters (GCs) and early-type galaxies. The intermediate nature of UCDs, which are larger and more massive than typical GCs but more compact than typical dwarf galaxies, has triggered hot debate on whether UCDs should be considered galactic in origin or merely the most extreme GCs. Previous studies of various scaling relations, stellar populations and internal dynamics did not give an unambiguous answer to the primary origin of UCDs. In this contribution, we present the first ever detailed study of global dynamics of 97 UCDs (rh ≳ 10 pc) associated with the central cD galaxy of the Virgo cluster, M87. We found that UCDs follow a different radial number density profile and different rotational properties from GCs. The orbital anisotropies of UCDs are tangentially-biased within ~ 40 kpc of M87 and become radially-biased with radius further out. In contrast, the blue GCs, which have similar median colors to our sample of UCDs, become more tangentially-biased at larger radii beyond ~ 40 kpc. Our analysis suggests that most UCDs in M87 are not consistent with being merely the most luminous and extended examples of otherwise normal GCs. The radially-biased orbital structure of UCDs at large radii is in general agreement with the scenario that most UCDs originated from the tidally threshed dwarf galaxies.

We made dynamical black hole mass measurements from nineteen Seyfert 2 galaxies which host sub-parsec H2O maser disks using the H2O megamaser technique. The nearly perfect Keplerian rotation curves in many of these maser systems guarantee the high accuracy and precision of the black hole mass measurements. With the stellar velocity dispersion (σ∗) of the galaxy bulges measured with the Dupont 2.5 m telescope at Las Campanas Observatory in the South and the Apache Point Observatory (APO) 3.5m telescope in the North, we found that H2O maser galaxies, most of which host pseudo bulges rather than classical bulges, do not all follow the MBH–σ∗ relation shown in the literature. This result is well consistent with the latest findings by Kormendy & Ho (2013) that only early type galaxies and galaxies with classical bulges follow a tight MBH–σ∗ relation. Such a tight correlation may not exist in pseudo bulge galaxies.

We aim to investigate the formation of sub-galactic structure in the Lambda cold dark matter (CDM) cosmology. To accomplish our research goal, we have added various baryonic physics on the existing cosmological hydrodynamic code, GADGET-2. We performed two test runs to check our new implementations. We show our preliminary results from these test runs.

Supermassive black hole binaries (SMBHBs) are the products of frequent galaxy mergers. It is very hard to be detected in quiescent galaxy. By using one million particle direct N-body simulations on special many-core hardware (GPU cluster), we study the dynamical co-evolution of SMBHB and its surrounding stars, specially focusing on the evolution of stellar tidal disruption event (TDE) rates before and after the coalescence of the SMBHB. We find a boosted TDE rate during the merger of the galaxies. After the coalescence of two supermassive black holes (SMBHs), the post-merger SMBH can get a kick velocity due to the anisotropic GW radiations. Our results about the recoiling SMBH, which oscillates around galactic center, show that most of TDEs are contributed by unbound stars when the SMBH passing through galactic center. In addition, the TDE light curve in SMBHB system is significantly different from the curve for single SMBH, which can be used to identify the SMBHB.

The capture and disruption of stars by supermassive black holes (SMBHs), and the formation and coalescence of binaries, are inevitable consequences of the presence of SMBHs at the cores of galaxies. Pairs of active galactic nuclei (AGN) and binary SMBHs are important stages in the evolution of galaxy mergers, and an intense search for these systems is currently ongoing. In the early and advanced stages of galaxy merging, observations of the triggering of accretion onto one or both BHs inform us about feedback processes and BH growth. Identification of the compact binary SMBHs at parsec and sub-parsec scales provides us with important constraints on the interaction processes that govern the shrinkage of the binary beyond the “final parsec”. Coalescing binary SMBHs are among the most powerful sources of gravitational waves (GWs) in the universe. Stellar tidal disruption events (TDEs) appear as luminous, transient, accretion flares when part of the stellar material is accreted by the SMBH. About 30 events have been identified by multi-wavelength observations by now, and they will be detected in the thousands in future ground-based or space-based transient surveys. The study of TDEs provides us with a variety of new astrophysical tools and applications, related to fundamental physics or astrophysics. Here, we provide a review of the current status of observations of SMBH pairs and binaries, and TDEs, and discuss astrophysical implications.

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

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.

We summarise the recent progress in understanding the formation and evolution of globular clusters (GCs) in the context of galaxy formation and evolution. It is discussed that an end-to-end model for GC formation and evolution should capture four different phases: (1) star and cluster formation in the high-pressure interstellar medium of high-redshift galaxies, (2) cluster disruption by tidal shocks in the gas-rich host galaxy disc, (3) cluster migration into the galaxy halo, and (4) the final evaporation-dominated evolution of GCs until the present day. Previous models have mainly focussed on phase 4. We present and discuss a simple model that includes each of these four steps – its key difference with respect to previous work is the simultaneous addition of the high-redshift formation and early evolution of young GCs, as well as their migration into galaxy haloes. The new model provides an excellent match to the observed GC mass spectrum and specific frequency, as well as the relations of GCs to the host dark matter halo mass and supermassive black hole mass. These results show (1) that the properties of present-day GCs are reproduced by assuming that they are the natural outcome of regular high-redshift star formation (i.e. they form according to same physical processes that govern massive cluster formation in the local Universe), and (2) that models only including GC evaporation strongly underestimate their integrated mass loss over a Hubble time.

There is a supermassive black hole, a gaseous accretion disk and compact star cluster in the center of active galactic nuclei, as known today. So the activity of AGN can be represented as the result of interaction of these three subsystems. In this work we investigate the dynamical interaction of a central star cluster surrounding a supermassive black hole and a central accretion disk. The dissipative force acting on stars in the disk leads to an asymmetry in the phase space distribution of the central star cluster due to the rotating accretion disk. In our work we present some results of Stardisk model, where we see some changes in density and phase space of central star cluster due to influence of rotating gaseous accretion disk.

Narrow-line Seyfert 1 (NLS1) galaxies are a sub-class of active galactic nuclei (AGN) with relatively low-mass black holes, accreting near the Eddington rate. A small fraction of them is radio-loud and harbors relativistic jets. As a class, these provide us with new insights into the cause(s) of radio-loudness, the blazar phenomenon at low black hole masses, and the operation of radio-mode feedback. The NLS1 galaxy RXJ2314.9+2243 is remarkable for its multi-wavelength properties. We present new radio observations taken at Effelsberg, and a summary of the recent results from our multi-wavelength study. RXJ2314.9+2243 is radio-loud, luminous in the infrared, has a flat X-ray spectrum and peculiar UV spectrum, and hosts an exceptionally broad and blueshifted [OIII]λ5007 emission line, indicating the presence of a strong outflow. RXJ2314.9+2243 likely represents an extreme case of AGN induced feedback in the local universe.

In this work we study the stellar-dynamical hardening of unequal mass massive black hole (MBH) binaries in the central regions of galactic nuclei. We present a comprehensive set of direct N-body simulations of the problem, varying both the total mass and the mass ratio of the MBH binary. Our initial model starts as an axisymmetric, rotating galactic nucleus, to describe the situation right after the galaxies have merged, but the black holes are still unbound to each other. We confirm that results presented in earlier works (Berczik et al. 2006; Khan et al. 2013; Wang et al. 2014) about the solution of the “last parsec problem” (sufficiently fast black hole coalescence for black hole growth in cosmological context) are robust for both for the case of unequal black hole masses and large particle numbers. The MBH binary hardening rate depends on the reduced mass ratio through a single parameter function, which quantitatively quite well agrees with standard 3 body scattering theory (see e.g., Hills 1983). Based on our results we conclude that MBH binaries at high redshifts are expected to merge with a factor of ~ 2 more efficiently, which is important to determine the possible overall gravitational wave signals. However, we have not yet fully covered all the possible parameter space, in particular with respect to the preceding of the galaxy mergers, which may lead to a wider variety of initial models, such as initially more oblate and / or even significantly triaxial galactic nuclei. Our N-body simulations were carried out on a new special supercomputers using the hardware acceleration with graphic processing units (GPUs).

There is a growing array of supermassive black hole and nuclear star cluster scaling relations with their host spheroid, including a bent (black hole mass)–(host spheroid mass) Mbh–Msph relation and a different (massive compact object mass)–(host spheroid velocity dispersion) Mmco–σ relations for black holes and nuclear star clusters. By combining the observed Mbh ∝ σ5.5 relation with the observed Mnc ∝ σ1.6–2.7 relation, we derive the expression Mbh ∝ Mnc2–3.4, which should hold until the nuclear star clusters are eventually destroyed in the larger core-Sérsic spheroids. This new mass scaling relation helps better quantify the rapid evolutionary growth of massive black holes in dense star clusters, and the relation is consistently recovered when coupling the observed Mnc ∝ Msph0.6–1.0 relation with the recently observed quadratic relation Mbh ∝ Msph2 for Sérsic spheroids.

Recent studies of active galactic nuclei (AGN) found a statistical inverse scaling between the X-ray normalized excess variance σrms2 (variability amplitude) and the black hole mass spanning over MBH = 106 − 109M⊙. We present a study of this relation by including AGN with MBH = 105 − 106M⊙. It is found that the relation is no longer a simple extrapolation of the known inverse proportion, but starts to flatten around 106M⊙. This behavior can be understood by the shape of the power spectrum density of AGN and its dependence on the black hole mass.

Triple supermassive black holes (SMBH) can form during the hierarchical mergers of massive galaxies with an existing binary. Perturbations by a third black hole may accelerate the merging process of an inner binary, for example through the Kozai mechanism. We analyze the evolution of simulated hierarchical triple SMBHs in galactic centers, and find resonances in the evolution of the semi-major axis, the eccentricity and the inclination, for both the inner and the outer orbits of the triple system, which are not only Kozai like. Through resonant oscillations, SMBH can trigger a significant increase of the inner SMBH binary eccentricity shortening the merger timescale expected from gravitational wave (GW) emission. As hierarchical triple SMBHs may be frequent in massive galaxies, the influence of orbital resonances is of great importance to our understanding of black hole coalescence and gravitational wave detection. Although Kozai mechanism is believed to play an important role in this process, detailed studies on the pattern of these resonances is necessary.