Several binary black holes (BBHs) have been observed using gravitational wave detectors. For the formation mechanism of BBHs, two main mechanisms, isolated binary evolution and dynamical formation in dense star clusters, have been suggested. Future observations are expected to provide more information about BBH distributions, and it will help us to distinguish the two formation mechanisms. For the star cluster channel, globular clusters have mainly been investigated. However, recent simulations have suggested that BBH formation in open clusters is not negligible. We estimate a local merger rate density of BBHs originated from open clusters using the results of our $N$-body simulations of open clusters with four different metallicities. We find that the merger rate per cluster is the highest for our 0.1 solar metallicity model. Assuming a cosmic star formation history and a metallicity evolution with dispersion, we estimate the local merger rate density of BBHs originated from open clusters to be $\sim 70~{\rm yr}^{-1} {\rm Gpc}^{-3}$. This value is comparable to the merger rate density expected from the first and second observation runs of LIGO and Virgo. In addition, we find that BBH mergers obtained from our simulations can reproduce the distribution of primary mass and mass ratio of merging BBHs estimated from the LIGO and Virgo observations.

Black hole-main sequence star (BH-MS) binaries are one of the targets of the future data releases of the astrometric satellite {\it Gaia}. They are supposed to be formed in two main sites: a galactic field and star clusters. However, previous work has never predicted the number of BH-MS binaries originating in the latter site. In this paper, we estimate the number of BH-MS binaries formed in open clusters and detectable with {\it Gaia} based on the results of {\it N}-body simulations. By considering interstellar extinction in the Milky Way (MW) and observational constraints, we predict $\sim 10$ BH-MS binaries are observable. We also find that chemical abundance patterns of companion MSs will help us to identify the origin of the binaries as star clusters. Such MSs are not polluted by outflows of the BH progenitors, such as stellar winds and supernova ejecta. Chemical anomalies might be a good test to confirm the origin of binaries with relatively less massive MSs ($\lesssim 5M_{\odot}$), orbital periods ($\sim 1.5\;$year) and higher eccentricities ($e \gtrsim 0.1$).

We estimate formation rates of LB-1-like systems through dynamical interactions in the framework of the theory of stellar evolution before the discovery of the LB-1 system. The LB-1 system contains $\sim 70M_\odot$ black hole (BH), so-called pair instability (PI)-gap BH, and B-type star with solar metallicity, and has nearly zero eccentricity. The most efficient formation mechanism is as follows. In an open cluster, a naked helium (He) star (with $\sim 20M_\odot$) collides with a heavy main-sequence (MS) star (with $\sim 50M_\odot$) which has a B-type companion. The collision results in a binary consisting of the collision product and B-type star with a high eccentricity. The binary can be circularized through the dynamical tide with radiative damping of the collision-product envelope. Finally, the collision product collapses to a PI-gap BH, avoiding pulsational pair instability and pair instability supernovae because its He core is as massive as the pre-colliding naked He star. We find that the number of LB-1-like systems in the Milky Way galaxy is $\sim 0.01 (\rho_{\rm oc} / 10^4 M_\odot \mbox{pc}^{-3})$, where $\rho_{\rm oc}$ is the initial mass densities of open clusters. If we take into account LB-1-like systems with O-type companion stars, the number increases to $\sim 0.03 (\rho_{\rm oc} / 10^4 M_\odot \mbox{pc}^{-3})$. This mechanism can form LB-1-like systems at least 10 times more efficiently than the other mechanisms: captures of B-type stars by PI-gap BHs, stellar collisions between other type stars, and stellar mergers in hierarchical triple systems. We conclude that no dynamical mechanism can explain the presence of the LB-1 system.

The origin of the Hercules stream, the most prominent velocity substructure in the Solar neighbour disc stars, is still under debate. Recent accurate measurements of position, velocity, and metallicity provided by Tycho Gaia Astrometric Solution (TGAS) and RAdial Velocity Experiments (RAVE) have revealed that the Hercules stream is most clearly seen in the metal-rich region ([Fe/H] > 0), while it is not clearly seen in lower metallicity region ([Fe/H] < -0.25). By using a large number of chemo-dynamical 2D test-particle simulations with a rotating bar and/or spiral arms, we find that the observed [Fe/H] dependence of the Hercules stream is a natural consequence of the inside-out formation of the stellar disc and the existence of highly non-closed orbits in the rotating frame of the bar or spiral arms. Our successful models that reproduce the observed properties of the Hercules stream include not only fast-bar-only and fast-bar+spiral models, but also slow-bar+spiral models. This indicates that it is very difficult to estimate the pattern speed of the bar or spiral arms based only on the observations of the Hercules stream in the Solar neighbourhood. As a by-product of our simulations, we make some predictions about the locations across the Galactic plane where we can observe velocity bimodality that is not associated with the Hercules stream. These predictions can be tested by the Gaia Data Release 2, and such a test will improve our understanding of the evolution of the Milky Way stellar disc.

In order to investigate the formation rate of binary black holes (BBHs) in stellar clusters with a mass comparable to open clusters, we performed a series of direct $N$-body simulations of open clusters with a mass of $2.5\times10^3$ (Model A) and $10^4 M_{\odot}$ (Model B). Since such low-mass clusters would have been more populous than globular clusters when they were born, low-mass clusters are also candidates as the origin of BBHs which are the source of the gravitational waves. In model A, most of BBHs merged within 10 Gyr formed via dynamically formed main-sequence binary stars and stable and unstable mass transfer between them since open clusters collapse within the main-sequence life-time of massive stars. These binaries, therefore, have little eccentricities. The fraction of such binaries among all merging BBHs increases as the cluster mass decrease due to the shorter relaxation time. In our simulations, $4.0\times10^{-5}$ and $1.7\times10^{-5}$ BBHs per solar mass merged within 10 Gyr for models A and B, respectively. These values correspond to $\sim$ 20-50% of the number of mergers per solar mass originated from globular clusters with a mass of $10^5$-$10^6M_{\odot}$. Thus, the contribution of BBHs originated from open clusters is not negligible. The estimated mergers rate density in the local universe is about 0.3 yr$^{-1}$ Gpc$^{-3}$ assuming a cluster mass function with a power of $-$2.

Observations of stars in the the solar vicinity show a clear tendency for old stars to have larger velocity dispersions. This relation is called the age-velocity dispersion relation (AVR) and it is believed to provide insight into the heating history of the Milky Way galaxy. Here, in order to investigate the origin of the AVR, we performed smoothed particle hydrodynamic simulations of the self-gravitating multiphase gas disks in the static disk-halo potentials. Star formation from cold and dense gas is taken into account, and we analyze the evolution of these star particles. We find that exponents of simulated AVR and the ratio of the radial to vertical velocity dispersion are close to the observed values. We also find that the simulated AVR is not a simple consequence of dynamical heating. The evolution tracks of stars with different epochs evolve gradually in the age-velocity dispersion plane as a result of: (1) the decrease in velocity dispersion in star forming regions, and (2) the decrease in the number of cold/dense/gas as scattering sources. These results suggest that the AVR involves not only the heating history of a stellar disk, but also the historical evolution of the ISM in a galaxy.

The formation and evolution of galactic spiral arms is not yet clearly understood despite many analytic and numerical work. Recently, a new idea has been proposed that local density enhancements (waklets) arising in the galactic disk connect with each other and make global spiral arms. However, the understanding of this mechanism is not yet sufficient. We analyze the interaction of wakelets by using N-body simulations including perturbing point masses, which are heavier than individual N-body particles and act as the seeds for wakelets. Our simulation facilitates more straightforward interpretation of numerical results than previous work by putting a certain number of perturbers in a well-motivated configuration. We detected a clear sign of non-linear interaction between wakelets, which make global spiral arms by connecting two adjacent wakelets. We found that the wave number of the strongest non-linear interaction depends on galactic disk mass and shear rate. This dependence is consistent with the prediction of swing amplification mechanism and other previous results. Our results provide unification of previous results which seemed not consistent with each other.