An analytical potential energy surface (PES) representation of the O+(S-4)+H-2(X (1)Sigma(g)(+)) system was developed by fitting around 600 CCSD(T)/cc-pVQZ ab initio points. Rate constant calculations for this reaction and its isotopic variants (D-2 and HD) were performed using the quasiclassical trajectory (QCT) method, obtaining a good agreement with experimental data. Calculations conducted to determine the cross section of the title reaction, considering collision energies (E-T) below 0.3 eV, also led to good accord with experiments. This PES appears to be suitable for kinetics and dynamics studies. Moreover, the QCT results show that, although the hypotheses of a widely used capture model are not satisfied, the resulting expression for the cross section can be applied within a suitable E-T interval, due to errors cancellation. This could be a general situation regarding the application of this simple model to ion-molecule processes. (C) 2004 American Institute of Physics.