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

24小时热门版块排行榜    

小木虫论坛-学术科研互动平台 » 计算模拟区 » 第一性原理 » Vasp&MedeA » VASP计算HOMO,LUMO轨道
292/3返回列表上一页123下一页
查看: 12156  |  回复: 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下一页
普通表情 龙 兔 虎 猫
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 高分子化学与物理调剂 +4 好好好1233 2026-02-28 8/400 2026-03-01 09:26 by 好好好1233
[考研] 290求调剂 +6 材料专硕调剂； 2026-02-28 7/350 2026-03-01 09:21 by L135790
[考研] 291分工科求调剂 +5 science饿饿 2026-03-01 5/250 2026-03-01 09:19 by L135790
[论文投稿] Optics letters投稿被拒求助 30+3 luckyry 2026-02-26 4/200 2026-03-01 09:06 by babero
[考研] 272求调剂 +4 田智友 2026-02-28 4/200 2026-03-01 06:43 by 刘兵
[考研] 285求调剂 +6 满头大汗的学生 2026-02-28 6/300 2026-03-01 06:29 by Trying]
[基金申请] 面上模板改不了页边距吧？ +5 ieewxg 2026-02-25 6/300 2026-03-01 00:10 by addressing
[考研] 304求调剂 +3 52hz~~ 2026-02-28 5/250 2026-03-01 00:00 by 52hz~~
[考研] 化工专硕348，一志愿985求调剂 +4 弗格个 2026-02-28 6/300 2026-02-28 22:00 by wang_dand
[考研] 材料学调剂 +5 提神豆沙包 2026-02-28 5/250 2026-02-28 21:34 by gaoxiaoniuma
[考研] 264求调剂 +3 巴拉巴拉根556 2026-02-28 3/150 2026-02-28 21:31 by gaoxiaoniuma
[考研] 311求调剂 +8 南迦720 2026-02-28 8/400 2026-02-28 21:30 by gaoxiaoniuma
[考研] 求调剂 +4 repeatt?t 2026-02-28 4/200 2026-02-28 21:16 by gaoxiaoniuma
[考研] 298求调剂 +8 人间唯你是清欢 2026-02-28 11/550 2026-02-28 20:26 by L135790
[考研] 276求调剂 +3 路lyh123 2026-02-28 4/200 2026-02-28 19:45 by 路lyh123
[考研] 0856材料求调剂 +10 hyf hyf hyf 2026-02-28 11/550 2026-02-28 18:50 by 无际的草原
[高分子] 求环氧树脂研发1名 +3 孙xc 2026-02-25 11/550 2026-02-28 16:57 by ichall
[考研] 寻找调剂 +3 LYidhsjabdj 2026-02-28 3/150 2026-02-28 12:59 by miniwendy
[硕博家园] 博士自荐 +6 科研狗111 2026-02-26 9/450 2026-02-28 12:32 by seaskyy
[硕博家园] 【博士招生】太原理工大学2026化工博士 +4 N1ce_try 2026-02-24 8/400 2026-02-26 08:40 by N1ce_try
信息提示
关闭
请填处理意见
关闭