The particle mass used in cosmology N-body simulations is close to $10^{8}M_{\odot}$, which is about $10^{65}$ times larger than the GeV scale expected in particle physics. However, self-gravity interacting particle systems made up of different particles mass have different statistical and dynamical properties. Here we demonstrate that, due to this particle mass difference, the nowaday cosmology N-body simulations can have introduced an excessive core collapse process, especially for the low mass halos at high redshift. Such dynamical effect introduces an excessive cuspy center for these small halos, and it implies a possible connection to the so called "small scale crisis" for CDM models. Our results show that there exist a physical limit in cosmological simulations, by using about $10^3$ particles to describe smallest halos, and we provide a simple suggestion based on it to relieve those effects from the bias.

We introduce a new 2-D hexagon technique to probe the topological structure of the universe, in which we map regions of the sky with high and low galaxy densities onto a 2-D lattice of hexagon unit cells, We define filled cells as corresponding to high density regions and empty cells as corresponding to low density region. The number of filled cells and empty cells are kept same by controlling the size of the cells. By analyzing the six neighbors of each hexagon we can get and compare statistical topological properties of the high density and low density regions in the universe, in order to have a better understanding of the evolution of the universe. We apply this hexagon method on 2MASS data and discover significant topological differences between the high density and low density regions. Both regions have significant (>5 Sigma) topological shifts from the binomial distribution or the random distribution.

Galaxy clusters are the most massive objects in the Universe and comprise a high-temperature intracluster medium of about 10^7 K, believed to offer a main foreground effect for cosmic microwave background (CMB) data in the form of the thermal Sunyaev-Zel'dovich (SZ) effect. This assumption has been confirmed by SZ signal detection in hundreds of clusters but, in comparison with the huge numbers of clusters within optically selected samples from Sloan Digital Sky Survey (SDSS) data, this only accounts for a few per cent of clusters. Here we introduce a model-independent new method to confirm the assumption that most galaxy clusters can offer the thermal SZ signal as their main foreground effect. For the Wilkinson Microwave Anisotropy Probe (WMAP) seven-year data (and a given galaxy cluster sample), we introduced a parameter d1 as the nearest-neighbour cluster angular distance of each pixel, then we classified data pixels as 'to be' (d1--> 0 case) or 'not to be' (d1 large enough) affected by the sample clusters. By comparing the statistical results of these two kinds of pixels, we can see how the sample clusters affect the CMB data directly. We find that the Planck Early Sunyaev-Zel'dovich (ESZ) sample and X-ray samples (~10^2 clusters) can lead to obvious temperature depression in the WMAP seven-year data, which confirms the SZ effect prediction. However, each optically selected sample (>10^4 clusters) shows an opposite result: the mean temperature rises to about 10uK. This unexpected qualitative scenario implies that the main foreground effect of most clusters is $not$ always the expected SZ effect. This may be the reason why the SZ signal detection result is lower than expected from the model.

We present the Nearest Neighbor Distance (NND) analysis of SDSS DR5 galaxies. We give NND results for observed, mock and random sample, and discuss the differences. We find that the observed sample gives us a significantly stronger aggregation characteristic than the random samples. Moreover, we investigate the direction of NND and find that the direction has close relation with the size of the NND for the observed sample.

We have calculated the ground state electronic structure of He under pressure from 0 to 1500 GPa using both all-electron full-potential and pseudopotential methods based on the density functional theory (DFT). We find that throughout this pressure range, pseudopotentials yield essentially the same energy-volume curve for all of bcc, fcc, and hcp configurations as does the full-potential method, a strong indication that pseudopotential approximation works well for He both as the common element in some giant planets and as detrimental impurities in fusion reactor materials. The hcp lattice is always the most stable structure and bcc the least stable one. Since the energy preference of hcp over fcc and bcc is within 0.01 eV below 100 GPa and about 0.1 eV at 1500 GPa, on the same order of the error bar in local or semi-local density approximations in DFT, phase transitions can only be discussed with more precise description of electron correlation in Quantum Monte Carlo or DFT-based GW methods.

An important concern in the application of gamma-ray bursts (GRBs) to cosmology is that the calibration of GRB luminosity/energy relations depends on the cosmological model, due to the lack of a sufficient low-redshift GRB sample. In this paper, we present a new method to calibrate GRB relations in a cosmology-independent way. Since objects at the same redshift should have the same luminosity distance and since the distance moduli of Type Ia supernovae (SNe Ia) obtained directly from observations are completely cosmology independent, we obtain the distance modulus of a GRB at a given redshift by interpolating from the Hubble diagram of SNe Ia. Then we calibrate seven GRB relations without assuming a particular cosmological model and construct a GRB Hubble diagram to constrain cosmological parameters. From the 42 GRBs at $1.4<z\le6.6$, we obtain $\Omega_{\rm M}=0.25_{-0.05}^{+0.04}$, $\Omega_{\Lambda}=0.75_{-0.04}^{+0.05}$ for the flat $\Lambda$CDM model, and for the dark energy model with a constant equation of state $w_0=-1.05_{-0.40}^{+0.27}$, which is consistent with the concordance model in a 1-$\sigma$ confidence region.

Numerical simulation is an important tool to help us understand the process of structure formation in the universe. However many simulation results of cold dark matter (CDM) halos on small scale are inconsistent with observations: the central density profile is too cuspy and there are too many substructures. Here we point out that these two problems may be connected with a hitherto unrecognized bias in the simulation halos. Although CDM halos in nature and in simulation are both virialized systems of collisionless CDM particles, gravitational encounter cannot be neglected in the simulation halos because they contain much less particles. We demonstrate this by two numerical experiments, showing that there is a difference on the microcosmic scale between the natural and simulation halos. The simulation halo is more akin to globular clusters where gravitational encounter is known to lead to such drastic phenomena as core collapse. And such artificial core collapse process appears to link the two problems together in the bottom-up scenario of structure formation in the $\Lambda$CDM universe. The discovery of this bias also has implications on the applicability of the Jeans Theorem in Galactic Dynamics.

Numerical simulation plays an important role in the study of structure formation of the universe. However, the mass resolution in current simulations is still poor. Due to technical difficulties, it is necessary to use both greatly reduced number density of particles and greatly raised unit particle mass. Consequently, the particle masses used in cosmological simulations are about $10^{70}$ times larger than the $Gev$ candidates in particle physics. This is a huge physical bias that cannot be neglected when interpreting the results of the simulations. Here we discuss how such a bias affects Cold Dark Matter (CDM) cosmological simulations. We find that the small scale properties of the CDM particle system are changed in two aspects. 1) An upper limit is imposed on the spatial resolution of the simulation results. 2) Most importantly, an unexpected short mean free path is produced, and the corresponding two body scattering cross section is close to the value expected in the Self-Interaction Dark Matter(SIDM) model. Since the mean free path of real CDM particle systems is much longer than that in the simulations, our results imply 1) that there is probably no 'cusp problem' in real CDM halos, and 2) that a much longer time is needed to form new virialized halos in real CDM particle systems than in the simulations. This last result can help us understand the 'substructure problem'. Our discussion can also explain why the massive halos in the simulations may have smaller concentration coefficients.