Optimal pair density functional for the description of nuclei with large neutron excess