计算化学公社

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

[新手求助] 如何生成原子排序一致的分子团簇PDB文件以实现批量电荷分配?

[复制链接 Copy URL]

2

帖子

0

威望

301

eV
积分
303

Level 3 能力者

本帖最后由 化学兔 于 2025-4-13 21:10 编辑

老师们好,我在使用分子力场计算外重组能时,给团簇的所有分子设定原子电荷需要使每个分子的原子排序保持一致(例如图1,indigo分子团簇中原子排序每30一循环,来自帖子第二部分:http://sobereva.com/332)。我将molecule.cif运用Multiwfn产生中心分子+相邻分子的团簇结构并根据键连关系重排原子序号,但是发现问题:在导出的 molecule-cluster-multiwfn.pdb 文件中,如图2所示,前52个原子排序和之后的52个原子排序不能保持一致(1个分子包含52原子);我又尝试使用Materials Studio打开ciif文件,导出原子团簇并使用Multiwfn重排原子序号,发现 molecule-cluster-ms-multiwfn.pdb 文件中同种分子的原子顺序仍不一致(图3)。由于团簇pdb文件中原子数很多(3900个原子)且排序不一致,我无法给原子1~52、53~104...3849~3900都依次赋上单分子的原子电荷,请问老师有什么好的办法可以生成原子排序保持一致的坐标文件呢?谢谢老师
图1-图3见下图:

molecule.cif

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

我的分子

indigo-cluster.pdb

353.98 KB, 下载次数 Times of downloads: 0

http://sobereva.com/332,参考帖子示例

molecule-cluster-multiwfn.pdb

69.12 KB, 下载次数 Times of downloads: 0

Multiwfn产生中心分子+相邻分子的团簇结构(884个原子)

molecule-cluster-ms-multiwfn.pdb

304.75 KB, 下载次数 Times of downloads: 0

6万

帖子

99

威望

5万

eV
积分
120137

管理员

公社社长

2#
发表于 Post on 2025-4-14 02:00:08 | 只看该作者 Only view this author
简单方法没有
原理上,可以写程序自动依次处理团簇里每个分子,对这个分子和参考分子做RMSD叠合确定原子对应关系并重新排序,再写入到原团簇文件里
北京科音自然科学研究中心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

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

GMT+8, 2025-8-15 17:44 , Processed in 0.162882 second(s), 24 queries , Gzip On.

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