diff --git a/samples/gibbs_ensemble/create_fits.py b/samples/gibbs_ensemble/create_fits.py
index c51b0c7a674..2aa553a8f23 100644
--- a/samples/gibbs_ensemble/create_fits.py
+++ b/samples/gibbs_ensemble/create_fits.py
@@ -122,7 +122,7 @@ def rectilinear_law_fit(kT, p_c, A, kT_c):
 
 # Fit using scaling law
 fit_scaling, p_cov_scaling = curve_fit(
-    scaling_law_fit, kTs, y_scaling, p0=(1., 1.1), maxfev=6000)
+    scaling_law_fit, kTs, y_scaling, p0=(1., 1.), maxfev=6000)
 
 # Print fit values
 print("Fits using scaling law: B, T_c: {}".format(fit_scaling))