计算化学公社

标题: QE计算能带需要读取自洽波函数吗？ [打印本页]

作者
Author: 13619327021 时间: 2025-3-26 23:33
标题: QE计算能带需要读取自洽波函数吗？
各位老师好：
在做完scf后，做bands，输入文件的电子部分是否需要加startingwfc=‘file’？默认的是atomic+random，QE会自动读取吗？
就在这个地方：
&electrons
conv_thr = 1e-8
startingwfc=‘file’ ！这个需要加吗？
/

作者
Author: Freeman 时间: 2025-3-27 00:57
本帖最后由 Freeman 于 2025-3-27 00:58 编辑

要startingpot = 'file'而不是startingwfc = 'file'。两者都相当于“读取前一步的波函数”，但有差别，适用于不同的场景。

前者是读取空间格点电子势能，后者是读取轨道。能带计算就是要使用现有的电子势能，获取在这个势能分布下，k-path上的能带，所以要用startingpot。为什么不能用startingwfc呢？因为你前一步做自洽时选取的K点，跟这一步做非自洽时选取的K点肯定是不一样的；但前一步的wfc文件只存了自洽计算的K点上的轨道，就不能用于这一步了。

那么有没有只能用startingwfc而不能用startingpot的情况呢？有的。读取势能，要求前一步的哈密顿必须和这一步匹配。例如，如果前一步没有考虑hubbard U，这一步考虑了，那么两步的哈密顿就不匹配了。此时就只能用startingwfc了，并且，如上所述，要保证两步的K点一样。

我在周期性体系计算领域不算是专家哈。以上只是我自己的经验见解。如有错误，欢迎专家斧正。

作者
Author: 13619327021 时间: 2025-3-27 08:19

Freeman 发表于 2025-3-27 00:57
要startingpot = 'file'而不是startingwfc = 'file'。两者都相当于“读取前一步的波函数”，但有差别，适用 ...

谢谢解惑，我初学QE，看网上教程能带计算这块很多都不写，也是稀里糊涂的。

作者
Author: 丁越 时间: 2025-3-27 09:44
本帖最后由丁越于 2025-3-27 09:52 编辑

这两个关键词的目的不一样。
startingwfc关键词的目的在于产生自洽计算的初猜波函数，当startingwfc=‘atomic’时对应原子轨道叠加，startingwfc=‘file’时对应从外部文件中读取波函数，这在restart_mode = 'restart'的续算中是默认使用的。
startingpot是给nscf和能带计算用的，专门用于获取电子密度信息。当做能带计算时默认使用startingpot=‘file’，所以输入文件中不用单独写明。

作者
Author: 13619327021 时间: 2025-3-27 11:16

丁越发表于 2025-3-27 09:44
这两个关键词的目的不一样。
startingwfc关键词的目的在于产生自洽计算的初猜波函数，当startingwfc=‘ato ...

嗯嗯，明白了，非常谢谢丁老师解惑。

欢迎光临计算化学公社 (http://ccc.keinsci.com/)