对pw.x的一些疑问，望解答

bushiyueliang

在网上看到有说“pw.x约定：在同一目录，并保持outdir、prefix一致时，先运行vc-relax、relax计算，在接着的scf、nscf、bands计算会读取之前弛豫后的结构，而忽略此时的结构设置。所以，此时做scf计算，修改的地方是：将calculation=’vc-relax’改成calculation=’scf’，其他部分与上一步输入文件相同。”对此我有以下几点疑问：
1.outdir、prefix一致时我做完自洽再做非自洽是不是就把自洽产生的文件覆盖了？如果我做完vc-relax、scf、nscf，发现我scf计算有问题，要重新修改scf.in，再次提交，覆盖的文件对scf计算会产生影响吗

2.比如说我在一个目录内做了vc-relax计算，但是scf计算我想测试不同参数会产生什么影响，然后我复制了这个做完vc-relax的文件夹产生一个新的文件夹，此时scf计算还会读取弛豫后的结构吗？


3.计算bands时要加入nbnd="",那在scf计算的时候就要加入这个参数吗？这个加的早晚有什么影响？


非常感谢，望您不吝赐教！

bushiyueliang

大佬救救

bushiyueliang

大佬看看这里

万里云

答案：
1. 会覆盖scf的波函数，但不会覆盖电荷密度。如果算到中途发现错误，应该删掉原来的文件重新计算。
2. 建议用脚本手动提取优化后的结构。
3. QE默认不开展宽，所以对于半导体和绝缘体不用加nbnds。对于金属不开展宽很难收敛，单独加 nbnds没有用。

建议：
1. 全自动的vc-relax会遇到严重的Pulay应力问题，最简单粗暴的解决方法就是增大动能截断。先根据文献估计一个大差不离的晶格常数，做scf收敛测试，确定动能截断、展宽和k点。根据收敛测试结果，把动能截断增大到1.2-1.5倍，做变胞优化。2. nscf是很久以前pw.x还不支持calculation="bands"时的功能，现在除了少数读取波函数的第三方程序需要nscf外，大部分场景可以用bands代替。从原理上来讲，做完scf电荷密度就确定了，没必要再搞nscf。

bushiyueliang

万里云 发表于 2024-1-30 11:33
答案：
1. 会覆盖scf的波函数，但不会覆盖电荷密度。如果算到中途发现错误，应该删掉原来的文件重新计算。 ...

谢谢大佬！！！！

bushiyueliang

万里云 发表于 2024-1-30 11:33
答案：
1. 会覆盖scf的波函数，但不会覆盖电荷密度。如果算到中途发现错误，应该删掉原来的文件重新计算。 ...

大佬我还有个问题，我用的模守恒赝势，之前用的80Ry的截断能，按您的意思是vc-relax计算时要加大到100左右吗？然后自洽和非自洽可以再调回80？不知道我的理解对不对

万里云

bushiyueliang 发表于 2024-1-31 14:35
大佬我还有个问题，我用的模守恒赝势，之前用的80Ry的截断能，按您的意思是vc-relax计算时要加大到100左 ...

是这个意思。如果截断大了算得慢，可以用固定体积的变胞，然后拟合E-V曲线把能量最低的V求出来，再在这个V下做一次固定体积的变胞。

bushiyueliang

万里云 发表于 2024-2-1 10:42
是这个意思。如果截断大了算得慢，可以用固定体积的变胞，然后拟合E-V曲线把能量最低的V求出来，再在这个 ...

好的，谢谢谢谢！！！

卡开发发

其他地方写了一遍，这里重新打一遍
1、如果重新提交SCF的话，理论上说restart_mode='from_scratch'应该没太大影响，但建议最好清空一下。我自己做非自洽（包括nscf和band）时候都会把之前的scf备份一下。
2、要在后续scf中用之前驰豫后的结构，得到之前的out文件里面去拷贝坐标和晶胞参数到新的inp当中才行。
3、平面波一般不会求解全部的本征值，nbnd就是控制求解多少本征值的参数，场景大概如下：
（1）绝缘体的情况一般occupations用'fixed'就行，这种情况不牵扯到设置nbnd。
（2）窄带半导体、金属需要加速BZ积分收敛，这时候会引入电子分数占据，需要求解更多本征值，以及磁性体系alpha或beta有一个会排到更高的能级，这时候有可能要适当调大nbnd，虽然默认的nbnd也适当做了放大，如果默认值出现对角化不收敛之类的问题再考虑放大也可以。
（3）计算非自洽（态密度或者能带计算）的时候如果要得到比较高能级的信息，那么可以适当放大一些nbnd。

卡开发发

万里云 发表于 2024-2-1 10:42
是这个意思。如果截断大了算得慢，可以用固定体积的变胞，然后拟合E-V曲线把能量最低的V求出来，再在这个 ...

我觉得其实这么做也不是最稳妥的，因为尽管固定体积，但实际上可能同一个体积的两个轴以上的长度发生不小的变化，其实实际格点间距还是会变的。可能可以这样，用这种方法初猜一个体积，然后在这个基础上再进行变胞优化，但是仍然要保证两个过程都具有各自足够的动能截断（尤其是后者变胞的情况尽量让stress的分量收敛），这样也许还凑合。

bushiyueliang

卡开发发 发表于 2024-2-1 17:01
我觉得其实这么做也不是最稳妥的，因为尽管固定体积，但实际上可能同一个体积的两个轴以上的长度发生不小 ...

谢谢老师！！！