Multistate adiabatic and diabatic PESs were computed for the H+ + O2 collision system in Jacobi coordinates, (R,r,γ) using the cc-pVTZ basis set and the ic-MRCI level of theory. In addition, all possible interaction potentials and nonadiabatic coupling matrix elements among those different electronic states were also computed. Comparisons with earlier computed interaction potentials were made wherever possible, and the differences between them is attributed to the multistate diabatization and the chosen level of theory and basis set. Focusing our attention on the ground-state (GS) and the first excited-state (ES) PES, quantum dynamics were performed using the 2 × 2 diabatic potential submatrix obtained from the multistate (four) diabatic potential matrix within the VCC-RIOSA scheme at two experimentally reported collision energies, E(cm) = 9.5 and 23 eV. The scattering quantities were computed for two experimentally observed collision processes, namely, the inelastic vibrational excitation (IVE), H+ + O2 (X3Σg(−),v = 0) → H+ + O2 (X3Σg(−),v′), and the vibrational charge transfer (VCT), H+ + O2 (X3Σg(−),v = 0) → H (2S) + O (X2Πg,v′′). Comparisons were made with experimental results and found an overall improvement relative to the earlier computed results, and the discrepancies, if any, could be brought down to minimum by further modification in employed ab initio PESs and the interaction potential.