Text this: Reaction dynamics of P(<sup>4</sup>S) + O<sub>2</sub>(X<sup>3</sup>Σ<sup>−</sup><sub>g</sub>)  →  O(<sup>3</sup>P) + PO(X<sup>2</sup>Π) on a global CHIPR potential energy surface of PO<sub>2</sub>(X<sup>2</sup>A<sub>1</sub>): implications for atmospheric modelling