Reaction dynamics of P(4S) + O2(X3Σ−g) → O(3P) + PO(X2Π) on a global CHIPR potential energy surface of PO2(X2A1): implications for atmospheric modelling

Abstract

Abstract. The reaction dynamics of P(4S) + O2(X3Σg-) → O(3P) + PO(X2Π) are thought to be important in atmospheric and interstellar chemistry. Based on the state-of-the-art ab initio energy points, we analytically constructed a global potential energy surface (PES) for the ground-state PO2(X2A1) using the combined-hyperbolic-inverse-power-representation (CHIPR) method. A total of 6471 energy points were computed by the multireference configuration interaction method with the Davidson correction and aug-cc-pV5Z basis set. The analytical CHIPR PES reproduces ab initio energies accurately with a root-mean-square deviation of 91.5 cm−1 (or 0.262 kcal mol−1). The strongly bound valence region of the PES has complicated topographical features with multiple potential wells and barriers. The attributes of the important intermediates are carefully validated with our geometry optimization results, as well as previous experimental and computational results. Finally, the reaction probability, integral cross sections, and rate constants for P(4S) + O2(X3Σg-) → O(3P) + PO(X2Π) are calculated using the quasi-classical trajectory and time-dependent wave packet methods. The trends of probability and integral cross section versus the collision energy can be divided into three stages, which are governed by the entrance barriers or exothermicity of the reaction. The rate constant demonstrates strong Arrhenius linear behaviour at relatively low temperatures but deviates from this pattern at high temperatures. The calculated cross sections and rate constants are helpful for modelling the phosphorus chemistry in atmospheric and interstellar media.