Monte Carlo simulations were used to calculate the yields for the primary species (e^–_aq, H^•, H₂, ^•OH and H₂O₂) formed from the radiolysis of neutral liquid water by mono-energetic 2 MeV neutrons at temperatures between 25–350°C. The 2 MeV neutron was taken as representative of a fast neutron flux in a reactor. For light water, the moderation of these neutrons generated elastically scattered recoil protons of ∼1.264, 0.465, 0.171 and 0.063 MeV, which at 25°C, had linear energy transfers (LETs) of ∼22, 43, 69 and 76 keV/μm, respectively. Neglecting the radiation effects due to oxygen ion recoils and assuming that the most significant contribution to the radiolysis came from these first four recoil protons, the fast neutron yields could be estimated as the sum of the yields for these protons after allowance was made for the appropriate weightings according to their energy. Yields were calculated at 10⁻⁷, 10⁻⁶ and 10⁻⁵ s after the ionization event at all temperatures, in accordance with the time range associated with the scavenging capacities generally used for fast neutron radiolysis experiments. The results of the simulations agreed reasonably well with the experimental data, taking into account the relatively large uncertainties in the experimental measurements, the relatively small number of reported radiolysis yields, and the simplifications included in the model. Compared with data obtained for low-LET radiation (⁶⁰Co γ rays or fast electrons), our computed yields for fast neutron radiation showed essentially similar temperature dependences over the range of temperatures studied, but with lower values for yields of free radicals and higher values for molecular yields. This general trend is a reflection of the high-LET character of fast neutrons. Although the results of the simulations were consistent with the experiment, more experimental data are required to better describe the dependence of radiolytic yields on temperature and to test more thoroughly our modeling calculations.