24小时热门版块排行榜     意得辑英文润色,师兄用过都说好

查看: 400  |  回复: 8

voleyes

银虫 (正式写手)

[求助] QE中计算态密度时,SCF计算完后,为什么还要再进行nscf计算?已有1人参与

能否scf计算完直接进行dos计算,nscf计算有什么意义吗?
回复此楼
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

漫天飘雪

版主 (职业作家)

认真做事,踏实做人

你可以一步成的
http://blog.sciencenet.cn/u/pfliu89
2楼2018-03-16 08:12:40
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

lyp_ontheway

金虫 (初入文坛)

【答案】应助回帖

★ ★ ★ ★ ★
感谢参与,应助指数 +1
voleyes: 金币+5, ★★★★★最佳答案 2018-03-17 03:18:11
计算态密度时需要一定的K点密度。为了节省时间,人们一般都会先在较稀疏的K网格下进行scf计算得到合理的电荷密度分布,再在较密的K网格下进行nscf计算得到态密度。当然也可以直接在较密的K网格下进行scf计算得到态密度,这样得到的结果没有什么区别。
3楼2018-03-16 11:14:26
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

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
4楼2018-03-17 03:21:46
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

lyp_ontheway

金虫 (初入文坛)

★ ★
漫天飘雪: 金币+2, 谢谢交流 2018-03-19 17:49:03
automatic就是Monkhost-Pack网格。我做nscf时一般把每个方向的k点都加倍。
5楼2018-03-19 16:51:51
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

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是之前输出过这个找不到吗?应该怎么解决啊?谢谢
6楼2018-06-08 20:06:05
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

sbingyi

金虫 (小有名气)

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

在重复Ni的费米面计算时(Examples-6.1/PP/example02),到态密度这一步,非自洽pw.x &lt;ni.dos.in&gt; 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
7楼2018-06-09 21:30:57
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

liuxinbao

铁虫 (初入文坛)

自洽计算选取稀疏的k点是为了减小计算量,算出相对准确的波函数等信息,而非自洽计算是为了插值得到其他k点的信息以便进行输出,插值计算量比自洽迭代计算量小很多,所以算的比较快,可以选取很多的k点
8楼2018-07-10 10:28:03
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

录过天晴

新虫 (初入文坛)

引用回帖:
7楼: Originally posted by sbingyi at 2018-06-09 21:30:57
我把自洽中k点改为automatic(12 12 12 0 0 0),在非自洽时为automatic(24 24 24 0 0 0),又出现了新的情况,dos.x &lt;ni.dos2.in&gt; ni.dos2.out运行时没有什么反应,但是也没有结果输出,不产生ni. ...

我想问一下你出现这样的问题是怎么解决的,我也遇到了类似的问题,计算dos没有输出,不知道怎么解决。
9楼2018-10-10 08:59:55
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 voleyes 的主题更新
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[论文投稿] 读博交流 +6 Lpab 2018-12-10 10/500 2018-12-10 20:51 by 韩暖宝
[公派出国] CSC公派出国及美签证申请过程 +15 zkunking 2018-12-10 21/1050 2018-12-10 20:28 by 苜蓿情缘
[仿真模拟] 某大型多晶硅企业仿真模拟岗位招聘信息 +6 猎猎头 2018-12-05 9/450 2018-12-10 19:45 by 林枫LKW
[物理] 关于波数和波长的关系 20+3 zhaobozhou 2018-12-05 6/300 2018-12-10 19:36 by zhufengjuan
[论文投稿] 多久给编辑发邮件 +6 材料人ing 2018-12-10 14/700 2018-12-10 17:05 by 429574511
[教师之家] 怎样高中晶胞3D模型制作 +3 静静是我19 2018-12-08 7/350 2018-12-10 15:00 by 卜卜卜卜卜
[高分子] PDMS固化问题,为什么一直是油状的,不能成型 50+3 在我辛中 2018-12-09 6/300 2018-12-10 13:49 by yh323
[找工作] 做个简单的调查,机械类或材料类的大家,工作几年的或现在刚签约的薪酬。 +7 yanwentao 2018-12-08 9/450 2018-12-10 12:54 by yanwentao
[留学DIY] 求定位!普通一本!EE专业均分83.3,paper较多,直接申请英国PHD的成功率大么? +5 鸟语丶 2018-12-08 5/250 2018-12-10 11:27 by sunoldmoon
[公派出国] csc工资延迟 +3 小熊样儿 2018-12-07 9/450 2018-12-10 04:22 by VPMHM
[论文投稿] 唉,不想上了 +52 xiaohedoudou 2018-12-06 94/4700 2018-12-09 23:02 by 南极大冒险
[硕博家园] 虽然你很美,梦也很美,但我应该清醒清醒了 (金币+5) +26 D调年华 2018-12-09 34/1700 2018-12-09 18:14 by abc-404
[基金申请] 标书会送到自己的导师手上审吗 +9 jiamen88 2018-12-04 9/450 2018-12-09 16:46 by 朱平_Hugh
[考博] 兰大考博 +4 草莓酱22 2018-12-07 4/200 2018-12-09 14:39 by fighting2016
[基金申请] 这样,青基能结题吗? +13 89509306 2018-12-07 15/750 2018-12-09 11:14 by 太初问花
[硕博家园] 兰州与北京,该如何选择 +8 陌上枫声 2018-12-08 18/900 2018-12-08 22:29 by 陌上枫声
[公派出国] 2019年有哪位同学想到意大利博洛尼亚大学读博(土木-结构健康监测方向)可以站内联系 +6 chengzb187 2018-12-07 13/650 2018-12-07 21:38 by chengzb187
[出国考试] 我的留学梦啊 +4 aienglish606 2018-12-05 9/450 2018-12-07 12:41 by aienglish606
[硕博家园] 心情极度压抑 +27 想想办法 2018-12-04 33/1650 2018-12-07 11:43 by 观千山暮雪
[硕博家园] 本人外科医生,目前想辞职做区域医学经理 +11 小青年1026 2018-12-04 17/850 2018-12-07 10:43 by davidrobin
信息提示
请填处理意见