24小时热门版块排行榜    

查看: 3361  |  回复: 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 的主题更新
普通表情 高级回复 (可上传附件)
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 本科郑州大学物理学院,一志愿华科070200学硕,346求调剂 +4 我不是一根葱 2026-03-18 4/200 2026-03-19 09:11 by 浮云166
[考研] 344求调剂 +6 knight344 2026-03-16 7/350 2026-03-18 20:13 by walc
[考研] 297求调剂 +8 戏精丹丹丹 2026-03-17 8/400 2026-03-18 14:30 by laoshidan
[考研] 085601材料工程专硕求调剂 +6 慕寒mio 2026-03-16 6/300 2026-03-18 14:26 by 007_lilei
[考研] 收复试调剂生 +4 雨后秋荷 2026-03-18 4/200 2026-03-18 14:16 by elevennnne
[考研] 304求调剂 +12 小熊joy 2026-03-14 13/650 2026-03-18 12:34 by Linda Hu
[考研] 工科材料085601 279求调剂 +6 困于星晨 2026-03-17 6/300 2026-03-18 10:21 by kkcoco25
[硕博家园] 湖北工业大学 生命科学与健康学院-课题组招收2026级食品/生物方向硕士 +3 1喜春8 2026-03-17 5/250 2026-03-17 17:18 by ber川cool子
[考研] 26考研求调剂 +6 丶宏Sir 2026-03-13 6/300 2026-03-17 16:13 by 醉在风里
[考研] 一志愿苏州大学材料工程(085601)专硕有科研经历三项国奖两个实用型专利一项省级立项 +6 大火山小火山 2026-03-16 8/400 2026-03-17 15:05 by 无懈可击111
[考研] 290求调剂 +6 孔志浩 2026-03-12 11/550 2026-03-17 14:41 by 周舟舟77
[考研] 一志愿,福州大学材料专硕339分求调剂 +3 木子momo青争 2026-03-15 3/150 2026-03-17 07:52 by laoshidan
[考研] 11408 一志愿西电,277分求调剂 +3 zhouzhen654 2026-03-16 3/150 2026-03-17 07:03 by laoshidan
[考研] 求老师收留调剂 +4 jiang姜66 2026-03-14 5/250 2026-03-15 20:11 by Winj1e
[考研] 070305求调剂 +3 mlpqaz03 2026-03-14 4/200 2026-03-15 11:04 by peike
[考研] 085601材料工程315分求调剂 +3 yang_0104 2026-03-15 3/150 2026-03-15 10:58 by peike
[考研] 中科大材料专硕319求调剂 +3 孟鑫材料 2026-03-13 3/150 2026-03-14 18:10 by houyaoxu
[考研] 266求调剂 +4 学员97LZgn 2026-03-13 4/200 2026-03-14 08:37 by zhukairuo
[考研] 304求调剂 +6 Mochaaaa 2026-03-12 7/350 2026-03-13 22:18 by 星空星月
[考研] 求调剂 +3 程雨杭 2026-03-12 3/150 2026-03-13 15:06 by JourneyLucky
信息提示
请填处理意见