35 |
|
* |
36 |
|
* [1] Meineke, et al., J. Comp. Chem. 26, 252-271 (2005). |
37 |
|
* [2] Fennell & Gezelter, J. Chem. Phys. 124, 234104 (2006). |
38 |
< |
* [3] Sun, Lin & Gezelter, J. Chem. Phys. 128, 24107 (2008). |
38 |
> |
* [3] Sun, Lin & Gezelter, J. Chem. Phys. 128, 234107 (2008). |
39 |
|
* [4] Kuang & Gezelter, J. Chem. Phys. 133, 164101 (2010). |
40 |
|
* [5] Vardeman, Stocker & Gezelter, J. Chem. Theory Comput. 7, 834 (2011). |
41 |
|
*/ |
174 |
|
myDerivC = 0.0; |
175 |
|
} else if (idat.shiftedForce) { |
176 |
|
myPotC = D_e * (expfnc2C - 2.0 * expfncC); |
177 |
< |
myDerivC = 2.0 * D_e * beta * (expfnc2C - expfnc2C); |
177 |
> |
myDerivC = 2.0 * D_e * beta * (expfncC - expfnc2C); |
178 |
|
myPotC += myDerivC * ( *(idat.rij) - *(idat.rcut) ); |
179 |
|
} else { |
180 |
|
myPotC = 0.0; |