  We study stellar models for Betelgeuse using the HR diagram and surface
abundances as observational constraints. Previous studies on Betelgeuse have
not systematically investigated the surface abundances, but we believe they can
be impacted by, and thus be used as an observational constraint for various
parameters such as initial mass, rotation, and overshoot scheme. We investigate
stellar models with varying initial mass as they evolve past the main sequence,
and we examine the red supergiant (RSG) properties in detail. For each mass, we
vary the initial rotation up to $\sim$ 300 $km\:s^{-1}$, and test two different
overshoot parameters. Overall, the acceptable initial mass range is 12 to 25
$M_\odot$, but for non-rotating models only, the range is decreased to 15 to 24
$M_\odot$. Also for rotating models, we find that $v/v_{\rm K} = 0.3$ is the
upper limit for initial rotation, as more rapidly rotating models are unable to
fit to Betelgeuse's surface abundances as an RSG. In addition, we report two
possibilities for the current stage of evolution, core helium burning or core
carbon burning and beyond. We find that certain 17 $M_\odot$ models could fit
to both stages. Finally, we discuss the implications of our results in the
context of merger scenarios which have been suggested as a mechanism to attain
the observed surface velocity of Betelgeuse.
