  Felsenstein's classical model for Gaussian distributions on a phylogenetic
tree is shown to be a toric variety in the space of concentration matrices. We
present an exact semialgebraic characterization of this model, and we
demonstrate how the toric structure leads to exact methods for maximum
likelihood estimation. Our results also give new insights into the geometry of
ultrametric matrices.
