Many simulations have been performed to elucidate the formation process of first stars. In first star formation, radiative feedback is a key process in determining stellar masses. However, previous simulations which follow the feedback process don’t resolve the small scale ( 10 AU) to realize long-term calculation, and the structure near massive protostars is still unknown. To clarify how the radiation from the protostar works, we need to resolve small scale and calculate the interaction between the radiation and the dense gas in such a region. As a first step towards understanding the phenomenon in this region, we perform the high-resolution simulation around the massive protostar without radiative transfer. We find that dense gas covers the protostar even in the polar direction and the HII region cannot expand. Solving the radiative transfer for getting accurate results is our future work. We are currently developing the new radiation hydrodynamics code for that.