We describe a rapid, accurate method for calculating rovibrational distributions in diatomic products from
elementary chemical reactions. The basis of the model is momentum interconversion at a critical configuration
defined in terms of molecular dimensions of the species involved. This approach shares common elements
with recent models of inelastic processes and the kinematic reactive model of Elsum and Gordon. We point
out that these and related approaches represent a development of Newtonian mechanics equivalent to that
followed in the conventional formulation of classical mechanics, but one in which motive force for change
at the molecular level is attributed to dp/dt rather than to dV/dq. This leads to a particularly transparent form
of mechanics that uses only familiar data such as bond length, mass, spectroscopic constants, and velocity,
yet may be applied to the highly resolved single collision experiments of molecular reaction dynamics. We
describe key aspects of the computational method, e.g., the definition of the critical configuration, the disposal
of reaction enthalpy, the manner of assigning product vibrational states, and the way in which conservation
of energy is ensured. Examples are chosen to illustrate the range of reactions to which the method may be
applied. Each would represent a challenge to conventional theory. We show that velocity-angular momentum
diagrams may be used to interpret data and to give physical insight into the origins of observed rotational
distributions. Good agreement is obtained between experimental and calculated (v,j) distributions for a wide
range of elementary reactions suggesting that our model, despite its simplicity, captures the principal physics
of chemical change at the molecular level