I have been wrestling with how to properly deal with bonded (1-2 1-3 and 1-4) particle neighbor interactions, or whether they have to be dealt with, in the long range electrostatic solver codes (PME, SPME, etc). For instance, in the short range calculation, these interactions are typically excluded by omitting them from the neighbor list, or are simply scaled by the appropriate factor (0.0, 0.5, 0.833 etc) in the short range calculation. However, I am confused as to where this would happen, or how it could happen, in the long range solver. In the long range calculation, we, to my understanding, perform the following steps.interpolate charges on the mesh, and obtain total charge densityFFT and solve for the electrostatic fieldFFT back and obtain electrostatic potentialinterpolate forces back to particlesI am not seeing where in these steps we could exclude bonded neighbor interactions. The electric field is obtained via the total charge density which contains information from all particles, whether or not they are bonded. Is it even necessary to exclude these interactions here?
Proper handling of 1-2 1-3 1-4 bonded neighbors in long range electrostatic solver
molecular dynamics