During optimization PCM (IEFPCM) solvent calculations of a molecule whose geometry was previously optimized in the gas phase, the error on total polarizations charges repeatedly came up in the output file. The calculations are being performed at the B3LYP/6-31G(d) level of theory, with the software Gaussian 16. Can someone give suggestions on how to overcome this error?