The structure and phase stability of a series of lithiated titania polymorphs were determined using energy minimizations of the periodic bulk crystal structures and both density functional theory (DFT) and a potential shell model. The DFT calculations were performed spin unrestricted following the linear combination of atomic orbital approach with the B3LYP exchange-correlation potential. For the potential shell model, a new set of force field parameters was derived, independently of the DFT calculations, to describe the lithium−lattice interactions. The eight polymorphs considered in this study are the rutile, anatase, brookite, TiO2−B, ramsdellite, hollandite, spinel, and hexagonal structures. The lithium to titanium ratio, x, of each lithiated titania polymorph was varied from 0.0 to 1.0 with 0.25 increments. The potential model predictions were found to be in good agreement with the structure and energetics of the lithiated titania polymorphs determined from DFT calculations, at all lithium contents. Both computational approaches indicate the following relationships. The naturally occurring titania polymorphs (i.e., rutile, anatase, brookite, and TiO2−B) were found to be the most stable of the eight phases in the absence of lithium. Anatase, brookite, and ramsdellite become energetically favored over rutile upon lithium insertion. The hexagonal and spinel polymorphs have stabilities approaching that of rutile with increasing lithium content and showed essentially equivalent total energies for x = 1.0. This prediction is consistent with numerous experimental studies, which have reported that rutile electrodes, in particular those that consist of nanomaterials, can undergo phase transformations to the hexagonal or spinel structure upon lithium insertion. The calculations indicate that the main factors controlling the relative stability of the lithiated titania polymorphs are the lithium bonding environment, the arrangement of LiO
x
and TiO6 polyhedra, and the extent of lattice deformation upon lithiation.