Proton oxygen processes are important in the field of proton therapy. In present work, the \(^{16}O(p,d)^{15}O\) , \(^{16}O(p,np)^{15}O\) reactions have been studied within the frame of Bertini and Liège’s intranuclear cascade models. The reactions have been simulated using the Geant4 program and the results were compared with experimental data from the EXFOR database. It has been found that both models predict a significantly higher probability for the \(^{16}O(p,np)^{15}O\) reaction channel compared with the \(^{16}O(p,np)^{15}O\) channel. It was shown that the theoritical models are consistent with the experimental data in energy range up to 80 MeV. The numerical results obtained within the Liège intranuclear cascade model reproduce the new data of T.Masuda quite accurately.
