Lattice Study of QCD Phase Structure by Canonical Approach - Towards determining the phase transition line

We compare higher moments of baryon numbers measured at the RHIC heavy ion collision experiments with those by the lattice QCD calculations. We employ the canonical approach, in which we can access the real chemical potential regions avoiding the sign problem. In the lattice QCD simulations, we study several fits of the number density in the pure imaginary chemical potential, and analyze how these fits affects behaviors at the real chemical potential. In the energy regions between $\sqrt{s}_{NN}$=19.6 and 200 GeV, the susceptibility calculated at $T/T_c=0.93$ is consistent with experimental data at $0 \le \mu_B/T < 1.5$, while the kurtosis shows similar behavior with that of the experimental data in the small $\mu_B/T$ regions $0 \le \mu_B/T < 0.3$. The experimental data at $\sqrt{s}_{NN}=$ 11.5 shows quite different behavior. The lattice result in the deconfinement region,$T/T_c=1.35$, is far from experimental data.