计算化学公社

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

[Quantum ESPRESSO] 新版本ONCV*SP生成的赝势与的老版本的有差别吗?

[复制链接 Copy URL]

300

帖子

6

威望

2711

eV
积分
3131

Level 5 (御坂)

本帖最后由 Freeman 于 2024-8-1 17:38 编辑

大家好。


我最近需要用ONCV结合DFT+U算个体系。Hubbard U需要原子轨道信息,可是现有的http://quantum-simulation.org/potentials/sg15_oncv/上的ONCV没有这个内容。查资料发现,那些赝势是老版本ONCV*SP生成的,而老版本ONCV*SP并没有加入原子轨道信息的功能,所以需要用户自己用新版本ONCV*SP生成。于是我下载了现有的sg15_oncv_2015-10-07.tar.gz中的输入文件,原封不动地喂给ONCV*SP让它重新生成一次。新的ONCV文件确实有原子轨道了(在<PP_PSWFC>栏目里),可是我看到了几点可疑之处。

<PP_HEADER>里
1、has_wfc="F"意思是没有原子轨道吗?我能否手动改成T?
2、total_psenergy和rho_cutoff变化很大。
<PP_MESH>下一行
3、size变大了。

请问这些是正常现象吗?

以下是Mn_ONCV_PBE-1.0的例子
Mn_ONCV_PBE-1.0.dat (711 Bytes, 下载次数 Times of downloads: 2)
Mn_ONCV_PBE-1.0.upf (112.7 KB, 下载次数 Times of downloads: 3)

mn.upf (260.77 KB, 下载次数 Times of downloads: 1)

O-N-C-V-P-S-P是违禁词吗?倒数第三个P自己变成*了。

3621

帖子

3

威望

1万

eV
积分
18429

Level 6 (一方通行)

第一原理惨品小作坊

2#
发表于 Post on 2024-8-1 22:38:18 | 只看该作者 Only view this author
1、不能且无用。这个关键字的解释可以看此处。UPF的格式我没具体研究过,但按我理解一般模守恒赝势似乎不太会用到全电子波函数信息,这个关键字很可能给PAW用的。
2、UPF一般会把生成或者转换赝势的输入文件放在PP_INFO里面,你可以看到你的两个输入文件其实就不一样,产生出来的赝势参数应该就不一样。或许你非要用和之前一样的方式产生就把PP_INFO的部分拷贝下来用这个进行生成(不排除以前那套参数有什么坑才导致这套参数被改了)。
3、同2。
4、有人针对这个问题已经产生了现成的UPF可以直接下崽,大可不必自己折腾。

评分 Rate

参与人数
Participants 1
eV +5 收起 理由
Reason
ghifi37 + 5 GJ!

查看全部评分 View all ratings

日常打哑谜&&探寻更多可能。
原理问题不公开讨论,非商业性质讨论欢迎私聊。
本周忙

300

帖子

6

威望

2711

eV
积分
3131

Level 5 (御坂)

3#
 楼主 Author| 发表于 Post on 2024-8-6 20:15:10 | 只看该作者 Only view this author
卡开发发 发表于 2024-8-1 22:38
1、不能且无用。这个关键字的解释可以看此处。UPF的格式我没具体研究过,但按我理解一般模守恒赝势似乎不太 ...

mn_new.upf (260.77 KB, 下载次数 Times of downloads: 0)

这是以Mn_ONCV_PBE-1.0.UPF的开头为输入文件生成的,结果total_psenergy, rho_cutoff和mesh_size还是不一样。

我之前就注意到这个了。另外,http://quantum-simulation.org/potentials/sg15_oncv/说生成过程中有些数据会更改,
In particular, the values of the cutoff radii are modified slightly to match a node of a logarithmic mesh used in oncV*Sp before they are printed. This introduces a small discrepancy that prevents exact reproduction of the potentials. If you wish to reproduce the potentials, make sure to use the actual input files (*.dat) that were used to generate the SG15 potentials.
导致输入文件和输出文件的开头可能会不一样。所以我当时才用这个网站上提供的dat文件。


3621

帖子

3

威望

1万

eV
积分
18429

Level 6 (一方通行)

第一原理惨品小作坊

4#
发表于 Post on 2024-8-17 10:11:21 | 只看该作者 Only view this author
本帖最后由 卡开发发 于 2024-8-17 10:16 编辑
Freeman 发表于 2024-8-6 20:15
这是以Mn_ONCV_PBE-1.0.UPF的开头为输入文件生成的,结果total_psenergy, rho_cutoff和mesh_size还是 ...

那你可能得看一下数据是怎么改动的以及generator的版本是否是一致的,必要的情况只能去翻generator的源码了。
既然有人折腾完(虽然数据不太一样),我个人看法是重新折腾一下必要性不大。
日常打哑谜&&探寻更多可能。
原理问题不公开讨论,非商业性质讨论欢迎私聊。
本周忙

300

帖子

6

威望

2711

eV
积分
3131

Level 5 (御坂)

5#
 楼主 Author| 发表于 Post on 2024-8-17 16:46:47 | 只看该作者 Only view this author
卡开发发 发表于 2024-8-17 10:11
那你可能得看一下数据是怎么改动的以及generator的版本是否是一致的,必要的情况只能去翻generator的源码 ...

确实,我也打算用别人做出来的得了。

本版积分规则 Credits rule

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

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

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