计算化学公社

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

[综合交流] 如果聚合物单体的原子数很多,那么应如何获取聚合物的原子电荷?

[复制链接 Copy URL]

198

帖子

0

威望

1154

eV
积分
1352

Level 4 (黑子)

如题,如果在进行涉及聚合物的动力学计算时,想要获得较准确的各项参数,其中也包括原子电荷,http://sobereva.com/soft/Sobtop中sob老师提到可以通过如下方法获得原子电荷,那么,如果聚合物的单体原子数就很多的话,这种情况下进行DFT计算可能会非常费劲,甚至算不动,有没有什么好的办法可以获得聚合度为n的聚合物的各项参数(包括原子电荷)呢?

0.png (428.18 KB, 下载次数 Times of downloads: 10)

0.png
社会浮躁

1480

帖子

0

威望

4548

eV
积分
6028

Level 6 (一方通行)

2#
发表于 Post on 2025-4-11 16:58:45 | 只看该作者 Only view this author
本帖最后由 牧生 于 2025-4-11 17:01 编辑

ztop可以。   http://bbs.keinsci.com/forum.php ... &highlight=ztop

但我一直在等那个完美的、便利的解决方案

评分 Rate

参与人数
Participants 1
eV +5 收起 理由
Reason
NGC626 + 5 谢谢

查看全部评分 View all ratings

又菜又爱玩

320

帖子

0

威望

1485

eV
积分
1805

Level 5 (御坂)

3#
发表于 Post on 2025-4-11 23:12:22 | 只看该作者 Only view this author

6万

帖子

99

威望

5万

eV
积分
120151

管理员

公社社长

4#
发表于 Post on 2025-4-12 01:21:33 | 只看该作者 Only view this author
单体原子数很多的情况,意味着你不需要把三个单元全都用完整单元,相邻的单元可以充分简化。例如的氨基酸残基用Multiwfn算RESP电荷,你只需要把一个残基两边用乙酰和甲胺封闭就行了(这就体现出了中间的残基在两个肽键之间),乙酰和甲胺也没多少原子,并不需要非得是靠三肽模型来用Multiwfn算中间的残基的RESP电荷。
北京科音自然科学研究中心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

6万

帖子

99

威望

5万

eV
积分
120151

管理员

公社社长

5#
发表于 Post on 2025-4-12 01:39:24 | 只看该作者 Only view this author
我专门改写了http://sobereva.com/soft/Sobtop网页里的内容,避免有人再有类似问题
一般来说,正确得到聚合物RESP电荷的做法是构造一个寡聚物(目的是让被计算的残基在寡聚物中所处的电子结构环境和实际体系足够相似。如果重复单元很小像聚乙烯,建议含五个单元;如果重复单元稍大如聚碳酸酯,三个单元就够;如果重复单元很大,如包含几十个原子,那就只保留它与相邻单元成键的一部分,如几个到十几个原子左右即可)

北京科音自然科学研究中心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

198

帖子

0

威望

1154

eV
积分
1352

Level 4 (黑子)

6#
 楼主 Author| 发表于 Post on 2025-4-12 17:23:28 | 只看该作者 Only view this author
sobereva 发表于 2025-4-12 01:39
我专门改写了http://sobereva.com/soft/Sobtop网页里的内容,避免有人再有类似问题

感谢sob老师
社会浮躁

198

帖子

0

威望

1154

eV
积分
1352

Level 4 (黑子)

7#
 楼主 Author| 发表于 Post on 2025-4-13 18:40:13 | 只看该作者 Only view this author
本帖最后由 NGC626 于 2025-4-13 18:41 编辑
sobereva 发表于 2025-4-12 01:39
我专门改写了http://sobereva.com/soft/Sobtop网页里的内容,避免有人再有类似问题

sob老师好,请问sobtop的介绍中提到的计算聚合物电荷的方法,"第一次算RESP电荷的时候把首端单元的总电荷约束为0",以此类推,得到首端、中间、尾端单元的原子电荷,请问为什么计算某部分的电荷时,是把这部分的总电荷约束为0,而不是约束除它以外的其他部分的电荷为0呢?
社会浮躁

6万

帖子

99

威望

5万

eV
积分
120151

管理员

公社社长

8#
发表于 Post on 2025-4-14 01:37:19 | 只看该作者 Only view this author
NGC626 发表于 2025-4-13 18:40
sob老师好,请问sobtop的介绍中提到的计算聚合物电荷的方法,"第一次算RESP电荷的时候把首端单元的总电荷 ...

这样就已经达到目的了,干嘛还用其它做法
北京科音自然科学研究中心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

198

帖子

0

威望

1154

eV
积分
1352

Level 4 (黑子)

9#
 楼主 Author| 发表于 Post on 2025-4-17 14:01:54 | 只看该作者 Only view this author
sobereva 发表于 2025-4-14 01:37
这样就已经达到目的了,干嘛还用其它做法

sob老师好,请问在按照sobtop的介绍中计算聚合物原子电荷的时候,假设三聚物的总电荷为0,首端、中间、末端单元的电荷分别为-2,1,1,那么是不是意味着四聚物的总电荷是1,五聚物的总电荷是2,以此类推呢?
社会浮躁

6万

帖子

99

威望

5万

eV
积分
120151

管理员

公社社长

10#
发表于 Post on 2025-4-17 14:48:53 | 只看该作者 Only view this author
NGC626 发表于 2025-4-17 14:01
sob老师好,请问在按照sobtop的介绍中计算聚合物原子电荷的时候,假设三聚物的总电荷为0,首端、中间、末 ...

我不晓得你的这个聚合物是什么情况。弄清楚为什么残基带净电荷,从逻辑和实际状态角度思考该是什么样。诸如中性水中的DNA当中每个重复单元都有去质子化的磷酸基,因此带负电荷,这会相应地由体系中的抗衡离子中和其电荷。
北京科音自然科学研究中心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

本版积分规则 Credits rule

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

GMT+8, 2025-8-16 02:57 , Processed in 0.174567 second(s), 24 queries , Gzip On.

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