Several neoclassical models for the calculation of poloidal rotation in tokamaks were rederived within a common framework, extended to include additional physics and numerically compared. The importance of new physics phenomena not usually included in poloidal rotation calculations (e.g., poloidal electric field, V×B force resulting from enhanced radial particle flow arising from the ionization of recycling neutrals) was examined. Extensions of the Hirshman–Sigmar, Kim–Diamond–Groebner, and Stacey–Sigmar poloidal rotation models are presented.