24小时热门版块排行榜    

查看: 12096  |  回复: 28

jpchou

专家顾问 (著名写手)

???????:
10?: Originally posted by xianglannihao at 2014-10-08 21:44:30
???? OUTCAR,??????????????? fermi level ?? state??лл?????????????????O2?????HOMO????????????6??????HOMO???????pi???????????????????????????????????????????ò??????? ...

???@????N?
PARCHG ?
???@???????? charge distribution
????]????????
charge ????????
???? pi orbital
??? wave function
wave function ?????????
??????N???
???? wave function
????? http://www.andrew.cmu.edu/user/feenstra/wavetrans/
11楼2014-10-09 02:37:03
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

zuocuiping

木虫 (职业作家)

引用回帖:
7楼: Originally posted by jpchou at 2014-10-08 18:11:12
spin up 和 spin down 的 occupancy 不一樣
就表示你計算的系統是具有磁性 (spin moment = spin up - spin down)
不同 k 點的 occupancy 當然有可能是不一樣 (ex. metal system)
要找所謂的 HOMO 當然就要找所有 ...

对于固体
k点很多
怎么找出k点呢
比如我的48个k点
有没有好点的方法去找呢
谢谢
12楼2014-10-10 08:51:17
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

jpchou

专家顾问 (著名写手)

【答案】应助回帖

????? fortran, shell script, or C

?????Щ?????????
??????????? HOMO ??? 333 ?? band ????
?? grep  "  333  "  OUTCAR  
????????? 48 ?? k ?c??? 333 band ???????
???????? 332 or 334
13楼2014-10-10 14:54:48
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

xianglannihao

铜虫 (小有名气)

???????:
11?: Originally posted by jpchou at 2014-10-09 02:37:03
???@????N?
PARCHG ?
???@???????? charge distribution
????]????????
charge ????????
???? pi orbital
??? wave function
wave function ?????????
??????N???
???? wave function
...

??PARCHG????????????????????wave function??????????????????????O2?????HOMO??????????????????2pi???????????????????????????????????CO??LUMO???????2pi???????????????????
VASP????HOMO,LUMO???
CO-LUMO-2pi.jpg

14楼2014-10-10 15:24:04
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

jpchou

专家顾问 (著名写手)

【答案】应助回帖

★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★
xianglannihao: 金币+15 2014-10-11 21:14:12
引用回帖:
14楼: Originally posted by xianglannihao at 2014-10-10 15:24:04
是PARCHG,也就是电荷密度,应该算是wave function的平方,所以没有正负。对O2分子的HOMO,其电荷密度形状应该要是2pi反键轨道形状才是对的,不知哪里出了问题。我算CO的LUMO就显示了2pi反键轨道形状,如下图。

C ...


你做的是 spin-singlet oxygen molecule (magnetization moment = 0)
如果你做 spin-triplet oxygen molecule (magnetization moment = 2)
就會看到 pi* 了
15楼2014-10-10 17:39:23
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

xianglannihao

铜虫 (小有名气)

引用回帖:
15楼: Originally posted by jpchou at 2014-10-10 17:39:23

你做的是 spin-singlet oxygen molecule (magnetization moment = 0)
如果你做 spin-triplet oxygen molecule (magnetization moment = 2)
就會看到 pi* 了...

我是设置了ISPIN=2,即考虑自旋极化,从OUTCAR文件可以看到自旋向上的有5个电子,向下的有7个电子,即有2个单电子,且算得mag=-2,应该计算的是三重态的O2。除了设置ISPIN=2,还要设置什么参数吗?如MAGMOM,NUPDOWN,你是怎么设置得到pi*?非常感谢。
16楼2014-10-10 18:28:45
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

jpchou

专家顾问 (著名写手)

引用回帖:
16楼: Originally posted by xianglannihao at 2014-10-10 18:28:45
我是设置了ISPIN=2,即考虑自旋极化,从OUTCAR文件可以看到自旋向上的有5个电子,向下的有7个电子,即有2个单电子,且算得mag=-2,应该计算的是三重态的O2。除了设置ISPIN=2,还要设置什么参数吗?如MAGMOM,NUPDO ...


那就怪了
我剛快速測了一下
mag moment=2 的第六個 band 會看到 anti-bonding pi orbital
mag moment=0 的第六個 band 的確會看到你畫的圖,一模一樣
所以才會認為你算成了 singlet oxygen molecule


如果你也是計算 triplet oxygen molecule,不懂為什麼你會看不到 anti-bonding pi ........

參數我是隨便給的 (只想單純看一下粗略結果)
所以請別太在意,如果有些參數很怪的話
#
NBANDS = 12
NELECT = 12
LREAL  = Auto
GGA    = PE
ISMEAR = -1
SIGMA  = 0.0005
ISPIN  = 2
ENCUT  = 400.00 eV
ENAUG  = 800.00 eV
NELM   = 60      
EDIFF  = 1E-04   
ALGO   = fast
#ALGO  = All
#TIME  = 0.4
NSW    = 1
EDIFFG = -0.02
IBRION = 2
ISIF   = 2
#ISYM  = 2            soft symmetrization
NWRITE = 2          most of the important information is written

NUPDOWN = 2



PBE exchange-correlation function with vasp-5.3.5

structure
   1.00000000000000     
    12.0000000000000000    0.0000000000000000    0.0000000000000000
     0.0000000000000000   11.0000000000000000    0.0000000000000000
     0.0000000000000000    0.0000000000000000   10.0000000000000000
   O
     2
Selective dynamics
Direct
  0.0514218175332555  0.0000000000000000  0.0000000000000000   T   T   T
  0.9485781824667445  0.0000000000000000  0.0000000000000000   T   T   T

  0.00000000E+00  0.00000000E+00  0.00000000E+00
  0.00000000E+00  0.00000000E+00  0.00000000E+00

symmetry is D2h
17楼2014-10-10 18:51:45
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

xianglannihao

铜虫 (小有名气)

引用回帖:
17楼: Originally posted by jpchou at 2014-10-10 18:51:45

那就怪了
我剛快速測了一下
mag moment=2 的第六個 band 會看到 anti-bonding pi orbital
mag moment=0 的第六個 band 的確會看到你畫的圖,一模一樣
所以才會認為你算成了 singlet oxygen molecule


...

我测试了下,只设置ISPIN= 2,从OUTCAR文件可以看到自旋向上的有5个电子,向下的有7个电子,即有2个单电子,且算得mag=-2,计算的是三重态的O2,但HOMO显示为反键的σ轨道;设置ISPIN  = 2,NUPDOWN = 2,算得mag=2, 计算的是三重态的O2,HOMO显示为反键的2pi*轨道。谢谢您的指点。
18楼2014-10-11 19:52:17
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

nk居士

铜虫 (初入文坛)

送红花一朵
引用回帖:
7楼: Originally posted by jpchou at 2014-10-08 18:11:12
spin up 和 spin down 的 occupancy 不一樣
就表示你計算的系統是具有磁性 (spin moment = spin up - spin down)
不同 k 點的 occupancy 當然有可能是不一樣 (ex. metal system)
要找所謂的 HOMO 當然就要找所有 ...

您好,我是vasp新手,现在在做这样的东西,但是查了很多资料,感觉还是没什么头绪,现在我想算一个有机分子在表面的吸附,看看这个分子在吸附前和吸附后的HOMO和LUMO有什么变化。这个分子拿vasp放在一个晶格里面是一个k点,依照您的方法,在OUTCAR里的确很容易找到HOMO和LUMO,但是slab确是一个周期性结构,有十几个k点。当然再算吸附体系的时候因为体系较大只有一个k点,这个时候的OUTCAR的occupation就很复杂了,用yang wang 大神的vaspmo也只能可视化分子轨道,请问这个时候我想看有机分子的HOMO和LUMO应该怎么看呢?又怎么把分子的HOMO LUMO 和 surface的“HOMO” “LUMO"区分开呢?这个问题纠结很久了,如能解答,感激不尽~
19楼2014-10-20 09:26:38
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

jpchou

专家顾问 (著名写手)

???????:
19?: Originally posted by nk??? at 2014-10-20 09:26:38
?????????vasp???????????????????????????????????????о???????????????????????????л????????????????????????????????????????????HOMO??LUMO?????仯???????????vasp????????????????? ...

?????????
?????????y??k ?c?????c?????????????? HOMO and LUMO
?L????y??k ?c???????c
??r??????M???? ??B ???I??]?? "HOMO" or "LUMO"
(????]?????^?? molecule orbital)
????????????? VBM (valence band maximum) and CBM (conduction band minimum)
??r?? VBM and CBM ???????????? k ?c????
???????? direct band gap & indirect band gap ?????

???ε??f
?ЙC???????????棬??? VBM (HOMO) & CBM (LUMO) ????
?c????ЙC??????? HOMO & LUMO ????~ ???????
??e????~ ????ЙC????~????????? k
???ЙC???????????棬????????е? k

???^??????c????? HOMO LUMO
?@?????????
??????????????Щ states ???? ????
??Щ???? surface states
???~????
??? occupancy ??????????????
20楼2014-10-27 19:05:57
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 xianglannihao 的主题更新
信息提示
请填处理意见