| 查看: 1695 | 回复: 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» 猜你喜欢
售SCI一区T0P文章，我:8.O.55.1.O.5.4,科目全,可+急
已经有16人回复
-
售SCI一区T0P文章，我:8.O.55.1.O.5.4,科目全,可+急
已经有14人回复
-
售SCI一区文章，我:8.O.55.1.O.54,科目齐全,可伽急
已经有15人回复
-
售SCI一区T0P文章，我:8O.55.1.O.54,科目全,可伽急
已经有17人回复
-
售SCI-T0P文章，我:8O.5.5.1.O.54,科目齐全,可+急
已经有21人回复
-
售SCI一区T0P文章，我:8.O.55.1.O.5.4,科目全,可+急
已经有21人回复
-
售SCI一区文章，我:8O5.5.1.O5.4,科目全,可伽急
已经有8人回复
-
售SCI一区文章，我:8.O.55.1.O.54,科目齐全,可伽急
已经有8人回复
-
售SCI一区T0P文章，我:8.O55.1.O.54,科目全,可十急
已经有8人回复
-
售SCI文章，我:8O5.5.1.O.54,科目齐全,可+急
已经有6人回复
Nater1ver
木虫 (正式写手)
- 应助: 1 (幼儿园)
- 金币: 3501.2
- 红花: 14
- 帖子: 772
- 在线: 48小时
- 虫号: 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
