24小时热门版块排行榜    

查看: 3528  |  回复: 9

小白如水

铁杆木虫 (小有名气)

[求助] opls下L-J参数转换成amber下L-J参数的问题

我想获得amber下乙烯的力场,查到了乙烯在opls下的参数,想转换到amber下。
查了相关文献,有的把甲烷opls下的参数用到amber力场中,可是我前后对比发现怎么数值不对应呢?比如甲烷中的碳opls下L-J参数是0.35nm和0.276kj/mol,文献中引入到amber下,怎么就变成0.3207和0.251了呢?
请大家帮帮忙。谢谢

amber下甲烷L-J参数



opls下甲烷L-J参数(表中第五六行为甲烷,单位换算后还是对不上)
回复此楼
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
祝福
2楼2012-01-08 19:19:48
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

小白如水

铁杆木虫 (小有名气)

谢谢,挺着急的,文章返修让加上模拟这一块儿,时间比较紧
3楼2012-01-08 19:29:53
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

043114076

木虫 (正式写手)

【答案】应助回帖

★ ★
感谢参与,应助指数 +1
chaizhm(金币+2): 谢谢~ 2012-01-09 09:26:47
小白如水(金币+1): 2012-01-09 21:44:27
能不能把文献附上
你是要把opls的参数转化为amber的形式呢, 还是要和amber力场混合使用? 一般是使用一套力场参数, 实在不行了再引用其他的参数.
4楼2012-01-09 08:48:02
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

小白如水

铁杆木虫 (小有名气)

引用回帖:
: Originally posted by 043114076 at 2012-01-09 08:48:02:
能不能把文献附上
你是要把opls的参数转化为amber的形式呢, 还是要和amber力场混合使用? 一般是使用一套力场参数, 实在不行了再引用其他的参数.

两篇文献分别是
分子模拟研究气体在室温离子液体中的溶解度,Acta Phys.-Chim.Sin.,2005,21(101):1138-1142
Development and Testing of the OPLS All-Atom Force Field
on Conformational Energetics and Properties of Organic Liquids,J. Am. Chem. Soc. 1996, 118, 11225-11236
    opls和amber的力场形式好像是一样的,只是算范德华作用的时候一个用的是L-J参数的算数平均一个用的是几何平均。我也不知道该怎么办。刚接触分子模拟不久,参数弄得我很混乱。。。。请大家指点。
5楼2012-01-09 10:15:31
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

043114076

木虫 (正式写手)

【答案】应助回帖

★ ★ ★
小白如水(金币+1): 2012-01-09 21:43:47
御剑江湖(金币+3): 2012-01-11 19:45:43
应该是中文文献的问题, 不知道从哪引用的参数, 根本对不上
力场的问题最好与比较好的文献中一致, 先解决比较特殊的原子类型
不过你这也太猛了, 完全不明白就敢往文章里写
6楼2012-01-09 13:50:48
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

小白如水

铁杆木虫 (小有名气)

引用回帖:
: Originally posted by 043114076 at 2012-01-09 13:50:48:
应该是中文文献的问题, 不知道从哪引用的参数, 根本对不上
力场的问题最好与比较好的文献中一致, 先解决比较特殊的原子类型
不过你这也太猛了, 完全不明白就敢往文章里写

我还没有改完文章呢。。。原来投的时候没有做分子动力学模拟啊,审稿人非让我加上我也很无奈,实验室倒是恰好有人刚开始做这方面的研究,我就用他们的方法计算,无奈gromacs里没找到乙烯的amber力场参数,搜完文献让我更纠结。。。我感觉应该不会写错的,汪文川老师课题组创建的这个全原子力场也发了很好的文章,这个虽然投到了中文期刊我想也不会有这种低级错误吧。

继续求教大家怎么转化力场参数
7楼2012-01-09 21:43:38
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

043114076

木虫 (正式写手)

【答案】应助回帖

★ ★ ★
chaizhm(金币+3): 谢谢~ 2012-01-10 10:30:23
小白如水(金币+8): 2012-01-10 18:21:33
汪老师的文献使用的甲烷力场确实引用了OPLS力场但引错了, 他们只开发了离子液体的力场, 这个要分清楚. 是不是低级错误, 这个仔细看就能看出.
乙烯的力场应该很常见, amber力场应该也包含了乙烯的力场, 从amber的原始文献或者官方网站上应该可以查到, 至于gromacs怎么加入这些参数, 手册中应该写的很清楚吧, 注意单位就行了. (没用过gromacs, 只用过dlpoly和namd)
8楼2012-01-09 22:43:44
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

小白如水

铁杆木虫 (小有名气)

引用回帖:
8楼: Originally posted by 043114076 at 2012-01-09 22:43:44:
汪老师的文献使用的甲烷力场确实引用了OPLS力场但引错了, 他们只开发了离子液体的力场, 这个要分清楚. 是不是低级错误, 这个仔细看就能看出.
乙烯的力场应该很常见, amber力场应该也包含了乙烯的力场, 从amber的 ...

谢谢你,我还是不要纠结于参数问题了,选opla-aa的参数先算起来再说,马上要放假了,争取年前搞定
9楼2012-01-10 18:21:00
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

夜月魂

新虫 (初入文坛)

请问amber的Bond Stretching potential parameters中 kb (kcal mol -1   )     与       OPLASS的bond 中kj/mol换算后数还是对不上,这个是如何换算的0.0
10楼2020-12-18 11:38:39
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 小白如水 的主题更新
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[文学芳草园] 梦想 +4 myrtle 2026-08-26 4/200 2026-08-31 08:43 by 逸域隐帆
[基金申请] 投票:  有多少人是今天查系统知道结果的? +17 爱看书的可乐 2026-08-26 19/950 2026-08-31 08:28 by xiangy672
[基金申请] 有没有仍没收到信息的 +7 德尚中行 2026-08-27 8/400 2026-08-30 20:52 by purplejack
[基金申请] 国自然面上复盘~欢迎讨论 (金币+15) +15 晴天加油 2026-08-26 16/800 2026-08-29 18:28 by symmetry
[基金申请] 国自然评审意见 +13 wangmingqi 2026-08-28 19/950 2026-08-29 10:22 by Poppy1104
[基金申请] 为什么资助数各大高校都创新高,自己申请怎么就这么难 +12 Kittylucky 2026-08-27 13/650 2026-08-29 00:04 by superceng
[基金申请] 怎么查啊 +6 huang1991js 2026-08-26 6/300 2026-08-28 08:42 by winsaint
[基金申请] 哪位高人中了,把查询到的截图贴出来让我看看,让我长长见识 +5 yuleib84 2026-08-26 6/300 2026-08-28 00:02 by yudaoqian88
[基金申请] 基金未中,这种答复是模板吗? +5 zhaosm1982 2026-08-27 6/300 2026-08-27 16:00 by lfy8008
[基金申请] 怎么看青基中了没有啊 +5 叶九微 2026-08-26 5/250 2026-08-27 10:35 by l_zh2008
[基金申请] 为什么 国际(地区)合作与交流项目 没有放榜? 10+3 majunge000 2026-08-26 11/550 2026-08-27 08:42 by 北京莱茵编辑
[基金申请] 我不理解! +15 Edward_pc 2026-08-26 23/1150 2026-08-26 20:34 by zzuzxg
[基金申请] 出来了 +9 trojank 2026-08-26 9/450 2026-08-26 14:25 by 宝贝虫子
[基金申请] 为什么国自然不能直接公布 +4 bjdxyxy 2026-08-26 4/200 2026-08-26 13:12 by qingmu1201
[基金申请] 国际合作可查了,中了面上 (EPI+1)(金币+50) +18 Ldrop2023 2026-08-26 18/900 2026-08-26 11:15 by cmrandy
[基金申请] 国合可查了 +3 paperzjh 2026-08-26 3/150 2026-08-26 10:41 by LemmonTr
[基金申请] 今天务委会开完了,明天出结果吗 +19 angus9576 2026-08-25 23/1150 2026-08-26 10:03 by zp519
[基金申请] 国合现在查不到了吗? +10 chengyan1220 2026-08-24 20/1000 2026-08-26 08:57 by peasantsprig
[基金申请] 牛来!米来!面来! +8 beefly 2026-08-26 8/400 2026-08-26 08:37 by xuzhipiao
[基金申请] 没有任何消息-是不是就凉了 +9 图啦图啦 2026-08-24 10/500 2026-08-25 11:59 by 南海小哥
信息提示
请填处理意见