diff --git a/src/GRANULAR/pair_gran_hertz_history.cpp b/src/GRANULAR/pair_gran_hertz_history.cpp
index 14cbc8eb0a20b7973c1342a6616e5c8519876990..e52aac10dbbb92ef15a8572115cd21756f35c9f9 100644
--- a/src/GRANULAR/pair_gran_hertz_history.cpp
+++ b/src/GRANULAR/pair_gran_hertz_history.cpp
@@ -183,7 +183,6 @@ void PairGranHertzHistory::compute(int eflag, int vflag)
         ccel = kn*(radsum-r)*rinv - damp;
         polyhertz = sqrt((radsum-r)*radi*radj / radsum);
         ccel *= polyhertz;
-        if (ccel < 0.0) ccel = 0.0;
 
         // relative velocities
 
@@ -391,7 +390,6 @@ double PairGranHertzHistory::single(int i, int j, int itype, int jtype,
   ccel = kn*(radsum-r)*rinv - damp;
   polyhertz = sqrt((radsum-r)*radi*radj / radsum);
   ccel *= polyhertz;
-  if (ccel < 0.0) ccel = 0.0;
 
   // relative velocities
 
diff --git a/src/GRANULAR/pair_gran_hooke.cpp b/src/GRANULAR/pair_gran_hooke.cpp
index 2d6d113ab7dbb8a7f707ab5ee11861f13e359040..66e1c3dd7a8bc09a2b8649e388fce6ff832b0cdd 100644
--- a/src/GRANULAR/pair_gran_hooke.cpp
+++ b/src/GRANULAR/pair_gran_hooke.cpp
@@ -161,7 +161,6 @@ void PairGranHooke::compute(int eflag, int vflag)
 
         damp = meff*gamman*vnnr*rsqinv;
         ccel = kn*(radsum-r)*rinv - damp;
-        if (ccel < 0.0) ccel = 0.0;
 
         // relative velocities
 
@@ -303,7 +302,6 @@ double PairGranHooke::single(int i, int j, int itype, int jtype, double rsq,
 
   damp = meff*gamman*vnnr*rsqinv;
   ccel = kn*(radsum-r)*rinv - damp;
-  if (ccel < 0.0) ccel = 0.0;
 
   // relative velocities
 
diff --git a/src/GRANULAR/pair_gran_hooke_history.cpp b/src/GRANULAR/pair_gran_hooke_history.cpp
index 1378709cb563406354c63fe30236091f3d1ce33f..e9662c9e7341ab77c6b76b69ac3ac773df501f4f 100644
--- a/src/GRANULAR/pair_gran_hooke_history.cpp
+++ b/src/GRANULAR/pair_gran_hooke_history.cpp
@@ -223,7 +223,6 @@ void PairGranHookeHistory::compute(int eflag, int vflag)
 
         damp = meff*gamman*vnnr*rsqinv;
         ccel = kn*(radsum-r)*rinv - damp;
-        if (ccel < 0.0) ccel = 0.0;
 
         // relative velocities
 
@@ -688,7 +687,6 @@ double PairGranHookeHistory::single(int i, int j, int itype, int jtype,
 
   damp = meff*gamman*vnnr*rsqinv;
   ccel = kn*(radsum-r)*rinv - damp;
-  if (ccel < 0.0) ccel = 0.0;
 
   // relative velocities