I am optimizing a system containing cations in the presence of iron porphyrin in Gaussian. Using the recommended basis set 6-311g+(d) is turning out to be very computationally expensive. Can anyone please suggest me how to specify different basis sets for different elements. Also, is it advisable to optimize the geometry at a primitive basis set and then re-optimize using that geometry and a higher basis set? Any help would be appreciated.