Updated default iteration paramter for Gd case
Added calculation of isolated 4f-spin system