The purpose of this study is to calculate atomic (H(0)) density distribution in JAEA 10 ampere negative ion source. A collisional radiative model is developed for the calculation of the H(0) density distribution. The non-equilibrium feature of the electron energy distribution function (EEDF), which mainly determines the H(0) production rate, is included by substituting the EEDF calculated from 3D electron transport analysis. In this paper, the H(0) production rate, the ionization rate, and the density distribution in the source chamber are calculated. In the region where high energy electrons exist, the H(0) production and the ionization are enhanced. The calculated H(0) density distribution without the effect of the H(0) transport is relatively small in the upper region. In the next step, the effect should be taken into account to obtain more realistic H(0) distribution.