We report a six-dimensional CCSD(T)/aug-cc-pVTZ dipole moment surface for the electronic ground state of PH3 computed ab initio on a large grid of 10 080 molecular geometries. Parameterized, analytical functions are fitted through the ab initio data, and the resulting dipole moment functions are used, together with a potential energy function determined by refining an existing ab initio surface in fittings to experimental wavenumber data, for simulating absorption spectra of the first three polyads of PH3, i.e., (ν2, ν4), (ν1, ν3, 2ν2, 2ν4, ν2 + ν4), and (ν1 + ν2, ν3 + ν2, ν1 + ν4, ν3 + ν4, 2ν2 + ν4, ν2 + 2ν4, 3ν2, 3ν4). The resulting theoretical transition moments show excellent agreement with experiment. A line-by-line comparison of the simulated intensities of the ν2/ν4 band system with 955 experimental intensity values reported by Brown et al. [L.R. Brown, R.L. Sams, I. Kleiner, C. Cottaz, L. Sagui, J. Mol. Spectrosc. 215 (2002) 178–203] gives an average absolute percentage deviation of 8.7% (and a root-mean-square deviation of 0.94 cm−1 for the transition wavenumbers). This is very remarkable since the calculations rely entirely on ab initio dipole moment surfaces and do not involve any adjustment of these surfaces to reproduce the experimental intensities. Finally, we predict the line strengths for transitions between so-called cluster levels (near-degenerate levels formed at high rotational excitation) for J up to 60.