diff --git a/neurolib/models/fhn/timeIntegration.py b/neurolib/models/fhn/timeIntegration.py
index 1a4bcbe8..f3be5874 100644
--- a/neurolib/models/fhn/timeIntegration.py
+++ b/neurolib/models/fhn/timeIntegration.py
@@ -57,7 +57,6 @@ def timeIntegration(params):
         # no self-feedback delay
         Dmat[np.eye(len(Dmat)) == 1] = np.zeros(len(Dmat))
     Dmat_ndt = np.around(Dmat / dt).astype(int)  # delay matrix in multiples of dt
-    params["Dmat_ndt"] = Dmat_ndt
 
     # Additive or diffusive coupling scheme
     coupling = params["coupling"]
diff --git a/neurolib/models/hopf/timeIntegration.py b/neurolib/models/hopf/timeIntegration.py
index 15c49980..20e4fe60 100644
--- a/neurolib/models/hopf/timeIntegration.py
+++ b/neurolib/models/hopf/timeIntegration.py
@@ -62,7 +62,6 @@ def timeIntegration(params):
         Dmat = dp.computeDelayMatrix(lengthMat, signalV)
         Dmat[np.eye(len(Dmat)) == 1] = np.zeros(len(Dmat))
     Dmat_ndt = np.around(Dmat / dt).astype(int)  # delay matrix in multiples of dt
-    params["Dmat_ndt"] = Dmat_ndt
     # ------------------------------------------------------------------------
 
     # Initialization
diff --git a/neurolib/models/wc/timeIntegration.py b/neurolib/models/wc/timeIntegration.py
index 5dcc830b..0c83d55f 100644
--- a/neurolib/models/wc/timeIntegration.py
+++ b/neurolib/models/wc/timeIntegration.py
@@ -62,7 +62,6 @@ def timeIntegration(params):
         Dmat = dp.computeDelayMatrix(lengthMat, signalV)
         Dmat[np.eye(len(Dmat)) == 1] = np.zeros(len(Dmat))
     Dmat_ndt = np.around(Dmat / dt).astype(int)  # delay matrix in multiples of dt
-    params["Dmat_ndt"] = Dmat_ndt
     # ------------------------------------------------------------------------
     # Initialization
     # Floating point issue in np.arange() workaraound: use integers in np.arange()
diff --git a/neurolib/models/ww/timeIntegration.py b/neurolib/models/ww/timeIntegration.py
index f60027c3..5142e965 100644
--- a/neurolib/models/ww/timeIntegration.py
+++ b/neurolib/models/ww/timeIntegration.py
@@ -67,7 +67,6 @@ def timeIntegration(params):
         # no self-feedback delay
         Dmat[np.eye(len(Dmat)) == 1] = np.zeros(len(Dmat))
     Dmat_ndt = np.around(Dmat / dt).astype(int)  # delay matrix in multiples of dt
-    params["Dmat_ndt"] = Dmat_ndt
 
     # # Additive or diffusive coupling scheme
     # version = params["version"]