24小时热门版块排行榜    

查看: 3467  |  回复: 7
当前只显示满足指定条件的回帖,点击这里查看本话题的所有回帖

043114076

木虫 (正式写手)

[交流] 使用COMPASS 力场的同学请注意 已有5人参与

在使用compass力场计算N2和H2混合时发现一个很不正常的现象跟大家讨论一下.
单独的H2和N2的参数分别是:
        r0        epsilon/K
H     1.421  10.869   0.000
N     3.801  30.092   0.000
混合以后
H-N  3.388   1.885   0.000
可以看到混合后的epsilon 竟然比H-H之间还要小. 个人认为这是不合理的. 究其原因是由于compass力场的混合法则是six power, 当两个元素r0 相差比较大时, 就会出现混合后的epsilon 比二者都小很多. 所以只要是与H相关的混合epsilon的值就变得非常奇怪. 对于几何法则和LB法则就没有这个问题. 查看了compass力场的原始文献, 发现拟合都是针对纯的气体. 如果使用这样的参数模拟就会低估H2气体与吸附质之间的相互作用.
尝试跟力场作者联系, 但是没有找到联系方式.
原始文献: http://pubs.acs.org/doi/abs/10.1021/jp992913p
回复此楼
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

HITxzh

新虫 (初入文坛)

★
小木虫: 金币+0.5, 给个红包,谢谢回帖
请问一下 为什吗在lammps里面我使用不了class2的相关命令例如 bond_style class2  dihedral_style class2
angle_style class2    pair_style    lj/class2  使用这些命令总是报错:invalid bond style 或 invalid angle style 或invalid pair_style   
我是做有机物吸附的 势函数参数使用compass力场参数
5楼2013-07-11 21:35:18
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
查看全部 8 个回答

huanghl2010

木虫 (著名写手)

写个comment给作者。
MOFs
2楼2012-05-02 22:06:13
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

piaoma

木虫 (小有名气)

★
小木虫: 金币+0.5, 给个红包,谢谢回帖
MS 里的所谓通用力场,本身是模拟固体的。你拿来算气体或液体,当然不准了。另,貌似MS里的MD模拟,慎用慎用。我曾经出现过,MS里的水的RDF都出错。我觉得MS里castep模块做的不错,结果挺准
3楼2012-05-03 19:23:05
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

043114076

木虫 (正式写手)

引用回帖:
3楼: Originally posted by piaoma at 2012-05-03 19:23:05:
MS 里的所谓通用力场,本身是模拟固体的。你拿来算气体或液体,当然不准了。另,貌似MS里的MD模拟,慎用慎用。我曾经出现过,MS里的水的RDF都出错。我觉得MS里castep模块做的不错,结果挺准

呵呵, 没那么夸张. MD 方法主要的决定因素还是力场吧, 时刻保持警惕, 商业软件也不能大意.
castep 历史悠久了, vasp 跟它还有渊源.
4楼2012-05-03 21:49:40
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
普通表情 龙 兔 虎 猫 高级回复 (可上传附件)
信息提示
请填处理意见