State-to-state inelastic rate coefficients of phosphine in collision with He at low to moderate temperature

Several phosphorus-bearing molecules, such as the phosphine of interest here, have been detected in astrophysical media. With the aim of satisfying the precision required by the astrophysical community, we report the rate coefficients of PH3 in collision with helium from low to moderate temperature. To this end, we constructed the first three-dimensional potential energy surface (3D-PES) of the PH3–He van der Waals complex, which governs the nuclear motions. The 3D-PES was worked out by means of the standard coupled cluster with single, double and non-iterative triple excitation approach, in conjunction with the aug-cc-pVQZ basis set and complemented by mid-bond functions. This 3D-PES presents a well of 34.92 cm−1 at {R, θ, Φ}  = {5.76 a0, 90°, 60°}. Afterwards, we incorporated this 3D-PES into time-independent close-coupling quantum dynamical computations to derive the inelastic cross-sections of rotational excitation of (ortho-) para-PH3 after collision with He up to (1000) 500 cm−1. Subsequently, we evaluated the rate coefficients for temperatures up to (100 K) 50 K populating the (41) 42 low-lying rotational levels of (ortho-) para-PH3. These data were derived by averaging the cross-sections thermally over the Maxwell–Boltzmann velocity distribution. No general propensity rules are found. We also performed a comparison with the rates for NH3–He. Differences are observed that invalidate the use of NH3 rates for deducing accurate abundances of phosphine in cold astrophysical media. Our results should be of great help in determining accurate PH3 abundances and, more generally, constraining the interstellar PH3 chemistry better.