| 查看: 536 | 回复: 4 | |||
| 当前主题已经存档。 | |||
[交流]
【求助】如何建立二面角的力场参数？
|
|||
|
最近想模拟一个体系，但是缺少部分力场参数，所以想自己创建，不过对二面角参数的一些概念比较模糊，所以向大家请教！文献中给出的二面角公式一般有两种，一种就是Kn/2[1+cos(n*phi-gamma)], 另一种就是傅里叶形式，V1/2（1 + cosΦi） +V2/2(1 - cos2Φi) +V3/2(1 + cos3Φi)。 问题1. 对于第一个公式在拟合二面角的能量曲线的时候是怎么确定力场参数里面的gamma（phase offset）和n（The periodicity of the torsion）的？我看力场里面的gamma一般是0度或180度，但是不知道这两个是怎么来的，不能是其它角度吗？ 问题2. 这两种形式间的二面角力场参数如何转化？例如有傅里叶形式的二面角参数0，0 ，1.4510，0，怎么将它转化成第一种形式的二面角参数？ 问题3. 这个问题是针对二面角参数的拟合。有文献提到公式拟合的二面角能量曲线应当是ab intio得到的能量曲线剪掉其中包含的1-4相互作用能量部分，在这里我有一个问题，如果有一个二面角参数未知，我将这个二面角的力场参数设为0，能得到他的1-4相互作用能量，但是当同一个分子中有两个以上的未知二面角参数的时候又应该怎么办呢？ 问题比较多，请高手们不吝赐教！ [ Last edited by wuli8 on 2009-11-13 at 22:53 ] |
回复此楼» 猜你喜欢
311求调剂
已经有21人回复
-
343求调剂085601
已经有7人回复
-
327求调剂
已经有9人回复
-
297求调剂
已经有5人回复
-
材料调剂
已经有7人回复
-
283分求调剂
已经有5人回复
-
336求调剂
已经有8人回复
-
329求调剂，一志愿西北工业大学，材料工程（085601）
已经有14人回复
-
兽医调剂
已经有3人回复
-
材料334求调剂
已经有4人回复
» 本主题相关商家推荐: (我也要在这里推广)
yyx19840628
木虫 (著名写手)
鹰的眼睛，狼的耳朵，豹的速 ...
- 应助: 4 (幼儿园)
- 贵宾: 0.41
- 金币: 3521.9
- 散金: 50
- 红花: 8
- 帖子: 2321
- 在线: 129小时
- 虫号: 508615
- 注册: 2008-02-21
- 性别: MM
- 专业: 环境污染化学
2楼2008-12-10 11:39:00
lianghust
铁杆木虫 (正式写手)
- 应助: 4 (幼儿园)
- 金币: 6353
- 散金: 1
- 帖子: 711
- 在线: 322.1小时
- 虫号: 186357
- 注册: 2006-02-16
- 性别: GG
- 专业: 理论和计算化学
3楼2008-12-10 11:47:38
|
我在分子模拟论坛上也问了这个问题，有个人说“傅立叶的1，3力常数都是0，实际上就是和第一个一样的，你确定n和gamma使余弦项与傅立叶形式相同就可以啦”！但是这里我不是很明白，对于傅里叶形式的力场参数1，3项不是0的时候就不能转化成第一种形式了吗？ 还有一个问题，对于第一个公式的中的n，molecular Modelling Principles And Applications 2Ed这本书上说是the multiplicity; its value gives the number of minimum points in the function as the bond is rotated through 360度. 而gamma (the phase factor) determines where the torsional angle passes through its minimum value.但是我有点不理解他所举的例子，书上说For example, the energy profile for rotation about the single bond between two sp3 carbon atoms could be represented by a single torsional term with n=3 and gamma=0. this would give a threefold rotational profile with minima at torsion angles of +60,-60 and 180, and maxima at +/-120 and 0. 它所给的图如下，  在图中n=3这个能够理解，但是对gamma值就不是很明白了！他明明说gamma应该取能量的最低点，可是是在gamma=0的时候分明是能量最高点啊！这是什么原因呢？ |
4楼2008-12-17 09:13:18
wgpig
木虫 (小有名气)
- 应助: 0 (幼儿园)
- 金币: 2504.6
- 帖子: 191
- 在线: 273.7小时
- 虫号: 35477
- 注册: 2004-02-03
- 性别: GG
- 专业: 生物大分子结构与功能
5楼2008-12-17 09:48:53
