• Journal of The Korean Astronomical Society
• /
• v.48 no.1
• /
• pp.75-82
• /
• 2015

We use a volume-limited sample of quasars in the Sloan Digital Sky Survey (SDSS) DR7 quasar catalog to identify quasar groups and address their statistical significance. This quasar sample has a uniform selection function on the sky and nearly a maximum possible contiguous volume that can be drawn from the DR7 catalog. Quasar groups are identified by using the Friend-of-Friend algorithm with a set of fixed comoving linking lengths. We find that the richness distribution of the richest 100 quasar groups or the size distribution of the largest 100 groups are statistically equivalent with those of randomly-distributed points with the same number density and sky coverage when groups are identified with the linking length of $70h^{-1}Mpc$. It is shown that the large-scale structures like the huge Large Quasar Group (U1.27) reported by Clowes et al. (2013) can be found with high probability even if quasars have no physical clustering, and does not challenge the initially homogeneous cosmological models. Our results are statistically more reliable than those of Nadathur (2013), where the test was made only for the largest quasar group. It is shown that the linking length should be smaller than $50h^{-1}Mpc$ in order for the quasar groups identified in the DR7 catalog not to be dominated by associations of quasars grouped by chance. We present 20 richest quasar groups identified with the linking length of $70h^{-1}Mpc$ for further analyses.

• Journal of The Korean Astronomical Society
• /
• v.47 no.3
• /
• pp.87-98
• /
• 2014

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 initial 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 assigned 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.

• Journal of The Korean Astronomical Society
• /
• v.36 no.1
• /
• pp.1-12
• /
• 2003

Nonthermal particles can be produced due to incomplete thermalization at collisionless shocks and further accelerated to very high energies via diffusive shock acceleration. In a previous study we explored the cosmic ray (CR) acceleration at cosmic shocks through numerical simulations of CR modified, quasi-parallel shocks in 1D plane-parallel geometry with the physical parameters relevant for the shocks emerging in the large scale structure formation of the universe (Kang & Jones 2002). Specifically we considered pancake shocks driven by accretion flows with $U_o = 1500 km\;s^{-l}$ and the preshock gas temperature of $T_o = 10^4 - 10^8K$. In order to consider the CR acceleration at shocks with a broader range of physical properties, in this contribution we present additional simulations with accretion flows with $U_o = 75 - 1500 km\;s^{-l}$ and $T_o = 10^4K$. We also compare the new simulation results with those reported in the previous study. For a given Mach number, shocks with higher speeds accelerate CRs faster with a greater number of particles, since the acceleration time scale is $t_{acc}\;{\propto}\;U_o^{-2}$. However, two shocks with a same Mach number but with different shock speeds evolve qualitatively similarly when the results are presented in terms of diffusion length and time scales. Therefore, the time asymptotic value for the fraction of shock kinetic energy transferred to CRs is mainly controlled by shock Mach number rather than shock speed. Although the CR acceleration efficiency depends weakly on a well-constrained injection parameter, $\epsilon$, and on shock speed for low shock Mach numbers, the dependence disappears for high shock Mach numbers. We present the 'CR energy ratio', ${\phi}(M_s)$, for a wide range of shock parameters and for $\epsilon$ = 0.2 - 0.3 at terminal time of our simulations. We suggest that these values can be considered as time-asymptotic values for the CR acceleration efficiency, since the time-dependent evolution of CR modified shocks has become approximately self-similar before the terminal time.

• Journal of The Korean Astronomical Society
• /
• v.35 no.3
• /
• pp.111-121
• /
• 2002

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 $\ge$ 70 (as defined by the APM team) is found to be $r_0 = 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 $\ge$ 0 and 1 are measured to be $r_0 = 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 $r_0= 0.40d_c + 3.2$ where $d_c$ 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 $L_x {\ge}0.65{\times} 10^{44}\;h^{-2}erg\;s^{-1}$ is $30.3_{-6.5}^{+8.2}\;h^{-1}$ Mpc, and that of the BCS clusters with $L_x {\ge}0.70{\times} 10^{44}\;h^{-2}erg\;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.

• Journal of The Korean Astronomical Society
• /
• v.36 no.3
• /
• pp.89-95
• /
• 2003

Observations of dark matter dominated dwarf and low surface brightness disk galaxies favor density profiles with a flat-density core, while cold dark matter (CDM) N-body simulations form halos with central cusps, instead. This apparent discrepancy has motivated a re-examination of the microscopic nature of the dark matter in order to explain the observed halo profiles, including the suggestion that CDM has a non-gravitational self-interaction. We study the formation and evolution of self-interacting dark matter (SIDM) halos. We find analytical, fully cosmological similarity solutions for their dynamics, which take proper account of the collisional interaction of SIDM particles, based on a fluid approximation derived from the Boltzmann equation. The SIDM particles scatter each other elastically, which results in an effective thermal conductivity that heats the halo core and flattens its density profile. These similarity solutions are relevant to galactic and cluster halo formation in the CDM model. We assume that the local density maximum which serves as the progenitor of the halo has an initial mass profile ${\delta}M / M {\propto} M^{-{\epsilon}$, as in the familiar secondary infall model. If $\epsilon$ = 1/6, SIDM halos will evolve self-similarly, with a cold, supersonic infall which is terminated by a strong accretion shock. Different solutions arise for different values of the dimensionless collisionality parameter, $Q {\equiv}{\sigma}p_br_s$, where $\sigma$ is the SIDM particle scattering cross section per unit mass, $p_b$ is the cosmic mean density, and $r_s$ is the shock radius. For all these solutions, a flat-density, isothermal core is present which grows in size as a fixed fraction of $r_s$. We find two different regimes for these solutions: 1) for $Q < Q_{th}({\simeq} 7.35{\times} 10^{-4}$), the core density decreases and core size increases as Q increases; 2) for $Q > Q_{th}$, the core density increases and core size decreases as Q increases. Our similarity solutions are in good agreement with previous results of N-body simulation of SIDM halos, which correspond to the low-Q regime, for which SIDM halo profiles match the observed galactic rotation curves if $Q {\~} [8.4 {\times}10^{-4} - 4.9 {\times} 10^{-2}]Q_{th}$, or ${\sigma}{\~} [0.56 - 5.6] cm^2g{-1}$. These similarity solutions also show that, as $Q {\to}{\infty}$, the central density acquires a singular profile, in agreement with some earlier simulation results which approximated the effects of SIDM collisionality by considering an ordinary fluid without conductivity, i.e. the limit of mean free path ${\lambda}_{mfp}{\to} 0$. The intermediate regime where $Q {\~} [18.6 - 231]Q_{th}$ or ${\sigma}{\~} [1.2{\times}10^4 - 2.7{\times}10^4] cm^2g{-1}$, for which we find flat-density cores comparable to those of the low-Q solutions preferred to make SIDM halos match halo observations, has not previously been identified. Further study of this regime is warranted.

• Journal of The Korean Astronomical Society
• /
• v.47 no.2
• /
• pp.77-86
• /
• 2014

Using a cosmological ${\Lambda}CDM$ simulation, we analyze the differences between the widely-used spin parameters suggested by Peebles and Bullock. The dimensionless spin parameter ${\lambda}$ 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 ${\lambda}^{\prime}$ 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 ${\lambda}$ and ${\lambda}^{\prime}$ 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, ${\lambda}^{\prime\prime}$, which is simply designed to roughly recover the value of ${\lambda}$ but to use the same halo quantities as used in ${\lambda}^{\prime}$. If the Plummer potential is adopted, the ${\lambda}^{\prime\prime}$ is related to the Bullock's definition as ${\lambda}^{\prime\prime}=0.80{\times}(1+z)^{-1/12}{\lambda}^{\prime}$. Hence, the new spin parameter ${\lambda}^{\prime\prime}$ distribution becomes consistent with a log-normal distribution frequently seen for the ${\lambda}^{\prime}$ while its mean value is much closer to that of ${\lambda}$. On the other hand, in case of the Newtonian potential model, we obtain the relation of ${\lambda}^{\prime\prime}=(1+z)^{-1/8}{\lambda}^{\prime}$; there is no significant difference at z = 0 as found by others but ${\lambda}^{\prime}$ becomes more overestimated than ${\lambda}$ or ${\lambda}^{\prime\prime}$ at higher redshifts. We also investigate the dependence of halo spin parameters on halo mass and redshift. We clearly show that although the ${\lambda}^{\prime}$ for small-mass halos with $M_h$ < $2{\times}10^{12}M_{\odot}$ 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.