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

24小时热门版块排行榜    

CyRhmU.jpeg
小木虫论坛-学术科研互动平台 » 计算模拟区 » 分子模拟 » Monte Carlo » 怎么自己做力场-及气体力场的选择
71/1返回列表
查看: 1719  |  回复: 6
只看楼主@他人 存档 新回复提醒 (忽略) 收藏    在APP中查看
本帖产生 1 个 模拟EPI ，点击这里进行查看

zhang_jaj

木虫 (正式写手)

[求助] 怎么自己做力场-及气体力场的选择

大家好，我想知道自己做力场的过程，还有做起来难么？比如我想做的是H2，D2，T2的力场还有他们要混合的力场，LJ势的。

还有现在有个力场用在上面的气体中，那个比较好，或者什么力场适合于气体，非常感谢！
回复此楼

» 猜你喜欢

» 本主题相关价值贴推荐，对您同样有帮助:
1楼2012-07-08 22:58:59
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
回帖置顶 ( 共有1个 )

ChemiAndy

木虫 (正式写手)

【答案】应助回帖

★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★
感谢参与，应助指数 +1
zhang_jaj: 金币+15, ★★★★★最佳答案, 谢谢指教 2012-07-09 07:27:16
chaizhm: 金币+5, 谢谢~ 2012-07-09 09:19:18
zh1987hs: 模拟EPI+1, 谢谢！ 2012-07-09 09:36:15
zhang_jaj: 金币+5, ★★★★★最佳答案, 谢谢 2012-07-09 19:40:23
你这个H2， D2， T2的气体力场应该很好做，可以通过气相量子化学计算拟合全部参数。

因为只有一个键，键参数可以通过量子力学力常数计算直接得到。LJ势可以通过做H2-H2, H2-D2, H2-T2, D2-T2的分子间势能面拟合得到。单质气体还无需考虑静电参数。所以十分容易开发。

一般的开发力场步骤的讨论可以看这里：http://muchong.com/bbs/viewthread.php?tid=3744217

储氢曾经很热，应该很容易找到大量的现有H2力场参数。你看看它们怎么做的参数，照样做个D2和T2的即可。你这个模拟的难点不在势参数，而在于量子效应，需要特殊的蒙特卡洛方法，对每一步模拟的坐标不确定性进行计算。
赞一下(2人) 回复此楼

» 本帖已获得的红花（最新10朵）

zhang_jaj
2楼2012-07-09 04:18:25
 已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
普通回帖

zhang_jaj

木虫 (正式写手)
送鲜花一朵
引用回帖:
2楼: Originally posted by ChemiAndy at 2012-07-09 04:18:25
你这个H2， D2， T2的气体力场应该很好做，可以通过气相量子化学计算拟合全部参数。

因为只有一个键，键参数可以通过量子力学力常数计算直接得到。LJ势可以通过做H2-H2, H2-D2, H2-T2, D2-T2的分子间势能面拟合得 ...

非常感谢你的指导。希望后面碰到问题可以继续向您请教
赞一下 回复此楼
3楼2012-07-09 07:28:33
 已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

zhang_jaj

木虫 (正式写手)
引用回帖:
2楼: Originally posted by ChemiAndy at 2012-07-09 04:18:25
你这个H2， D2， T2的气体力场应该很好做，可以通过气相量子化学计算拟合全部参数。

因为只有一个键，键参数可以通过量子力学力常数计算直接得到。LJ势可以通过做H2-H2, H2-D2, H2-T2, D2-T2的分子间势能面拟合得 ...

看来再没人回答了，所以这些金币都是你的，谢谢！
赞一下 回复此楼
4楼2012-07-09 19:40:50
 已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

zhang_jaj

木虫 (正式写手)
引用回帖:
2楼: Originally posted by ChemiAndy at 2012-07-09 04:18:25
你这个H2， D2， T2的气体力场应该很好做，可以通过气相量子化学计算拟合全部参数。

因为只有一个键，键参数可以通过量子力学力常数计算直接得到。LJ势可以通过做H2-H2, H2-D2, H2-T2, D2-T2的分子间势能面拟合得 ...

你好，我想问一下。因为之前用vasp的时候，有人说可以把H的POTCAR里面的质量修改当做D来用，所以想在用MC我涉及到LJ选力场，就像把力场文件改质量来用，但是不知道到LJ的力场最后构建的原子间的对势有没有涉及到质量因素。也就是这样修改算同位素合理吗，能把同位素区分出来吗？

非常感谢！
赞一下 回复此楼
5楼2012-07-13 15:16:43
 已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

ChemiAndy

木虫 (正式写手)
★ ★ ★
zh1987hs: 金币+3, 谢谢 2012-07-14 09:10:15
力场不能这样简单改。因为通常Bond项的能量是这样的：
E=k(b-b0)^2
这个参数k就是振动力常数，由振动频率决定；很明显，H2和D2由于质量不同，振动力常数不同。

VASP中可以直接改，是因为它是要通过电子结构计算出来两原子间受力和振动频率的；而在力场中，这个东西被固定了，参数化了。
赞一下 回复此楼

» 本帖已获得的红花（最新10朵）

zhang_jaj
6楼2012-07-13 23:38:58
 已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

zhang_jaj

木虫 (正式写手)
送鲜花一朵
引用回帖:
6楼: Originally posted by ChemiAndy at 2012-07-13 23:38:58
力场不能这样简单改。因为通常Bond项的能量是这样的：
E=k(b-b0)^2
这个参数k就是振动力常数，由振动频率决定；很明显，H2和D2由于质量不同，振动力常数不同。

VASP中可以直接改，是因为它是要通过电子结构计算 ...

哦，谢谢指导啊。如果这样，在MD和MC中力场都不能这么改了。只能自己做，呵呵
赞一下 回复此楼
7楼2012-07-14 08:39:52
 已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 zhang_jaj 的主题更新
71/1返回列表
普通表情 龙 兔 虎 猫
信息提示
关闭
请填处理意见
关闭