Hi everyone,
I've been working on implementing external electric fields into AMBER the last few weeks.
Using a force for Fx, Fy, Fz to be the charge of the atom times the charge of the electric field in the x,y,z direction.
I was wondering if anyone knew whether or not there should be a virial component in this.
It seems from here
https://mailman-1.sys.kth.se/pipermail/gromacs.org_gmx-users/2004-July/011473.html that external forces should not be included, but I think since the forces effect the internal movements of the system, the virial should still be added (or would be subtracted because its external?) similar to how it is done in the nonbonded code (pairs_calc.i).
So virial vxx would be something like force generated by electric field on x * (x coodinate of molecule - x origin of periodic cell).
Any thoughts?
Thanks,
Charlie Lin
Graduate Student
Walker Lab
San Diego Supercomputer Center
_______________________________________________
AMBER-Developers mailing list
AMBER-Developers.ambermd.org
http://lists.ambermd.org/mailman/listinfo/amber-developers
Received on Fri Feb 19 2016 - 15:30:02 PST
