  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.
