应《网络安全法》要求，自2017年10月1日起，未进行实名认证将不得使用互联网跟帖服务。为保障您的帐号能够正常使用，请尽快对帐号进行手机号验证，感谢您的理解与支持！

24小时热门版块排行榜    

小木虫论坛-学术科研互动平台 » 计算模拟区 » 分子模拟 » 其他 » 使用COMPASS 力场的同学请注意
81/1返回列表
查看: 2809  |  回复: 7
只看楼主 @他人 存档 新回复提醒 (忽略) 收藏    在APP中查看

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
回复此楼

» 猜你喜欢

» 本主题相关价值贴推荐，对您同样有帮助:
1楼 2012-05-02 20:46:47
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

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的回帖

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的回帖

043114076

木虫 (正式写手)
引用回帖:
5楼: Originally posted by HITxzh at 2013-07-11 21:35:18
请问一下 为什吗在lammps里面我使用不了class2的相关命令例如 bond_style class2  dihedral_style class2
angle_style class2    pair_style    lj/class2  使用这些命令总是报错：invalid bond style 或 invalid  ...

具体lammps的东西我没用过，不懂，还是向lammps用户求教吧。好好看看手册也许自己也能找到答案。
赞一下 回复此楼
6楼2013-07-12 09:25:04
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

pigrass

专家顾问 (著名写手)

伊卡斯特的小猪
★
小木虫: 金币+0.5, 给个红包，谢谢回帖
虽然对各力场的细节不怎么了解，但是曾经一个晶体用不同力场优化，晶胞参数明显不同
赞一下 回复此楼
#有问题发论坛，不接受私聊#爱晶体，爱MaterialsStudio#MS专辑点击右边#http://emuch.net/taotie/41673.html#物性数据#http://webbook.nist.gov/chemistry/#
7楼2013-07-15 10:18:56
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

ChemiAndy

木虫 (正式写手)
★
小木虫: 金币+0.5, 给个红包，谢谢回帖
分子力场参数的状态可迁移性历来就是个大问题。纯态下拟合获取的参数迁移到混合态体系必须做验证。这些验证至少要包括混合体系的密度，蒸发焓。一般相对实验误差不超过10%才能算是勉强可接受的。
赞一下 回复此楼
8楼2013-07-17 05:54:44
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 043114076 的主题更新
81/1返回列表
普通表情 龙 兔 虎 猫 高级回复 (可上传附件)
信息提示
关闭
请填处理意见
关闭