查看: 3617  |  回复: 12
当前只显示满足指定条件的回帖,点击这里查看本话题的所有回帖

[交流] 【解决】讨论,siesta中作HOMO 与LUMO图的方法

第一种是用denchar, 但是存在问题是denchar会给出你所有的波函数,需要一个个查看
第二种是用LDOS, 计算在Fermi level 附近的能量的电子分布,转换成*。XSF 查看,但是不能确定具体的LUMO和HOMO的具体位置
大家都用哪种?

[ Last edited by minmin_0082003 on 2011-3-2 at 09:20 ]
回复此楼

» 收录本帖的淘帖专辑推荐

第一性原理相关文档 siesta问题

» 猜你喜欢

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

» 抢金币啦!回帖就可以得到:

查看全部散金贴

已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
★ ★ ★
小木虫(金币+0.5):给个红包,谢谢回帖交流
zzy870720z(金币+2): 谢谢指点 2011-03-01 14:11:30
引用回帖:
Originally posted by minmin_0082003 at 2011-02-25 02:49:47:
还有,如果你不只是找homo跟Lumo的情况下怎么办呢,比如你想找个O2的π轨道,排除一个个筛选的办法

WriteWaveFunctions        T                 
WaveFuncKPointsScale     ReciprocalLatticeVectors
%block WaveFuncKPoints
0.000 0.000 0.000 from    1 to 14     
%endblock WaveFuncKPoints
这个是在siesta计算的fdf文件中设置的。然后计算的结果再用denchar处理。denchar有自己的输入文件。WriteWaveFunctions        T 你是不是忘记了写这一个选项啊。
具体的关于波函数的设置请参照使用手册3.0-b来说是60页。也就是6.13小节。里面有很细节的讲述。另外注意波函数的输出只有用diagon的时候才可以用。OderN是不行的。
7楼2011-03-01 13:41:59
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
查看全部 13 个回答
★
minmin_0082003(金币+10): 2011-02-24 18:19:44
zzy870720z(金币+1): 谢谢交流 2011-02-25 12:07:08
引用回帖:
Originally posted by minmin_0082003 at 2011-02-24 07:30:31:
第一种是用denchar, 但是存在问题是denchar会给出你所有的波函数,需要一个个查看
第二种是用LDOS, 计算在Fermi level 附近的能量的电子分布,转换成*。XSF 查看,但是不能确定具体的LUMO和HOMO的具体位置
大 ...

谁说denchar会给出所有的波函数?你可以在输入文件中设置要计算那些轨道的。提前数好了,你的HOMO,LUMO是哪几条,不就OK了!
请看一下WaveFuncKPoints的用法。
2楼2011-02-24 15:26:57
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
引用回帖:
Originally posted by zhangguangping at 2011-02-24 15:26:57:
谁说denchar会给出所有的波函数?你可以在输入文件中设置要计算那些轨道的。提前数好了,你的HOMO,LUMO是哪几条,不就OK了!
请看一下WaveFuncKPoints的用法。

你说的这个定义,我是已经定义过了
不过也有可能是数错了
我首先是做一下优化,看的*.EIG 里面的能量,因为我只要Gamma点,所以没有考虑其他的,我认为的HOMO是小于费米能级的所有数中最大的一个,而HOMO是大于费米能级中最小的一个数。然后我就开始数这两个数的位置,假设这两个数排在第19位与23位,那我该这样定义
%block WaveFuncKPoints
0.000 0.000 0.000 19 23 # Gamma wavefuncs 3 4
%endblock WaveFuncKPoints
对吧?嘿嘿
3楼2011-02-24 18:19:27
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
★ ★ ★ ★
小木虫(金币+0.5):给个红包,谢谢回帖交流
youzhizhe(金币+3): 谢谢交流。 2011-02-24 19:52:03
引用回帖:
Originally posted by minmin_0082003 at 2011-02-24 11:19:27:
你说的这个定义,我是已经定义过了
不过也有可能是数错了
我首先是做一下优化,看的*.EIG 里面的能量,因为我只要Gamma点,所以没有考虑其他的,我认为的HOMO是小于费米能级的所有数中最大的一个,而HOMO是大 ...

用你这个方法也行。不过怎么可能HOMO是19,LUMO是23呢?计算一下你的体系里面有多少电子,按照你体系的自旋简并情况,就知道HOMO,LUMO的位置了。
我不清楚你是不是你的格式不对的原因。
%block WaveFuncKPoints
0.000 0.000 0.000 19 23 # Gamma wavefuncs 3 4
%endblock WaveFuncKPoints
我写的话,我会这样写:
%block WaveFuncKPoints
0.000 0.000 0.000 from 19 to 23    # Gamma wavefuncs 3 4
%endblock WaveFuncKPoints

加上一个from to
手册上有这个,所以要按照手册来好些!我不清楚你的那个写法管不管用,起码我没用过。我这样写的,都可以的。
4楼2011-02-24 19:24:17
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
普通表情 龙 兔 虎 猫 高级回复 (可上传附件)
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[硕博家园] 极限拉扯,如何抉择 +6 Isak0320 2026-10-12 8/400 2026-10-12 20:09 by caixiaoye1
[基金申请] 入职后第一次管课题组的专利缴费,对账才发现官费比别人多掏了一倍 +6 13108017953 2026-10-08 6/300 2026-10-12 20:07 by w3w3y3
[有机交流] OpenAI又攻克数学难题了,离AI做有机合成还有多远? +4 wangyl_123 2026-10-10 8/400 2026-10-12 18:09 by wangyl_123
[教师之家] 课题组会帮你通过考核么? +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一区T0P文章,我:8.O55.1.O.54,科目全,可十急 +5 qR2AdCWTVtkw 2026-10-09 8/400 2026-10-11 22:26 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 9/450 2026-10-11 22:15 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,科目齐全,可+急 +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 jing816y 2026-10-08 9/450 2026-10-10 09:05 by zapen
信息提示
请填处理意见