diff --git a/system.jl b/system.jl
index 7b6a4c7..1c860ef 100644
--- a/system.jl
+++ b/system.jl
@@ -97,7 +97,7 @@ function solve_system(s::system; reinitialize_densities=true, monitor_print=true
         monitor_print && println("Total binding E per nucleon = $E_per_A")
 
         # check convergence
-        abs(previous_E_per_A - E_per_A) < 0.1 && break
+        abs(previous_E_per_A - E_per_A) < 0.0001 && break
         previous_E_per_A = E_per_A
     end
 end