24小时热门版块排行榜

返回列表

当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖

tonglijia

木虫 (正式写手)

应助: 12 (小学生)
金币: 1150.6
散金: 4
红花: 13
帖子: 519
在线: 305.9小时
虫号: 1776353
注册: 2012-04-24
性别: GG
专业: 金属材料的制备科学与跨学

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

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

39HU]L$Y3]N2T507S~X]]HA.png

回复此楼

» 猜你喜欢

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

VASP计算HOMO,LUMO轨道已经有28人回复
STM原理及其第一性原理模拟的物理意义已经有18人回复
【求助】vasp计算Partial Charge，能带交叉的如何确定IBAND？已经有13人回复

人，要靠自己！！！

1楼 2015-05-22 17:14:34

已阅回复此楼关注TA 给TA发消息送TA红花 TA的回帖

tonglijia

木虫 (正式写手)

应助: 12 (小学生)
金币: 1150.6
散金: 4
红花: 13
帖子: 519
在线: 305.9小时
虫号: 1776353
注册: 2012-04-24
性别: GG
专业: 金属材料的制备科学与跨学

引用回帖:

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.
这个是我的ＩＮＣＡＲ文件，您看有什么问题吗？以前也做过一些非自旋极化的ｐｄｃ没有遇到问题。
还有您提到的“因为VASP里面是每条能带一个自旋”这句话怎么理解？假设我设置IBAND＝24那么ｖａｓｐ应该计算的是第24条自旋向上的能带和第24条自旋向下的能带的总得电荷密度吗？也就是说PARCHG是由第24条自旋向上的能带和第24条自旋向下的能带的总电荷密度组成。
谢谢您的回复！