Abstract: The Monte Carlo method is usually used to simulate the feature of the external radiotherapy of proton beam to optimize the clinical therapy scenario. The selection of a suitable physical model is very critical to ensure the correction and the efficiency of the simulation. This work studied nine Geant4 physical models, which are often used in the simulation of radiotherapy of proton beam with energy of 60-250 MeV. The dosage features and the microcosmic secondary yields in different materials were compared. The result shows that, although the standard EM, the low energy Penelope model and the low energy Livermore model can output the right dose distribution, they can’t simulate the baryon yields. The parametrization driven model LHEP and the Quark gluon string pre-compound model both can’t simulate the complex ion transportation. Thus the above five models are not enough for the research calculation in the radiation protection and the radiation damage. The QGSP_BIC_EMY model in Geant4 hadron therapy example is very bad for underestimating the dose ratio of Bragg peak to entrance and absolute dose. The QGSP_BERT model, QGSP_BIC model and FTFP_BERT model are suitable for the radiotherapy proton beam with energy of 60-250 MeV.