Ludendorff coronal flattening index of the Total Solar Eclipse (TSE) on March 9, 2016, was calculated at various distances in solar radius. As a result, we obtained the coronal flattening index $\left(\epsilon =a+b\right)$ at a distance of 2 solar radii is 0.16. The $24^{th}$ solar cycle phase based on the 2016 TSE event obtained -0.64 which showed the corona is pre-minimum type. Resulted coronal flattening index value gives a predicted maximum amplitude of the monthly sunspot number ($W_{max}$) for the $25^{th}$ solar cycle to be $70\pm65$. Therefore, the solar activity for $25^{th}$ solar cycle predicted to be lower than the current solar cycle, which has a maximum sunspot number value of 146 in February 2014

The global structure of solar corona observed in optical window is governed by the global magnetic field with different characteristics over solar activity cycle. Ludendorff flattening index becomes a popular measure of the global structure of solar corona as observed during eclipse. In this study, 15 digital images of solar corona from 1991 to 2016 were analyzed in order to construct the coronal flattening profiles as a function of radius. In most of the cases, the profile can be modeled with 2nd order polynomial function so that the radius with maximum flattening index ($R_{\text{max}}$) can be determined. Along with this value, Ludendorff index ($a+b$) was also calculated. Both Ludendorff index and $R_{\text{max}}$ show anti-correlation with monthly sunspot number, though the $R_{\text{max}}$ values are more scattered. The variation of $R_{\text{max}}$ can be regarded as the impact of changing coronal brightness profile over equator.

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.

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.

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.

Optical observation through sub-meter telescope equipped with CCD camera becomes alternative method for increasing orbital debris detection and surveillance. This observational mode is expected to eye medium-sized objects in higher orbits (e.g. MEO, GTO, GSO \& GEO), beyond the reach of usual radar system. However, such observation of fast-moving objects demands special treatment and analysis technique. In this study, we performed photometric analysis of the satellite track images photographed using rehabilitated Schmidt Bima Sakti telescope in Bosscha Observatory. The Hough transformation was implemented to automatically detect linear streak from the images. From this analysis and comparison to USSPACECOM catalog, two satellites were identified and associated with inactive Thuraya-3 satellite and Satcom-3 debris which are located at geostationary orbit. Further aperture photometry analysis revealed the periodicity of tumbling Satcom-3 debris. In the near future, it is not impossible to apply similar scheme to establish an analysis pipeline for optical space surveillance system hosted in Indonesia.

We present geometrical and kinematical analysis of Moreton waves and EIT waves observed on June 3, 2012 and Moreton waves observed on July 6, 2012. The Moreton waves were recorded in H$\alpha$ images of Global Oscillation Network Group (GONG) archive and EIT waves obtained from SDO/AIA observations, especially in 193 nm channel. The observed wave of June 3 has angular span of about $70^{\circ}$ with a broad wave front associated to NOAA active region 11496. It was found that the speed of the wave that started propagating at 17.53 UT is between 950 to 1500 km/s. Related to this wave occurrence, there was solar type II and III radio bursts. The speed of the EIT in this respect about 247 km/sec. On the other hand, the wave of July 6 may be associated to X1.1 class flare that occurred at 23.01 UT around the 11514 active region. From the kinematical analysis, the wave propagated with the initial velocity of about 1180 km/s which is in agreement with coronal shock velocity derived from type II radio burst observation, $v\sim1200$ km/s.

In this study, we present geometrical and kinematical analysis of Moreton wave observed in 2012 June 3rd and July 6th, recorded in H-$\alpha$ images of Global Oscillation Network Group (GONG) archive. These large-scale waves exhibit different features compared to each other. The observed wave of June 3rd has angular span of about $70^{\circ}$ with a diffuse wave front associated to NOAA active region 11496. It was found that the speed of the wave that started propagating at 17.53 UT is about $931\pm80$ km/s. The broadness nature of this Moreton wave can be interpreted as the vertical extension of the wave over the chromosphere. On the other hand, the wave of July 6th may be associated to X1.1 class flare that occurred at 23.01 UT around the 11515 active region. From the kinematical analysis, the wave propagated with the initial velocity of about $994\pm70$ km/s which is in agreement with coronal shock velocity derived from type II radio burst observation, $v\sim1100$ km/s. These two identified waves add the inventory of the large-scale waves observed in 24th Solar Cycle.

Sundial is the ancient or classic timekeeper device, especially prior to the invention of mechanical clock. In the classical Islamic civilization, the daily movement of the Sun becomes main indicator of praying time, which can be deduced using sundial. This kind of device probably permeated to Indonesia during the Islamic acculturation. Since then, the development of astronomical knowledge, technology, art and architectural in classical Indonesia are partially reflected into sundial. These historical attractions of sundial demand comprehensive documentation and investigation of Indonesian sundial which are rarely found in the current literatures. The required spatial and temporal information regarding Indonesian sundial can be collected by general public through citizen science scheme. This concept may answer scientific curiosity of a research and also educate the people, expose them with science. In this article, general scheme of citizen science are discussed, its application for sundial study in Indonesia is proposed as Suryakala-Nusantara program.

BM Orionis, eclipsing binary system that located in the center of Orion Nebula Cluster posses several enigmatic problems. Its intrinsic nature and nebular environment make it harder to measure the physical parameters of the system, but it is believed as Algol type binary where secondary component is pre-main sequence star with larger radius. To assure this, several stellar models ($M_1=5.9$ {M_{\odot}} and $M_2=2.0$ {M_{\odot}}) are created and simulated using MESA. Models with rigid rotation of $\omega=10^{-5}$ rad/s exhibit considerable similar properties during pre-main sequence stage, but 2.0 {\msun} at assumed age of $\sim10^6$ is $6.46$ times dimmer than observed secondary star. There must be an external mechanism to fill this luminosity gap. Then, simulated post-main sequence binary evolution of BM Ori that involves mass transfer shows that primary star will reach helium sequence with the mass of $\sim0.8$ {\msun} before second stage mass transfer.

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.

Variable stars, or particularly eclipsing binaries, are very essential astronomical occurrence. Surveys are the backbone of astronomy, and many discoveries of variable stars are the results of surveys. All-Sky Automated Survey (ASAS) is one of the observing projects whose ultimate goal is photometric monitoring of variable stars. Since its first light in 1997, ASAS has collected 50,099 variable stars, with 11,076 eclipsing binaries among them. In the present work we focus on the period determination of the eclipsing binaries. Since the number of data points in each ASAS eclipsing binary light curve is sparse, period determination of any system is a not straightforward process. For 30 samples of such systems we compare the implementation of Lomb-Scargle algorithm which is an Fast Fourier Transform (FFT) basis and Phase Dispersion Minimization (PDM) method which is non-FFT basis to determine their period. It is demonstrated that PDM gives better performance at handling eclipsing detached (ED) systems whose variability are non-sinusoidal. More over, using semi-automatic recipes, we get better period solution and satisfactorily improve 53% of the selected object's light curves, but failed against another 7% of selected objects. In addition, we also highlight 4 interesting objects for further investigation.