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.

We investigate the dependence of the extended X-ray emission from the halos of optically luminous early-type galaxies on the small-scale (the nearest neighbor distance) and large-scale (the average density inside the 20 nearest galaxies) environments. We cross-match the 3rd Data Release of the Second XMMNewton Serendipitous Source Catalog (2XMMi-DR3) to a volume-limited sample of the Sloan Digital Sky Survey (SDSS) Data Release 7 with Mr < −19.5 and 0.020 < z < 0.085, and find 20 early-type galaxies that have extended X-ray detections. The X-ray luminosity of the galaxies is found to have a tighter correlation with the optical and near infrared luminosities when the galaxy is situated in the low large-scale density region than in the high large-scale density region. Furthermore, the X-ray to optical (r-band) luminosity ratio, LX/Lr, shows a clear correlation with the distance to the nearest neighbor and with large-scale density environment only where the galaxies in pair interact hydrodynamically with seperations of rp < rvir. These findings indicate that the galaxies in the high local density region have other mechanisms that are responsible for their halo X-ray luminosities than the current presence of a close encounter, or alternatively, in the high local density region the cooling time of the heated gas halo is longer than the typical time between the subsequent encounters.

We have analyzed H and Ks -band images of the Arches cluster obtained using the NIRC2 instrument on Keck with the laser guide star adaptive optics (LGS AO) system. With the help of the LGS AO system, we were able to obtain the deepest ever photometry for this cluster and its neighborhood, and derive the background-subtracted present-day mass function (PDMF) down to 1.3M⊙ for the 5"-9" annulus of the cluster. We find that the previously reported turnover at 6M⊙ is simply due to a local bump in the mass function (MF), and that the MF continues to increase down to our 50 % completeness limit (1.3M⊙) with a power-law exponent of Γ= -0.91 for the mass range of 1.3 < M/M⊙ < 50. Our numerical calculations for the evolution of the Arches cluster show that the Γ values for our annulus increase by 0.1-0.2 during the lifetime of the cluster, and thus suggest that the Arches cluster initially had Γ of -1.0~-1.1 which is only slightly shallower than the Salpeter value.

The Fokker-Planck (FP) model is one of the commonly used methods for studies of the dynamical evolution of dense spherical stellar systems such as globular clusters and galactic nuclei. The FP model is numerically stable in most cases, but we find that it encounters numerical difficulties rather often when the effects of tidal shocks are included in two-dimensional (energy and angular momentum space) version of the FP model or when the initial condition is extreme (e.g., a very large cluster mass and a small cluster radius). To avoid such a problem, we have developed a new integration scheme for a two-dimensional FP equation by adopting an Alternating Direction Implicit (ADI) method given in the Douglas-Rachford split form. We find that our ADI method reduces the computing time by a factor of ~2 compared to the fully implicit method, and resolves problems of numerical instability.

We have developed a control electronics system for an infrared detector array of KASINICS (KASI Near Infrared Camera System), which is a new ground-based instrument of the Korea Astronomy and Space science Institute (KASI). Equipped with a 512×512 InSb array (ALADDIN III Quadrant, manufactured by Raytheon) sensitive from 1 to 5μm, KASINICS will be used at J, H, Ks, and L-bands. The controller consists of DSP(Digital Signal Processor), Bias, Clock, and Video boards which are installed on a single VME-bus backplane. TMS320C6713DSP, FPGA(Field Programmable Gate Array), and 384-MB SDRAM(Synchronous Dynamic Random Access Memory) are included in the DSP board. DSP board manages entire electronics system, generates digital clock patterns and communicates with a PC using USB 2.0 interface. The clock patterns are downloaded from a PC and stored on the FPGA. UART is used for the communication with peripherals. Video board has 4 channel ADC which converts video signal into 16-bit digital numbers. Two video boards are installed on the controller for ALADDIN array. The Bias board provides 16 dc bias voltages and the Clock board has 15 clock channels. We have also coded a DSP firmware and a test version of control software in C-language. The controller is flexible enough to operate a wide range of IR array and CCD. Operational tests of the controller have been successfully finished using a test ROIC (Read-Out Integrated Circuit).

We examine the corecollapse times of isolated, two-mass-component star clusters using Fokker-Planck models. With initial condition of Plummer models, we find that the corecollapse times of clusters with M1/M2 >> 1 are well correlated with (N1/N2)^0.5 (m1/m2)^2 Trh, where (M1/M2) and (m1/m2) are the light to heavy component total and individual mass ratios, respectively, N1/N2 is the number ratio, and Trh is the initial half-mass relaxation time scale. We also find two-component cluster parameters that best match multi-component (thus more realistic) clusters with power-law mass functions.

Two-component models (normal star and degenerate star components) are the simplest realization of clusters with a mass spectrum because the high mass stars quickly evolve off leaving degenerate stars behind, while low mass stars survive for a long time as main-sequence stars. In the present study we examine the post-collapse evolution of globular clusters using two-component Fokker-Planck models that include three-body binary heating. We confirm that a simple parameter є ≡ (Etot/trh)/(Ec/trc) well describes the occurrence of gravothermal oscillations of two-component clusters. Also, we find that the degree of instability depends on the steepness of the mass function such that clusters with a steeper mass function are less exposed to instability.