The Korean Astronomical Society (KAS) Education & Public Outreach Committee has provided education services for children and school teachers in Cambodia over the past three years from 2016 to 2018. In the first year, 2016, one KAS member visited Pusat to teach astronomy to about 50 children, and in the following two years of 2017 and 2018, three and six KAS members, respectively, executed education workshops for ∼ 20 (per each year) local school teachers in Sisophon. It turned out that it is desirable to include both teaching of astronomical knowledge and making experiments and observations in the education in order for the program to be more effective. Language barrier was the main obstacle in conveying concepts and knowledge, and having a good interpreter was very important. It happens that some languages, such as the Khmer of Cambodia, do not have astronomical terminologies, so that lecturers and even the education participants together are needed to communicate and create appropriate words. Handout hardcopies of the education materials (presentation files, lecture/experiment summaries, terminologies, etc.) are extremely helpful for the participants. Actual performing of assembling and using astronomical telescopes for night sky observations has been lifetime experience for some of the participants, which might promote zeal for knowledge and education. It is hoped that these education services for developing countries like Cambodia can be regularly continued in the future, and further extended to other countries such as Laos and Myanmar.

The National Foundation Day of Korea (개천절, 開天節) is currently celebrated on October 3 in Gregorian calendar. We review the history of dating the National Foundation Day of Korea and make a suggestion that it be celebrated on October 3 in the lunar calendar. We present numerous historical records on heaven-worship rites supporting the date October 3 in the lunar calendar. It is pointed out that October 3 in the solar calendar has been adopted in 1949 by the National Assembly with the thought that the lunar calendar is inferior and behind the times. The thought originates from misunderstanding on the value of the lunar calendar and from the ignorance of importance of history and tradition. Since there are now many national holidays that follow the lunar calendar, the logic of the National Assembly in 1949 also makes no sense. We emphasize that the lunar calendar should be followed for the National Foundation Day of Korea for its historical and symbolic characteristics restoration. We also investigate the year of the foundation of the first country of Korea, Dangun Joseon. It is found that even though the majority of the literature before late 15th century recorded the beginning year of Dangun Joseon dynasty to be equal to that of Liao Dynasty (堯), it was accidentally changed to the 25th year of Liao Dynasty in 1484 through a misinterpretation of the previous records. We claim that the beginning year of Dangun Joseon should be set to that of Liao Dynasty as recorded in the original literature in the earlier days. According to the two main opinions accepted by Korea, the beginning year of Liao Dynasty was 2357 B.C. or 2333 B.C., which correspond to the year of Gap-Jin (the 41st year of the sexagenary cycle) or Mu-Jin (the 4th year of the sexagenary cycle), respectively.

We develop a parallel cosmological hydrodynamic simulation code designed for the study of formation and evolution of cosmological structures. The gravitational force is calculated using the TreePM method and the hydrodynamics is implemented based on the smoothed particle hydrodynamics. The ini- tial displacement and velocity of simulation particles are calculated according to second-order Lagrangian perturbation theory using the power spectra of dark matter and baryonic matter. The initial background temperature is given by Recfast and the temperature uctuations at the initial particle position are as- signed according to the adiabatic model. We use a time-limiter scheme over the individual time steps to capture shock-fronts and to ease the time-step tension between the shock and preshock particles. We also include the astrophysical gas processes of radiative heating/cooling, star formation, metal enrichment, and supernova feedback. We test the code in several standard cases such as one-dimensional Riemann problems, Kelvin-Helmholtz, and Sedov blast wave instability. Star formation on the galactic disk is investigated to check whether the Schmidt-Kennicutt relation is properly recovered. We also study global star formation history at different simulation resolutions and compare them with observations.

A novel method to characterize the topology of the early-universe intergalactic medium during the epoch of cosmic reionization is presented. The 21-cm radiation background from high redshift is analyzed through calculation of the 2-dimensional (2D) genus. The radiative transfer of hydrogen- ionizing photons and ionization-rate equations are calculated in a suite of numerical simulations under various input parameters. The 2D genus is calculated from the mock 21-cm images of high-redshift universe. We construct the 2D genus curve by varying the threshold differential brightness temperature, and compare this to the 2D genus curve of the underlying density field. We find that (1) the 2D genus curve reflects the evolutionary track of cosmic reionization and (2) the 2D genus curve can discriminate between certain reionization scenarios and thus indirectly probe the properties of radiation-sources. Choosing the right beam shape of a radio antenna is found crucial for this analysis. Square Kilometre Array (SKA) is found to be a suitable apparatus for this analysis in terms of sensitivity, even though some deterioration of the data for this purpose is unavoidable under the planned size of the antenna core.

We present the relation between the genus in cosmology and the Betti numbers for excursion sets of three- and two-dimensional smooth Gaussian random fields, and numerically investigate the Betti numbers as a function of threshold level. Betti numbers are topological invariants of figures that can be used to distinguish topological spaces. In the case of the excursion sets of a three-dimensional field there are three possibly non-zero Betti numbers; β0 is the number of connected regions, β1 is the number of circular holes (i.e., complement of solid tori), and β2 is the number of three-dimensional voids (i.e., complement of three-dimensional excursion regions). Their sum with alternating signs is the genus of the surface of excursion regions. It is found that each Betti number has a dominant contribution to the genus in a specific threshold range. β0 dominates the high-threshold part of the genus curve measuring the abundance of high density regions (clusters). β1 dominates the genus near the median thresholds which measures the topology of negatively curved iso-density surfaces, and β2 corresponds to the low-threshold part measuring the void abundance. We average the Betti number curves (the Betti numbers as a function of the threshold level) over many realizations of Gaussian fields and find that both the amplitude and shape of the Betti number curves depend on the slope of the power spectrum n in such a way that their shape becomes broader and their amplitude drops less steeply than the genus as n decreases. This behaviour contrasts with the fact that the shape of the genus curve is fixed for all Gaussian fields regardless of the power spectrum. Even though the Gaussian Betti number curves should be calculated for each given power spectrum, we propose to use the Betti numbers for better specification of the topology of large scale structures in the universe.

We construct several Milky Way-like galaxy models containing a gas halo (as well as gaseous and stellar disks, a dark matter halo, and a stellar bulge) following either an isothermal or an NFW density profile with varying mass and initial spin. In addition, galactic winds associated with star formation are tested in some of the simulations. We evolve these isolated galaxy models using the GADGET-3 N-body/hydrodynamic simulation code, paying particular attention to the effects of the gaseous halo on the evolution. We find that the evolution of the models is strongly affected by the adopted gas halo component, particularly in the gas dissipation and the star formation activity in the disk. The model without a gas halo shows an increasing star formation rate (SFR) at the beginning of the simulation for some hundreds of millions of years and then a continuously decreasing rate to the end of the run at 3 Gyr. Whereas the SFRs in the models with a gas halo, depending on the density profile and the total mass of the gas halo, emerge to be either relatively flat throughout the simulations or increasing until the middle of the run (over a gigayear) and then decreasing to the end. The models with the more centrally concentrated NFW gas halo show overall higher SFRs than those with the isothermal gas halo of the equal mass. The gas accretion from the halo onto the disk also occurs more in the models with the NFW gas halo, however, this is shown to take place mostly in the inner part of the disk and not to contribute significantly to the star formation unless the gas halo has very high density at the central part. The rotation of a gas halo is found to make SFR lower in the model. The SFRs in the runs including galactic winds are found to be lower than those in the same runs but without winds. We conclude that the effects of a hot gaseous halo on the evolution of galaxies are generally too significant to be simply ignored. We also expect that more hydrodynamical processes in galaxies could be understood through numerical simulations employing both gas disk and gas halo components.

We have measured the correlation functions of the optically selected clusters of galaxies in the Abell and the APM catalogs, and of the X-ray clusters in the X-ray-Brightest Abell-type Clusters of galaxies (XBACs) catalog and the Brightest Clusters Sample (BCS). The same analysis method and the same method of characterizing the resulting correlation functions are applied to all observational samples. We have found that the amplitude of the correlation function of the APM clusters is much higher than what has been previously claimed, in particular for richer subsamples. The correlation length of the APM clusters with the richness R ≥ 70 (as defined by the APM team) is found to be r0 = 25.4 -3.0 +3.1 h-1 Mpc. The amplitude of correlation function is about 2.4 times higher than that of Croft et al. (1997). The correlation lengths of the Abell clusters with the richness class RC ≥ 0 and 1 are measured to be r0 = 17.4 -1.1 +1.2 and 21.0 -2.8 +2.8 h-1 Mpc, respectively, which is consistent with our results for the APM sample at the similar level of richness. The richness dependence of cluster correlations is found to be r0= 0.40dc + 3.2 where dc is the mean intercluster separation. This is identical in slope with the Bahcall & West (1992)'s estimate, but is inconsistent with the weak dependence of Croft et al. (1997). The X-ray bright Abell clusters in the XBACs catalog and the X-ray selected clusters in the BCS catalog show strong clustering. The correlation length of the XBACs clusters with Lx ≥ 0.65×10 44 h-2 erg s-1 is 30.3 -6.5 +8.2 h-1 Mpc, and that of the BCS clusters with Lx ≥ 0.70×10 44 h-2erg s-1 is 30.2 -8.9 +9.8 h-1 Mpc. The clustering strength of the X-ray clusters is much weaker than what is expected from the optical clusters.

We have made a topological study of cosmic microwave background (CMB) polarization maps by simulating the AMiBA experiment results. A ACDM CMB sky is adopted to make mock interferometric observations designed for the AMiBA experiment. CMB polarization fields are reconstructed from the AMiBA mock visibility data using the maximum entropy method. We have also considered effects of Galactic foregrounds on the CMB polarization fields. The genus statistic is calculated from the simulated Q and U polarization maps, where Q and U are Stokes parameters. Our study shows that the Galactic foreground emission, even at low Galactic latitude, is expected to have small effects on the CMB polarization field. Increasing survey area and integration time is essential to detect non-Gaussian signals of cosmological origin through genus measurement.

We have developed a cosmological N-body code which can simulate unprecedently large number of massive particles. This code is based on the Particle-Mesh scheme, and utilize the recent fast I/O devices to swap all variables. Using the new code we have simulated the formation and evolution of structures at high redshifts in the standard Cold Dark Matter (CDM) cosmogony. A simulation evolving 1024^3 particles on a 2048^3 mesh with the initial standard CDM power spectrum is being made. This is the first billion particle cosmological simulation with initial conditions representing the theoretical model over the widest range of space. A smaller, but still very large CDM simulation with 512^3 particles on a 1024^3 mesh has been completed. We have found that the galaxy-scale CDM halos with diameters of tens of kpcs undergo complete collapse before redshift 10. Our results clearly indicate that galactic and subgalactic structures have formed far before redshift 5 which is the present upper limit to the epoch of observed structures. We emphasize that the non-linear evolution of the galactic and subgalactic-scale structures starts as early as z ~ 50, and that cosmological simulations must start at such high redshifts. A high mass resolution is also indispensable to accurately represent the theoretical model in the initial conditions down to subgalactic scales, and to correctly study the subsequent formation and evolution of structures through hierarchical clustering.

We have analyzed the content of the Korean stone star chart. Ch'on-Sang-Yul-Cha-Bun-Ya-Ji-Do(here-after Ch'on-Sang-Do). In the star map we have found 1468 stars, 4 more than the Chinese star catalog Bo-Chun-Ga. The four extra stars form a constellation, Jong Dae Boo. The map projection law used in the star chart is found to be the polar equtorial and equidistance projection. The linear distance of an object on Ch'on-Sang-Do from the center is linearly proportional to the north polar angular distance. We have found from a statistical analysis that most stars with declination lower than 50 are at positions representing the epoch of around the first century. On the other hand, stars near the north pole with declination higher than 50 are at the epoch of about 1300, which is close to the time the chart was engraved. This implies that the original Ko-Gu-Rye Dynasty's star chart has been revised by astronomers of Cho-Sun Dynasty. We have also shown that stars on Ch'on-Sang-Do are engraved in such a way that their area is linearly proportional to the visual magnitude.