It is difficult to find or identifying the binary black holes in parsec scales, since the dual AGN may be merged quickly. It is required to explore more possibilities to identifying binary black holes in parsec scales, we give some discussions, especially with the VLBI methods.

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.

The circum-nuclear region in an active galaxy is often complex with presence of high excitation gas, collimated radio outflow, and star formation activities, besides the actively accreting supermassive black hole. The unique spatial resolving power of Chandra X-ray imaging spectroscopy enables more investigations to disentangle the active galactic nuclei and starburst activities. For galaxies in the throes of a violent merging event such as NGC6240, we were able to resolve the high temperature gas surrounding its binary active black holes and discovered a large scale soft X-ray halo.

I review the evolution of binary supermassive black holes and focus on the stellar-dynamical mechanisms that may help to overcome the final-parsec problem – the possible stalling of the binary at a separation much larger than is required for an efficient gravitational wave emission. Recent N-body simulations have suggested that a departure from spherical symmetry in the nucleus of the galaxy may keep the rate of interaction of stars with the binary at a high enough level so that the binary continues to shrink rather rapidly. However, a major problem of all these simulations is that they do not probe the regime where collisionless effects are dominant – in other words, the number of particles in the simulation is still not sufficient to reach the asymptotic behavior of the system. I present a novel Monte Carlo method for simulating both collisional and collisionless evolution of non-spherical stellar systems, and apply it for the problem of binary supermassive black hole evolution. I show that in triaxial galaxies the final-parsec problem is largely non-existent, while in the axisymmetric case it seems to still exist in the limit of purely collisionless regime relevant for real galaxies, but disappears in the N-body simulations where the feasible values of N are still too low to get rid of collisional effects.

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.

Almost all galaxies along the Hubble sequence host a compact massive object (CMO) in their center. The CMO can be either a supermassive black hole (SMBH) or a very dense stellar cluster, also known as nuclear star cluster (NSC). Generally, heavier galaxies (mass ≳ 1011 M⊙) host a central SMBH while lighter show a central NSC. Intermediate mass hosts, instead, contain both a NSC and a SMBH. One possible formation mechanisms of a NSC relies on the dry-merger (migratory) scenario, in which globular clusters (GCs) decay toward the center of the host galaxy and merge. In this framework, the absence of NSCs in high-mass galaxies can be imputed to destruction of the infalling GCs by the intense tidal field of the central SMBH. In this work, we report preliminary results of N-body simulations performed using our high-resolution, direct, code HiGPUs, to investigate the effects of a central SMBH on a single GC orbiting around it. By varying either the mass of the SMBH and the mass of the host galaxy, we derived an upper limit to the mass of the central SMBH, and thus to the mass of the host, above which the formation of a NSC is suppressed.

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.

Nuclear stellar clusters (NSCs) are dense stellar systems known to exist at the center of most of the galaxies. Some of them host a central massive black hole (MBH). They are though to form through in-situ star formation following the infall of gas to the galactic center and/or because of the infall and merger of several stellar clusters. Here we explore the latter scenario by means of detailed self-consistent N-body simulations, proving that a NSC built by the infall and following merger of stellar clusters shows many of the observed features of the Milky Way NSC. We also explore the possibility that the infalling clusters host intermediate mass black holes (IMBHs). Once decayed to the center, the IMBHs act as massive-perturbers accelerating the relaxation of the NSC, filling the loss-cone and boosting the tidal disruption rate of stars up to a value larger than the observational estimates, therefore providing a cumulative constraint on the existence of IMBHs in NSCs. Studying how the properties of the infalling clusters map to the properties of the resulting NSC, we find that, in the IMBHs-free case, the infall mechanism is able to produce many different observational signatures in the form of age segregation.

Galactic white dwarf binaries will be an abundant source of gravitational waves in the mHz frequency band of space-based detectors such as eLISA. A few thousand to a few tens of thousands of these systems will be individually resolvable by eLISA, depending on the final detector configuration. The remaining tens of millions of close white dwarf binaries will create an unresolvable anisotropic foreground of gravitational waves that will be comparable to the instrument noise of eLISA at frequencies below about a mHz. Both the resolvable binaries and the foreground can be used to better understand this population. Careful choice of the initial orientation of eLISA can mitigate this foreground in searches for other sources.

We present a set of, large scale direct N-body simulations of the galaxy collision with the central Supermassive Black Hole Binary (SMBHB) system. Based on our simulations which include the accurate Post Newtonian (PN) relativistic dynamical corrections we can estimated the merging time for the real astrophysical object. Each galaxy initially was represented as a set of particles (up to N=500k) with Plummer distribution. The SMBHBs system is described using the two special high mass, i.e. “relativistic”, particles. The interaction between these two particles have an extra PN correction terms (up to 3.5PN). Merging time upper limit was obtained for the closely interacting galaxy system NGC 6240.

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.