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.

Using a cosmological CDM simulation, we analyze the differences between the widely-used spin pa- rameters suggested by Peebles and Bullock. The dimensionless spin parameter λ proposed by Peebles is theoretically well-justified but includes an annoying term, the potential energy, which cannot be directly obtained from observations and is computationally expensive to calculate in numerical simulations. The Bullock’s spin parameter λ′ avoids this problem assuming the isothermal density profile of a virialized halo in the Newtonian potential model. However, we find that there exists a substantial discrepancy between λ and λ′ depending on the adopted potential model (Newtonian or Plummer) to calculate the halo total energy and that their redshift evolutions differ to each other significantly. Therefore, we introduce a new spin parameter, λ′′, which is simply designed to roughly recover the value of λ but to use the same halo quantities as used in λ′. If the Plummer potential is adopted, the λ′′ is related to the Bullock’s definition as λ′′ = 0.80 × (1 + z)−1/12λ′. Hence, the new spin parameter λ′′ distribution becomes consistent with a log-normal distribution frequently seen for the λ′ while its mean value is much closer to that of λ. On the other hand, in case of the Newtonian potential model, we obtain the relation of λ′′ = (1 + z)−1/8λ′; there is no significant difference at z = 0 as found by others but λ′ becomes more overestimated than λ or λ′′ at higher redshifts. We also investigate the dependence of halo spin parameters on halo mass and redshift. We clearly show that although the λ′ for small-mass halos with Mh < 2× 1012M⊙ seems redshift independent after z = 1, all the spin parameters explored, on the whole, show a stronger correlation with the increasing halo mass at higher redshifts.

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.