Contracted Gaussian-type function sets are proposed for polarization functi
ons of the atoms Ga-Kr and In-Xe. We also report polarization functions for
Ca and Sr. A segmented contraction scheme is used for its compactness and
computational efficiency. The contraction coefficients and orbital exponent
s are fully optimized to minimize the deviation from accurate atomic natura
l orbitals. The present polarization functions yield more than 99% of atomi
c correlation energies predicted by accurate natural orbitals of the same s
ize.