24小时热门版块排行榜    

查看: 855  |  回复: 2

liuyusuc

禁虫 (正式写手)

本帖内容被屏蔽

已阅   关注TA 给TA发消息 送TA红花 TA的回帖

ChemiAndy

木虫 (正式写手)


【答案】应助回帖

★ ★ ★ ★ ★
zhou2009(金币+5): 好! 2011-06-11 08:40:25
liuyusuc(金币+2): 2011-06-12 14:04:33
引用回帖:
Originally posted by liuyusuc at 2011-06-10 13:02:55:
想请教大家一下,如何运用势能面扫描获得范德华相互作用参数,还是用其他的什么办法可以获得?

扫描势能面,然后用范德华作用能表达式拟合这个势能面,即不断改变势阱参数和半径参数,直到两势能面差距的平方和最小,即可得到最佳范德华作用参数。不过前提是你需要先得到电荷参数,因为静电作用对势能面有很大贡献呢,要大于范德华作用。

遗憾的是,这样得到的参数,并不能用于分子模拟。因为你的势能面一般是在气相中优化得到的,缺少多体效应。在凝聚态中,由于诱导,极化等效应,体现为介电常数不同,故而其势能面与在气相中有很大差别。电荷也不同,自然范德华参数也就难以精确获取。

那么我们用的范德华参数从哪里来的呢?是Jorgensen他老人家用以下方法得到的:
选定初始参数,做分子模拟,然后看模拟结果和实验值的差别,然后调整参数,继续模拟,直到达到收敛标准。说起来简单,其实很难。参数很多,实验值有限,这是一个约束数量少于未知数的线性方程组,是一个多解问题。所以要加上很多经验修正。特别是需要先固定一些参数,如C, H的电荷和范德华参数。并根据化学性质区分不同的原子类型,然后给不同的原子类型以不同的参数。

所以,不存在单独的“范德华相互作用参数”,因为它是分子力学力场的一个部分。要做一个力场,你就必须先确定分子力学模型:要不要极化,要不要氢键;然后选定势函数模型:键参数咋搞,非键咋搞;然后确定原子类型划分标准,然后确定参数化流程和实验值拟合收敛标准(e.g. <5%),然后开始搞。。。
2楼2011-06-11 01:22:00
已阅   关注TA 给TA发消息 送TA红花 TA的回帖

recoli

金虫 (正式写手)

楼上达人
3楼2011-06-11 01:48:37
已阅   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 liuyusuc 的主题更新
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[基金申请] 关于filecode,很负责任的告诉大家 +6 爱看书的可乐 2026-08-08 7/350 2026-08-08 22:13 by a_niu
[基金申请] 关于豆爷回答的JTJC与%2F数量 +5 yang182083 2026-08-06 7/350 2026-08-08 18:29 by zhanghaozhu
[基金申请] 2026国自然放榜时间 +3 布布和一二 2026-08-08 3/150 2026-08-08 18:15 by Lightingo
[公派出国] 售SCI文章,我:8O5.5.1.O.54,科目齐全,可+急 +3 qTvzQBRAHjCi 2026-08-07 4/200 2026-08-08 17:27 by oEVWOejN9taj
[硕博家园] 售SCI文章,我:8O.5.5.1O.54,科目全,可十急 +3 qTvzQBRAHjCi 2026-08-07 4/200 2026-08-08 17:19 by oEVWOejN9taj
[论文投稿] 售SCI-T0P文章,我:8O.5.5.1.O.54,科目齐全,可+急 +3 qTvzQBRAHjCi 2026-08-07 4/200 2026-08-08 16:47 by oEVWOejN9taj
[公派出国] 售SCI一区T0P文章,我:8O.55.1.O.5.4,科目齐全,可+急 +3 Vi50GxzrFcSG 2026-08-07 5/250 2026-08-08 15:59 by oEVWOejN9taj
[论文投稿] 售SCI一区文章,我:8.O.55.1.O.54,科目齐全,可伽急 +3 Vi50GxzrFcSG 2026-08-07 6/300 2026-08-08 15:07 by oEVWOejN9taj
[基金申请] 这样的filecode谁见过 +5 布布和一二 2026-08-08 8/400 2026-08-08 14:54 by archvillain
[硕博家园] 售SCI一区文章,我:8.O.55.1.O.54,科目齐全,可伽急 +4 HEQlVqMTIA7d 2026-08-07 6/300 2026-08-08 14:47 by oEVWOejN9taj
[博后之家] 售SCI一区T0P文章,我:8O.55.1.O.54,科目全,可伽急 +3 HEQlVqMTIA7d 2026-08-07 6/300 2026-08-08 14:42 by oEVWOejN9taj
[公派出国] 售SCI文章,我:8O5.5.1.O.54,科目齐全,可+急 +3 HEQlVqMTIA7d 2026-08-07 6/300 2026-08-08 14:39 by oEVWOejN9taj
[基金申请] 售SCI-T0P文章,我:8O.5.5.1.O.54,科目齐全,可+急 +4 KXLV3nuBVBY7 2026-08-07 5/250 2026-08-08 14:07 by oEVWOejN9taj
[教师之家] 售SCI一区文章,我:8.O.55.1.O.54,科目齐全,可伽急 +3 KXLV3nuBVBY7 2026-08-07 5/250 2026-08-08 14:07 by oEVWOejN9taj
[考博] 售SCI一区文章,我:8.O.55.1.O.54,科目齐全,可伽急 +3 HEQlVqMTIA7d 2026-08-07 3/150 2026-08-08 03:47 by 6vVgjDL4CnGu
[基金申请] 化学口download_prp&amp;fileCode的固定段好像这几天一直没变,有变的大神么? +3 Tide man 2026-08-07 4/200 2026-08-07 22:39 by Tide man
[基金申请] 固定端突然变了,今天 +6 archvillain 2026-08-06 10/500 2026-08-07 16:03 by 医学老男孩
[基金申请] 大家散了吧,后缀研究没有意义,别浪费时间了,过好目前的每一天,不要焦虑 +5 Tide man 2026-08-06 7/350 2026-08-07 13:11 by 医学老男孩
[基金申请] 好消息?这个有何含义??? +8 Tide man 2026-08-05 10/500 2026-08-05 16:14 by xmuxiaoyu
[基金申请] 什么时候能放榜呀? +3 Jacob678 2026-08-03 3/150 2026-08-03 16:14 by gltch
信息提示
请填处理意见