diff --git a/src/ext/symmfunc.f90 b/src/ext/symmfunc.f90
index dc50b31..700e9db 100644
--- a/src/ext/symmfunc.f90
+++ b/src/ext/symmfunc.f90
@@ -638,11 +638,11 @@ subroutine sf_F1_ijk(vecRij, vecRik, Rij, Rik, cost, lambda, &
 
     Fijk  = arg**zeta
 
-    dFijk_dRj = prefactor*( -cost*( vecRij/Rij ) &
-                            + vecRik/Rij )
+    dFijk_dRj = prefactor*( -cost*( vecRij/Rij/Rij ) &
+                            + vecRik/Rij/Rik )
 
-    dFijk_dRk = prefactor*( -cost*( vecRik/Rik ) &
-                            + vecRij/Rik )
+    dFijk_dRk = prefactor*( -cost*( vecRik/Rik/Rik ) &
+                            + vecRij/Rik/Rij )
 
   end subroutine sf_F1_ijk