24小时热门版块排行榜    

查看: 1997  |  回复: 3
当前只显示满足指定条件的回帖,点击这里查看本话题的所有回帖

tonglijia

木虫 (正式写手)

[求助] 考虑自旋(ISPIN=2)Band Decomposed Chargedensity生成的PARCHG.*.ALL文件 已有1人参与

如题我在进行考虑自旋的某一条能带(全K点)的电荷密度时发现生成的PARCHG.*.ALL文件中两部分电荷密度的值完全相同,小木虫中很多虫友都说PARCHG.*.ALL文件和CHGCAR文件的格式基本相同(除了augmentation occupaacies有区别外),如果这样理解的话应该是第一部分自旋相加,第二部分自旋相减,也就是说两部分的数值一定不相同才对,但是我在计算某条自旋劈裂很大的能带时(也就是说自旋向上和自旋向下的能带相差很大)发现两部分是一样的;请问PARCHG.*.ALL文件中的格式到底是怎么样的?

考虑自旋(ISPIN=2)Band Decomposed Chargedensity生成的PARCHG.*.ALL文件
39HU]L$Y3]N2T507S~X]]HA.png
回复此楼
人,要靠自己!!!
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

tonglijia

木虫 (正式写手)

引用回帖:
2楼: Originally posted by pangrui1985 at 2015-05-23 11:07:56
这取决于你的INCAR吧,因为VASP里面是每条能带一个自旋,是不是这里你弄错了。我做过固定能量范围的PARCHG,两部分不一样。

SYSTEM = aln

ISTART = 1
ICHARG = 11  

ENCUT = 700   
EDIFF = 1e-5   
EDIFFG = -0.02

ISMEAR = 0
SIGMA = 0.01

LREAL = Auto
PREC = Normal
ROPT = 1E-4 1E-4 1E-4
# ISYM = 0

ISPIN = 2
MAGMOM = 5*1 5*1 2*1
# NUPDOWN = 2

IBRION  = -1
# ISIF    = 2
ALGO   = Normal
# IALGO= 48
  
NSW = 0
# NELMIN = 4         
NELMDL = -5

NPAR = 2   
# NCORE = 12
# LPLANE = .TURE.
# LSCALU = .FALSE.   
# NSIM = 4

NBANDS = 70

# LCHARG = .FALSE.
# LWAVE = .FALSE.
# NWRITE  = 1  
# NBLOCK  = 1

LPARD = .TRUE.
IBAND = 24 25 26 27
# KPUSE = 60
LSEPB = .TRUE.
# LSEPK = .TRUE.
这个是我的INCAR文件,您看有什么问题吗?以前也做过一些非自旋极化的pdc没有遇到问题。
还有您提到的“因为VASP里面是每条能带一个自旋”这句话怎么理解?假设我设置IBAND=24那么vasp应该计算的是第24条自旋向上的能带和第24条自旋向下的能带的总得电荷密度吗?也就是说PARCHG是由第24条自旋向上的能带和第24条自旋向下的能带的总电荷密度组成。
谢谢您的回复!
人,要靠自己!!!
3楼2015-05-24 10:11:04
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
查看全部 4 个回答

pangrui1985

铜虫 (小有名气)

【答案】应助回帖

★ ★
感谢参与,应助指数 +1
ljw4010: 金币+2, 谢谢交流! 2015-05-23 20:47:13
这取决于你的INCAR吧,因为VASP里面是每条能带一个自旋,是不是这里你弄错了。我做过固定能量范围的PARCHG,两部分不一样。
2楼2015-05-23 11:07:56
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

韩D小希

铁虫 (正式写手)

您好,请问PARCHG的自旋问题,您解决了么?我的PARCHG也是上下两部分一样,请教下是怎么回事呢?
谢谢
4楼2016-02-20 10:50:21
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[博后之家] 售SCI一区T0P文章,我:8O.55.1.O.5.4,科目齐全,可+急 +3 TZqyNTBa07NP 2026-11-13 3/150 2026-11-15 11:58 by xSrQNSSjcgoX
[博后之家] 售SCI一区文章,我:8.O.55.1.O.54,科目齐全,可伽急 +3 SX5mJmPWDfOO 2026-11-13 3/150 2026-11-15 08:40 by xSrQNSSjcgoX
[硕博家园] 售一区SCI文章T0P,我:8O.551.O54,科目全,可十急 +3 uKOSCocCb6Fm 2026-11-13 3/150 2026-11-15 05:00 by xSrQNSSjcgoX
[考研] 售SCI一区T0P文章,我:8O.55.1.O.54,科目全,可伽急 +3 uKOSCocCb6Fm 2026-11-13 3/150 2026-11-15 04:05 by xSrQNSSjcgoX
[公派出国] 售一区SCI文章T0P,我:8O.551.O54,科目全,可十急 +3 uKOSCocCb6Fm 2026-11-12 7/350 2026-11-15 02:48 by xSrQNSSjcgoX
[考研] 售SCI一区T0P文章,我:8.O.55.1.O.5.4,科目全,可+急 +3 uKOSCocCb6Fm 2026-11-12 8/400 2026-11-15 02:26 by xSrQNSSjcgoX
[硕博家园] 售SCI文章,我:8O.5.5.1O.54,科目全,可十急 +3 uKOSCocCb6Fm 2026-11-12 9/450 2026-11-15 02:04 by xSrQNSSjcgoX
[硕博家园] 售SCI一区T0P文章,我:8.O.55.1.O.5.4,科目全,可+急 +6 FChKX8meu1ds 2026-11-09 10/500 2026-11-13 11:08 by Equinoxhua
[硕博家园] 售SCI一区T0P文章,我:8O.55.1.O.54,科目全,可伽急 +5 C6CRn66RYYXS 2026-11-10 9/450 2026-11-13 11:02 by Equinoxhua
[基金申请] 售SCI一区T0P文章,我:8O.55.1.O.54,科目全,可伽急 +6 C6CRn66RYYXS 2026-11-10 7/350 2026-11-13 10:55 by Equinoxhua
[公派出国] 售SCI一区T0P文章,我:8.O55.1.O.54,科目全,可十急 +3 TgvkbZZl9tqB 2026-11-11 4/200 2026-11-12 12:51 by ECLTXnJnRSXf
[考博] 售一区SCI文章T0P,我:8O.551.O54,科目全,可十急 +4 cdqGD6ftHmbl 2026-11-11 5/250 2026-11-12 12:28 by ECLTXnJnRSXf
[公派出国] 售SCI一区T0P文章,我:8.O.55.1.O.54,科目齐全,可+急 +3 Z7qfCsdIiD1T 2026-11-10 7/350 2026-11-12 11:31 by HajghKOMynmz
[硕博家园] 售SCI一区T0P文章,我:8O.55.1.O.5.4,科目齐全,可+急 +4 Z7qfCsdIiD1T 2026-11-10 6/300 2026-11-12 11:22 by HajghKOMynmz
[教师之家] 售SCI一区T0P文章,我:8O.55.1.O.54,科目全,可伽急 +4 Z7qfCsdIiD1T 2026-11-10 7/350 2026-11-12 11:22 by HajghKOMynmz
[公派出国] 售SCI一区文章,我:8.O.551.O.5.4,科目全,可伽急 +3 czAiolzR0IJC 2026-11-10 8/400 2026-11-12 11:11 by HajghKOMynmz
[考研] 售SCI一区T0P文章,我:8O.55.1.O.5.4,科目齐全,可+急 +7 czAiolzR0IJC 2026-11-10 9/450 2026-11-12 11:00 by HajghKOMynmz
[考研] 售SCI一区T0P文章,我:8.O.55.1.O.5.4,科目全,可+急 +5 C6CRn66RYYXS 2026-11-10 7/350 2026-11-12 10:16 by HajghKOMynmz
[考研] 售SCI文章,我:8O.5.5.1O.54,科目全,可十急 +7 0TEMgxdT1eIw 2026-11-09 9/450 2026-11-12 09:24 by HajghKOMynmz
[硕博家园] 售SCI一区T0P文章,我:8.O.55.1.O.5.4,科目全,可+急 +3 C6CRn66RYYXS 2026-11-10 5/250 2026-11-12 06:03 by 9fOTqj074tIg
信息提示
请填处理意见