| 查看: 2092 | 回复: 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的时候没有考虑到这一点。 不知道大家以前有没有注意到这个问题。 |
回复此楼» 猜你喜欢
材料科学与工程339求调剂
已经有9人回复
-
材料与化工调剂一志愿大连海事085600，349
已经有9人回复
-
292求调剂
已经有8人回复
-
275求调剂
已经有7人回复
-
材料调剂
已经有11人回复
-
288求调剂 一志愿哈工大 材料与化工
已经有25人回复
-
085602化学工程268分蹲调剂
已经有4人回复
-
282求调剂
已经有14人回复
-
301求调剂
已经有3人回复
-
总分343，求生物学调剂
已经有3人回复
» 本主题相关价值贴推荐，对您同样有帮助:
- 用VASP计算Mg的PDOS的时候出现d电子是什么原因？
已经有6人回复
- 分析TDOS中间隙态组成时PDOS的计算和讨论
已经有10人回复
- VASP计算出来的PDOS怎么画图啊？求指教
已经有5人回复
- 【求助成功 】PWscf 进行pdos计算报错 error while reading from file
已经有13人回复
- 【求助】vasp怎么计算PDOS?
已经有10人回复
- 【求助】Abinit中计算PDOS时怎样确定投影原子半径？
已经有4人回复
- 【求助】MS５.０ Castep 计算不出部分态密度（PDOS）怎么办？（站内有人提问MS4.3，�
已经有19人回复
- 【求助】请教siesta3.0-beta如何计算PDOS、overlap population及coop
已经有12人回复
hn19870519
实习版主
- 1ST强帖: 1
- 应助: 45 (小学生)
- 金币: 1225.2
- 散金: 1861
- 红花: 11
- 帖子: 379
- 在线: 203.1小时
- 虫号: 867945
- 注册: 2009-10-10
- 性别: GG
- 专业: 凝聚态物性 II ：电子结构
4楼2013-08-17 10:08:06
wyxhenu
兑换贵宾
- 应助: 13 (小学生)
- 金币: 3075.5
- 散金: 772
- 红花: 5
- 帖子: 279
- 在线: 250.6小时
- 虫号: 700545
- 注册: 2009-02-13
- 专业: 凝聚态物性I:结构、力学和
2楼2013-08-08 18:27:43
chenweiguang
管理员
- 应助: 11 (小学生)
- 金币: 576.2
- 红花: 4
- 帖子: 281
- 在线: 153.6小时
- 虫号: 441663
- 注册: 2007-10-27
- 性别: GG
- 专业: 凝聚态物性 II ：电子结构
3楼2013-08-16 20:51:13
chenweiguang
管理员
- 应助: 11 (小学生)
- 金币: 576.2
- 红花: 4
- 帖子: 281
- 在线: 153.6小时
- 虫号: 441663
- 注册: 2007-10-27
- 性别: GG
- 专业: 凝聚态物性 II ：电子结构
|
你比如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
10