The "Event Generator Mode" has been introduced into the PHITS code in the treatment of low energy transport phenomena. In this mode, the evaluated nuclear data for neutrons and a special evaporation model are combined so as to trace all correlations of ejectiles keeping the energy and the momentum conservation in a collision. We have calculated the double differential cross section, the kerma coefficient and the displacement cross section by the event generator mode of PHITS, and compared them with the evaluated nuclear data for the validation of the mode in neutron-induced reactions under 30 MeV. As results, all calculated results by the event generator mode reproduce the evaluated nuclear data well below 20 MeV. As an application of this mode, we indicate a new approach to the deposit energy distribution for the thermal neutron transport and the averaged quality factor for neutrons.