计算化学公社

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

[GROMACS] Gromacs可以模拟计算硝酸钠硝酸钾类熔盐体系的热物理性质吗?

[复制链接 Copy URL]

18

帖子

0

威望

73

eV
积分
91

Level 2 能力者

老师好:
       我做的方向为熔盐热物理性质模拟计算,导师让用Gromacs模拟计算硝酸钠、硝酸钾类硝酸盐的二元至多元盐体系的密度、熔点、粘度、热导率等物性。在文献看到研究者都用的是Lammps,想求问一下Gromacs可是实现熔盐物性的计算吗?如果可以有没有详细的操作手册呢?
多谢回复

6万

帖子

99

威望

6万

eV
积分
125155

管理员

公社社长

2#
发表于 Post on 2022-4-30 07:42:28 | 只看该作者 Only view this author
你看文献里模拟这些体系常用、被普遍接受的参数是哪来的,形式是否是GROMACS能支持的,如果能支持,就能弄到GROMACS里面用。GROMACS支持的势函数形式在手册里都说了。原理上GROMACS支持的形式对于模拟这类体系没什么问题。

具体GROMACS怎么用于算这些体系不是几句话能说明白的。只能告诉你硝酸根之类阴离子的拓扑文件可以用sobtop(http://sobereva.com/soft/Sobtop)构建,文献里的参数也都可以往sobtop的力场文件里添加,参考网页里二氧化硅的例子中添加参数的过程。
如果是基本没有GROMACS基础知识的话,应当从头系统性好好学学,免得鼓捣很久也不得要领还可能犯低级错误。可以考虑购买北京科音分子动力学与GROMACS培训班(http://www.keinsci.com/workshop/KGMX_content.html)的往届资料或参加下一届。

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

18

帖子

0

威望

73

eV
积分
91

Level 2 能力者

3#
 楼主 Author| 发表于 Post on 2022-5-2 16:26:03 | 只看该作者 Only view this author
sobereva 发表于 2022-4-30 07:42
你看文献里模拟这些体系常用、被普遍接受的参数是哪来的,形式是否是GROMACS能支持的,如果能支持,就能弄 ...

老师您好,模拟硝酸盐一般用什么力场呢?

6万

帖子

99

威望

6万

eV
积分
125155

管理员

公社社长

4#
发表于 Post on 2022-5-3 11:35:38 | 只看该作者 Only view this author
peixing 发表于 2022-5-2 16:26
老师您好,模拟硝酸盐一般用什么力场呢?

可以尝试GAFF结合RESP2电荷,这对于溶液下的模拟也完全没问题
但对于熔盐的情况,如果去搜文献,可能能搜到更适合的参数
北京科音自然科学研究中心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

18

帖子

0

威望

73

eV
积分
91

Level 2 能力者

5#
 楼主 Author| 发表于 Post on 2022-5-4 16:23:25 | 只看该作者 Only view this author
sobereva 发表于 2022-5-3 11:35
可以尝试GAFF结合RESP2电荷,这对于溶液下的模拟也完全没问题
但对于熔盐的情况,如果去搜文献,可能能 ...

      老师您好,我用GAFF+RESP模拟了NaNO3-KNO3-Ca(NO3)2三元熔盐体系,RESP电荷采用0.78的修正因子,粒子数设置为Na:K:Ca:NO3=100:540:360:1360,计算423.15K密度误差8%,粘度2%,523.15K密度误差9.16%,粘度6%;粒子数整体放大一倍,423.15K下密度误差8.2%,粘度27.7%。
      这个误差正常吗?Gromacs模拟熔盐体系的熔点也是通过计算温度-扩散系数数据来判别吗?
      可以计算热导率和生成焓吗?

6万

帖子

99

威望

6万

eV
积分
125155

管理员

公社社长

6#
发表于 Post on 2022-5-5 09:07:28 | 只看该作者 Only view this author
peixing 发表于 2022-5-4 16:23
老师您好,我用GAFF+RESP模拟了NaNO3-KNO3-Ca(NO3)2三元熔盐体系,RESP电荷采用0.78的修正因子,粒 ...

密度属于比较容易模拟准确的量,当前偏得太大
可以尝试调节NO3的原子电荷。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

18

帖子

0

威望

73

eV
积分
91

Level 2 能力者

7#
 楼主 Author| 发表于 Post on 2022-5-5 09:51:00 | 只看该作者 Only view this author
本帖最后由 peixing 于 2022-5-5 17:28 编辑
sobereva 发表于 2022-5-5 09:07
密度属于比较容易模拟准确的量,当前偏得太大
可以尝试调节NO3的原子电荷。LJ参数也可能有优化的余地。 ...

    老师您好,我给NO3和其他阳离子的原电荷分别乘了0.77、0.78、0.79、0.8系数,并使每个粒子的电荷和为系数值.密度和粘度分别如图片所示,标为黄色,密度偏差为8.5~7%标为绿色,粘度偏差3~50%。密度和粘度值都随电荷系数增大而增大,目前最优系数是0.78,就是密度8%,粘度3%。我曾用0.9的系数计算密度偏差0.7%,但粘度就偏差的离谱了,所以是不是存在当密度偏差小粘度就会偏差大或粘度偏差小密度大的情况存在呢?
    LJ参数该怎么优化呢?有一篇用lammps采用LJ势函数模拟计算太阳盐(NaNO3-KNO3)的文章,说明LJ势是可用的,但我目前不知道怎么调节LJ参数。
    熔盐方面的文章都是基于lammps的,与gromacs相通吗?
    模拟熔点没有别的办法了吗?怎么计算热导率和生成焓呢,计算这两个性质需要调节什么参数?
谢谢老师!

202205051716197059..png (21.39 KB, 下载次数 Times of downloads: 20)

202205051716197059..png

239

帖子

0

威望

5249

eV
积分
5488

Level 6 (一方通行)

8#
发表于 Post on 2022-5-5 18:50:03 | 只看该作者 Only view this author
熔点是很难算的。直接加热固体会得到过热态,熔点会被高估。
相关方法可以看这篇的introduction: https://link.aps.org/doi/10.1103/PhysRevE.91.033310

6万

帖子

99

威望

6万

eV
积分
125155

管理员

公社社长

9#
发表于 Post on 2022-5-5 19:39:11 | 只看该作者 Only view this author
peixing 发表于 2022-5-5 09:51
老师您好,我给NO3和其他阳离子的原电荷分别乘了0.77、0.78、0.79、0.8系数,并使每个粒子的电荷和为 ...

既然文献里都有参数,就直接拿过来用
若是gmx支持的势函数形式,拿什么程序跑原理上都一样
LJ参数的优化是个非线性优化问题,比较折腾,可以参考比如OPC水模型原文的例子
熔点相对比较简单的模拟方法是8L提到的那篇文章里introduction部分的Two-phases approach

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

18

帖子

0

威望

73

eV
积分
91

Level 2 能力者

10#
 楼主 Author| 发表于 Post on 2022-5-5 20:04:07 | 只看该作者 Only view this author
sobereva 发表于 2022-5-5 19:39
既然文献里都有参数,就直接拿过来用
若是gmx支持的势函数形式,拿什么程序跑原理上都一样
LJ参数的优 ...

  老师,文献中的LJ参数如图所示,跟我用的形式不太一样,我用的是amber生成的拓扑文件,怎么修改呢?
   热导率和生成焓该用什么具体的指令来计算呢?在科音的培训里没有发现相关内容。
  谢谢老师!

202205051958245089..png (45.93 KB, 下载次数 Times of downloads: 41)

202205051958245089..png

本版积分规则 Credits rule

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

GMT+8, 2026-2-24 04:33 , Processed in 0.172175 second(s), 23 queries , Gzip On.

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