| 查看: 2116 | 回复: 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的时候没有考虑到这一点。 不知道大家以前有没有注意到这个问题。 |
回复此楼» 猜你喜欢
材料科学与工程320求调剂，080500
已经有11人回复
-
368化学求调剂
已经有6人回复
-
0703化学调剂 348分
已经有14人回复
-
347材料专硕求调剂
已经有17人回复
-
一志愿哈工大，初试329，求环境科学与工程调剂！
已经有11人回复
-
一志愿211，0703化学305分求调剂
已经有20人回复
-
0703化学
已经有30人回复
-
277求调剂
已经有14人回复
-
273求调剂
已经有48人回复
-
一志愿华东理工085601材料工程303分求调剂
已经有13人回复
» 本主题相关价值贴推荐，对您同样有帮助:
- 用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人回复
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
hn19870519
木虫 (正式写手)
- 1ST强帖: 1
- 应助: 45 (小学生)
- 金币: 1225.2
- 散金: 1861
- 红花: 11
- 帖子: 379
- 在线: 203.1小时
- 虫号: 867945
- 注册: 2009-10-10
- 性别: GG
- 专业: 凝聚态物性 II ：电子结构
4楼2013-08-17 10:08:06
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
6楼2013-12-30 13:55:17
7楼2013-12-30 13:56:01
8楼2013-12-31 00:28:29
highpeter
金虫 (初入文坛)
- 应助: 0 (幼儿园)
- 金币: 952.1
- 帖子: 43
- 在线: 113小时
- 虫号: 898205
- 注册: 2009-11-09
- 性别: GG
- 专业: 凝聚态物性 II ：电子结构
9楼2015-09-14 20:02:16
10楼2015-09-15 19:54:31
