版块导航: 正在加载中...

登录注册

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

24小时热门版块排行榜

返回列表

hexin1982

新虫 (小有名气)

应助: 3 (幼儿园)
金币: 546
帖子: 148
在线: 259.6小时
虫号: 1835611
注册: 2012-05-26
专业: 凝聚态物性 II ：电子结构

[求助] PDOS中如何确定各分轨道已有1人参与

下面是我计算得到的Fe的3d电子的PDOS
# E (eV) ldos(E) pdos(E) pdos(E) pdos(E) pdos(E) pdos(E)
-5.000  0.000E+00  0.000E+00  0.000E+00  0.000E+00  0.000E+00  0.000E+00
-4.900  0.000E+00  0.000E+00  0.000E+00  0.000E+00  0.000E+00  0.000E+00
-4.800  0.000E+00  0.000E+00  0.000E+00  0.000E+00  0.000E+00  0.000E+00
-4.700  0.000E+00  0.000E+00  0.000E+00  0.000E+00  0.000E+00  0.000E+00
-4.600  0.000E+00  0.000E+00  0.000E+00  0.000E+00  0.000E+00  0.000E+00
-4.500  0.000E+00  0.000E+00  0.000E+00  0.000E+00  0.000E+00  0.000E+00
-4.400  0.000E+00  0.000E+00  0.000E+00  0.000E+00  0.000E+00  0.000E+00
-4.300  0.000E+00  0.000E+00  0.000E+00  0.000E+00  0.000E+00  0.000E+00
-4.200  0.000E+00  0.000E+00  0.000E+00  0.000E+00  0.000E+00  0.000E+00
-4.100  0.000E+00  0.000E+00  0.000E+00  0.000E+00  0.000E+00  0.000E+00
-4.000  0.000E+00  0.000E+00  0.000E+00  0.000E+00  0.000E+00  0.000E+00
-3.900  0.000E+00  0.000E+00  0.000E+00  0.000E+00  0.000E+00  0.000E+00
-3.800  0.000E+00  0.000E+00  0.000E+00  0.000E+00  0.000E+00  0.000E+00
-3.700  0.000E+00  0.000E+00  0.000E+00  0.000E+00  0.000E+00  0.000E+00
-3.600  0.000E+00  0.000E+00  0.000E+00  0.000E+00  0.000E+00  0.000E+00
-3.500  0.000E+00  0.000E+00  0.000E+00  0.000E+00  0.000E+00  0.000E+00
-3.400  0.000E+00  0.000E+00  0.000E+00  0.000E+00  0.000E+00  0.000E+00
-3.300  0.000E+00  0.000E+00  0.000E+00  0.000E+00  0.000E+00  0.000E+00
-3.200  0.000E+00  0.000E+00  0.000E+00  0.000E+00  0.000E+00  0.000E+00
-3.100  0.000E+00  0.000E+00  0.000E+00  0.000E+00  0.000E+00  0.000E+00
-3.000  0.000E+00  0.000E+00  0.000E+00  0.000E+00  0.000E+00  0.000E+00
-2.900  0.000E+00  0.000E+00  0.000E+00  0.000E+00  0.000E+00  0.000E+00
-2.800  0.000E+00  0.000E+00  0.000E+00  0.000E+00  0.000E+00  0.000E+00
-2.700  0.000E+00  0.000E+00  0.000E+00  0.000E+00  0.000E+00  0.000E+00
。。。。。。。。。。。。。。。

请问如何分辨这5个PDOS分别属于哪个轨道？这里的Fe原子处于八面体晶场中。

回复此楼

» 猜你喜欢

1楼 2016-07-24 21:30:30

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

hexin1982

新虫 (小有名气)

应助: 3 (幼儿园)
金币: 546
帖子: 148
在线: 259.6小时
虫号: 1835611
注册: 2012-05-26
专业: 凝聚态物性 II ：电子结构

怪我粗心，这里的给出的数据全为0，后面的不为0的

赞一下

回复此楼

2楼2016-07-24 21:33:17

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

jpchou

专家顾问 (著名写手)

专家经验: +20
1ST强帖: 5
应助: 231 (大学生)
金币: 2342
散金: 4367
红花: 103
帖子: 1061
在线: 757.3小时
虫号: 585915
注册: 2008-08-08
性别: GG
专业: 凝聚态物性 II ：电子结构
管辖: 第一性原理

???@???n???????????????? VASP output files (DOSCAR, EIGENVAL, PROCAR....)
?????????????^??n???

赞一下

回复此楼

3楼2016-07-25 17:28:33

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

souledge

专家顾问 (著名写手)

专家经验: +174
1ST强帖: 19
应助: 303 (大学生)
金币: 11306.4
散金: 1123
红花: 108
帖子: 2814
在线: 916.1小时
虫号: 941486
注册: 2010-01-12
专业: 结构陶瓷
管辖: 第一性原理

【答案】应助回帖

感谢参与，应助指数 +1

在QE的源代码里有提示：

CODE:

  CHARACTER (len=1)  :: l_label(0:3)=(/'s','p','d','f'/)

  CHARACTER (len=7)  :: lm_label(1:7,1:3)=reshape( (/ &

    'z      ','x      ','y      ','       ','       ','       ','       ', &

    'z2     ','xz     ','yz     ','x2-y2  ','xy     ','       ','       ', &

    'z3     ','xz2    ','yz2    ','zx2-zy2','xyz    ','x3-3xy2','3yx2-y3' /), (/7,3/) )

projwfc.x的输出中也有这样的提示：

CODE:

     state # 129: atom  33 (Ti ), wfc  1 (l=0 m= 1)

     state # 130: atom  33 (Ti ), wfc  2 (l=1 m= 1)

     state # 131: atom  33 (Ti ), wfc  2 (l=1 m= 2)

     state # 132: atom  33 (Ti ), wfc  2 (l=1 m= 3)

     state # 133: atom  33 (Ti ), wfc  3 (l=0 m= 1)

     state # 134: atom  33 (Ti ), wfc  4 (l=2 m= 1)

     state # 135: atom  33 (Ti ), wfc  4 (l=2 m= 2)

     state # 136: atom  33 (Ti ), wfc  4 (l=2 m= 3)

     state # 137: atom  33 (Ti ), wfc  4 (l=2 m= 4)

     state # 138: atom  33 (Ti ), wfc  4 (l=2 m= 5)

在PP/Doc中，还有更加详细的文档说明：

CODE:

   ::: Orbital Order

      Order of m-components for each l in the output:

          1, cos(phi), sin(phi), cos(2*phi), sin(2*phi), .., cos(l*phi), sin(l*phi)

      where phi is the polar angle:x=r cos(theta)cos(phi), y=r cos(theta)sin(phi)

      This is determined in file flib/ylmr2.f90 that calculates spherical harmonics.

      for l=1:

        1 pz     (m=0)

        2 px     (real combination of m=+/-1 with cosine)

        3 py     (real combination of m=+/-1 with sine)

      for l=2:

        1 dz2    (m=0)

        2 dzx    (real combination of m=+/-1 with cosine)

        3 dzy    (real combination of m=+/-1 with sine)

        4 dx2-y2 (real combination of m=+/-2 with cosine)

        5 dxy    (real combination of m=+/-2 with sine)

赞一下(1人)

回复此楼

思想重于技巧，内涵重于表象

4楼2016-07-25 21:45:05

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

souledge

专家顾问 (著名写手)

专家经验: +174
1ST强帖: 19
应助: 303 (大学生)
金币: 11306.4
散金: 1123
红花: 108
帖子: 2814
在线: 916.1小时
虫号: 941486
注册: 2010-01-12
专业: 结构陶瓷
管辖: 第一性原理

???????:

3?: Originally posted by jpchou at 2016-07-25 17:28:33
???@???n???????????????? VASP output files (DOSCAR, EIGENVAL, PROCAR....)
?????????????^??n???

?????QE (Quantum-Espresso) ???????????VASP??

赞一下

回复此楼

思想重于技巧，内涵重于表象

5楼2016-07-25 21:46:56

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

hexin1982

新虫 (小有名气)

应助: 3 (幼儿园)
金币: 546
帖子: 148
在线: 259.6小时
虫号: 1835611
注册: 2012-05-26
专业: 凝聚态物性 II ：电子结构

引用回帖:

5楼: Originally posted by souledge at 2016-07-25 21:46:56

这个是QE (Quantum-Espresso) 的输出，不是VASP的...

谢谢回复，我弄清楚了，就是被这个m困扰。
for l=2:
      1 dz2 (m=0)
      2 dzx (real combination of m=+/-1 with cosine)
      3 dzy (real combination of m=+/-1 with sine)
      4 dx2-y2 (real combination of m=+/-2 with cosine)
      5 dxy (real combination of m=+/-2 with sine)
这前面的12345用m来表示，与磁量子数弄混了。

赞一下

回复此楼

6楼2016-07-25 21:59:57

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