Current and future neutrino oscillation experiments utilize information of hadronic final states to improve sensitivities on oscillation parameters measurements. Among the physics of hadronic systems in neutrino interactions, the hadronization model controls multiplicities and kinematics of final state hadrons from the primary interaction vertex. For relatively high invariant mass events, neutrino interaction generators rely on the PYTHIA6 hadronization program. Here, we show a possible improvement of this process in neutrino event generators, by utilizing expertise from the HERMES experiment. Next, we discuss the possibility to implement the PYTHIA8 program in neutrino interaction generators, including GENIE and NEUT. Finally, we show preliminary comparisons of PYTHIA8 predictions with neutrino hadron multiplicity data from bubble chamber experiments within the GENIE hadronization validation tool.