diff --git a/jaxquantum/quantum/solvers.py b/jaxquantum/quantum/solvers.py
index 560cd4e..53bf514 100644
--- a/jaxquantum/quantum/solvers.py
+++ b/jaxquantum/quantum/solvers.py
@@ -103,7 +103,7 @@ def f(
     term = ODETerm(f)
     solver = Dopri5()
     saveat = SaveAt(ts=t_list)
-    stepsize_controller = PIDController(rtol=1e-6, atol=1e-6)
+    stepsize_controller = PIDController(rtol=1e-6, atol=1e-8)
 
     sol = diffeqsolve(
         term,