查看: 4037  |  回复: 11

neweroica

木虫 (著名写手)

[交流] 【求助】请教一下PROCAR的格式 已有1人参与

当我用LORBIT=2或12,计算得到的PROCAR中,有两处原子轨道的信息,如:

band   1 # energy  -20.57931908 # occ.  2.00000000

ion      s     py     pz     px    dxy    dyz    dz2    dxz    dx2    tot
  1  0.031  0.001  0.001  0.001  0.000  0.000  0.000  0.000  0.000  0.035
  2  0.029  0.003  0.001  0.000  0.000  0.000  0.000  0.000  0.000  0.033
  3  0.099  0.000  0.000  0.000  0.000  0.000  0.000  0.000  0.000  0.099
……    ……    ……    ……    ……    ……    ……    ……    ……    ……    ……    ……
tot  0.621  0.012  0.029  0.011  0.000  0.000  0.000  0.000  0.000  0.672

ion      s     py     pz     px    dxy    dyz    dz2    dxz    dx2
  1 -0.442  0.084  0.079 -0.093  0.000  0.000  0.000  0.000  0.000
  1 -0.065  0.012  0.012 -0.014  0.000  0.000  0.000  0.000  0.000
  2 -0.429  0.133 -0.056 -0.035  0.000  0.000  0.000  0.000  0.000
  2 -0.063  0.020 -0.008 -0.005  0.000  0.000  0.000  0.000  0.000
  3 -0.791 -0.028 -0.011  0.020  0.000  0.000  0.000  0.000  0.000
  3 -0.117 -0.004 -0.002  0.003  0.000  0.000  0.000  0.000  0.000
……    ……    ……    ……    ……    ……    ……    ……    ……    ……    ……    ……

band   2 # energy  -19.85477436 # occ.  2.00000000
……    ……    ……    ……    ……    ……    ……    ……    ……    ……    ……    ……


请问:
1. 第二段的数值(红色和蓝色)按手册说,是相位系数(phase factor),但为什么每个原子有两行数值? 我起先怀疑是实部(红色)和虚部(蓝色),但发现对于Gamma点,这两行数值都不为零(Gamma点只有实部)。 我不明白。


2. 第一段的数值(绿色)是什么?  phase factor是不就是原子轨道的系数?那为什么又要第一行呢?

请大家指教!  万分感谢!!!

[ Last edited by neweroica on 2010-1-15 at 18:55 ]
回复此楼

» 猜你喜欢

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

Simpler, stronger, more tolerant
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

calos818

木虫 (著名写手)

新人菜鸟

我也很想知道,帮你顶下吧
2楼2010-01-15 21:17:49
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

njuxyh

木虫 (著名写手)

★ ★ ★
fegg7502(金币+1,VIP+0):thank you very much! 1-17 18:03
neweroica(金币+2,VIP+0):谢谢。那为什么spin up 和spin down的数据对双占据的轨道数值也会不一样呢? 1-19 00:41
是 spin up 和四spin down 的数据。楼上设了Lorbit, 是算磁性的参数阿
3楼2010-01-16 17:17:20
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

neweroica

木虫 (著名写手)

那为什么spin up 和spin down的数据对双占据的轨道数值也会不一样呢?
而且对于非Gamma点,虚数部分的轨道系数在哪里找?
Simpler, stronger, more tolerant
4楼2010-01-19 00:42:01
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

wangyujia

木虫 (正式写手)

★ ★ ★ ★ ★ ★ ★
neweroica(金币+1):谢谢参与
neweroica(金币+5,VIP+0):非常感谢!待我进一步研究一下。 1-20 17:10
ice_rain(金币+1,VIP+0):谢谢分享 1-20 22:45
不太懂你的问题,不过我以前下过一个处理PROCAR文件的fortran程序,可能会对你有帮助。
5楼2010-01-20 14:47:25
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

hchfox

木虫 (正式写手)

★ ★
neweroica(金币+1):谢谢参与
neweroica(金币+1,VIP+0):请说明一下你的理由。 1-20 17:11
虽然我不太明白,但我知道你的说法肯定是错的,呵。。。
引用回帖:
Originally posted by njuxyh at 2010-1-16 17:17:
是 spin up 和四spin down 的数据。楼上设了Lorbit, 是算磁性的参数阿

6楼2010-01-20 16:23:01
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

hchfox

木虫 (正式写手)

★ ★
neweroica(金币+1,VIP+0):虽然不是很可信,但我只是想在研究的时候弄一些定性的图像,并不是用来发表的。而且很多搞vasp的都说他们可以做coop。所以,我一定要成功! 1-20 17:06
aylayl08(金币+1,VIP+0):谢谢讨论 1-20 20:12
想搞COOP啊,经过我的验证vasp算的COOP不是很可信,劝你还是放弃。。。
引用回帖:
Originally posted by neweroica at 2010-1-15 18:51:
当我用LORBIT=2或12,计算得到的PROCAR中,有两处原子轨道的信息,如:

band   1 # energy  -20.57931908 # occ.  2.00000000

ion      s     py     pz     px    dxy    dyz    dz2    dxz    dx2    to ...

7楼2010-01-20 16:26:38
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

njuxyh

木虫 (著名写手)

★ ★
neweroica(金币+1,VIP+0):我的理由前面说过了:为什么spin up 和spin down的数据对双占据的轨道数值也会不一样呢? 而且对于非Gamma点,虚数部分的轨道系数在哪里找? 1-20 17:11
aylayl08(金币+1,VIP+0):谢谢讨论 1-20 20:12
呵呵,
我也是看到设的参数,
才认为和磁性有关。
既然楼上那么肯定我是错的,
那就请说说理由吧。
毕竟我也想知道那两行是什么意思。
8楼2010-01-20 16:41:44
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

hchfox

木虫 (正式写手)

★ ★
aylayl08(金币+2,VIP+0):谢谢指点 1-20 20:11
别的我不知道,我就知道LORBIT跟磁性无任何关系,呵。。。
以前搞过,在我的印象中,好像有实虚部(可能翻译不太对)的说法,自己再深抠下吧。。。
引用回帖:
Originally posted by njuxyh at 2010-1-20 16:41:
呵呵,
我也是看到设的参数,
才认为和磁性有关。
既然楼上那么肯定我是错的,
那就请说说理由吧。
毕竟我也想知道那两行是什么意思。

9楼2010-01-20 17:16:09
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

westmonster

金虫 (正式写手)

★
neweroica(金币+1):谢谢参与
一直没搞清楚这个PROCAR的格式问题,楼主有进展吗? 赐教啊~~
蝉,于其脱壳展翅,则蛰居地下,似无声无息,实则以备有所为。
10楼2010-04-27 19:27:22
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 neweroica 的主题更新
普通表情 龙 兔 虎 猫 高级回复 (可上传附件)
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[有机交流] OpenAI又攻克数学难题了,离AI做有机合成还有多远? +4 wangyl_123 2026-10-10 8/400 2026-10-12 18:09 by wangyl_123
[硕博家园] 极限拉扯,如何抉择 +5 Isak0320 2026-10-12 7/350 2026-10-12 18:03 by 有人要害朕
[教师之家] 课题组会帮你通过考核么? +3 好运来~9527 2026-10-10 6/300 2026-10-12 14:55 by 简单化xn
[公派出国] 售一区SCI文章T0P,我:8O.551.O54,科目全,可十急 +5 qR2AdCWTVtkw 2026-10-09 9/450 2026-10-11 22:34 by XQNyDCa1Tk5n
[公派出国] 售SCI一区文章,我:8O5.5.1.O5.4,科目全,可伽急 +5 qR2AdCWTVtkw 2026-10-09 10/500 2026-10-11 22:26 by XQNyDCa1Tk5n
[找工作] 售SCI一区T0P文章,我:8.O55.1.O.54,科目全,可十急 +5 qR2AdCWTVtkw 2026-10-09 8/400 2026-10-11 22:26 by XQNyDCa1Tk5n
[考博] 售SCI一区文章,我:8.O.55.1.O.54,科目齐全,可伽急 +3 qR2AdCWTVtkw 2026-10-09 7/350 2026-10-11 22:11 by XQNyDCa1Tk5n
[教师之家] 售SCI文章,我:8O5.5.1.O.54,科目齐全,可+急 (EPI+-1)(金币-50) +5 qR2AdCWTVtkw 2026-10-09 10/500 2026-10-11 22:06 by XQNyDCa1Tk5n
[基金申请] 售SCI一区T0P文章,我:8.O.55.1.O54,科目全,可伽急 +5 qR2AdCWTVtkw 2026-10-09 9/450 2026-10-11 22:06 by XQNyDCa1Tk5n
[教师之家] 售SCI一区T0P文章,我:8.O.55.1.O.5.4,科目全,可+急 (EPI+-1)(金币-50) +5 qR2AdCWTVtkw 2026-10-09 10/500 2026-10-11 14:25 by VOMexwnDzOFZ
[教师之家] 售SCI文章,我:8O.5.5.1O.54,科目全,可十急 (EPI+-1)(金币-50) +4 qR2AdCWTVtkw 2026-10-09 7/350 2026-10-11 14:05 by VOMexwnDzOFZ
[考研] 售SCI文章,我:8O.5.5.1O.54,科目全,可十急 +5 qR2AdCWTVtkw 2026-10-09 8/400 2026-10-11 13:51 by VOMexwnDzOFZ
[硕博家园] 售SCI-T0P文章,我:8O.5.5.1.O.54,科目齐全,可+急 +5 qR2AdCWTVtkw 2026-10-09 8/400 2026-10-11 13:45 by VOMexwnDzOFZ
[找工作] 售SCI一区T0P文章,我:8O.55.1.O.5.4,科目齐全,可+急 +3 qR2AdCWTVtkw 2026-10-09 7/350 2026-10-11 10:54 by VOMexwnDzOFZ
[找工作] 售SCI一区T0P文章,我:8O.55.1.O.5.4,科目齐全,可+急 +4 qR2AdCWTVtkw 2026-10-09 7/350 2026-10-11 10:34 by VOMexwnDzOFZ
[考博] 招收博士生一名 +3 mandolin 2026-10-10 3/150 2026-10-11 10:21 by 不是山谷!
[公派出国] 售SCI一区T0P文章,我:8O.55.1.O.54,科目全,可伽急 +4 qR2AdCWTVtkw 2026-10-09 4/200 2026-10-11 03:14 by M0zrYVSAnxpE
[有机交流] 各位大佬,路线设计求助 +6 zapen 2026-10-08 8/400 2026-10-10 15:59 by zapen
[基金申请] 入职后第一次管课题组的专利缴费,对账才发现官费比别人多掏了一倍 +5 13108017953 2026-10-08 5/250 2026-10-10 11:15 by wdwd123
[有机交流] 我应该选哪条合成路线让实验成功率高一些? +5 jing816y 2026-10-08 9/450 2026-10-10 09:05 by zapen
信息提示
请填处理意见