| 查看: 1370 | 回复: 2 | |||
乐石||木虫 (正式写手)
|
[交流]
LAMMPS, CHARMM-gen force field, topology file builder, PDB to DATA已有1人参与
|
|
可能有时候大家会觉得有机物的data结构文件搭建比较费时，所以会用一些软件进行从pdb到data的自动转换。但是，MS也好，其他软件也好，总有找不到信息的时候，因为ff数据库不可能涵盖所有原子种类。Amber会稍好一些, 因为它只关心原子和几个原子相连，不太关心和谁相连，但是这样可能会降低准确度。如果你用charmm ff的话，可以参照我写的python代码。力场参数文件为par_all36_cgenff.prm。 https://github.com/yangl9/topobuilder/ 该代码也只能够实现从pdb到data转换的半自动化。因为任何时候，par_all36_cgenff.prm是不能给出所有原子种类信息的，这种参数测试工作太耗时了。这时我们需要做一些比较靠谱的近似，cgenff就可以给出和你想要的结构差不多的而且存在于par_all36_cgenff.prm文件中的立场参数。但是问题在于cgenff有时给出的信息也非常不靠谱，比如penalty超过50.（ https://cgenff.paramchem.org/commonFiles/newsEvents.php）这时我们需要手动更正一些信息，这就是半自动化。 我的代码能够自动从par_all36_cgenff.prm搜取已存在的信息，而对于那些不存在的，允许手动输入。当你输入键的信息之后，代码会根据所输入键（1-2）的信息去找其所对应的角的力场参数（1-2-3）；如果你输入了角（1-2-3）的信息，代码会根据所输入角的信息去找相关dihedral的力场参数（1-2-3-4）。（如果存在一个dihedral对应多行信息的情况，也需要去手动更正。） 有任何问题欢迎随时反馈，邮箱请见代码首行。 18.03.2019 |
回复此楼» 本帖已获得的红花（最新10朵）
Nater1ver» 猜你喜欢
孩子确诊有中度注意力缺陷
已经有14人回复
-
三甲基碘化亚砜的氧化反应
已经有4人回复
-
请问下大家为什么这个铃木偶联几乎不反应呢
已经有5人回复
-
请问有评职称，把科研教学业绩算分排序的高校吗
已经有5人回复
-
2025冷门绝学什么时候出结果
已经有3人回复
-
天津工业大学郑柳春团队欢迎化学化工、高分子化学或有机合成方向的博士生和硕士生加入
已经有4人回复
-
康复大学泰山学者周祺惠团队招收博士研究生
已经有6人回复
-
AI论文写作工具：是科研加速器还是学术作弊器？
已经有3人回复
-
论文投稿，期刊推荐
已经有4人回复
-
请问2026国家基金面上项目会启动申2停1吗
已经有5人回复
Nater1ver
木虫 (正式写手)
- 应助: 1 (幼儿园)
- 金币: 3253.6
- 红花: 14
- 帖子: 718
- 在线: 46.8小时
- 虫号: 1526557
- 注册: 2011-12-06
- 专业: 遥感机理与方法
2楼2019-05-15 21:12:36
乐石||
木虫 (正式写手)
- 应助: 26 (小学生)
- 金币: 2416.4
- 散金: 538
- 红花: 9
- 帖子: 977
- 在线: 807.5小时
- 虫号: 1474843
- 注册: 2011-11-03
- 性别: GG
- 专业: 半导体材料
3楼2019-05-23 10:37:04