implemented the use of lmfit package. Still, fail to reproduce the PNAS result for SDF