计算化学公社

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

[GROMACS] 求助rdf曲线的波动和收敛的问题

[复制链接 Copy URL]

141

帖子

0

威望

507

eV
积分
648

Level 4 (黑子)

跳转到指定楼层 Go to specific reply
楼主
请问一下各位老师,建立了一个含有Zn离子和两种混合溶剂的5*5*5的溶剂盒子,做完分子模拟后得到的rdf曲线不是尖锐的峰并且在2nm以内的范围还有较大的波动,没有收敛到1.出现这种情况是正常的吗?是不是不能直接计算分子和分子间的rdf,应该选择分子中特定的配位原子计算呢?

202302230127596659..png (55.4 KB, 下载次数 Times of downloads: 7)

202302230127596659..png

6万

帖子

99

威望

5万

eV
积分
120132

管理员

公社社长

2#
发表于 Post on 2023-2-23 04:39:41 | 只看该作者 Only view this author
什么峰型和具体体系、计算的对象有关,没法一概而论。要从体系特征(也结合轨迹图考察)、计算设置、rdf计算原理上理解结果

明显已经往1收敛了,再增加统计距离自然就会基本精确收敛到1

“分子和分子间”没有确切定义。也不存在没法计算,诸如GROMACS里可以计算质心之间、几何中心之间的rdf等等
北京科音自然科学研究中心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

141

帖子

0

威望

507

eV
积分
648

Level 4 (黑子)

3#
 楼主 Author| 发表于 Post on 2023-2-23 15:57:16 | 只看该作者 Only view this author
sobereva 发表于 2023-2-23 04:39
什么峰型和具体体系、计算的对象有关,没法一概而论。要从体系特征(也结合轨迹图考察)、计算设置、rdf计 ...

谢谢sob老师,还想请问您一下,如果设置的rdf曲线是计算针对Zn和resname ACN的话,他是会默认计算的是质心间的距离还是最近原子的距离啊?如果是最近原子的话,要怎么设置计算质心位置的距离呢?

6万

帖子

99

威望

5万

eV
积分
120132

管理员

公社社长

4#
发表于 Post on 2023-2-24 05:54:34 | 只看该作者 Only view this author
sxh3231796 发表于 2023-2-23 15:57
谢谢sob老师,还想请问您一下,如果设置的rdf曲线是计算针对Zn和resname ACN的话,他是会默认计算的是质 ...

是Zn与ACN各个原子间的结果取平均

参考我的GROMACS培训班里的例子


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

141

帖子

0

威望

507

eV
积分
648

Level 4 (黑子)

5#
 楼主 Author| 发表于 Post on 2023-2-24 19:33:05 | 只看该作者 Only view this author
sobereva 发表于 2023-2-24 05:54
是Zn与ACN各个原子间的结果取平均

参考我的GROMACS培训班里的例子

谢谢sob老师

141

帖子

0

威望

507

eV
积分
648

Level 4 (黑子)

6#
 楼主 Author| 发表于 Post on 2023-6-25 19:51:26 | 只看该作者 Only view this author
sobereva 发表于 2023-2-24 05:54
是Zn与ACN各个原子间的结果取平均

参考我的GROMACS培训班里的例子

想在问一下sob老师,在计算溶剂盒子的md模拟时,算出来的rdf曲线所有有机分子到Zn离子的距离都大于0.6nm,这种结果是合理的吗?我附上了我的top文件,请问是计算的设定有什么问题吗?

rdf.docx

183.16 KB, 下载次数 Times of downloads: 2

box.top

512 Bytes, 下载次数 Times of downloads: 4

6万

帖子

99

威望

5万

eV
积分
120132

管理员

公社社长

7#
发表于 Post on 2023-6-26 04:06:30 | 只看该作者 Only view this author
sxh3231796 发表于 2023-6-25 19:51
想在问一下sob老师,在计算溶剂盒子的md模拟时,算出来的rdf曲线所有有机分子到Zn离子的距离都大于0.6nm, ...

取决于你怎么统计的
参数合理的话,Zn到有机分与之最近的原子间距离到不了6埃的程度,如果是与小分子质心间的距离倒是正常
北京科音自然科学研究中心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

141

帖子

0

威望

507

eV
积分
648

Level 4 (黑子)

8#
 楼主 Author| 发表于 Post on 2023-6-26 14:58:19 | 只看该作者 Only view this author
sobereva 发表于 2023-6-26 04:06
取决于你怎么统计的
参数合理的话,Zn到有机分与之最近的原子间距离到不了6埃的程度,如果是与小分子质 ...

统计的确实是到小分子质心的距离,但是计算最近的原子也有5埃左右差不多,想问下有可能是分子数量填充太少了的原因吗?我是按照0.5mol/L的密度配比填进5*5*5的盒子里,需要把盒子填满吗?在盒子里已经包含有机溶剂分子的情况下还需要设定其他的溶剂化效应的参数吗?

6万

帖子

99

威望

5万

eV
积分
120132

管理员

公社社长

9#
发表于 Post on 2023-6-26 21:54:33 | 只看该作者 Only view this author
sxh3231796 发表于 2023-6-26 14:58
统计的确实是到小分子质心的距离,但是计算最近的原子也有5埃左右差不多,想问下有可能是分子数量填充太 ...

做NPT就没必要管初始密度。如果你是指构成盒子的比例,设成和实际研究的对象一致。
北京科音自然科学研究中心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

141

帖子

0

威望

507

eV
积分
648

Level 4 (黑子)

10#
 楼主 Author| 发表于 Post on 2023-7-2 19:45:37 | 只看该作者 Only view this author
sobereva 发表于 2023-6-26 21:54
做NPT就没必要管初始密度。如果你是指构成盒子的比例,设成和实际研究的对象一致。

想再麻烦问下sob老师,用oplss/aa力场计算,更改参数得到的结果还是距离在5-6埃左右,但是换了GAFF力场在做计算就变成很标准的2埃配位,请问要怎么解释这种结果啊?是Zn---tms 和乙腈的溶剂体系

6万

帖子

99

威望

5万

eV
积分
120132

管理员

公社社长

11#
发表于 Post on 2023-7-3 01:32:06 | 只看该作者 Only view this author
sxh3231796 发表于 2023-7-2 19:45
想再麻烦问下sob老师,用oplss/aa力场计算,更改参数得到的结果还是距离在5-6埃左右,但是换了GAFF力场在 ...

用GAFF就完了,我没看出有任何理由非得用OPLS-AA。更何况GAFF还兼容Merz的单原子离子力场参数,里面还有专门给重现好配位数而优化的LJ参数。
北京科音自然科学研究中心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

141

帖子

0

威望

507

eV
积分
648

Level 4 (黑子)

12#
 楼主 Author| 发表于 Post on 2023-7-3 14:58:28 | 只看该作者 Only view this author
sobereva 发表于 2023-7-3 01:32
用GAFF就完了,我没看出有任何理由非得用OPLS-AA。更何况GAFF还兼容Merz的单原子离子力场参数,里面还有 ...

谢谢sob老师,因为之前给过合作方opls/aa的结果,有什么可以解释为什么算出来结果会差距这么大的理论支持吗?

6万

帖子

99

威望

5万

eV
积分
120132

管理员

公社社长

13#
发表于 Post on 2023-7-4 02:41:04 | 只看该作者 Only view this author
sxh3231796 发表于 2023-7-3 14:58
谢谢sob老师,因为之前给过合作方opls/aa的结果,有什么可以解释为什么算出来结果会差距这么大的理论支持 ...

着重对比金属离子的LJ参数
北京科音自然科学研究中心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

141

帖子

0

威望

507

eV
积分
648

Level 4 (黑子)

14#
 楼主 Author| 发表于 Post on 2023-7-4 14:55:32 | 只看该作者 Only view this author
sobereva 发表于 2023-7-4 02:41
着重对比金属离子的LJ参数

好的好的,谢谢sob老师

本版积分规则 Credits rule

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

GMT+8, 2025-8-15 01:28 , Processed in 0.173327 second(s), 23 queries , Gzip On.

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