Gallium oxide ([Formula: see text]) attracts considerable technological interest because of its high Baliga’s figure-of-merit and high breakdown voltages. As the models for the breakdown behavior of n-doped [Formula: see text] that consider soft (barrier lowering) and hard (avalanche effect) breakdowns are still lacking, in this study, we model the breakdown operations in <001> oriented Schottky barrier diodes considering both the soft- and hard-breakdown phenomena. The completion of the impact ionization model of [Formula: see text]-[Formula: see text] in <001> orientation is proposed by determining the hole impact ionization coefficient, thereby reproducing hard breakdown operations. Moreover, a barrier lowering model is determined for reproducing soft breakdown operations. The outcomes of the proposed modeling investigation are expected to be crucial for predicting the reverse-biased operations of [Formula: see text]-[Formula: see text] in <001> orientation to facilitate further technological development and applications of [Formula: see text].