Re: [lammps-users] Tabulated potential and long range interaction
Meral Sharkas <eng_meral.sh10@...16...>
Wed, 30 Aug 2017 20:59:20 +0000 (UTC)
Thank you very much for your reply however I didn't include the Colombic range in my table but when I used the pair write command I found that my table alone(without hybrid overlay with long range potential) evaluates a potential similar to the potential added to colombic range using hybrid overlay command
So do you have any indication what may be the reason for that?
I created my tabulated potential for born Mayer potential and ZBL potential , then I used "pair_style hybrid/overlay table linear 10000 ewald coul/long 10" command to combine the long range potential. I got very high temperature and pressure system.
I troubled shoot the problem by once involving the colombic range and once removing it, and I used the "pair_write" command, I found that the potential curve is nearly the same in both cases as if the tabulated potential contains a long range..
I read on the mailing list some of the colleagues were talking about excluding the colombic part from the potential, however my potentials are both for short ranges
Vb=Ae^-r/rho (A,rho are constants)
Vzbl= coloumbic potential* screening function
Can someone help me what am I missing here, is it a problem of cut off distances, should i exclude some terms from the potential or another thing?
it is much safer to not include coulomb in the tabulation as there are many, many ways to make mistakes, but rather add it with hybrid/overlay to the tabulated potential during the simulation.