Direct determination of pair natural orbitals