计算化学公社

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

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

  [复制链接 Copy URL]

156

帖子

0

威望

3743

eV
积分
3899

Level 5 (御坂)

本帖最后由 superrice 于 2025-9-22 21:02 编辑

[楼主本人已经从鸿之微离职了,请不要再私信了。有问题请直接跟帖,会有工程师为您解答]
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 15
eV +63 收起 理由
Reason
梁慷 + 5
wwy223 + 5 好物!
樊樊樊 + 4
iamback2011 + 2 赞!
MUMBER + 3 赞!
PLwang + 3
耿斌 + 4 赞!
Huschein + 5 赞!
七尺贱 + 5 初学者的福音
XYP + 5
波波波 + 4 好物!
honest + 5 好物!
zjulfp + 3 好物!
funok + 5 好物!
五月雨 + 5 牛!

查看全部评分 View all ratings

2

帖子

0

威望

99

eV
积分
101

Level 2 能力者

193#
发表于 Post on 前天 09:45 | 只看该作者 Only view this author
本帖最后由 magicrao 于 2026-1-21 09:48 编辑

请教一个聚合物建模的问题,我的目标是构建一个4臂PLGA聚合物分子,但在操作过程中发现当添加的臂数超过2臂时(即在第3根臂上添加重复单元时),把对应的MRV字符串复制到AutoFF中试图生成3d结构就会报错(2臂时是可以正常生成3d结构的),请问这个有办法解决吗?

202601210944325443..png (22.9 KB, 下载次数 Times of downloads: 1)

202601210944325443..png

7

帖子

0

威望

301

eV
积分
308

Level 3 能力者

192#
发表于 Post on 2026-1-11 14:00:32 | 只看该作者 Only view this author
牧生 发表于 2026-1-11 13:44
没有解决,没有途径,不做粗粒化。

哈哈好吧~

1557

帖子

0

威望

4947

eV
积分
6504

Level 6 (一方通行)

191#
发表于 Post on 2026-1-11 13:44:23 | 只看该作者 Only view this author
wwy223 发表于 2026-1-11 13:27
好吧,那您后面是通过啥途径得到粗粒化分子的模型和拓扑文件的呢

没有解决,没有途径,不做粗粒化。
又菜又爱玩

7

帖子

0

威望

301

eV
积分
308

Level 3 能力者

190#
发表于 Post on 2026-1-11 13:27:55 | 只看该作者 Only view this author
牧生 发表于 2026-1-10 17:19
我用粗粒化,从来没下载成功过。你好歹还有成功的

好吧,那您后面是通过啥途径得到粗粒化分子的模型和拓扑文件的呢

1557

帖子

0

威望

4947

eV
积分
6504

Level 6 (一方通行)

189#
发表于 Post on 2026-1-10 17:19:54 | 只看该作者 Only view this author
wwy223 发表于 2026-1-10 16:50
您好,请问这个问题您最后解决了吗,我今天在粗粒化SDS和BS-12的时候也碰到了同样的问题,点下载没反应, ...

我用粗粒化,从来没下载成功过。你好歹还有成功的
又菜又爱玩

7

帖子

0

威望

301

eV
积分
308

Level 3 能力者

188#
发表于 Post on 2026-1-10 16:50:11 | 只看该作者 Only view this author
牧生 发表于 2024-8-20 18:14
反馈一下,粗粒化一个分子后,比如十六烷基三甲基氯化铵的阳离子,点下载没反应,换了几个浏览器都不行。。 ...

您好,请问这个问题您最后解决了吗,我今天在粗粒化SDS和BS-12的时候也碰到了同样的问题,点下载没反应,但是我粗粒化另一个比较复杂的分子(包含一百多个原子)却可以下载,不过这个分子不带电,我感觉不能下载跟粗粒化的分子带电荷有关系。

17

帖子

0

威望

67

eV
积分
84

Level 2 能力者

187#
发表于 Post on 2026-1-5 14:47:14 | 只看该作者 Only view this author
AUTOFF不支持生成不对称的含氮阴离子的CL&P参数吗,想要获得如图所示的离子的CL&P参数,n原子一直显示DU,选了别的力场作为补充(oplsaa/l、oplsaa/m),n原子参数有了,但二面角参数有的缺失

342

帖子

1

威望

1736

eV
积分
2098

Level 5 (御坂)

186#
发表于 Post on 2026-1-5 11:52:20 | 只看该作者 Only view this author
downpuor 发表于 2026-1-5 10:58
运用AUTOFF自带的绘图工具 创建分子结构后 生成的拓补函数忽略了H是为什么呀?

用的联合原子力场吧

1

帖子

0

威望

105

eV
积分
106

Level 2 能力者

185#
发表于 Post on 2026-1-5 10:58:20 | 只看该作者 Only view this author
运用AUTOFF自带的绘图工具 创建分子结构后 生成的拓补函数忽略了H是为什么呀?

3

帖子

0

威望

85

eV
积分
88

Level 2 能力者

184#
发表于 Post on 2025-9-10 01:25:53 | 只看该作者 Only view this author
您好,请问为什么IRMOF1中的C原子的 lennard-jones  值都是一样的呢,视乎好像区分不开苯环上的C、CH的区别

202509100123212933..png (44.53 KB, 下载次数 Times of downloads: 1)

202509100123212933..png

71

帖子

0

威望

6801

eV
积分
6872

Level 6 (一方通行)

183#
发表于 Post on 2025-5-31 18:49:29 | 只看该作者 Only view this author
xiaowu759 发表于 2025-5-30 16:59
生成的文件中重命名了原子名称,由于原子数较多,名称字符溢出,导致同名,请问如何解决?

能否不重命名,这样才方便核对修改

71

帖子

0

威望

6801

eV
积分
6872

Level 6 (一方通行)

182#
发表于 Post on 2025-5-30 16:59:02 | 只看该作者 Only view this author
生成的文件中重命名了原子名称,由于原子数较多,名称字符溢出,导致同名,请问如何解决?

62

帖子

0

威望

605

eV
积分
667

Level 4 (黑子)

181#
发表于 Post on 2025-5-23 00:11:29 | 只看该作者 Only view this author
slxc920113 发表于 2025-5-22 10:14
超过100个原子的体系,不会做结构优化,直接计算电荷,如果你上传的结构本身很离谱,那么半经验方法计算 ...

学到了,谢谢老师

342

帖子

1

威望

1736

eV
积分
2098

Level 5 (御坂)

180#
发表于 Post on 2025-5-22 10:14:43 | 只看该作者 Only view this author
sarphuart 发表于 2025-5-22 09:31
不是为不为零的问题,而是生成的电荷基本都是瞎编,明显不在合理的范围之内,试了三四个200左右的原子都 ...

超过100个原子的体系,不会做结构优化,直接计算电荷,如果你上传的结构本身很离谱,那么半经验方法计算出来的电荷也会很离谱。所以才有专门的聚合物模块,根据单体来计算电荷。

评分 Rate

参与人数
Participants 1
eV +2 收起 理由
Reason
sarphuart + 2 谢谢

查看全部评分 View all ratings

62

帖子

0

威望

605

eV
积分
667

Level 4 (黑子)

179#
发表于 Post on 2025-5-22 09:31:24 | 只看该作者 Only view this author
slxc920113 发表于 2025-5-21 18:44
OPLS自带的电荷不能保证总和为零的。

不是为不为零的问题,而是生成的电荷基本都是瞎编,明显不在合理的范围之内,试了三四个200左右的原子都这样

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

GMT+8, 2026-1-23 19:39 , Processed in 0.246569 second(s), 26 queries , Gzip On.

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