计算化学公社

 找回密码 Forget password
 注册 Register
Views: 17062|回复 Reply: 157
打印 Print 上一主题 Last thread 下一主题 Next thread

[程序/脚本开发] 分子模拟可视化力场辅助工具AuToFF发布

  [复制链接 Copy URL]

151

帖子

0

威望

3279

eV
积分
3430

Level 5 (御坂)

跳转到指定楼层 Go to specific reply
楼主
本帖最后由 superrice 于 2023-5-17 17:10 编辑

AuToFF (Auxiliary Tools of Force Field) 是一款完全自主开发的可视化力场辅助工具,对外免费开放,采用网页端图形界面设计风格,用户即可点击鼠标实现在线使用,轻松实现分子建模和力场拓扑文件生成,即可一站式获得分子动力学计算软件的力场拓扑文件或建模工具的输入文件,继而直接进行分子动力学模拟,可以极大的降低分子动力学模拟用户使用门槛。
软件免费使用链接https://cloud.hzwtech.com/web/pe ... e-tool/auto-ff-home
软件使用手册链接https://autoff.readthedocs.io/en/latest/index.html
AuToFF工具解决了现有的力场拓扑文件生成工具的诸多局限性,如:无法适应所有的力场形式和多款分子动力学模拟计算软件,缺乏复杂体系建模功能,不能正确分配体系的力场参数,力场参数缺失以及难以计算分配体系中的原子电荷等。

主要功能:
  • 实现2D分子建模及其3D分子显示
  • 支持复杂体系建模,包括碳材料(富勒烯、石墨烯、碳纳米管)和聚合物体系(均聚物、嵌段共聚物、无规共聚物)
  • 提供晶体材料结构库和金属/共轭有机框架结构库
  • 支持现有的主流力场(GAFF、CVFF、OPLS、CGenFF)、二代力场(CFF、COMPASS)、通用力场形式(UFF、Dreiding)、联合原子力场(TraPPE-UA、OPLS-UA、GROMOS)、13种水模型力场(OPC, OPC3, SPC, SPC/E, SPC/Eb, TIP3P, TIPS3P, TIP3P-FB, TIP4P, TIP4P/Ew, TIP4P/2005, TIP4P/ICE, TIP4P/ε )
  • 针对特定体系支持使用力场,如离子液体体系(CL&P、OPLS-2009IL、Merz离子力场、Haiyang Zhang离子力场),界面体系(interface力场),MOF体系(UFF4MOF、UFF4MOFII)、黏土体系(ClayFF)
  • 支持九种类型原子电荷的计算,包括基于机器学习的第一性原理精度的RESP原子电荷、MMFF94、AM1-BCC、Qeq、RESP(GFN2-XTB精度)、CM5、1.2*CM5、1.14*CM1A和1.14*CM1A-LBCC
  • 针对MOF/COF等周期性体系支持多种电荷计算方法,包括GMP-QEq、MEPO-QEq、m-CBAC、ffp4mof
  • 通过选择杂化力场进行单个力场参数缺失的补齐
  • 支持多种结构数据文件导入和分子结构文件下载
  • 支持生成简单分子体系、复杂分子体系、金属/共轭有机框架体系、晶体材料体系、混合体系的力场拓扑文件
  • 支持生成多款分子动力学计算软件的力场拓扑文件或建模工具的输入文件,包括LAMMPS, GROMACS, AMBER, Moltemplate, OpenMM, Tinker, CHARMM

AuToFF未来还有很多新功能会逐步更新,如下:
  • mSeminario方法拟合力常数(已开发完,近期上线)
  • 集成modified ClayFF,OPLS-DES等一些其他的LJ和Morse形式的小众力场
  • 集成可极化力场
  • 非混合规则的原子对LJ参数检索
  • 类似GaussView的3D分子建模
  • 基于空间群的晶体建模
  • ReaxFF反应力场和多体势的数据库
  • 支持生成NAMD和RASPA计算软件输入文件
  • 周期性体系的AM1-BCC、EQEq、I-QEq和MMFF94电荷计算
  • 碳材料添加官能团的建模功能
  • 类似ChemDraw的二维聚合物建模
  • ANI机器学习力场
  • 蛋白质、核酸、多糖体系的残基库、磷脂双分子层建模
  • 全原子—Martini粗粒化力场映射

敬请期待~~~

感谢广大用户使用AuToFF工具,欢迎您提出需求及改进建议,我们都将会积极处理与回复~




常见问题:

  • 上传的结构文件无法显示?
    通常原因都是文件格式不规范导致的,AuToFF开发的时候严格基于文件的官方格式文档,对于一些不规范的格式兼容性较差,多见于pdb格式(symbol列没有写)和cif格式(结构未精修,存在分数占据)。用户可以检查文件格式,或者用其他软件,如VESTA打开后重新保存为其他格式。

    还有一种可能的原因就体系的原子数目太多,目前为了满足大部分用户需求以及保证用户体验,限制上传的结构不超过1000个原子。对于聚合物等大体系可以在聚合物模块让AuToFF来搭建结构,而不要建完模型再上传,因为通过AuToFF的聚合物模块来搭建结构,程序可以充分利用标度算法来减少原子类型识别和原子电荷计算的计算量。对于晶体材料,可以先上传结构较小的晶胞,待最后一步生成对应分子动力学模拟软件的拓扑格式的时候,通过扩胞来创建更大的体系。


    2. 识别完原子类型之后,无法点击下一步操作?

    如果当前的结构中有原子是选择力场中不支持的,程序自动会禁止下一步操作。用户需要返回上一步,更换其他的力场,或者添加备选的杂化力场,如UFF力场来达到能够匹配所有原子的目的。另外需要说明的是,AuToFF是基于力场规则定义好的SMART语法和当前结构中原子的连接关系来判断原子类型的,如果存在不合理的键连关系,也可能导致无法识别原子类型。用户可以通过第三方建模软件修改成键关系来进一步调整,然后重新上传结构。未来AuToFF会直接支持网页端的3d分子建模,用户编辑结构会更加容易。


    3. 体系总电荷不为零?

    若选择力场参数有自带的电荷,如OPLS-AA,ClayFF等,且没有额外选择原子电荷的计算,那么程序不能保证所有原子的总电荷是用户期望的净电荷数。用户可以自行编辑更改电荷数达到目的,或者可以选择AuToFF支持的多种电荷计算方法,实现电荷平衡。


    4. 导入优化后的分子结构进行力场参数生成,可以不改变键长,键角参数吗?

    可以实现的,考虑到用户可能已经用量子化学计算方法优化了分子结构,用户想要将DFT优化得到的键长、键角替换力场中的平衡键长和键角,用户可以在“根据力场选择原子类型”功能区的3D显示的下方,设置使用当前结构的键长,键角的操作。






    5. 力场参数检索方式的模糊和精确按钮如何区分?

    如果使用“模糊检索”,且当用户选择了备选力场进行杂化的时候,程序会按照力场的排序依次检索,如果找到了参数,就不会继续往下检索。因此当前体系不是一套力场能够完全匹配,必须使用杂化力场的时候,必须选择“模糊检索”。
    如果只选择了一种力场,但是这个力场中存在多个原子类型可以匹配某一个原子的时候,如果使用“精确检索”,程序则之后使用匹配到的第一个原子类型进行参数检索(在前端界面上,用户可以点击对应的原子更换优先级)。如果使用“模糊检索”,则程序会按照匹配到的原子类型的优先级,依次进行参数的检索,指直到找到参数。比如,对于羰基C=O,力场中存在精确的C=O定义和更广泛的C=定义,两种原子类型都可以匹配当前碳原子,如果选择“精确检索”,则程序只会按照C=O进行参数的检索。如果选择了“模糊检索”,那么如果当前力场库中找不到C=O相关的参数,程序会进一步检索C=的参数来将就用。
    一般默认使用模糊检索即可。


    6. 可以修改力场吗?

    “生成拓扑文件”功能区可以进行编辑修改力场参数,对于较大的体系,一一修改可能会比较繁琐,建议用户下载参数后用文本编辑器自行批量修改。


    7. 是否支持力场参数的拟合?

    近期会上线mSeminario方法拟合力常数,但是需要用户自己上传Gaussian、ORCA等量化软件计算的Hessian文件,且保证量化计算的原子序号和上传到AuToFF的结构文件的序号保持一致。但是该方法只支持拟合bond和angle的参数。对于多重二面角的参数拟合,AuToFF计划会支持,但是由于需要调用超算资源做量化的势能面扫描,因此这个的开发计划比较后。


    8. CgenFF力场没有像官方那样提供电荷?

    CgenFF的官方电荷模型是基于bond charge increasement 规则来生成了,参见文献:10.1021/ci3003649,由于文献中没有全部公开相关的BCI参数,因此用户只能通过CgenFFF官网上传结构来生成,大多数时候也可以用其他电荷模型将就使用,比如MMFF94和RESP电荷。




评分 Rate

参与人数
Participants 9
eV +41 收起 理由
Reason
耿斌 + 4 赞!
Huschein + 5 赞!
七尺贱 + 5 初学者的福音
XYP + 5
波波波 + 4 好物!
honest + 5 好物!
zjulfp + 3 好物!
funok + 5 好物!
五月雨 + 5 牛!

查看全部评分 View all ratings

28

帖子

0

威望

92

eV
积分
120

Level 2 能力者

159#
发表于 Post on 2024-10-17 10:33:53 | 只看该作者 Only view this author
我使用了Autoff生成了拓扑文件,如何在论文中去引用?

63

帖子

0

威望

550

eV
积分
613

Level 4 (黑子)

158#
发表于 Post on 2024-9-22 16:28:00 | 只看该作者 Only view this author
非常感谢大佬提供如此强力的工具,在下有一个小问题。
在构建COF模型的时候 如果每层是一个分子,构建4x4x4的结构相当于四个分子,但是得到的拓扑文件.top里边分子数只是1,这样会有影响吗?还是我需要先生成4x4x1的,后续自己添加多层?
感谢回复。

1379

帖子

0

威望

3992

eV
积分
5371

Level 6 (一方通行)

157#
发表于 Post on 2024-8-20 18:14:26 | 只看该作者 Only view this author
反馈一下,粗粒化一个分子后,点下载没反应,换了几个浏览器都不行。。但是用全原子模型,就可以正常下载得到的文件。
又菜又爱玩

101

帖子

0

威望

451

eV
积分
552

Level 4 (黑子)

156#
发表于 Post on 2024-8-14 11:05:20 | 只看该作者 Only view this author
superrice 发表于 2023-5-23 09:55
EC结构不属于离子液体,属于溶剂,请移步全原子模块使用。类似于 FSI等离子液体分子才可在此模块正确使用

我想请问一下在拟合力场参数的时候,对于缺失的参数会使用一个较大的力常数进行替代,这样的作法有什么依据吗

312

帖子

2

威望

2306

eV
积分
2658

Level 5 (御坂)

155#
发表于 Post on 2024-8-13 18:06:44 | 只看该作者 Only view this author
slxc920113 发表于 2023-9-26 16:59
GNN的原始训练集的精度很高,而且是RESP2电荷,所以如果元素支持的时候,它的精度更高。GFN2的精度和def2 ...

“GNN的原始训练集的精度很高”,请问这是哪一篇文献了?

265

帖子

0

威望

1209

eV
积分
1474

Level 4 (黑子)

154#
发表于 Post on 2024-6-24 20:39:10 | 只看该作者 Only view this author
davi 发表于 2024-6-20 22:25
@slxc920113 老师,您好。我往上浏览看到了moltemplate版本没法下载,成功下载了lammps版,和文献有点不太 ...

不是published的版本,而是2018年的MS软件内的版本,improper的参数如果为0,AuToFF会自动省略掉的,因为没有意义,浪费计算量而已。

105

帖子

0

威望

502

eV
积分
607

Level 4 (黑子)

153#
发表于 Post on 2024-6-20 22:25:27 | 只看该作者 Only view this author
@slxc920113 老师,您好。我往上浏览看到了moltemplate版本没法下载,成功下载了lammps版,和文献有点不太一样。文献做的是2,2,4-3甲基己烷,其补充材料给出了compass参数。我使用Autoff结果有些不太一样。希望您能够看看
这是文献链接
https://ars.els-cdn.com/content/ ... 1219302808-mmc1.pdf

下面是我的产生的力场文件,其中标注了一些不一样的地方。

UNK.data

12.21 KB, 下载次数 Times of downloads: 1

105

帖子

0

威望

502

eV
积分
607

Level 4 (黑子)

152#
发表于 Post on 2024-6-20 21:25:30 | 只看该作者 Only view this author
slxc920113 发表于 2024-6-20 13:44
晶体模块有,你可以创建一个盒子,放入你的分子,保存为cif文件上传

老师,是没有办法下载下来吗?还是我哪里操作不对呢

509

帖子

1

威望

4249

eV
积分
4778

Level 6 (一方通行)

151#
发表于 Post on 2024-6-20 14:22:09 | 只看该作者 Only view this author
slxc920113 发表于 2023-12-14 17:31
在判断成键的时候,考虑了周期性边界条件,具体实现的算法类似于MD中的邻居列表法+镜像盒子。具体来说就 ...

请问您描述的这部分判断成键的算法是只适用于周期性体系的吗?

105

帖子

0

威望

502

eV
积分
607

Level 4 (黑子)

150#
发表于 Post on 2024-6-20 13:58:30 | 只看该作者 Only view this author
slxc920113 发表于 2024-6-20 13:44
晶体模块有,你可以创建一个盒子,放入你的分子,保存为cif文件上传

确实可以,老师。但是我最后下载有个必填扩胞XYZ,不知道怎么填。貌似怎么填都是下载失败昂?不太理解扩胞含义

265

帖子

0

威望

1209

eV
积分
1474

Level 4 (黑子)

149#
发表于 Post on 2024-6-20 13:44:28 | 只看该作者 Only view this author
davi 发表于 2024-6-20 13:35
@slxc920113 老师,您好。我想使用Autoff实现异丙基环己烷的Compass力场,是不是没有办法啊。貌似Autoff没 ...

晶体模块有,你可以创建一个盒子,放入你的分子,保存为cif文件上传

105

帖子

0

威望

502

eV
积分
607

Level 4 (黑子)

148#
发表于 Post on 2024-6-20 13:35:04 | 只看该作者 Only view this author
@slxc920113 老师,您好。我想使用Autoff实现异丙基环己烷的Compass力场,是不是没有办法啊。貌似Autoff没法支持已经公开的Compass参数,Moltemplate中是有,感觉没有Autoff方便。

2

帖子

0

威望

111

eV
积分
113

Level 2 能力者

147#
发表于 Post on 2024-6-4 13:04:37 | 只看该作者 Only view this author

我试了,在晶体材料模块,基本上所有的材料导出moltemplate都会报错,包括平台预设的材料。

265

帖子

0

威望

1209

eV
积分
1474

Level 4 (黑子)

146#
发表于 Post on 2024-6-3 11:09:55 | 只看该作者 Only view this author
reinhold 发表于 2024-5-31 23:10
您好,晶体材料下,导出moltemplate格式文件总提示下载错误,而lammps格式正常

什么晶体?

手机版 Mobile version|北京科音自然科学研究中心 Beijing Kein Research Center for Natural Sciences|京公网安备 11010502035419号|计算化学公社 — 北京科音旗下高水平计算化学交流论坛 ( 京ICP备14038949号-1 )|网站地图

GMT+8, 2024-11-27 06:08 , Processed in 0.235614 second(s), 26 queries , Gzip On.

快速回复 返回顶部 返回列表 Return to list