计算化学公社

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

[GROMACS] 求助生成TFSI离子的gaff参数或者oplsaa参数

[复制链接 Copy URL]

17

帖子

0

威望

67

eV
积分
84

Level 2 能力者

我想模拟阳离子为li,阴离子为tfsi的一个聚合物体系,但是用sobtop生成的gaff力场的 itp文件会缺失二面角参数,我尝试使用ligpargen生成oplsaa参数也有缺失,想求助一下各位大佬还有什么方法可以得到呢?

tfsi1.mol2

940 Bytes, 下载次数 Times of downloads: 1

1557

帖子

0

威望

4947

eV
积分
6504

Level 6 (一方通行)

2#
发表于 Post on 2025-12-14 20:52:01 | 只看该作者 Only view this author
本帖最后由 牧生 于 2025-12-15 11:07 编辑

sobtop生成的物质的itp一般是合理的。。。如果自己很介意那个参数,就去看看sobtop的例1




又菜又爱玩

6万

帖子

99

威望

5万

eV
积分
124667

管理员

公社社长

3#
发表于 Post on 2025-12-15 01:10:07 | 只看该作者 Only view this author
把你给sobtop用的输入文件上传,大概率输入文件本来就不合理(诸如mol2里是复合物而不是单个离子,或者连接关系本来就不对,这些问题在http://sobereva.com/soft/Sobtop教程里都明确强调了)
sobtop产生TFSI的拓扑文件毫无问题,产生的有问题只可能是用法不对。遇到问题先检查自己的用法试图解决,别上来就想着换程序,否则纯粹是绕弯路
北京科音自然科学研究中心http://www.keinsci.com)致力于计算化学的发展和传播,长期开办极高质量的各种计算化学类培训:初级量子化学培训班中级量子化学培训班高级量子化学培训班量子化学波函数分析与Multiwfn程序培训班分子动力学与GROMACS培训班CP2K第一性原理计算培训班,内容介绍以及往届资料购买请点击相应链接查看。这些培训是计算化学从零快速入门以及进一步全面系统性提升研究水平的高速路!培训各种常见问题见《北京科音办的培训班FAQ》
欢迎加入北京科音微信公众号获取北京科音培训的最新消息,并避免错过网上有价值的计算化学文章!
欢迎加入人气极高、专业性特别强的理论与计算化学综合交流群思想家公社QQ群(群号见此链接),合计达一万多人。北京科音培训班的学员在群中可申请VIP头衔,提问将得到群主Sobereva的最优先解答。
思想家公社的门口Blog:http://sobereva.com(发布大量原创计算化学相关博文)
Multiwfn主页:http://sobereva.com/multiwfn(十分强大、极为流行的量子化学波函数分析程序)
Google Scholar:https://scholar.google.com/citations?user=tiKE0qkAAAAJ
ResearchGate:https://www.researchgate.net/profile/Tian_Lu

17

帖子

0

威望

67

eV
积分
84

Level 2 能力者

4#
 楼主 Author| 发表于 Post on 2025-12-15 09:58:43 | 只看该作者 Only view this author
本帖最后由 soyolalla 于 2025-12-15 10:03 编辑
sobereva 发表于 2025-12-15 01:10
把你给sobtop用的输入文件上传,大概率输入文件本来就不合理(诸如mol2里是复合物而不是单个离子,或者连接 ...

老师我已修改帖子上传了mol2文件,我的是有的二面角参数缺失,我可以使用m2seminario的方式生成参数吗

202512150958503777..png (71.21 KB, 下载次数 Times of downloads: 0)

202512150958503777..png

17

帖子

0

威望

67

eV
积分
84

Level 2 能力者

5#
 楼主 Author| 发表于 Post on 2025-12-15 10:08:53 | 只看该作者 Only view this author
牧生 发表于 2025-12-14 20:52
sobtop生成的物质的itp一般是合理的。。。如果自己很介意那个参数,就去看看sobtop的例1

老师想问一下您是哪种方法计算的?我运用的是 opt freq=noraman B3LYP/6-31+g(d,p) em=GD3BJ,直接生成参数有的二面角参数缺失,目前想用sob老师的m2seminario方法

1557

帖子

0

威望

4947

eV
积分
6504

Level 6 (一方通行)

6#
发表于 Post on 2025-12-15 10:15:20 | 只看该作者 Only view this author
soyolalla 发表于 2025-12-15 10:08
老师想问一下您是哪种方法计算的?我运用的是 opt freq=noraman B3LYP/6-31+g(d,p) em=GD3BJ,直接生成参 ...

看这个例子

http://sobereva.com/soft/Sobtop/#ex1
又菜又爱玩

1208

帖子

6

威望

2430

eV
积分
3758

Level 5 (御坂)

傻傻的木瓜

7#
发表于 Post on 2025-12-15 10:38:40 | 只看该作者 Only view this author
本帖最后由 Uus/pMeC6H4-/キ 于 2025-12-15 11:11 编辑
牧生 发表于 2025-12-14 20:52
sobtop生成的物质的itp一般是合理的。。。如果自己很介意那个参数,就去看看sobtop的例1

可以上传一下结构优化的输入输出文件以及给sobtop用的mol2文件么?这个拓扑文件感觉不太对劲,按对称性的话(F3CSO2)2N-应该每个元素都是一种原子吧,然而这里面O5和其他氧原子不是一种;本来N原子只与两个S原子成键的,但这键角项里有个N2-O5,又产生了由mSeminario确定的几个键角项O5-N2-S15、N2-O5-S10及几个二面角项C3-S15-N2-O5、O5-N2-S15-O6、O5-N2-S15-O14。

编辑:从chg文件看一开始几何优化就跑到别的不是TFSI-的结构去了,难怪会有N-O键(汗)这也提醒所有懒人脚本的用户,人工检查结构合理性很有必要。

√546=23.36664289109

1557

帖子

0

威望

4947

eV
积分
6504

Level 6 (一方通行)

8#
发表于 Post on 2025-12-15 10:54:50 | 只看该作者 Only view this author
本帖最后由 牧生 于 2025-12-15 11:25 编辑
Uus/pMeC6H4-/キ 发表于 2025-12-15 10:38
可以上传一下结构优化的输入输出文件以及给sobtop用的mol2文件么?这个拓扑文件感觉不太对劲,按对称性的 ...

果然有问题。。我改改

应该是初始文件十分不合理导致的,而我当时没有认真去看。


这次好了



90076-65-6.itp (9.43 KB, 下载次数 Times of downloads: 9)

90076-65-6.mol2 (917 Bytes, 下载次数 Times of downloads: 4)

又菜又爱玩

1208

帖子

6

威望

2430

eV
积分
3758

Level 5 (御坂)

傻傻的木瓜

9#
发表于 Post on 2025-12-15 14:47:52 | 只看该作者 Only view this author
牧生 发表于 2025-12-15 10:54
果然有问题。。我改改

应该是初始文件十分不合理导致的,而我当时没有认真去看。

那些C-S-N-S或者O-S-N-S二面角是刚性的还是柔性的?总担心用mSeminario方法得到的谐振势描述会被质疑,可能还是得柔性扫描、拟合周期性旋转的势能曲线吧,这就比较复杂了。由于这个极小点的构象不对称,两个Tf基团上对应的原子电荷不同,这也不好解释。当然Cs对称性的构象在气相都有破坏对称性的振动虚频,不太确定能否拿非极小点的结构算原子电荷,但在溶液中跑MD并没局限于极小点附近,还是应该能采样到对称的瞬时构象才对,这样还是需要用柔性二面角描述旋转与用等价约束拟合原子电荷。
√546=23.36664289109

17

帖子

0

威望

67

eV
积分
84

Level 2 能力者

10#
 楼主 Author| 发表于 Post on 2025-12-15 15:49:11 | 只看该作者 Only view this author
Uus/pMeC6H4-/キ 发表于 2025-12-15 14:47
那些C-S-N-S或者O-S-N-S二面角是刚性的还是柔性的?总担心用mSeminario方法得到的谐振势描述会被质疑,可 ...

老师可以讲解一下如何进行柔性扫描、拟合周期性旋转的势能曲线吗,我有点不太理解

6万

帖子

99

威望

5万

eV
积分
124667

管理员

公社社长

11#
发表于 Post on 2025-12-15 23:43:44 | 只看该作者 Only view this author
soyolalla 发表于 2025-12-15 15:49
老师可以讲解一下如何进行柔性扫描、拟合周期性旋转的势能曲线吗,我有点不太理解

详谈使用Gaussian做势能面扫描
http://sobereva.com/474http://bbs.keinsci.com/thread-12660-1-1.html
对得到的扫描曲线,自己用Origin或其它方式拟合成扭转势
北京科音自然科学研究中心http://www.keinsci.com)致力于计算化学的发展和传播,长期开办极高质量的各种计算化学类培训:初级量子化学培训班中级量子化学培训班高级量子化学培训班量子化学波函数分析与Multiwfn程序培训班分子动力学与GROMACS培训班CP2K第一性原理计算培训班,内容介绍以及往届资料购买请点击相应链接查看。这些培训是计算化学从零快速入门以及进一步全面系统性提升研究水平的高速路!培训各种常见问题见《北京科音办的培训班FAQ》
欢迎加入北京科音微信公众号获取北京科音培训的最新消息,并避免错过网上有价值的计算化学文章!
欢迎加入人气极高、专业性特别强的理论与计算化学综合交流群思想家公社QQ群(群号见此链接),合计达一万多人。北京科音培训班的学员在群中可申请VIP头衔,提问将得到群主Sobereva的最优先解答。
思想家公社的门口Blog:http://sobereva.com(发布大量原创计算化学相关博文)
Multiwfn主页:http://sobereva.com/multiwfn(十分强大、极为流行的量子化学波函数分析程序)
Google Scholar:https://scholar.google.com/citations?user=tiKE0qkAAAAJ
ResearchGate:https://www.researchgate.net/profile/Tian_Lu

1208

帖子

6

威望

2430

eV
积分
3758

Level 5 (御坂)

傻傻的木瓜

12#
发表于 Post on 2025-12-16 00:39:28 | 只看该作者 Only view this author
本帖最后由 Uus/pMeC6H4-/キ 于 2025-12-16 01:11 编辑

社长,J. Phys. Chem. B 2021, 125, 26, 7231–7240的附录列出了用antechamber生成的GAFF力场的TFSI-拓扑文件,但是其中二面角的functype全是3,这样的拓扑文件可以和sobtop生成的其他functype为2或9的拓扑文件共用于同一体系吗?如何把这种Ryckaert-Bellemans形式的参数等价地换算成普通周期性形式呢?

编辑:该作者的另一篇文章Mater. Adv., 2021,2, 7691-7700的附录除NTf2-外还有OTf-的拓扑文件
√546=23.36664289109

6万

帖子

99

威望

5万

eV
积分
124667

管理员

公社社长

13#
发表于 Post on 2025-12-16 00:46:47 | 只看该作者 Only view this author
Uus/pMeC6H4-/キ 发表于 2025-12-16 00:39
社长,J. Phys. Chem. B 2021, 125, 26, 7231–7240的附录列出了用antechamber生成的GAFF力场的TFSI-拓扑 ...

不同形式的二面角项可以混用,不用转换

评分 Rate

参与人数
Participants 1
eV +1 收起 理由
Reason
Uus/pMeC6H4-/キ + 1 谢谢

查看全部评分 View all ratings

北京科音自然科学研究中心http://www.keinsci.com)致力于计算化学的发展和传播,长期开办极高质量的各种计算化学类培训:初级量子化学培训班中级量子化学培训班高级量子化学培训班量子化学波函数分析与Multiwfn程序培训班分子动力学与GROMACS培训班CP2K第一性原理计算培训班,内容介绍以及往届资料购买请点击相应链接查看。这些培训是计算化学从零快速入门以及进一步全面系统性提升研究水平的高速路!培训各种常见问题见《北京科音办的培训班FAQ》
欢迎加入北京科音微信公众号获取北京科音培训的最新消息,并避免错过网上有价值的计算化学文章!
欢迎加入人气极高、专业性特别强的理论与计算化学综合交流群思想家公社QQ群(群号见此链接),合计达一万多人。北京科音培训班的学员在群中可申请VIP头衔,提问将得到群主Sobereva的最优先解答。
思想家公社的门口Blog:http://sobereva.com(发布大量原创计算化学相关博文)
Multiwfn主页:http://sobereva.com/multiwfn(十分强大、极为流行的量子化学波函数分析程序)
Google Scholar:https://scholar.google.com/citations?user=tiKE0qkAAAAJ
ResearchGate:https://www.researchgate.net/profile/Tian_Lu

17

帖子

0

威望

67

eV
积分
84

Level 2 能力者

14#
 楼主 Author| 发表于 Post on 2025-12-16 14:31:06 | 只看该作者 Only view this author
Uus/pMeC6H4-/キ 发表于 2025-12-16 00:39
社长,J. Phys. Chem. B 2021, 125, 26, 7231–7240的附录列出了用antechamber生成的GAFF力场的TFSI-拓扑 ...

老师我刚看了这个参数文件,tfsi离子不是一个对称的离子吗,怎么硫的原子类型还能出现两种啊?这种结果可靠吗,而且比如说s-c,s-n键长也不相等

1208

帖子

6

威望

2430

eV
积分
3758

Level 5 (御坂)

傻傻的木瓜

15#
发表于 Post on 2025-12-16 15:07:37 | 只看该作者 Only view this author
soyolalla 发表于 2025-12-16 14:31
老师我刚看了这个参数文件,tfsi离子不是一个对称的离子吗,怎么硫的原子类型还能出现两种啊?这种结果可 ...

Tf2N-结构里两边的三氟甲磺酰基确实会有排斥,用量子化学方法优化出来的能量最低构象并不是镜面对称的,而对称的构象并不是能量极小点;如果原文作者是基于优化后的结构制作的拓扑文件,当然不能保证完全对称嘛,这种细节就不是问我了。在补充材料里有完整拓扑文件的文献不见得很容易找,实在要获得对称的参数还得多花些功夫(当然找到了也欢迎贴出doi供其他人参考),要么就仍然自己制作。

虽然sy和s6看起来是两种原子类型,但具体有什么差别还得检查文件里的原子电荷、质量以及力场定义里的L-J参数等。视研究目的与精度,直接把比较相近的参数值取平均也未尝不可接受。
√546=23.36664289109

本版积分规则 Credits rule

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

GMT+8, 2026-1-23 15:14 , Processed in 0.337933 second(s), 24 queries , Gzip On.

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