计算化学公社

 找回密码 Forget password
 注册 Register
楼主 Author: 不会扣篮的后卫
打印 Print 上一主题 Last thread 下一主题 Next thread

[CP2K] 求助,在CP2K中使用DFTB方法时,出现了报错

[复制链接 Copy URL]

282

帖子

0

威望

2873

eV
积分
3155

Level 5 (御坂)

16#
 楼主 Author| 发表于 Post on 2022-2-24 07:08:00 | 只看该作者 Only view this author
sobereva 发表于 2022-2-12 03:07
建议先考虑直接支持的GFN1-xTB,省得折腾参数

卢老师,您好。我在使用GFN1-xTB方法计算结构的能量时(该结构由 6个chemical unit的CuFeO2组成),
出现了以下的warning,
Cu和O原子的信息都被忽略了。
请问卢老师,这个warning该怎么解决呢?
(附上相关的输入和输出文件)


*** WARNING in qs_kind_types.F:2043 :: Information provided in the input ***
*** file regarding POTENTIAL for KIND <Cu> will be ignored!              ***


*** WARNING in qs_kind_types.F:2049 :: Information provided in the input ***
*** file regarding BASIS for KIND <Cu> will be ignored!                  ***


*** WARNING in qs_kind_types.F:2043 :: Information provided in the input ***
*** file regarding POTENTIAL for KIND <O> will be ignored!               ***


*** WARNING in qs_kind_types.F:2049 :: Information provided in the input ***
*** file regarding BASIS for KIND <O> will be ignored!                   ***





files.zip

81.06 KB, 下载次数 Times of downloads: 4

6万

帖子

99

威望

6万

eV
积分
125131

管理员

公社社长

17#
发表于 Post on 2022-2-24 11:27:48 | 只看该作者 Only view this author
不会扣篮的后卫 发表于 2022-2-24 07:08
卢老师,您好。我在使用GFN1-xTB方法计算结构的能量时(该结构由 6个chemical unit的CuFeO2组成),
出 ...

不用管
GFN-xTB方法自身定义了基函数
北京科音自然科学研究中心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

282

帖子

0

威望

2873

eV
积分
3155

Level 5 (御坂)

18#
 楼主 Author| 发表于 Post on 2022-2-24 13:17:04 | 只看该作者 Only view this author
sobereva 发表于 2022-2-24 11:27
不用管
GFN-xTB方法自身定义了基函数

好的,谢谢卢老师!

282

帖子

0

威望

2873

eV
积分
3155

Level 5 (御坂)

19#
 楼主 Author| 发表于 Post on 2022-3-4 14:06:06 | 只看该作者 Only view this author
sobereva 发表于 2022-2-24 11:27
不用管
GFN-xTB方法自身定义了基函数

请教卢老师,我在用CP2K跑6个chemical unit CuFeO2组成的优化结构输入文件里,选择用GFN1-xTB方法run_type:cell_opt,自旋多重度分别设为1(3个Fe1 alpha 电子数为5,3个Fe2 beta 电子数为5)
31(6个Fe1 alpha 电子数为 5).
在跑自旋多重度分别设为1的输入文件时,
得到理想的输出结果,优化结构成功收敛,

但在跑自旋多重度分别设为31的输入文件时,
优化结构没有成功收敛。

尤其是压力变的很不稳定。
第一个问题,请问卢老师为什么在高自旋的情况下,压力变化这么大?有什么输入文件的参数设置,可以使压力稳定和优化结构完成收敛吗?
(我尝试了选择run_type:Optimizing structure(cell is fixed)跑自旋多重度分别设为1(3个Fe1 alpha 电子数为5,3个Fe2 beta 电子数为5)
31(6个Fe1 alpha 电子数为 5)的输入文件 ,这样可以避免产生关于压力
Pressure Deviation [bar]   =     38288.6656706373
  Pressure Tolerance [bar]   =       100.0000000000
  Conv. for  PRESSURE        =                   NO       的信息报告,但是在run_type:Optimizing structure(cell is fixed)这种情况下,自旋多重度分别设为131都没有完成收敛。)


第二个问题

请问卢老师,这个&OUTER_SCF参数设置有什么作用吗?
(附上相关的文件)

All_files.zip

4.81 MB, 下载次数 Times of downloads: 0

由于文件尺寸超额,先不上传Fix_cell type的输入输出文件了

210

帖子

3

威望

2110

eV
积分
2380

Level 5 (御坂)

20#
发表于 Post on 2022-3-4 14:23:39 | 只看该作者 Only view this author
不会扣篮的后卫 发表于 2022-3-4 14:06
请教卢老师,我在用CP2K跑6个chemical unit CuFeO2组成的优化结构输入文件里,选择用GFN1-xTB方法[/backc ...

你使用OT方法时,普通SCF迭代到步数上限未收敛会进入outer SCF迭代,这个是outer SCF最大步数和收敛限设置的标准

6万

帖子

99

威望

6万

eV
积分
125131

管理员

公社社长

21#
发表于 Post on 2022-3-4 18:50:05 | 只看该作者 Only view this author
不会扣篮的后卫 发表于 2022-3-4 14:06
请教卢老师,我在用CP2K跑6个chemical unit CuFeO2组成的优化结构输入文件里,选择用GFN1-xTB方法[/backc ...

你这个体系也不大,完全没必要用GFN-xTB
本来GFN-xTB是对于分子体系训练的参数,用于这种过渡金属氧化物,又牵扯不同自旋的复杂情况,可靠度很有限。
北京科音自然科学研究中心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

282

帖子

0

威望

2873

eV
积分
3155

Level 5 (御坂)

22#
 楼主 Author| 发表于 Post on 2022-3-5 03:41:08 | 只看该作者 Only view this author
Eudaimonia 发表于 2022-3-4 14:23
你使用OT方法时,普通SCF迭代到步数上限未收敛会进入outer SCF迭代,这个是outer SCF最大步数和收敛限设 ...

谢谢老师!

282

帖子

0

威望

2873

eV
积分
3155

Level 5 (御坂)

23#
 楼主 Author| 发表于 Post on 2022-3-5 03:43:41 | 只看该作者 Only view this author
本帖最后由 不会扣篮的后卫 于 2022-3-5 03:44 编辑
sobereva 发表于 2022-3-4 18:50
你这个体系也不大,完全没必要用GFN-xTB
本来GFN-xTB是对于分子体系训练的参数,用于这种过渡金属氧化物 ...

谢谢卢老师!我尝试一下其他方法去。

173

帖子

0

威望

1800

eV
积分
1973

Level 5 (御坂)

24#
发表于 Post on 2022-5-21 19:33:21 | 只看该作者 Only view this author
你好,我的体系里面有FE
出现报错: Failure in assigning KINDS <FE> and <FE> to a DFTB interaction pair!

想请教一下,您是怎么解决的这个问题。
独立之精神 自由之思想

本版积分规则 Credits rule

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

GMT+8, 2026-2-20 06:42 , Processed in 0.192670 second(s), 22 queries , Gzip On.

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