当前位置: 首页 > 第一原理 >QE中计算态密度时,SCF计算完后,为什么还要再进行nscf计算?

QE中计算态密度时,SCF计算完后,为什么还要再进行nscf计算?

作者 voleyes
来源: 小木虫 400 8 举报帖子
+关注

能否scf计算完直接进行dos计算,nscf计算有什么意义吗? 返回小木虫查看更多

今日热帖
  • 精华评论
  • 漫天飘雪

    你可以一步成的

  • lyp_ontheway

    计算态密度时需要一定的K点密度。为了节省时间,人们一般都会先在较稀疏的K网格下进行scf计算得到合理的电荷密度分布,再在较密的K网格下进行nscf计算得到态密度。当然也可以直接在较密的K网格下进行scf计算得到态密度,这样得到的结果没有什么区别。

  • voleyes

    引用回帖:
    3楼: Originally posted by lyp_ontheway at 2018-03-16 11:14:26
    计算态密度时需要一定的K点密度。为了节省时间,人们一般都会先在较稀疏的K网格下进行scf计算得到合理的电荷密度分布,再在较密的K网格下进行nscf计算得到态密度。当然也可以直接在较密的K网格下进行scf计算得到态密 ...

    谢谢!较密的k网格一般都设为多少呢?
    比如我这个cell是2x2的,
    scf计算k点是这样设置的
    K_POINTS AUTOMATIC
    8 8 1 0 0 0
    如果nscf计算应该设置为多少比较合适呢?
    另外,也是可以选择自动生成吗AUTOMATIC

  • lyp_ontheway

    automatic就是Monkhost-Pack网格。我做nscf时一般把每个方向的k点都加倍。

  • sbingyi

    我把我的求助帖搬过来了,麻烦帮我看一下。

    在重复Ni的费米面计算时(Examples-6.1/PP/example02),到态密度这一步,非自洽pw.x <ni.dos.in> ni.dos.out结束后,进行总态密度计算时,输入dos.x <ni.dos2,in> ni.dos2.out时提示:
    Program DOS v.6.1 (svn rev. 13369) starts on  7Jun2018 at  4:35:17

         This program is part of the open-source Quantum ESPRESSO suite
         for quantum simulation of materials; please cite
             "P. Giannozzi et al., J. Phys.:Condens. Matter 21 395502 (2009);
              URL http://www.quantum-espresso.org",
         in publications or presentations arising from this work. More details at
         http://www.quantum-espresso.org/quote

         Parallel version (MPI), running on     1 processors

       Info: using nr1, nr2, nr3 values from input

       Info: using nr1, nr2, nr3 values from input
    ni.dos.in:
    &control
        calculation='nscf'
        prefix='ni',
        pseudo_dir = './',
        outdir='./tmp'
    /
    &system
        ibrav=2, celldm(1) =6.48, nat=1, ntyp=1,
        nspin = 2,  starting_magnetization(1)=0.7,
        ecutwfc = 24.0, ecutrho = 288.0, nbnd=8,
        occupations='tetrahedra'
    /
    &electrons
        conv_thr = 1.0e-10
        mixing_beta = 0.7
    /
    ATOMIC_SPECIES
    Ni 58.69 Ni.pz-n-rrkjus_psl.1.0.0.UPF
    ATOMIC_POSITIONS
    Ni 0.0 0.0 0.0
    K_POINTS {automatic}
    12 12 12 0 0 0
    ni.dos2,in:
    &dos
       outdir='./tmp'
       prefix='ni',
       fildos='ni.dos',
       Emin=5.0, Emax=25.0, DeltaE=0.1
    /
    请问nr1, nr2, nr3是之前输出过这个找不到吗?应该怎么解决啊?谢谢,

  • sbingyi

    引用回帖:
    6楼: Originally posted by sbingyi at 2018-06-08 20:06:05
    我把我的求助帖搬过来了,麻烦帮我看一下。

    在重复Ni的费米面计算时(Examples-6.1/PP/example02),到态密度这一步,非自洽pw.x <ni.dos.in> ni.dos.out结束后,进行总态密度计算时,输入dos.x & ...

    我把自洽中k点改为automatic(12 12 12 0 0 0),在非自洽时为automatic(24 24 24 0 0 0),又出现了新的情况,dos.x <ni.dos2.in> ni.dos2.out运行时没有什么反应,但是也没有结果输出,不产生ni.dos

  • liuxinbao

    自洽计算选取稀疏的k点是为了减小计算量,算出相对准确的波函数等信息,而非自洽计算是为了插值得到其他k点的信息以便进行输出,插值计算量比自洽迭代计算量小很多,所以算的比较快,可以选取很多的k点

猜你喜欢
下载小木虫APP
与700万科研达人随时交流
  • 二维码
  • IOS
  • 安卓