24小时热门版块排行榜    

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

hjlyyc

木虫 (知名作家)

无门派

[求助] 如何利用GULP拟合分子间势能曲线得出力场参数 已有1人参与

想利用两个分子间相互作用力做出的势能曲线拟合得出各个原子之间力场参数,不知GULP能否拟合出结果?

看到Gulp的observable中关于能量的energy,这个应该是total energy,能否通过设置力场参数使得分子间作用力 (这个应该可能通过设置力场参数使之为0)、静电作用力(设置电荷为0)、原子本身能量为0,从而使得energy的值就是分子间作用力?但原子本身能量怎样才能设置为0?

Gulp中有一个reaction energy,不知能否把两分子间作用力理解为reaction energy,来进行拟合呢?

看到mannual中有一句话:Unless these values have been converted to a binding or lattice energy with reference to the dissociated state of the species within the system then it is inappropriate to fit the absolute values of the energies. Consequently, the easiest solution is to include an additive energy shift parameter in the fit, such that only relative energies are actually fitted. 不知这个如何实现?

多谢高人指点!

[ Last edited by hjlyyc on 2013-7-29 at 05:28 ]
回复此楼
天道酬勤!
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

gx730

金虫 (正式写手)

★
月只蓝: 金币+1, 欢迎交流 2013-12-13 18:31:28
引用回帖:
6楼: Originally posted by hjlyyc at 2013-12-13 17:03:38
不好意思,这个问题我还没解决。
理论上你做出了势能曲线,势能是包含分子中各个原子的势能总和,你将每个原子之间的势能按照势能公式写出,再应用到势能曲线,就可以计算出或者拟合出各个原子之间的势能参数。
...

不好意思,我没有理解清楚你所说的。希望你能详细给我讲一下,非常感谢。
我现在用的是TraPPE联合原子力场,势能函数L-J 12-6,只需要拟合两个参数:范德华直径和势阱
我所研究的分子内有六个基团,把每个原子之间的势能按照势能公式写出,再应用到势能曲线,这时候也就是说会有12个拟合参数,这样该如何操作呢?
7楼2013-12-13 17:39:33
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
查看全部 11 个回答

hjlyyc

木虫 (知名作家)

无门派

高人快出现吧
天道酬勤!
3楼2013-07-31 20:54:39
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

gx730

金虫 (正式写手)

★
月只蓝: 金币+1, 欢迎交流提问 2013-12-13 18:30:36
你好
我也想知道你是怎么解决的。有机分子力场参数的拟合用GULP可以实现吗?如何进行?
我现在主要是想拟合分子内某个基团的非键L-J参数,我用MS计算能量的时候是两个分子之间的能量,最后拟合出来的应该是整个分子的非键L-J参数,如何得到某个基团的呢?
5楼2013-12-13 09:53:40
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

hjlyyc

木虫 (知名作家)

无门派

★
月只蓝: 金币+1, 感谢分享心得。 2013-12-14 08:54:26
引用回帖:
5楼: Originally posted by gx730 at 2013-12-13 09:53:40
你好
我也想知道你是怎么解决的。有机分子力场参数的拟合用GULP可以实现吗?如何进行?
我现在主要是想拟合分子内某个基团的非键L-J参数,我用MS计算能量的时候是两个分子之间的能量,最后拟合出来的应该是整个分 ...

不好意思,这个问题我还没解决。
理论上你做出了势能曲线,势能是包含分子中各个原子的势能总和,你将每个原子之间的势能按照势能公式写出,再应用到势能曲线,就可以计算出或者拟合出各个原子之间的势能参数。
我现在可以选几个点计算出势能参数,但参数多了,的需要拟合,这就不会啦。
天道酬勤!
6楼2013-12-13 17:03:38
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[公派出国] 售一区SCI文章T0P,我:8O.551.O54,科目全,可十急 +5 qR2AdCWTVtkw 2026-10-09 9/450 2026-10-11 22:34 by XQNyDCa1Tk5n
[找工作] 售SCI一区T0P文章,我:8.O55.1.O.54,科目全,可十急 +5 qR2AdCWTVtkw 2026-10-09 8/400 2026-10-11 22:26 by XQNyDCa1Tk5n
[公派出国] 售SCI一区文章,我:8O5.5.1.O5.4,科目全,可伽急 +5 qR2AdCWTVtkw 2026-10-09 10/500 2026-10-11 22:26 by XQNyDCa1Tk5n
[考博] 售SCI一区T0P文章,我:8.O55.1.O.54,科目全,可十急 +5 qR2AdCWTVtkw 2026-10-09 9/450 2026-10-11 22:15 by XQNyDCa1Tk5n
[考博] 售SCI一区文章,我:8.O.55.1.O.54,科目齐全,可伽急 +3 qR2AdCWTVtkw 2026-10-09 7/350 2026-10-11 22:11 by XQNyDCa1Tk5n
[基金申请] 售SCI一区T0P文章,我:8.O.55.1.O54,科目全,可伽急 +5 qR2AdCWTVtkw 2026-10-09 9/450 2026-10-11 22:06 by XQNyDCa1Tk5n
[教师之家] 售SCI文章,我:8O5.5.1.O.54,科目齐全,可+急 (EPI+-1)(金币-50) +5 qR2AdCWTVtkw 2026-10-09 10/500 2026-10-11 22:06 by XQNyDCa1Tk5n
[有机交流] OpenAI又攻克数学难题了,离AI做有机合成还有多远? +3 wangyl_123 2026-10-10 6/300 2026-10-11 15:04 by Abin0406
[教师之家] 售SCI一区T0P文章,我:8.O.55.1.O.5.4,科目全,可+急 (EPI+-1)(金币-50) +5 qR2AdCWTVtkw 2026-10-09 10/500 2026-10-11 14:25 by VOMexwnDzOFZ
[教师之家] 售SCI文章,我:8O.5.5.1O.54,科目全,可十急 (EPI+-1)(金币-50) +4 qR2AdCWTVtkw 2026-10-09 7/350 2026-10-11 14:05 by VOMexwnDzOFZ
[考研] 售SCI文章,我:8O.5.5.1O.54,科目全,可十急 +5 qR2AdCWTVtkw 2026-10-09 8/400 2026-10-11 13:51 by VOMexwnDzOFZ
[硕博家园] 售SCI-T0P文章,我:8O.5.5.1.O.54,科目齐全,可+急 +5 qR2AdCWTVtkw 2026-10-09 8/400 2026-10-11 13:45 by VOMexwnDzOFZ
[找工作] 售SCI一区T0P文章,我:8O.55.1.O.5.4,科目齐全,可+急 +3 qR2AdCWTVtkw 2026-10-09 7/350 2026-10-11 10:54 by VOMexwnDzOFZ
[找工作] 售SCI一区T0P文章,我:8O.55.1.O.5.4,科目齐全,可+急 +4 qR2AdCWTVtkw 2026-10-09 7/350 2026-10-11 10:34 by VOMexwnDzOFZ
[考博] 招收博士生一名 +3 mandolin 2026-10-10 3/150 2026-10-11 10:21 by 不是山谷!
[教师之家] 课题组会帮你通过考核么? +3 好运来~9527 2026-10-10 5/250 2026-10-11 09:00 by 好运来~9527
[公派出国] 售SCI一区T0P文章,我:8O.55.1.O.54,科目全,可伽急 +4 qR2AdCWTVtkw 2026-10-09 4/200 2026-10-11 03:14 by M0zrYVSAnxpE
[有机交流] 各位大佬,路线设计求助 +6 zapen 2026-10-08 8/400 2026-10-10 15:59 by zapen
[基金申请] 入职后第一次管课题组的专利缴费,对账才发现官费比别人多掏了一倍 +5 13108017953 2026-10-08 5/250 2026-10-10 11:15 by wdwd123
[有机交流] 我应该选哪条合成路线让实验成功率高一些? +5 jing816y 2026-10-08 9/450 2026-10-10 09:05 by zapen
信息提示
请填处理意见