diff --git a/src/USER-MISC/pair_coul_slater_cut.cpp b/src/USER-MISC/pair_coul_slater_cut.cpp index 7d7e727d0d..33363b869a 100644 --- a/src/USER-MISC/pair_coul_slater_cut.cpp +++ b/src/USER-MISC/pair_coul_slater_cut.cpp @@ -103,7 +103,7 @@ void PairCoulSlaterCut::compute(int eflag, int vflag) } if (eflag) ecoul = factor_coul * qqrd2e * - scale[itype][jtype] * qtmp*q[j] * rinv * + scale[itype][jtype] * qtmp*q[j] * rinv * (1 - (1 + r/lamda)*exp(-2*r/lamda)); if (evflag) ev_tally(i,j,nlocal,newton_pair,