The National Ignition Facility has recently achieved successful burning plasma and ignition using the inertial confinement fusion (ICF) approach. However, there are still many fundamental physics phenomena that are not well understood, including the kinetic processes in the hohlraum. Shan et al. [Phys. Rev. Lett. 120, 195001 (2018)0031-900710.1103/PhysRevLett.120.195001] utilized the energy spectra of neutrons to investigate the kinetic colliding plasma in a hohlraum of indirect drive ICF. However, due to the typical large spatial-temporal scales, this experiment could not be well simulated by using available codes at that time. Utilizing our advanced high-order implicit PIC code, LAPINS, we were able to successfully reproduce the experiment on a large scale of both spatial and temporal dimensions, in which the original computational scale was increased by approximately seven to eight orders of magnitude. Not only is the validity of the explanation of the experiment confirmed by our simulations, i.e., the abnormally large width of neutron spectra comes from beam-target nuclear fusions, but also a different physical insight into the source of energetic deuterium ions is provided. The acceleration of deuterium ions can be categorized into two components: one is propelled by a sheath electric field created by the charge separation at the onset, while the other is a result of the reflection of the potential of the shock wave. The robustness of the acceleration mechanism is analyzed with varying initial conditions, e.g., temperatures, drifting velocity, and ion components. This paper might serve as a reference for benchmark simulations of upcoming simulation codes and may be relevant for future research on mixtures and entropy increments at plasma interfaces.