| 1049 |
|
// indirect reaction field terms. |
| 1050 |
|
|
| 1051 |
|
*(idat.vpair) += indirect_vpair; |
| 1052 |
+ |
(*(idat.excludedPot))[ELECTROSTATIC_FAMILY] += epot; |
| 1053 |
|
(*(idat.pot))[ELECTROSTATIC_FAMILY] += indirect_Pot; |
| 1054 |
|
*(idat.f1) += indirect_dVdr; |
| 1055 |
|
|
| 1078 |
|
chg1 += *(sdat.flucQ); |
| 1079 |
|
// dVdFQ is really a force, so this is negative the derivative |
| 1080 |
|
*(sdat.dVdFQ) -= *(sdat.flucQ) * data.hardness + data.electronegativity; |
| 1081 |
+ |
cerr << "dVdFQ harmonic part = " << *(sdat.dVdFQ) << "\n"; |
| 1082 |
|
} |
| 1083 |
|
|
| 1084 |
|
if (summationMethod_ == esm_REACTION_FIELD) { |