with 6-31G basis for p
not working 6-31G !