24小时热门版块排行榜    

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

chenweiguang

金虫 (小有名气)

[交流] 有人注意过PDOS计算跟ISYM的关系吗? 已有6人参与

这个是vasp forum里面的一个帖子,注意admin的回帖:
http://cms.mpi.univie.ac.at/vasp-forum/forum_viewtopic.php?3.306
“the reason of this behaviour is that the
partial (and ONLY the partial ) magnetizations and
charge densities are not symmetrized before they are written to
OUTCAR. (the total charge density is, of course)
Therefore, the difference in the output stems from the difference
in the k-meshes which are generated:
if ISYM is switched off, the irreducible part of the Brillouin
zone is not determined, hence the k-points are spread all over
the 1st BZ, and the partial densities look symmetric, as they in fact are.
if ISYM is switched on, only the irreducible part of the BZ is filled
with k-points, but for the patrial charges no a-posteriori symmetrisation
is done (as should be and is in fact done for the charge densities).
Therefore, these partial quantities may differ slightly,
even for symmetry-equivalent atoms in a unit cell. ”

我也是在计算一个体系的px,py和pz的DOS的时候发现对称的体系不一样,觉得很奇怪。
刚开始以为是坐标系的原因,变换了坐标系还是不同。后来就搜到了上面这个帖子。
然后就做了测试,对于同一个体系,设置了ISYM=0和1,分别计算并比较。
发现总DOS完全一样,但是对于分波态密度就不一样了。
设置了ISYM=0,也就是说不不让布里渊区的K点进行简约,而是均匀撒点,这个结果px,py和pz完全一样,跟预想的一样。
这个问题在vasp的官方手册里没有看到,所以一直做PDOS的时候没有考虑到这一点。
不知道大家以前有没有注意到这个问题。
回复此楼

» 猜你喜欢

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

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

highpeter

金虫 (初入文坛)

★
小木虫: 金币+0.5, 给个红包,谢谢回帖
我也刚注意到这个问题。ISYM默认值算出来的PDOS对称性不对。应该是用ISYM=0。但很少看到有这个提示。
9楼2015-09-14 20:02:16
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
查看全部 10 个回答

wyxhenu

新虫 (小有名气)

★
小木虫: 金币+0.5, 给个红包,谢谢回帖
真没注意到,以后要试试

[ 发自手机版 http://muchong.com/3g ]
2楼2013-08-08 18:27:43
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

chenweiguang

金虫 (小有名气)

一个星期了,一直在关注这个问题,可是很可惜好像大家都不怎么关注。
我发现这个px,py,pz等与坐标系也是有很大关系的。
大家计算单个原子的PDOS不考虑这两点的话,差别有时候是很大的。
希望大家可以讨论交流一下。
3楼2013-08-16 20:51:13
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

hn19870519

木虫 (正式写手)

★
小木虫: 金币+0.5, 给个红包,谢谢回帖
碰到过楼主说的这个问题。关闭对称性之后px py pz就与坐标对称性一致了。楼主说的和坐标系有很大关系具体指的是?
4楼2013-08-17 10:08:06
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
普通表情 龙 兔 虎 猫 高级回复 (可上传附件)
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考博] 售SCI一区文章,我:8.O.55.1.O.54,科目齐全,可伽急 +10 ncTVBM6AOJHF 2026-11-04 22/1100 2026-11-06 17:18 by 5sJFDKp2kwuI
[公派出国] 售一区SCI文章T0P,我:8O.551.O54,科目全,可十急 +4 ELBZ6WwEaRj6 2026-11-05 5/250 2026-11-06 16:43 by 5sJFDKp2kwuI
[考博] 售SCI一区文章,我:8O5.5.1.O5.4,科目全,可伽急 +6 ELBZ6WwEaRj6 2026-11-05 8/400 2026-11-06 16:34 by 5sJFDKp2kwuI
[考博] 售SCI一区文章,我:8.O.55.1.O.54,科目齐全,可伽急 +6 ELBZ6WwEaRj6 2026-11-04 14/700 2026-11-06 16:32 by 5sJFDKp2kwuI
[公派出国] 售SCI-T0P文章,我:8O.5.5.1.O.54,科目齐全,可+急 +6 ELBZ6WwEaRj6 2026-11-04 11/550 2026-11-06 16:32 by 5sJFDKp2kwuI
[考研] 售SCI一区T0P文章,我:8.O.55.1.O.54,科目齐全,可+急 +5 ELBZ6WwEaRj6 2026-11-04 10/500 2026-11-06 16:23 by 5sJFDKp2kwuI
[公派出国] 售SCI一区T0P文章,我:8.O.55.1.O.54,科目齐全,可+急 +5 ttH3f1oA6ace 2026-11-04 10/500 2026-11-06 16:12 by 5sJFDKp2kwuI
[硕博家园] 售SCI一区文章,我:8O5.5.1.O5.4,科目全,可伽急 +4 ttH3f1oA6ace 2026-11-04 11/550 2026-11-06 16:10 by 5sJFDKp2kwuI
[公派出国] 售SCI一区文章,我:8.O.551.O.5.4,科目全,可伽急 +6 ttH3f1oA6ace 2026-11-04 12/600 2026-11-06 16:10 by 5sJFDKp2kwuI
[公派出国] 售SCI一区T0P文章,我:8O.55.1.O.5.4,科目齐全,可+急 +6 ncTVBM6AOJHF 2026-11-04 15/750 2026-11-06 13:49 by BkA9Bi5cMhkD
[考博] 售SCI一区T0P文章,我:8.O55.1.O.54,科目全,可十急 +7 ncTVBM6AOJHF 2026-11-04 17/850 2026-11-06 13:47 by BkA9Bi5cMhkD
[公派出国] 售SCI-T0P文章,我:8O.5.5.1.O.54,科目齐全,可+急 +8 ncTVBM6AOJHF 2026-11-04 12/600 2026-11-06 13:47 by BkA9Bi5cMhkD
[硕博家园] 售SCI一区T0P文章,我:8O.55.1.O.54,科目全,可伽急 +8 ncTVBM6AOJHF 2026-11-04 14/700 2026-11-06 13:36 by BkA9Bi5cMhkD
[硕博家园] 售SCI一区文章,我:8.O.55.1.O.54,科目齐全,可伽急 +3 ELBZ6WwEaRj6 2026-11-04 12/600 2026-11-06 05:32 by qL1cU46YIdK0
[考研] 售SCI-T0P文章,我:8O.5.5.1.O.54,科目齐全,可+急 +4 ttH3f1oA6ace 2026-11-04 14/700 2026-11-06 05:11 by qL1cU46YIdK0
[公派出国] 售SCI一区T0P文章,我:8.O.55.1.O54,科目全,可伽急 +5 kvewwGtKktqB 2026-11-04 14/700 2026-11-06 05:10 by qL1cU46YIdK0
[硕博家园] 售SCI一区T0P文章,我:8O.55.1.O.5.4,科目齐全,可+急 +6 kvewwGtKktqB 2026-11-04 16/800 2026-11-06 05:00 by qL1cU46YIdK0
[硕博家园] 售一区SCI文章T0P,我:8O.551.O54,科目全,可十急 +4 ELBZ6WwEaRj6 2026-11-04 8/400 2026-11-06 02:24 by qL1cU46YIdK0
[考研] 售SCI一区T0P文章,我:8.O55.1.O.54,科目全,可十急 +4 ttH3f1oA6ace 2026-11-04 10/500 2026-11-06 02:04 by qL1cU46YIdK0
[教师之家] 售SCI一区T0P文章,我:8O.55.1.O.54,科目全,可伽急 +4 ncTVBM6AOJHF 2026-11-04 12/600 2026-11-05 11:10 by PYZmcHr5kSaq
信息提示
请填处理意见