计算化学公社

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

[CP2K] cp2k结构优化是否需要设置温度

[复制链接 Copy URL]

66

帖子

0

威望

1089

eV
积分
1155

Level 4 (黑子)

cp2k结构优化是否需要设置温度,找到的做结构优化的例子为啥都没有设置温度?望指教

5万

帖子

99

威望

5万

eV
积分
112356

管理员

公社社长

2#
发表于 Post on 2021-3-8 23:02:06 | 只看该作者 Only view this author
仔细看
谈谈温度、压力、同位素设定对量子化学计算结果产生的影响
http://sobereva.com/423http://bbs.keinsci.com/thread-10159-1-1.html
北京科音自然科学研究中心http://www.keinsci.com)致力于计算化学的发展和传播,长期开办高质量的各种计算化学类培训:初级量子化学培训班中级量子化学培训班高级量子化学培训班量子化学波函数分析与Multiwfn程序培训班分子动力学与GROMACS培训班CP2K第一性原理计算培训班,内容介绍以及往届资料购买请点击相应链接查看。这些培训是计算化学从零快速入门以及进一步全面系统性提升研究水平的高速路!培训各种常见问题见《北京科音办的培训班FAQ》
欢迎加入“北京科音”微信公众号获取北京科音培训的最新消息、避免错过网上有价值的计算化学文章!
欢迎加入人气非常高、专业性特别强的综合性理论与计算化学交流QQ群“思想家公社QQ群”:1号:18616395,2号:466017436,3号:764390338,搜索群号能搜到哪个说明目前哪个能加,合计9000人。北京科音培训班的学员在群中可申请VIP头衔,提问将得到群主Sobereva的最优先解答。
思想家公社的门口Blog:http://sobereva.com(发布大量原创计算化学相关博文)
Multiwfn主页:http://sobereva.com/multiwfn(十分强大的量子化学波函数分析程序)
ResearchGate:https://www.researchgate.net/profile/Tian_Lu
Money and papers are rubbish, get a real life!

66

帖子

0

威望

1089

eV
积分
1155

Level 4 (黑子)

3#
 楼主 Author| 发表于 Post on 2021-3-9 09:33:10 | 只看该作者 Only view this author
谢谢答复

153

帖子

0

威望

1418

eV
积分
1571

Level 5 (御坂)

4#
发表于 Post on 2023-9-14 11:20:40 | 只看该作者 Only view this author
sobereva 发表于 2021-3-8 23:02
仔细看
谈谈温度、压力、同位素设定对量子化学计算结果产生的影响
http://sobereva.com/423(http://bbs. ...

卢老师好,看到一篇文献中使用vasp优化了晶体结构后讨论与实验晶体结构参数的表格时特意标注了0K和298K,请问老师是否文献这么做是没意义的?因为读了博文后讨论 结构优化 时没有温度这个概念。文献中是认为了结构优化是在0K,晶体结构是在298K。


Lure Over Love

5万

帖子

99

威望

5万

eV
积分
112356

管理员

公社社长

5#
发表于 Post on 2023-9-14 19:08:48 | 只看该作者 Only view this author
Gzh_NJ 发表于 2023-9-14 11:20
卢老师好,看到一篇文献中使用vasp优化了晶体结构后讨论与实验晶体结构参数的表格时特意标注了0K和298K, ...

标注成0K是不当的
北京科音自然科学研究中心http://www.keinsci.com)致力于计算化学的发展和传播,长期开办高质量的各种计算化学类培训:初级量子化学培训班中级量子化学培训班高级量子化学培训班量子化学波函数分析与Multiwfn程序培训班分子动力学与GROMACS培训班CP2K第一性原理计算培训班,内容介绍以及往届资料购买请点击相应链接查看。这些培训是计算化学从零快速入门以及进一步全面系统性提升研究水平的高速路!培训各种常见问题见《北京科音办的培训班FAQ》
欢迎加入“北京科音”微信公众号获取北京科音培训的最新消息、避免错过网上有价值的计算化学文章!
欢迎加入人气非常高、专业性特别强的综合性理论与计算化学交流QQ群“思想家公社QQ群”:1号:18616395,2号:466017436,3号:764390338,搜索群号能搜到哪个说明目前哪个能加,合计9000人。北京科音培训班的学员在群中可申请VIP头衔,提问将得到群主Sobereva的最优先解答。
思想家公社的门口Blog:http://sobereva.com(发布大量原创计算化学相关博文)
Multiwfn主页:http://sobereva.com/multiwfn(十分强大的量子化学波函数分析程序)
ResearchGate:https://www.researchgate.net/profile/Tian_Lu
Money and papers are rubbish, get a real life!

153

帖子

0

威望

1418

eV
积分
1571

Level 5 (御坂)

6#
发表于 Post on 2023-9-16 10:21:42 | 只看该作者 Only view this author
sobereva 发表于 2023-9-14 19:08
标注成0K是不当的

谢谢卢老师解惑
Lure Over Love

103

帖子

0

威望

1198

eV
积分
1301

Level 4 (黑子)

7#
发表于 Post on 2023-9-18 17:28:18 | 只看该作者 Only view this author
sobereva 发表于 2021-3-8 23:02
仔细看
谈谈温度、压力、同位素设定对量子化学计算结果产生的影响
http://sobereva.com/423(http://bbs. ...

看了老师的博文不禁有一个疑惑,既然压力对热力学数据以外的任何量子化学程序计算结果(比如优化的结构、振动频率、电子激发能、NMR、偶极矩等等)都不产生丝毫影响。

那为啥cp2k会在集合优化CELL_OPT下放一个EXTERNAL_PRESSURE的关键词,而且采用Multiwfn生成CELL_OPT优化的输入文件也有该关键词,在计算过程中会把压力当作一个收敛判据。

而且感觉CELL_OPT优化出来的结构比较不合理,晶格容易发生较大的变化。

77

帖子

1

威望

1938

eV
积分
2035

Level 5 (御坂)

8#
发表于 Post on 2023-9-18 18:18:04 | 只看该作者 Only view this author
1138711019 发表于 2023-9-18 17:28
看了老师的博文不禁有一个疑惑,既然压力对热力学数据以外的任何量子化学程序计算结果(比如优化的结构、 ...

注意看sob老师博文的第一部分最后一段,cp2k一般是用来算周期性体系的,EXTERNAL_PRESSURE是在外压下优化。

5万

帖子

99

威望

5万

eV
积分
112356

管理员

公社社长

9#
发表于 Post on 2023-9-18 18:46:29 | 只看该作者 Only view this author
1138711019 发表于 2023-9-18 17:28
看了老师的博文不禁有一个疑惑,既然压力对热力学数据以外的任何量子化学程序计算结果(比如优化的结构、 ...
周期性体系的DFT计算属于第一性原理计算范畴,不叫量子化学计算,下文说了
辨析计算化学中的任务类型和理论方法
http://sobereva.com/680http://bbs.keinsci.com/thread-38982-1-1.html


“CELL_OPT优化出来的结构比较不合理,晶格容易发生较大的变化” 完全没这回事。该怎么变怎么变,初始结构和晶胞参数烂的情况,不做变胞优化才明显不合理
北京科音自然科学研究中心http://www.keinsci.com)致力于计算化学的发展和传播,长期开办高质量的各种计算化学类培训:初级量子化学培训班中级量子化学培训班高级量子化学培训班量子化学波函数分析与Multiwfn程序培训班分子动力学与GROMACS培训班CP2K第一性原理计算培训班,内容介绍以及往届资料购买请点击相应链接查看。这些培训是计算化学从零快速入门以及进一步全面系统性提升研究水平的高速路!培训各种常见问题见《北京科音办的培训班FAQ》
欢迎加入“北京科音”微信公众号获取北京科音培训的最新消息、避免错过网上有价值的计算化学文章!
欢迎加入人气非常高、专业性特别强的综合性理论与计算化学交流QQ群“思想家公社QQ群”:1号:18616395,2号:466017436,3号:764390338,搜索群号能搜到哪个说明目前哪个能加,合计9000人。北京科音培训班的学员在群中可申请VIP头衔,提问将得到群主Sobereva的最优先解答。
思想家公社的门口Blog:http://sobereva.com(发布大量原创计算化学相关博文)
Multiwfn主页:http://sobereva.com/multiwfn(十分强大的量子化学波函数分析程序)
ResearchGate:https://www.researchgate.net/profile/Tian_Lu
Money and papers are rubbish, get a real life!

103

帖子

0

威望

1198

eV
积分
1301

Level 4 (黑子)

10#
发表于 Post on 2023-9-19 10:01:42 | 只看该作者 Only view this author
sobereva 发表于 2023-9-18 18:46
周期性体系的DFT计算属于第一性原理计算范畴,不叫量子化学计算,下文说了
辨析计算化学中的任务类型和理 ...

原来如此,以前一直搞不清楚量子化学计算和量子计算的区别,谢谢老师

CELL_OPT优化不合理是因为在Materials Project上下载了一个CuCl2练手,发现cp2k优化了接近2h,不知道是不是自己哪里设置出了问题,而且得到的晶格常数是a=8.3,同样的计算用vasp就很快就算完了,晶格常数为a=7.3,可能是之前一直用vasp感觉很可靠,也没有设置过压力相关参数,所以就感觉cp2k的计算有点不太对劲。

103

帖子

0

威望

1198

eV
积分
1301

Level 4 (黑子)

11#
发表于 Post on 2023-9-19 10:02:47 | 只看该作者 Only view this author
wsz 发表于 2023-9-18 18:18
注意看sob老师博文的第一部分最后一段,cp2k一般是用来算周期性体系的,EXTERNAL_PRESSURE是在外压下优化 ...

谢谢,我看到了,也看到了sob老师下面发的博文,现在搞得更清楚了~

5万

帖子

99

威望

5万

eV
积分
112356

管理员

公社社长

12#
发表于 Post on 2023-9-20 00:19:48 | 只看该作者 Only view this author
1138711019 发表于 2023-9-19 10:01
原来如此,以前一直搞不清楚量子化学计算和量子计算的区别,谢谢老师

CELL_OPT优化不合理是因为在Mate ...

CP2K做这种体系优化完全没有问题,结果异常肯定是哪里有硬伤,没输入文件别人也无从判断
北京科音自然科学研究中心http://www.keinsci.com)致力于计算化学的发展和传播,长期开办高质量的各种计算化学类培训:初级量子化学培训班中级量子化学培训班高级量子化学培训班量子化学波函数分析与Multiwfn程序培训班分子动力学与GROMACS培训班CP2K第一性原理计算培训班,内容介绍以及往届资料购买请点击相应链接查看。这些培训是计算化学从零快速入门以及进一步全面系统性提升研究水平的高速路!培训各种常见问题见《北京科音办的培训班FAQ》
欢迎加入“北京科音”微信公众号获取北京科音培训的最新消息、避免错过网上有价值的计算化学文章!
欢迎加入人气非常高、专业性特别强的综合性理论与计算化学交流QQ群“思想家公社QQ群”:1号:18616395,2号:466017436,3号:764390338,搜索群号能搜到哪个说明目前哪个能加,合计9000人。北京科音培训班的学员在群中可申请VIP头衔,提问将得到群主Sobereva的最优先解答。
思想家公社的门口Blog:http://sobereva.com(发布大量原创计算化学相关博文)
Multiwfn主页:http://sobereva.com/multiwfn(十分强大的量子化学波函数分析程序)
ResearchGate:https://www.researchgate.net/profile/Tian_Lu
Money and papers are rubbish, get a real life!

103

帖子

0

威望

1198

eV
积分
1301

Level 4 (黑子)

13#
发表于 Post on 2023-9-20 08:49:39 | 只看该作者 Only view this author
sobereva 发表于 2023-9-20 00:19
CP2K做这种体系优化完全没有问题,结果异常肯定是哪里有硬伤,没输入文件别人也无从判断

我是用Multiwfn生成的,附件已上传,请卢老师看看呢,应该是设置哪里有问题。

CuCl2.inp

4.7 KB, 下载次数 Times of downloads: 13

CuCl2-1.restart

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

5万

帖子

99

威望

5万

eV
积分
112356

管理员

公社社长

14#
发表于 Post on 2023-9-20 08:56:32 | 只看该作者 Only view this author
本来就是闭壳层的,设UKS毫无意义
没必要开smearing
PBEsol算晶格常数比PBE强
北京科音自然科学研究中心http://www.keinsci.com)致力于计算化学的发展和传播,长期开办高质量的各种计算化学类培训:初级量子化学培训班中级量子化学培训班高级量子化学培训班量子化学波函数分析与Multiwfn程序培训班分子动力学与GROMACS培训班CP2K第一性原理计算培训班,内容介绍以及往届资料购买请点击相应链接查看。这些培训是计算化学从零快速入门以及进一步全面系统性提升研究水平的高速路!培训各种常见问题见《北京科音办的培训班FAQ》
欢迎加入“北京科音”微信公众号获取北京科音培训的最新消息、避免错过网上有价值的计算化学文章!
欢迎加入人气非常高、专业性特别强的综合性理论与计算化学交流QQ群“思想家公社QQ群”:1号:18616395,2号:466017436,3号:764390338,搜索群号能搜到哪个说明目前哪个能加,合计9000人。北京科音培训班的学员在群中可申请VIP头衔,提问将得到群主Sobereva的最优先解答。
思想家公社的门口Blog:http://sobereva.com(发布大量原创计算化学相关博文)
Multiwfn主页:http://sobereva.com/multiwfn(十分强大的量子化学波函数分析程序)
ResearchGate:https://www.researchgate.net/profile/Tian_Lu
Money and papers are rubbish, get a real life!

103

帖子

0

威望

1198

eV
积分
1301

Level 4 (黑子)

15#
发表于 Post on 2023-9-21 13:22:23 | 只看该作者 Only view this author
sobereva 发表于 2023-9-20 08:56
本来就是闭壳层的,设UKS毫无意义
没必要开smearing
PBEsol算晶格常数比PBE强

卢老师,我把smearing关了,然后也关了UKS,分别计算了一下PBE(a=7.868)和PBEsol(a=6.829)的几何优化,跑了3h也都收敛了,得到得晶格常数相差也比较大,感觉对称性消失了,我是不是应该把KEEP_SYMMETRY 打开啊?我上传了restart文件,麻烦老师看看,谢谢您~

CuCl2-1.restart

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

CuCl2-pbesol-1.restart

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

本版积分规则 Credits rule

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

GMT+8, 2024-11-24 13:51 , Processed in 0.177291 second(s), 23 queries , Gzip On.

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