setSat_p loosen convergence criterion
add helium reference