Skip to content

call bounds in exfield to maintain PBC#177

Open
momokchung wants to merge 1 commit intomainfrom
dev
Open

call bounds in exfield to maintain PBC#177
momokchung wants to merge 1 commit intomainfrom
dev

Conversation

@momokchung
Copy link

Bounds should be called in exfield to maintain PBC in exfield.f.
Since external field electrostatic energy is a function of real space atomic coordinates and atoms can drift during MD simulation, any drifted molecules should be brought back to the real space prior to exfield routine.

(Calling bounds does not affect the force since force is independent of coordinates. However, if an atom is allowed to drift for too long, the energy might get too large/unwieldy.)

Sign up for free to join this conversation on GitHub. Already have an account? Sign in to comment

Labels

None yet

Projects

None yet

Development

Successfully merging this pull request may close these issues.

1 participant