Fragmentation and fission of giant molecular clouds occasionally results in a pair of gravitationally bound star clusters that orbit their mutual centre of mass for some time, under the influence of internal and external perturbations. We investigate the evolution of binary star clusters with different orbital configurations, with a particular focus on the Galactic tidal field. We carry out $N$-body simulations of evolving binary star clusters and compare our results with estimates from our semi-analytic model. The latter accounts for mass loss due to stellar evolution and two-body relaxation, and for evolution due to external tides. Using the semi-analytic model we predict the long-term evolution for a wide range of initial conditions. It accurately describes the global evolution of such systems, until the moment when a cluster merger is imminent. $N$-body simulations are used to test our semi-analytic model and also to study additional features of evolving binary clusters, such as the kinematics of stars, global cluster rotation, evaporation rates, and the cluster merger process. We find that the initial orientation of a binary star cluster with respect to the Galactic field, and also the initial orbital phase, are crucial for its fate. Depending on these properties, the binaries may experience orbital reversal, spiral-in, or vertical oscillation about the Galactic plane before they actually merge at $t\approx100$~Myr, and produce rotating star clusters with slightly higher evaporation rates. The merger process of a binary cluster induces an outburst that ejects $\sim10\%$ of the stellar members into the Galactic field.

The bulk motion of star clusters can be determined after careful membership analysis using parametric or non-parametric approaches. This study aims to implement non-parametric membership analysis based on Binned Kernel Density Estimator which accounts measurements errors (simply called BKDE-e) and to determine the average proper motion of each cluster. This method is applied to 178 selected star clusters with angular diameter less than 20 arc minutes. Proper motion data from UCAC4 are used for membership determination. Non-parametric analysis using BKDE-e successfully determined the average proper motion of 129 clusters, with good accuracy. Compared to COCD and NCOVOCC, there are 79 clusters with less than $3\sigma$ difference. Moreover, we are able to analyse distribution of the member stars in vector point diagram which is not always normal distribution.

Stellar membership determination of an open cluster is an important process to do before further analysis. Basically, there are two classes of membership determination method: parametric and non-parametric. In this study, an alternative of non-parametric method based on Binned Kernel Density Estimation that accounts measurements errors (simply called BKDE-e) is proposed. This method is applied upon proper motions data to determine cluster's membership kinematically and estimate the average proper motions of the cluster. Monte Carlo simulations show that the average proper motions determination using this proposed method is statistically more accurate than ordinary Kernel Density Estimator (KDE). By including measurement errors in the calculation, the mode location from the resulting density estimate is less sensitive to non-physical or stochastic fluctuation as compared to ordinary KDE that excludes measurement errors. For the typical mean measurement error of 7 mas/yr, BKDE-e suppresses the potential of miscalculation by a factor of two compared to KDE. With median accuracy of about 93%, BKDE-e method has comparable accuracy with respect to parametric method (modified Sanders algorithm). Application to real data from The Fourth USNO CCD Astrograph Catalog (UCAC4), especially to NGC 2682 is also performed. The mode of member stars distribution on Vector Point Diagram is located at $\mu_{\alpha}\cos\delta=-9.94\pm0.85$ mas/yr and $\mu_{\delta}=-4.92\pm0.88$ mas/yr. Although the BKDE-e performance does not overtake parametric approach, it serves a new view of doing membership analysis, expandable to astrometric and photometric data or even in binary cluster search.

Binarity among stellar clusters in galaxy is such a reality which has been realized for a long time, but still hides several questions and problems to be solved. Some of binary star clusters are formed by close encounter, but the others are formed together from similar womb. Some of them undergo separation process, while the others are in the middle of merger toward common future. The products of merger binary star cluster have typical characteristics which differ from solo clusters, especially in their spatial distribution and their stellar members kinematics. On the other hand, these merger products still have to face dissolving processes triggered by both internal and external factors. In this study, we performed N-body simulations of merger binary clusters with different initial conditions. After merging, these clusters dissolve with greater mass-loss rate because of their angular momentum. These rotating clusters also experience more deceleration caused by external tidal field.

We analysed the fine structure of the phase space distribution function of nearby subdwarfs using data extracted from various catalogues. Applying a new search strategy based on Dekker's theory of galactic orbits, we found four overdensely populated regions in phase space. Three of them were correlated with previously known star streams: the Hyades-Pleiades and Hercules streams in the thin disk of the Milky Way and the Arcturus stream in the thick disk. In addition we find evidence for another stream in the thick disk, which resembles closely the Arcturus stream and probably has the same extragalactic origin.

We present an analysis of the space motions of 742 subdwarf stars based on the sample of Carney et al. (1994, CLLA). Hipparcos parallaxes, TYC2+HIP proper motions and Tycho2 proper motions were combined with radial velocities and metallicities from CLLA. The kinematical behavior is discussed in particular in relation to their metallicities. The majority of these sample stars have metal abundances of [Fe/H] >-1 and represent the thick disk population. The halo component, with [Fe/H] <-1.6, is characterized by a low mean rotation velocity and a radially elongated velocity ellipsoid. In the intermediate metallicity range (-1.6 < [Fe/H] <-1), we find a significant number of subdwarfs with disklike kinematics. We interpret this population of stars as a metal-weak thick disk population.