应《网络安全法》要求，自2017年10月1日起，未进行实名认证将不得使用互联网跟帖服务。为保障您的帐号能够正常使用，请尽快对帐号进行手机号验证，感谢您的理解与支持！

24小时热门版块排行榜    

小木虫论坛-学术科研互动平台 » 计算模拟区 » 第一性原理 » Vasp&MedeA » 有人注意过PDOS计算跟ISYM的关系吗？
101/1返回列表
查看: 1931  |  回复: 9
只看楼主 @他人 存档 新回复提醒 (忽略) 收藏    在APP中查看

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的时候没有考虑到这一点。
不知道大家以前有没有注意到这个问题。
回复此楼

» 猜你喜欢

» 本主题相关价值贴推荐，对您同样有帮助:
1楼 2013-08-08 16:28:44
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

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的回帖

chenweiguang

金虫 (小有名气)
引用回帖:
4楼: Originally posted by hn19870519 at 2013-08-17 10:08:06
碰到过楼主说的这个问题。关闭对称性之后px py pz就与坐标对称性一致了。楼主说的和坐标系有很大关系具体指的是？

你比如fcc的（111）面
对于这种：
1.0
        3.5638182163         0.0000000000         0.0000000000
       -1.7819091082         3.0863571098         0.0000000000
        0.0000000000         0.0000000000        34.9181442261
和
   5.040
    0.5000000000000000 -0.5000000000000000  0.0000000000000000
    0.0000000000000000  0.5000000000000000 -0.5000000000000000
    4.0000000000000000  4.0000000000000000  4.0000000000000000
这两种是一样的，但是xyz的方向不一样，你得到的DOS的px,py,pz就不一样。
赞一下 回复此楼
5楼2013-08-22 00:21:13
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

foolcatwy

铁虫 (小有名气)
★
小木虫: 金币+0.5, 给个红包，谢谢回帖
这个问题我才关注到，我也发现论坛里的同学们好像不太注意这个问题，我是在提取能带成分的时候发现不对的，不知道楼主最后怎么解决的
赞一下 回复此楼
6楼2013-12-30 13:55:17
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

foolcatwy

铁虫 (小有名气)
引用回帖:
4楼: Originally posted by hn19870519 at 2013-08-17 10:08:06
碰到过楼主说的这个问题。关闭对称性之后px py pz就与坐标对称性一致了。楼主说的和坐标系有很大关系具体指的是？

你是指在计算的时候设置ISYM=0?就可以了吗
赞一下 回复此楼
7楼2013-12-30 13:56:01
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

gump_813

金虫 (小有名气)
★
小木虫: 金币+0.5, 给个红包，谢谢回帖
有道理 没注意过这个 话说为什么vasp不顺便把partial charge density也给对称一下呢
赞一下 回复此楼
8楼2013-12-31 00:28:29
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

highpeter

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

JJ20

金虫 (小有名气)
★
小木虫: 金币+0.5, 给个红包，谢谢回帖
这个问题我遇到过，也是算pdos不简并，所以一直是使用isym=-1,完全取消对称性去计算，具体没有做过太多测试，
赞一下 回复此楼
选择放弃只为保留尊严
10楼2015-09-15 19:54:31
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 chenweiguang 的主题更新
101/1返回列表
普通表情 龙 兔 虎 猫 高级回复 (可上传附件)
信息提示
关闭
请填处理意见
关闭