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

24小时热门版块排行榜    

小木虫论坛-学术科研互动平台 » 计算模拟区 » 第一性原理 » Vasp&MedeA » VASP计算HOMO,LUMO轨道
292/3返回列表上一页123下一页
查看: 12715  |  回复: 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下一页
普通表情 龙 兔 虎 猫
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[找工作] 售SCI-T0P文章，我:8O.5.5.1.O.54,科目齐全,可+急 +3 Vi50GxzrFcSG 2026-08-07 4/200 2026-08-08 10:39 by 3OOjAIS77qg2
[硕博家园] 售SCI-T0P文章，我:8O.5.5.1.O.54,科目齐全,可+急 +3 Vi50GxzrFcSG 2026-08-07 6/300 2026-08-08 10:19 by 3OOjAIS77qg2
[硕博家园] 售SCI一区文章，我:8.O.55.1.O.54,科目齐全,可伽急 +3 HEQlVqMTIA7d 2026-08-07 5/250 2026-08-08 08:59 by 6vVgjDL4CnGu
[基金申请] 售SCI-T0P文章，我:8O.5.5.1.O.54,科目齐全,可+急 +3 KXLV3nuBVBY7 2026-08-07 4/200 2026-08-08 08:21 by 6vVgjDL4CnGu
[考博] 售SCI一区文章，我:8.O.55.1.O.54,科目齐全,可伽急 +3 HEQlVqMTIA7d 2026-08-07 3/150 2026-08-08 03:47 by 6vVgjDL4CnGu
[基金申请] 基金中了 +14 laoda193707 2026-08-06 14/700 2026-08-08 00:23 by 实验小白ha
[基金申请] 娱乐 +6 Tide man 2026-08-03 6/300 2026-08-07 22:40 by 铁帽子农民
[基金申请] 关于代码变化问题，想知道的进来 +9 且听虎啸 2026-08-07 12/600 2026-08-07 21:53 by 理想三旬2025
[基金申请] 听说今天filecode变了 +24 布布和一二 2026-08-06 47/2350 2026-08-07 16:02 by zhiyanjiang
[基金申请] filecode变化情况 +6 布布和一二 2026-08-07 22/1100 2026-08-07 14:45 by 且听虎啸
[基金申请] 关于豆爷回答的JTJC与%2F数量 +4 yang182083 2026-08-06 6/300 2026-08-07 12:43 by 虫友是什么虫
[基金申请] filecode +14 等待解的谜 2026-08-06 19/950 2026-08-07 12:20 by wlwhappy
[基金申请] 求各位大神看下 100+6 hpkpkpkp 2026-08-05 33/1650 2026-08-06 14:49 by zhiyanjiang
[基金申请] 影响面上的因素 +8 布布和一二 2026-08-05 11/550 2026-08-06 10:41 by 宝贝虫子
[基金申请] 8月时间戳变的，举个手。玩一下，释放压力 +9 archvillain 2026-08-04 11/550 2026-08-05 20:06 by wlwhappy
[考博] 【2027博士申请】纳米药物递送方向 20+3 13586093586 2026-08-03 4/200 2026-08-05 09:59 by lfy8008
[基金申请] 有没有H口的？有收到消息的吗？ +3 超级海虾 2026-08-04 3/150 2026-08-04 17:26 by 学教育滴
[基金申请] 纯娱乐，不喜欢勿喷 +7 Tide man 2026-08-04 10/500 2026-08-04 15:10 by loufangrui
[基金申请] 面上提前没消息，有中的吗 +14 archvillain 2026-08-02 18/900 2026-08-04 14:42 by archvillain
[基金申请] 什么时候能放榜呀？ +3 Jacob678 2026-08-03 3/150 2026-08-03 16:14 by gltch
信息提示
关闭
请填处理意见
关闭