If I have a list of potential energies of a pair particles r1 E1 r2 E2 r3 E3 how can I obtain the interaction for the particle pairs in all the space ? sometime this is refered to as a look-up table in MD simulations. Any body know somebody do this calculations and programs before ?? Many thanks.
l*y
3 楼
Read any text book for MD. Like "Computer simulation of liquids". I guess you're talking about how to fit a potential. If you know the physical meaning, use some analytical potential for fitting. For example LJ potential. If you have no idear about the analytical form, you can try Sum(An/r^n) if the potential has the boundary condition U->Inf as r->0 and U->0 as r->Inf. If neither of above method is good, just use the tabulated potential as it is. In this case maybe you need to do interpolation
【在 w****g 的大作中提到】 : If I have a list of potential energies of a pair particles : r1 E1 : r2 E2 : r3 E3 : how can I obtain the interaction for the particle pairs in all the space ? : sometime this is refered to as a look-up table in MD simulations. : Any body know somebody do this calculations and programs before ?? : Many thanks.
w*g
4 楼
Thanks. In my case it would be a tabulated potential as it is.Do you know any free code can do such work ? I remember somebody claimed that this will be faster than the normal function calls. Please let me know if you know somebody did such work and have such programs thanks
【在 l*******y 的大作中提到】 : Read any text book for MD. : Like "Computer simulation of liquids". : I guess you're talking about how to fit a potential. : If you know the physical meaning, use some analytical : potential for fitting. For example LJ potential. : If you have no idear about the analytical form, you : can try Sum(An/r^n) if the potential has the boundary : condition U->Inf as r->0 and U->0 as r->Inf. : If neither of above method is good, just use the : tabulated potential as it is. In this case maybe you
l*y
5 楼
You can use GSL, to fit your data in a cubic spline form. The algorithm is also in the book "the numerical recipes". I think many MD code support tabulated potentials, like DL_POLY, GROMACS etc.
will
【在 w****g 的大作中提到】 : Thanks. In my case it would be a tabulated potential as it is.Do you know : any free code can do such work ? I remember somebody claimed that this will : be faster than the normal function calls. Please let me know if you know : somebody did such work and have such programs : thanks