I am calculating ligand-protein type interaction and noted that the set of bases available in Dmol for the GGA / LDA theory level overestimates the energy interaction among residues containing sulfur atoms and the ligand. Even for very large distances the interaction energy doesn't go to zero.