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

24小时热门版块排行榜    

小木虫论坛-学术科研互动平台 » 计算模拟区 » 第一性原理 » Vasp&MedeA » VASP计算HOMO,LUMO轨道
292/3返回列表上一页123下一页
查看: 12096  |  回复: 28
只看楼主 @他人 存档 新回复提醒 (忽略) 收藏    在APP中查看

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/
赞一下(2人) 回复此楼
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 ??????????????
赞一下(4人) 回复此楼
20楼2014-10-27 19:05:57
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 xianglannihao 的主题更新
292/3返回列表上一页123下一页
普通表情 龙 兔 虎 猫
信息提示
关闭
请填处理意见
关闭