Rate constants for the reaction O(P-3) + O-3((1)A(1)) --> 2O(2)(X(3) Sigma(g)(-)) are reported using variational transition state theory. A semi-empirical potential energy surface based on a double many-body expansion method due to Varandas and Pais is used to obtain transition state information. The computed rate constants are found to be in good to satisfactory agreement with available theoretical and experimental results. Preliminary studies of the reverse reaction, with one O-2 molecule in upsilon = 0 and the other in 27 less than or equal to upsilon less than or equal to 30 using a semiclassical wavepacket method reported recently show no reactivity on the potential energy surface used, contrary to some indirect evidence of reaction observed in recent experiments.