应《网络安全法》要求，自2017年10月1日起，未进行实名认证将不得使用互联网跟帖服务。为保障您的帐号能够正常使用，请尽快对帐号进行手机号验证，感谢您的理解与支持！

24小时热门版块排行榜    

小木虫论坛-学术科研互动平台 » 计算模拟区 » 第一性原理 » QE(Pwscf) » QE中计算态密度时，SCF计算完后，为什么还要再进行nscf计算？
91/1返回列表
查看: 5225  |  回复: 8
只看楼主 @他人 存档 新回复提醒 (忽略) 收藏    在APP中查看

voleyes

银虫 (正式写手)

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

能否scf计算完直接进行dos计算，nscf计算有什么意义吗？
回复此楼

» 猜你喜欢
1楼 2018-03-16 04:11:23
已阅   回复此楼   关注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计算得到态密度，这样得到的结果没有什么区别。
赞一下(4人) 回复此楼
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 的主题更新
91/1返回列表
普通表情 龙 兔 虎 猫
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 一志愿华中科技大学，080502，354分求调剂 +4 守候夕阳CF 2026-03-18 4/200 2026-03-18 22:16 by li123456789.
[考研] 一志愿天津大学化学工艺专业（081702）315分求调剂 +10 yangfz 2026-03-17 10/500 2026-03-18 20:14 by walc
[考研] 【同济软件】软件（085405）考研求调剂 +3 2026eternal 2026-03-18 3/150 2026-03-18 19:09 by 搏击518
[考研] 314求调剂 +8 无懈可击的巨人 2026-03-12 8/400 2026-03-18 14:50 by haxia
[考研] 一志愿西南交大，求调剂 +4 材化逐梦人 2026-03-18 4/200 2026-03-18 14:22 by 007_lilei
[考研] 070300化学319求调剂 +6 锦鲤0909 2026-03-17 6/300 2026-03-18 13:22 by Iveryant
[考研] 302求调剂 +10 呼呼呼。。。。 2026-03-17 10/500 2026-03-18 12:45 by Linda Hu
[考研] 304求调剂 +12 小熊joy 2026-03-14 13/650 2026-03-18 12:34 by Linda Hu
[考研] 0854，计算机类招收调剂 +3 胡辣汤放糖 2026-03-15 6/300 2026-03-18 12:09 by 上岸上岸……..
[考研] 301求调剂 +9 yy要上岸呀 2026-03-17 9/450 2026-03-18 08:58 by 无际的草原
[基金申请] 被我言中：新模板不强调格式了，假专家开始管格式了 +4 beefly 2026-03-14 4/200 2026-03-17 22:04 by 黄鸟于飞Chao
[考研] 085600材料与化工求调剂 +5 绪幸与子 2026-03-17 5/250 2026-03-17 16:40 by laoshidan
[考研] 一志愿，福州大学材料专硕339分求调剂 +3 木子momo青争 2026-03-15 3/150 2026-03-17 07:52 by laoshidan
[考研] 304求调剂 +4 ahbd 2026-03-14 4/200 2026-03-16 16:48 by 我的船我的海
[考研] 中科大材料与化工319求调剂 +3 孟鑫材料 2026-03-14 3/150 2026-03-14 20:10 by ms629
[考研] 297一志愿上交085600求调剂 +5 指尖八千里 2026-03-14 5/250 2026-03-14 17:26 by a不易
[考研] 265求调剂 +4 威化饼07 2026-03-12 4/200 2026-03-14 17:23 by userper
[考研] 290求调剂 +7 ADT 2026-03-12 7/350 2026-03-13 15:17 by JourneyLucky
[考研] 土木第一志愿276求调剂，科研和技能十分丰富，求新兴方向的导师收留 +3 土木小天才 2026-03-12 3/150 2026-03-13 15:01 by JourneyLucky
[考研] 290求调剂 +3 ADT 2026-03-13 3/150 2026-03-13 10:19 by peike
信息提示
关闭
请填处理意见
关闭