计算化学公社

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

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

[复制链接 Copy URL]

66

帖子

0

威望

1089

eV
积分
1155

Level 4 (黑子)

跳转到指定楼层 Go to specific reply
#
cp2k结构优化是否需要设置温度,找到的做结构优化的例子为啥都没有设置温度?望指教

170

帖子

0

威望

835

eV
积分
1005

Level 4 (黑子)

22#
发表于 Post on 2024-5-17 10:43:33 | 只看该作者 Only view this author
sobereva 发表于 2024-5-17 03:41
相对论效应的描述不是基组层面的事,是赝势或理论方法层面的事
CP2K中常用的MOLOPT系列基组搭配的GTH赝 ...

明白了,感谢卢老师解惑。

5万

帖子

99

威望

5万

eV
积分
112356

管理员

公社社长

21#
发表于 Post on 2024-5-17 03:41:19 | 只看该作者 Only view this author
Shana 发表于 2024-5-16 19:39
感谢卢老师的回复。
另外还想问一个问题,对于Pb、Bi这种重原子,采用DZVP,TZVP这种基组能够描述相对论 ...

相对论效应的描述不是基组层面的事,是赝势或理论方法层面的事
CP2K中常用的MOLOPT系列基组搭配的GTH赝势对这些元素体现了标量相对论效应
北京科音自然科学研究中心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!

170

帖子

0

威望

835

eV
积分
1005

Level 4 (黑子)

20#
发表于 Post on 2024-5-16 19:39:13 | 只看该作者 Only view this author
sobereva 发表于 2024-5-16 00:28
跟这个没关系
尝试用更好的基组
当前的基组只是保底水平

感谢卢老师的回复。
另外还想问一个问题,对于Pb、Bi这种重原子,采用DZVP,TZVP这种基组能够描述相对论效应吗?

5万

帖子

99

威望

5万

eV
积分
112356

管理员

公社社长

19#
发表于 Post on 2024-5-16 00:28:35 | 只看该作者 Only view this author
Shana 发表于 2024-5-15 10:24
卢老师,我最近在用cp2k算一个石榴石结构的材料,PBEsol泛函优化的晶格常数是12.341,vasp优化的和实验的 ...

跟这个没关系
尝试用更好的基组
当前的基组只是保底水平
北京科音自然科学研究中心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!

170

帖子

0

威望

835

eV
积分
1005

Level 4 (黑子)

18#
发表于 Post on 2024-5-15 10:24:06 | 只看该作者 Only view this author
sobereva 发表于 2023-9-22 09:01
PBEsol/DZVP-MOLOPT-SR-GTH对CuCl2优化完的结果,和实验晶胞参数吻合很好

卢老师,我最近在用cp2k算一个石榴石结构的材料,PBEsol泛函优化的晶格常数是12.341,vasp优化的和实验的晶格常数都是12.256,在cp2k中Ge的价电子是4,vasp中Ge的价电子是14,会不会是价电子影响的呢?

base2.inp

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

103

帖子

0

威望

1198

eV
积分
1301

Level 4 (黑子)

17#
发表于 Post on 2023-9-22 14:56:06 | 只看该作者 Only view this author
sobereva 发表于 2023-9-22 09:01
PBEsol/DZVP-MOLOPT-SR-GTH对CuCl2优化完的结果,和实验晶胞参数吻合很好

感谢卢老师,我好好研究研究这个测试算例。

5万

帖子

99

威望

5万

eV
积分
112356

管理员

公社社长

16#
发表于 Post on 2023-9-22 09:01:50 | 只看该作者 Only view this author
1138711019 发表于 2023-9-21 13:22
卢老师,我把smearing关了,然后也关了UKS,分别计算了一下PBE(a=7.868)和PBEsol(a=6.829)的几何优化 ...

PBEsol/DZVP-MOLOPT-SR-GTH对CuCl2优化完的结果,和实验晶胞参数吻合很好

CuCl2.inp (4.71 KB, 下载次数 Times of downloads: 11)

CuCl2-1.restart (5.03 KB, 下载次数 Times of downloads: 1)
北京科音自然科学研究中心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!

5万

帖子

99

威望

5万

eV
积分
112356

管理员

公社社长

15#
发表于 Post on 2023-9-22 05:16:56 | 只看该作者 Only view this author
1138711019 发表于 2023-9-21 13:22
卢老师,我把smearing关了,然后也关了UKS,分别计算了一下PBE(a=7.868)和PBEsol(a=6.829)的几何优化 ...

之前看错了,以为你算的是CuCl。CuCl2得用UKS,当顺磁性体系算
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 (黑子)

14#
发表于 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

5万

帖子

99

威望

5万

eV
积分
112356

管理员

公社社长

13#
发表于 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 (黑子)

12#
发表于 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

管理员

公社社长

11#
发表于 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 (黑子)

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

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

103

帖子

0

威望

1198

eV
积分
1301

Level 4 (黑子)

9#
发表于 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的计算有点不太对劲。

5万

帖子

99

威望

5万

eV
积分
112356

管理员

公社社长

8#
发表于 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!

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

GMT+8, 2024-11-24 17:04 , Processed in 0.193726 second(s), 25 queries , Gzip On.

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