--- branches/development/src/nonbonded/MAW.cpp 2013/02/20 15:39:39 1850 +++ branches/development/src/nonbonded/MAW.cpp 2013/05/15 15:09:35 1874 @@ -174,7 +174,7 @@ namespace OpenMD { myDerivC = 0.0; } else if (idat.shiftedForce) { myPotC = D_e * (expfnc2C - 2.0 * expfncC); - myDerivC = 2.0 * D_e * beta * (expfnc2C - expfnc2C); + myDerivC = 2.0 * D_e * beta * (expfncC - expfnc2C); myPotC += myDerivC * ( *(idat.rij) - *(idat.rcut) ); } else { myPotC = 0.0;