24小时热门版块排行榜    

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

frank_zhan

金虫 (正式写手)

[求助] 怎样得到HOMO和LUMO轨道,并把它显示出来?

请教一下,怎样设置才能指定输出的轨道呢?
比如我需要得到 LUMO或者HOMO,我需要怎样设置inputpp文件呢?或者怎样从输出文件中识别出HOMO,LUMO?
谢谢.
回复此楼
frank_zhan,zy2zhan@gmail.com
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

043114076

木虫 (正式写手)

【答案】应助回帖

★ ★ ★ ★ ★ ★ ★ ★ ★ ★
感谢参与,应助指数 +1
frank_zhan: 金币+10, ★★★很有帮助, 谢谢,请问你用这个做过么?能否给个例子看看? 2012-07-26 10:12:07
看PP 模块中的选项,plot_num=7 的选项就是你想要的homo和lumo了,稍微的不同是pwscf输出的是/psi/^2。lsign 设置为.true. 。kpoint 选择gamma点。
8楼2012-07-26 09:26:48
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

043114076

木虫 (正式写手)

【答案】应助回帖

★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★
frank_zhan: 金币+20, ★★★很有帮助, 谢谢。还有为什么kpoint要设置为Gamma点呢?你有考虑过有自旋的情况么? 2012-07-26 13:41:12
&inputpp
    prefix  = 'yourprefix'
    outdir = './tmp/'
    filplot = 'potfile'
    plot_num= 7
    kpoint= 1
    kband = 45
    lsign = .true.
/
&plot
    nfile = 1
    filepp(1) = 'potfile'
    iflag = 3
    output_format = 6
    fileout = 'out.cub'
/
9楼2012-07-26 12:49:19
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

043114076

木虫 (正式写手)

引用回帖:
9楼: Originally posted by 043114076 at 2012-07-26 12:49:19
&inputpp
    prefix  = 'yourprefix'
    outdir = './tmp/'
    filplot = 'potfile'
    plot_num= 7
    kpoint= 1
    kband = 45
    lsign = .true.
/
&plot
    nfile = 1
    filepp ...

kpoint 是指你的第几个k点,一般设为gamma点
kband 是你的第几个能级,你应该看一下你的pw输出找出第几个为homo和lumo
最后生成的是cube文件,可以使用gaussview或者其他的做图软件打开
10楼2012-07-26 12:58:40
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

043114076

木虫 (正式写手)

【答案】应助回帖

★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★
frank_zhan: 金币+20, ★★★★★最佳答案, 谢谢你的指导。我做这个时间不是很长,还有很多要学习啊,呵呵。 2012-07-26 19:18:52
注意,在pw中设置的k点,不一定是gamma点。在pp中的kpoint意思是选择pw中第几个k点来做图。我的理解是:gamma点的波函数对应于实空间,其他k点在实空间的意义不详。
自选体系没有做过pp,不过应该是类似的。我看到了你在pw_forum 上的帖子,O2确实是有磁性的,你需要看看结构化学书。对于磁性体系的homo和lumo,和你要区分清楚。
12楼2012-07-26 15:17:18
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 frank_zhan 的主题更新
信息提示
请填处理意见