24小时热门版块排行榜    

查看: 153  |  回复: 1
当前主题已经存档。

txf8378

金虫 (小有名气)

[交流] 【求助】关于DL-POLY的Table

看了n遍说明书,可能是天生愚钝吧,愣是没看懂Table文件是怎么回事,哪位好心的大侠帮忙解释一下这些参数的意义啊
1DL_POLY TEST CASE 1: Silica Potentials (Vessal) in internal units
  1.5212170E-03  7.6000000E+00      5000
Si4+    O_2-     
  9.5140002E+06  9.4702042E+06  9.4266099E+06  9.3832162E+06
  9.3400223E+06  9.2970273E+06  9.2542301E+06  9.2116300E+06
  9.1692260E+06  9.1270171E+06  9.0850026E+06  9.0431815E+06

[ Last edited by zdhlover on 2009-12-13 at 23:03 ]
回复此楼

» 猜你喜欢

已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

老虎大王

木虫 (著名写手)

★ ★ ★ ★ ★
lei0736(金币+5,VIP+0):谢谢 希望能赏脸开个MD 专题
我来回答你:
A. 第一行是title,也就是Manual上说的file header。
B。第二行有三个量,分别是mesh resolution(单位是ang),cutoff,以及table中的grid points 的数目。在你给的文件中,相当于dr=1.521217e-3,rc=7.6,ngrid=5000
C。第三行写的是原子类别。你的文件里分别是si4+和O2-,这个名字是可以自己起的,但要跟其它输入文件(如FIELD)一致,顺序也要跟FIELD里的一致。
D. 从第四行开始就是势函数的值了。每行四个数,一共是ngrid个数,则共有Int((ngrid+3)/4)行。
E然后是力,也就是势函数在每个点上对r的负导数,同上。

对于其它的原子间势函数,重复以上的C.D.E。

[ Last edited by 老虎大王 on 2008-12-25 at 13:06 ]
2楼2008-12-25 13:05:22
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 txf8378 的主题更新
普通表情 高级回复 (可上传附件)
信息提示
请填处理意见