ABSTRACTWe have developed a novel theoretical method to obtain independent fission product yields and prompt neutron observables by connecting mass and total kinetic energy (TKE) distributions calculated by a four-dimensional Langevin dynamical model to a Hauser-Feshbach statistical decay model. In the Langevin calculations, mass distributions corresponding to the standard I and II modes were obtained separately and superposed to reproduce the fission fragment yield of spontaneous fission of ${ ^{238,240,242}}{\rm{Pu}}$ 238,240,242Puand thermal neutron-induced fission of ${ ^{239}}{\rm{Pu}}$ 239Pu. This was achieved by using different neck parameters for these two modes in the two-center shell model shape parametrization, and a systematics of the superposing ratio was obtained as a function of $\left({N - Z} \right)/A$N−Z/Aof the fissioning nuclei. The Hauser-Feshbach calculations were performed using a nuclear reaction code TALYS for ${ ^{239}}{\rm{Pu}}\left({{\rm{n}},{\rm{f}}} \right)$ 239Pun,freaction in the incident energy range from thermal up to $5\ {\rm{MeV}}$5 MeV, and the calculated prompt fission observables were compared with experimental and evaluated data. Although further improvements are needed for the most important nuclides, it turned out that the present methodology has the capability to prepare fission-related nuclear data for nuclides for which measurements are difficult.