版块导航: 正在加载中...

登录注册

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

24小时热门版块排行榜

南方科技大学公共卫生及应急管理学院2025级博士研究生招生报考通知

返回列表

当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖

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

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

[ Last edited by minmin_0082003 on 2011-3-2 at 09:20 ]

回复此楼

» 抢金币啦！回帖就可以得到:
查看全部散金贴

1楼2011-02-24 14:30:31

已阅回复此楼关注TA 给TA发消息送TA红花 TA的回帖

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的回帖

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的回帖

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的回帖

zhangguangping

木虫 (著名写手)

★
小木虫(金币+0.5):给个红包，谢谢回帖交流
minmin_0082003(金币+10): 谢谢 2011-03-01 14:52:52

引用回帖:

Originally posted by minmin_0082003 at 2011-03-01 07:23:49:
估计问题出在这里吧
WriteWaveFunctions T
与%block WaveFuncKPoints
0.000 0.000 0.000 from 1 to 14
%endblock WaveFuncKPoints
一起设置的情况下，会出现所有的波函数。
不设置WriteW ...

这个和版本没有关系。
建议最好WriteWaveFunctions和COOP.write混合用，如果你只想得到一部分波函数的话。

赞一下(1人)

回复此楼

9楼2011-03-01 14:27:52

已阅回复此楼关注TA 给TA发消息送TA红花 TA的回帖