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

24小时热门版块排行榜    

小木虫论坛-学术科研互动平台 » 计算模拟区 » 第一性原理 » Siesta&Smeagol&Atk » 【解决】讨论，siesta中作HOMO 与LUMO图的方法
北京石油化工学院2026年研究生招生接收调剂公告
51/1返回列表
查看: 3060  |  回复: 12
只看楼主 @他人 存档 新回复提醒 (忽略) 收藏    在APP中查看
当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖

minmin_0082003

金虫 (正式写手)

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

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

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

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

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

» 猜你喜欢

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

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

查看全部散金贴
1楼 2011-02-24 14:30:31
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

zhangguangping

木虫 (著名写手)
★ ★ ★
小木虫(金币+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是不行的。
赞一下(2人) 回复此楼
7楼2011-03-01 13:41:59
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
查看全部 13 个回答

zhangguangping

木虫 (著名写手)
★
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的用法。
赞一下(1人) 回复此楼
2楼2011-02-24 15:26:57
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

minmin_0082003

金虫 (正式写手)
引用回帖:
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的回帖

zhangguangping

木虫 (著名写手)
★ ★ ★ ★
小木虫(金币+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
手册上有这个，所以要按照手册来好些！我不清楚你的那个写法管不管用，起码我没用过。我这样写的，都可以的。
赞一下(2人) 回复此楼
4楼2011-02-24 19:24:17
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
查看全部 13 个回答
普通表情 龙 兔 虎 猫 高级回复 (可上传附件)
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 315求调剂 +5 顺理成张 2026-04-03 6/300 2026-04-03 13:21 by 来看流星雨10
[考研] 315分 085602 求调剂 +15 26考研上岸版26 2026-04-02 15/750 2026-04-03 12:45 by xingguangj
[考研] 求材料调剂 一志愿南昌大学 328分 +4 yyy..... 2026-04-03 4/200 2026-04-03 12:31 by wxiongid
[考研] 309求调剂 +4 刘刘刘1231 2026-04-02 5/250 2026-04-03 12:04 by 1753564080
[考研] 08工科275分求调剂 +14 AaAa7420 2026-03-31 14/700 2026-04-03 11:13 by cocolv
[考研] 一志愿北京科技大学材料工程085601，求调剂 +12 cdyw 2026-04-02 12/600 2026-04-03 10:18 by 蓝云思雨
[考研] 复试调剂 +4 春日来信- 2026-04-03 4/200 2026-04-03 08:36 by JourneyLucky
[考研] 求调剂！生物与医药专硕 +4 逆转陆先生 2026-04-01 5/250 2026-04-03 08:33 by Jaylen.
[考研] 330分求调剂 +15 qzenlc 2026-03-29 15/750 2026-04-02 21:41 by 百灵童888
[考研] 326分求调剂 +3 于是乎呢 2026-04-01 3/150 2026-04-02 21:30 by jing_jenna
[考研] 285求调剂 +8 AZMK 2026-04-02 11/550 2026-04-02 20:16 by yulian1987
[考研] 298求B区调剂 +4 zzz，，r 2026-04-02 5/250 2026-04-02 12:17 by 土木硕士招生
[考研] 一志愿北交大材料工程总分358 +3 cs0106 2026-04-02 5/250 2026-04-02 11:37 by olim
[考研] 生物学327，求调剂 +5 书上的梅子 2026-04-01 6/300 2026-04-02 06:47 by ilovexiaobin
[考研] 求调剂 +4 图鉴212 2026-03-30 5/250 2026-04-01 15:32 by 图鉴212
[考研] 材料与化工调剂一志愿大连海事085600，349 +9 吃的不少 2026-03-30 9/450 2026-04-01 11:24 by wangjy2002
[考研] 288资源与环境专硕求调剂，不限专业，有学上就行 +25 lllllos 2026-03-30 26/1300 2026-04-01 09:52 by 一只好果子?
[考研] 080500-315分复试调剂 +9 上岸3821 2026-03-31 9/450 2026-03-31 17:29 by 唐沐儿
[考研] 求调剂 +10 家佳佳佳佳佳 2026-03-29 10/500 2026-03-30 18:34 by 544594351
[考研] 求佛 +7 迷人的哈哈 2026-03-28 7/350 2026-03-28 16:47 by 催化大白
信息提示
关闭
请填处理意见
关闭