新版本ONCV*SP生成的赝势与的老版本的有差别吗？

Freeman

大家好。


我最近需要用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: 4)

2024-8-1 17:31 上传 Uploaded

点击下载
Click to download
现成的输入文件

Mn_ONCV_PBE-1.0.upf (112.7 KB, 下载次数 Times of downloads: 4)

2024-8-1 17:31 上传 Uploaded

点击下载
Click to download
现成的赝势文件

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

2024-8-1 17:31 上传 Uploaded

点击下载
Click to download
重新生成的赝势文件

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

卡开发发

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

Freeman

卡开发发 发表于 2024-8-1 22:38
1、不能且无用。这个关键字的解释可以看此处。UPF的格式我没具体研究过，但按我理解一般模守恒赝势似乎不太 ...

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

2024-8-6 20:04 上传 Uploaded

点击下载
Click to download

这是以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文件。

卡开发发

Freeman 发表于 2024-8-6 20:15
这是以Mn_ONCV_PBE-1.0.UPF的开头为输入文件生成的，结果total_psenergy, rho_cutoff和mesh_size还是 ...

那你可能得看一下数据是怎么改动的以及generator的版本是否是一致的，必要的情况只能去翻generator的源码了。
既然有人折腾完（虽然数据不太一样），我个人看法是重新折腾一下必要性不大。

Freeman

卡开发发 发表于 2024-8-17 10:11
那你可能得看一下数据是怎么改动的以及generator的版本是否是一致的，必要的情况只能去翻generator的源码 ...

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