Non-convergence of Geant4 hadronic models for 10 and 30 MeV protons in 18O and 14N

We investigated Geant4 Monte Carlo simulation of protons interactions in 18O and 14N. With proton incident energies of 10 and 30 MeV, we compared the different Geant4 models available for proton inelastic collisions: the low-energy parameterisation, Bertini cascade and precompound models. The models produced different (1) combinations of secondary particles and (2) final states of secondary particles. The non-convergence suggests that the hadronic models, yet to be benchmarked for this energy range, are not ready for production-mode problem-solving.