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

24小时热门版块排行榜    

小木虫论坛-学术科研互动平台 » 计算模拟区 » 分子模拟 » Monte Carlo » 怎么自己做力场-及气体力场的选择
北京石油化工学院2026年研究生招生接收调剂公告
71/1返回列表
查看: 1898  |  回复: 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返回列表
普通表情 龙 兔 虎 猫
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 求调剂 +21 111623 2026-04-04 23/1150 2026-04-07 02:16 by BruceLiu320
[考研] 0703化学 +20 妮妮ninicgb 2026-04-04 23/1150 2026-04-06 23:06 by chenzhimin
[考研] 287求调剂 +3 通信学硕081000 2026-04-03 4/200 2026-04-06 21:03 by going home
[考研] 307求调剂 +3 所念及所望 2026-04-06 3/150 2026-04-06 17:30 by 土木硕士招生
[考研] 一志愿安徽某211 0703化学总分339求调剂 +7 晚风不晚 2026-04-04 7/350 2026-04-06 14:06 by houyaoxu
[考研] 专硕0854初试考材科基，求调剂 +6 3220548044 2026-04-06 9/450 2026-04-06 10:26 by barlinike
[考研] 0817化学工程与技术求调剂，一志愿中海洋319 +14 lv945 2026-04-04 14/700 2026-04-06 10:20 by 蓝云思雨
[考研] 324求调剂 +3 k可乐 2026-04-05 4/200 2026-04-06 09:54 by 蓝云思雨
[考研] 一志愿C9的化学工程（085602） 340分，感觉校内调剂无望，求调剂 +12 万事宜臻 2026-04-04 12/600 2026-04-06 07:46 by 无际的草原
[考研] 22408 总分320，一篇论文二作，两个国三，求调剂 +3 Leomulufu 2026-04-04 5/250 2026-04-05 19:04 by chongya
[考研] 0703求调剂 +6 zizimo 2026-03-31 6/300 2026-04-04 14:16 by 无际的草原
[考研] 085701求调剂 +7 龚禹铭 2026-04-04 8/400 2026-04-04 13:49 by 小小树2024
[考研] 305求调剂 +3 77Qi 2026-04-03 3/150 2026-04-03 23:01 by qzxyhcsy
[考研] 081200-11408-276学硕求调剂 +5 崔wj 2026-04-03 5/250 2026-04-03 15:06 by arrow8852
[考研] 交通运输考试264分求工科调剂 +4 jike777 2026-04-02 4/200 2026-04-02 21:53 by zllcz
[考研] 0856材料与化工调剂，339 +14 10213207 2026-03-31 14/700 2026-04-02 21:01 by 1104338198
[考研] 生物学327，求调剂 +5 书上的梅子 2026-04-01 6/300 2026-04-02 06:47 by ilovexiaobin
[考研] 省双一流重点一本大学招收调剂 +4 wwwwffffff 2026-03-31 7/350 2026-04-01 15:23 by wwwwffffff
[考研] 一志愿北京科技大学085601材料工程英一数二初试总分335求调剂 +5 双马尾痞老板2 2026-03-31 5/250 2026-04-01 09:04 by oooqiao
[考研] 070300化学专业279调剂 +10 哈哈哈^_^ 2026-03-31 10/500 2026-03-31 23:13 by liu823948201
信息提示
关闭
请填处理意见
关闭