The diatomic systems, PO and PO− are studied, using numerical Hartree-Fock (NHF) and coupled-cluster calculations. The latter employs a hybrid NHF and Slater orbital basis set. Highly accurate CCSD methods predict bond lengths accurate to <0.004 A and frequencies to ∼60 cm−1. In addition the electron affinity of PO is computed to be 0.89 eV compared to an experimental value of 1.09±0.01. Comparisons are made with SCF and MBPT(2) results for PO+ using conventional basis sets.