24小时热门版块排行榜    

查看: 3340  |  回复: 16

Prettyswan

木虫 (小有名气)

[交流] 【求助成功】DOSCAR中的信息不全,没有site-projected DOS 已有9人参与

这是我以前算的一个结果,DOSCAR:
引用回帖:
10  10   1   0
  0.1786177E+02  0.3912680E-09  0.3912680E-09  0.1166745E-08  0.5000000E-15
  1.00000000000000005E-004
  CAR
CaFe2As2                                
     20.00000000    -20.00000000 3501      3.96332741      1.00000000
    -20.000  0.0000E+00  0.0000E+00 #总的态密度信息,下同
    -19.989  0.0000E+00  0.0000E+00
    -19.977  0.0000E+00  0.0000E+00
…… ……
…… ……
  19.977  0.0000E+00  0.6800E+02
     19.989  0.0000E+00  0.6800E+02
     20.000  0.0000E+00  0.6800E+02 #总的态密度信息到此为止
     20.00000000    -20.00000000 3501      3.96332741      1.00000000 #分的态密度信息,下同
    -20.000  0.0000E+00  0.0000E+00  0.0000E+00  0.0000E+00  0.0000E+00  0.0000E+00  0.0000E+00  0.0000E+00  0.0000E+00
    -19.989  0.0000E+00  0.0000E+00  0.0000E+00  0.0000E+00  0.0000E+00  0.0000E+00  0.0000E+00  0.0000E+00  0.0000E+00
    -19.977  0.0000E+00  0.0000E+00  0.0000E+00  0.0000E+00  0.0000E+00  0.0000E+00  0.0000E+00  0.0000E+00  0.0000E+00
…… ……
…… ……

这是这次的结果,就只有开始的那三列总的,下面的那些分的态密度就没有。这是INCAR:
引用回帖:
System = CaFe2As2 ThCr2Si2 DOS
PREC = Accurate
ISPIN = 1
ISMEAR = -5 #for DOS,tethraedron method with Blochl corrections
SGIMA = 0.02
ISTART = 0
EDIFF = 10E-5
ICHARG = 2 #SC-calculation
#ISIF = 2 #stress and relaxation
NSW = 0 #for DOS, the steps of ions relaxtions
IBRION = -1 #for DOS, ions no update, not moved

LREAL = Auto #for supercell or largecell
LORBIT = 11
EMIN = -10 #energy-range for DOS
EMAX = 15
NEDOS = 3501 #number of grid kpoints DOS
ENLMDL = -5 #no update charge for 3steps
NELMIN = 2 #min of two steps
NELM = 60

LPLANE = .TRUE. #parallelisation for linux cluster
NPAR = 4 #number of nodes
LSCALU = .FALSE.
NSIM = 4

请问大家支招。

[ Last edited by zzy870720z on 2010-6-5 at 13:11 ]
回复此楼

» 收录本帖的淘帖专辑推荐

VASP 第一性

» 猜你喜欢

» 本主题相关价值贴推荐,对您同样有帮助:

已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
回帖支持 ( 显示支持度最高的前 50 名 )

Prettyswan

木虫 (小有名气)

引用回帖:
Originally posted by fanchen021 at 2010-06-03 23:04:34:
NPAR = 1 试试呢

这个NPAR是并行计算的控制参数,在不设置NPAR时,默认为1,我开始算的时候没有设置NAPR不过结果也是一样的,说明书上的建议是对于linux cluster来说是NPAR=number of nodes;而且以前有人这么说过,我也试过,不过确实不行。
6楼2010-06-04 10:05:04
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
普通回帖

fanchen021

铁杆木虫 (正式写手)


cenwanglai(金币+1):呵呵 2010-06-03 23:39:11
Prettyswan(金币+3): 2010-06-05 09:48:57
NPAR = 1 试试呢
2楼2010-06-03 23:04:34
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

锐利的碎片

木虫 (正式写手)

star watcher


cenwanglai(金币+1):RWIGS好象也不是非常必须的.因为对于含多种原子以及有真空层的结构,RWIGS的物理意义好象并不是很有意义. 2010-06-03 23:40:52
RWIGS
3楼2010-06-03 23:35:43
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

liu-jw09

铜虫 (初入文坛)

zzy870720z(金币+1):xiexie 2010-06-04 00:08:20
zzy870720z(金币-1):不好意思,昨天明明不是这个解释?为什么要修改呢?呵呵 2010-06-04 10:28:11
不懂,帮顶

[ Last edited by liu-jw09 on 2010-6-4 at 00:04 ]
4楼2010-06-04 00:01:28
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

Prettyswan

木虫 (小有名气)

引用回帖:
Originally posted by liu-jw09 at 2010-06-04 00:01:28:
不懂,帮顶

[ Last edited by liu-jw09 on 2010-6-4 at 00:04 ]

在LORBIT=11时,RWIGS是不被读取的,所以设置也没有用。除了LORBIT=5时不往DOSCAR中写信息,其他时候都是写的。
引用回帖:
VASP.4.6 behaviour:
integer RWIGS line in INCAR  files written   
0         line required             DOSCAR and PROCAR file   
1         line required             DOSCAR and lm decomposed PROCAR file   
2         line required             DOSCAR and lm decomposed PROCAR file + phase factors   
5         line required             PROOUT file   
10       not read                  DOSCAR and PROCAR file   
11       not read                  DOSCAR and lm decomposed PROCAR file   
12       not read                  DOSCAR and lm decomposed PROCAR file + phase factors   

这是我贴的那个完整的DOSCAR的INCAR:
引用回帖:
System = CaFe2As2
PREC = accurate
ISPIN = 1
ISMEAR = -5
SGIMA = 0.2
ISTART = 0
EDIFF = 10E-4
ICHARG = 2
#ISIF = 2
#NSW = 0
#IBRION = 2
NELMDL = -6
LORBIT = 11
EMIN = -20
EMAX = 20
NEDOS = 3501

[ Last edited by Prettyswan on 2010-6-4 at 10:15 ]
5楼2010-06-04 10:00:26
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

Prettyswan

木虫 (小有名气)

引用回帖:
Originally posted by 锐利的碎片 at 2010-06-03 23:35:43:
RWIGS

在LORBIT=11的时候,RWIGS是不被读取的,VASP会在POTCAR中寻找相关参数,只有当LORBIT为0、1、2时才被读取。4楼是我在第一次计算时的INCAR的信息,没有问题。
但是我在以前的计算中发现,如果当LORBIT为1和2时,使用doscar.f编译的程序读取分态密度时候会提示PROCAR中的数据错误;但是当11和12时,就不会有这个问题。

[ Last edited by Prettyswan on 2010-6-4 at 10:11 ]
7楼2010-06-04 10:09:25
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

liu-jw09

铜虫 (初入文坛)

昨天看错了,我以为是没有设置NEDOS的问题呢,可是发现楼主的INCAR里是设置了的,本想删帖的,可是不知道怎么删,就只有修改了
引用回帖:
Originally posted by liu-jw09 at 2010-06-04 00:01:28:
不懂,帮顶

[ Last edited by liu-jw09 on 2010-6-4 at 00:04 ]

8楼2010-06-04 12:44:59
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

gleerat

木虫 (正式写手)

引用回帖:
Originally posted by liu-jw09 at 2010-06-04 12:44:59:
昨天看错了,我以为是没有设置NEDOS的问题呢,可是发现楼主的INCAR里是设置了的,本想删帖的,可是不知道怎么删,就只有修改了

如果没有设置的话,会使用默认值NEDOS=300的。
明察、慎思、笃行
9楼2010-06-04 13:44:57
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

gleerat

木虫 (正式写手)

Prettyswan(金币+7): 2010-06-05 09:49:19
引用回帖:
The parallel version has some restrictions: The site projected DOS is not evaluated in the parallel version in the following cases:
VASP.4.5, NPAR 6 =1 no site projected DOS
VASP.4.6, NPAR 6 =1, LORBIT=0-5 no site projected DOS

引用回帖:
The wave function character is calculated, either by projecting the wavefunctions onto spherical harmonics that are non zero within spheres of a radius RWIGS around each ion (LORIT=1, 2), or using a quick projection scheme relying that works only for the PAW method (LORBIT=10,11,12, see below).

引用回帖:
If the projector augmented wave method is used, LORBIT can also be set to 10, 11 or 12. This alternative setting selects a quick method for the determination of the spd- and site projected wave function character and does not require the specification of a Wigner-Seitz radius in the INCAR file (the RWIGS line is neglected in this case). The method works only for PAW POTCAR files and not for ultrasoft or norm conserving pseudopotentials.

必须要NPAR=1,LORBIT=11或者12,必须要使用PAW方法,采用PAW势。试着改动一下。

[ Last edited by gleerat on 2010-6-4 at 18:32 ]
明察、慎思、笃行
10楼2010-06-04 18:28:01
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 Prettyswan 的主题更新
普通表情 高级回复 (可上传附件)
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考博] 东华理工大学化材专业26届硕士博士申请 +6 zlingli 2026-03-13 6/300 2026-03-15 20:00 by ryzcf
[考研] 290求调剂 +3 孔志浩 2026-03-12 8/400 2026-03-15 15:30 by 孔志浩
[考研] 080500,材料学硕302分求调剂学校 +4 初识可乐 2026-03-14 5/250 2026-03-14 21:08 by peike
[考研] 289求调剂 +4 这么名字咋样 2026-03-14 6/300 2026-03-14 18:58 by userper
[基金申请] 面上和青基一样限30页不合理 +5 wowsunflower 2026-03-10 7/350 2026-03-14 17:21 by kingkocxr
[考研] 【0703化学调剂】-一志愿华中师范大学-六级475 +5 Becho359 2026-03-11 5/250 2026-03-14 11:35 by 哦哦123
[考研] 308 085701 四六级已过求调剂 +7 温乔乔乔乔 2026-03-12 14/700 2026-03-14 10:49 by JourneyLucky
[考研] 332分材料工程调剂 +3 莓好时光海苔 2026-03-09 3/150 2026-03-14 02:03 by JourneyLucky
[考研] 一志愿北京化工大学材料与化工296分求调剂 +16 稻妻小编 2026-03-09 18/900 2026-03-14 02:00 by JourneyLucky
[考研] 328,0703考生求调剂,一志愿为东北师范大学 +4 观素律 2026-03-09 5/250 2026-03-14 01:24 by JourneyLucky
[基金申请] 有必要更换申报口吗 20+3 fannyamoy 2026-03-11 3/150 2026-03-14 00:52 by zhanghaozhu
[考研] 312求调剂 +6 陌宸希 2026-03-10 6/300 2026-03-14 00:40 by JourneyLucky
[考研] 285化工学硕求调剂(081700) +6 柴郡猫_ 2026-03-12 6/300 2026-03-13 20:46 by hmn_wj
[考研] 0703化学求调剂 +7 绿豆芹菜汤 2026-03-12 7/350 2026-03-13 17:25 by njzyff
[考研] 307求调剂 +5 超级伊昂大王 2026-03-12 5/250 2026-03-13 15:56 by 棒棒球手
[考研] 一志愿山大07化学 332分 四六级已过 本科山东双非 求调剂! +3 不想理你 2026-03-12 3/150 2026-03-13 14:18 by JourneyLucky
[论文投稿] 投稿问题 5+4 星光灿烂xt 2026-03-12 6/300 2026-03-13 14:17 by god_tian
[考博] 读博申请 +5 感dd 2026-03-10 7/350 2026-03-11 17:02 by QGZDSYS
[考研] 求调剂材料专硕293 +6 段_(:з」∠)_ 2026-03-10 6/300 2026-03-10 18:22 by ms629
[考研] 家人们 调剂不迷路 看这里 +8 likeihood 2026-03-09 13/650 2026-03-10 08:09 by likeihood
信息提示
请填处理意见