diff --git a/samples/electrophoresis.py b/samples/electrophoresis.py
index e69391a8fa3..d82ac8831ca 100644
--- a/samples/electrophoresis.py
+++ b/samples/electrophoresis.py
@@ -32,6 +32,8 @@
 
 logging.basicConfig(level=logging.INFO)
 
+# Use a fixed int for a deterministic behavior
+np.random.seed()
 
 required_features = ["P3M", "EXTERNAL_FORCES", "WCA"]
 espressomd.assert_features(required_features)
diff --git a/testsuite/scripts/samples/test_electrophoresis.py b/testsuite/scripts/samples/test_electrophoresis.py
index f1532a2b871..5302680b288 100644
--- a/testsuite/scripts/samples/test_electrophoresis.py
+++ b/testsuite/scripts/samples/test_electrophoresis.py
@@ -19,8 +19,28 @@
 import importlib_wrapper
 import numpy as np
 
+
+def set_seed(code):
+    np_seed = "np.random.seed()"
+    assert np_seed in code
+    return code.replace(np_seed, "np.random.seed(42)")
+
+
+def set_p3m_params(code):
+    p3m = "electrostatics.P3M(prefactor=1.0, accuracy=1e-2)"
+    assert p3m in code
+    return code.replace(
+        p3m, "electrostatics.P3M(prefactor=1, mesh=[16, 16, 16], cao=1, accuracy=1e-2, r_cut=3.7, alpha=0.2, tune=False)")
+
+
+def make_deterministic(code):
+    code = set_seed(code)
+    code = set_p3m_params(code)
+    return code
+
+
 sample, skipIfMissingFeatures = importlib_wrapper.configure_and_import(
-    "@SAMPLES_DIR@/electrophoresis.py", N_SAMPLES=400)
+    "@SAMPLES_DIR@/electrophoresis.py", N_SAMPLES=400, substitutions=make_deterministic, random_seeds=False)
 
 
 @skipIfMissingFeatures