24小时热门版块排行榜    

查看: 3646  |  回复: 13
本帖产生 1 个 1ST强帖 ,点击这里进行查看

mazuju028

木虫 (正式写手)

[求助] vasp5.2光学计算中怎么加入scissor operation?

如题,有没有人知道vasp5.2光学计算中怎么加入scissor operation? 有参数吗?

[ Last edited by mazuju028 on 2011-11-17 at 21:54 ]
回复此楼

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

第一性原理 vasp光学 计算化学-光学性质 纳米及第一性原理方向
催化

» 猜你喜欢

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

锐气藏于胸,和气浮于脸,才气见于事,义气施与人
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
回帖支持 ( 显示支持度最高的前 50 名 )

buct2010

金虫 (正式写手)

【答案】应助回帖

★ ★ ★
mazuju028(金币+5): 哈哈,太感谢你了。 2011-11-18 11:08:25
zzy870720z(金币+3): 谢谢侯老师指点 2011-11-18 12:11:19
zzy870720z(1ST强帖+1): 很详细 2011-11-19 09:41:56
引用回帖:
1楼: Originally posted by mazuju028 at 2011-11-17 21:51:33:
如题,有没有人知道vasp5.2光学计算中怎么加入scissor operation? 有参数吗?
[ Last edited by mazuju028 on 2011-11-17 at 21:54 ]

这个scissor operation就是把VASP计算得到的禁带值修正到实验值,比如计算的值为1.68 实验值3.0,那就自己采用scissor operation 1.32 ,把计算得到的禁带值修正到实验值(它只对结果修正,没有控制的输入参数)。还有要得到好的光吸收谱,首先要保证你的能带是准确的,用1.68的能带算出来的吸收谱肯定与实验值不符,用了scissor operation也只是把结果人为地修正一下,建议你采用准确的能带结构计算方法,对于过渡金属采用LDA+U(或GGA+U),其他体系可以用杂化泛函计算(GW,HSE06)计算。现在计算的文章也不好发,要用新方法来算才能被认可,我今年的一篇文章用LDA算的,用scissor operation做了人为修正,审稿人说这个scissor operation不合理,我们就用LDA+U把所有的能带重算了一遍。
2楼2011-11-18 08:24:22
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
普通回帖

fzx2008

荣誉版主 (著名写手)

优秀版主优秀版主

【答案】应助回帖

mazuju028(金币+2): 嗯啊,不过abinit的是基于DFPT算的,计算量很大啊。 2011-11-18 11:09:04
马哥,又开始用vasp啦?ABINIT中有个sciss的参数
3楼2011-11-18 10:06:42
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

mazuju028

木虫 (正式写手)

引用回帖:
3楼: Originally posted by fzx2008 at 2011-11-18 10:06:42:
马哥,又开始用vasp啦?ABINIT中有个sciss的参数

哈哈,abinit太慢了啊。尤其是光学性质
锐气藏于胸,和气浮于脸,才气见于事,义气施与人
4楼2011-11-18 11:07:19
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

mazuju028

木虫 (正式写手)

引用回帖:
2楼: Originally posted by buct2010 at 2011-11-18 08:24:22:
这个scissor operation就是把VASP计算得到的禁带值修正到实验值,比如计算的值为1.68 实验值3.0,那就自己采用scissor operation 1.32 ,把计算得到的禁带值修正到实验值(它只对结果修正,没有控制的输入参数) ...

谢谢啊,不过我还有个疑问,就是如果一定要用剪刀算符,在vasp5.2中是怎么实现的?(vasp4.6我知道要另外的一个输入文件OPTCRT)。在vasp5.2中不需要OPTCRT,直接自洽计算完,用前面的波函数(WAVCAR)和电荷密度(CHG*),进行光学计算就可以得到介电函数,如果要人为调整禁带宽度,应该在哪一步对什么文件进行调整,调整之后怎么用到光学计算中?呵呵,不知道你明白我的意思没有?
锐气藏于胸,和气浮于脸,才气见于事,义气施与人
5楼2011-11-18 11:19:15
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

buct2010

金虫 (正式写手)

【答案】应助回帖

★ ★
zzy870720z(金币+2): 谢谢指点 2011-11-19 09:42:08
mazuju028(金币+3): 谢谢指点 2011-11-19 17:00:40
引用回帖:
5楼: Originally posted by mazuju028 at 2011-11-18 11:19:15:
谢谢啊,不过我还有个疑问,就是如果一定要用剪刀算符,在vasp5.2中是怎么实现的?(vasp4.6我知道要另外的一个输入文件OPTCRT)。在vasp5.2中不需要OPTCRT,直接自洽计算完,用前面的波函数(WAVCAR)和电荷密度 ...

光学性质得到的是光子能量(eV)和介电常数虚部(吸收系数)的对应关系,你可以吧吸收带边右移剪刀算符,比如你没修正前,光子能量在1.68 eV开始吸收,修正量是1.32eV ,那你就把吸收谱整个右移1.32eV (光子能量,即横坐标减去1.32eV),使吸收边带移到3 eV。注意剪刀算符的大小是用能带结构确定的,吸收带的带边一般会比禁带宽度大一些(光学带宽大于禁带宽)。
6楼2011-11-19 09:34:46
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

mazuju028

木虫 (正式写手)

引用回帖:
6楼: Originally posted by buct2010 at 2011-11-19 09:34:46:
光学性质得到的是光子能量(eV)和介电常数虚部(吸收系数)的对应关系,你可以吧吸收带边右移剪刀算符,比如你没修正前,光子能量在1.68 eV开始吸收,修正量是1.32eV ,那你就把吸收谱整个右移1.32eV (光子能量 ...

用HSE折腾了一阵子,发现实在耗不起那时间。因为我的体系中没有3d,4f过渡金属元素,所以不能使用+U的方法。可能还得回到用剪刀差的方法上来。我还有一个疑问,如果我用剪刀差的话,对于介电函数实部需不需要做些处理?对于虚部,简单的把吸收边右移,可是实部呢?因为我要算折射率,所以同时需要虚部和实部的数据。希望指点一下。
锐气藏于胸,和气浮于脸,才气见于事,义气施与人
7楼2011-12-14 15:32:37
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

buct2010

金虫 (正式写手)

★ ★ ★
sunyang1988(金币+1): 谢谢提示 2011-12-14 17:13:10
mazuju028(金币+2): 谢谢交流 2011-12-14 17:29:14
引用回帖:
7楼: Originally posted by mazuju028 at 2011-12-14 15:32:37:
用HSE折腾了一阵子,发现实在耗不起那时间。因为我的体系中没有3d,4f过渡金属元素,所以不能使用+U的方法。可能还得回到用剪刀差的方法上来。我还有一个疑问,如果我用剪刀差的话,对于介电函数实部需不需要 ...

你所有的数据只要把对应的频率平移就可以了,不用再处理介电常数的实部和虚部
8楼2011-12-14 16:31:27
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

mazuju028

木虫 (正式写手)

引用回帖:
8楼: Originally posted by buct2010 at 2011-12-14 16:31:27:
你所有的数据只要把对应的频率平移就可以了,不用再处理介电常数的实部和虚部

我用MS做了个简单的对比,发现除了介电函数虚部是简单平移外,其它光学参数都不是简单平移。
根据Kramers–Kronig relations,由介电函数虚部得到实部,

从上式可以看出,介电函数实部,以及与介电函数实部有关的光学量都不是简单的平移。
锐气藏于胸,和气浮于脸,才气见于事,义气施与人
9楼2011-12-14 17:28:44
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

buct2010

金虫 (正式写手)


fzx2008(金币+1): 谢谢交流 2011-12-15 19:13:54
引用回帖:
9楼: Originally posted by mazuju028 at 2011-12-14 17:28:44:
我用MS做了个简单的对比,发现除了介电函数虚部是简单平移外,其它光学参数都不是简单平移。
根据Kramers–Kronig relations,由介电函数虚部得到实部,

从上式 ...

那就是不同情况用不同方法处理?我看的人家算光吸收谱的文献,明确说:“ here, we have used the scissor approximations to fit the calculated absorption edges to the experimental values of 3.2 and 3.37 eV for TiO2 and...”确实是算完吸收谱再平移的,请看 Appl. Phys. Lett. 92, 252104 (2008) 第三幅图,他也是用CASTEP算的光吸收谱。
10楼2011-12-15 09:30:38
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 mazuju028 的主题更新
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考博] 请问大家,我有希望进浙大普博么 +5 高艺文 2024-09-28 8/400 2024-09-29 00:08 by 小红豆
[教师之家] 东南大学电气学院电力电子方向招收硕士博士 +4 蜡笔小鑫1989 2024-09-26 7/350 2024-09-28 21:16 by 伍文童
[论文投稿] 求推荐医学类杂志 35+4 wshxtim1 2024-09-25 4/200 2024-09-28 13:56 by Wormaciae
[基金申请] 山东省基金公示结束 +6 ishow19 2024-09-24 14/700 2024-09-28 12:20 by ishow19
[电化学] 刚装好的电池电压掉很快怎么解决?? +3 花花你最棒 2024-09-26 6/300 2024-09-28 08:20 by cnofc
[基金申请] 广东省基金 +5 yingyueabc 2024-09-22 5/250 2024-09-27 22:57 by FengjiangLiu
[考博] 材料/电信/生物-2025普博生自荐-985本双非硕一区一作 +4 enowei0127 2024-09-23 8/400 2024-09-27 22:40 by enowei0127
[育儿交流] 熊孩子睡觉前很多事,啥原因? +7 quan2153 2024-09-25 14/700 2024-09-27 16:46 by quan2153
[基金申请] 请问大家的计划书填写列表中状态更新了吗? +7 Laker610 2024-09-25 9/450 2024-09-27 16:36 by 田田hj
[硕博家园] 当前读博士还有性价比吗? +8 苏东坡二世 2024-09-23 8/400 2024-09-27 11:39 by 半生梦君
[考博] 数学博导 +4 学术霸王 2024-09-25 6/300 2024-09-27 11:14 by 青古
[有机交流] 二氯甲烷的去除 +3 cgsa吧 2024-09-24 8/400 2024-09-27 10:55 by bear2007
[考博] 电催化,电池方向学生自荐 30+3 Freelancedd 2024-09-24 10/500 2024-09-27 10:27 by 莱茵润色
[论文投稿] 为什么审稿人一直拒绝审稿啊? +7 潇湘雪月 2024-09-24 13/650 2024-09-27 07:53 by weizhi111
[考博] 有机转材料申博难度 30+3 昕散 2024-09-22 10/500 2024-09-26 22:42 by 824282658
[有机交流] 请问胺的盐酸盐中氯化氢的氢会在核磁氢谱中出峰吗? +3 rommel1975 2024-09-25 4/200 2024-09-26 13:07 by 091602
[论文投稿] laser physics期刊投稿 5+3 mengxiangcz 2024-09-23 6/300 2024-09-25 15:08 by mengxiangcz
[有机交流] chemdraw结构式复制到word变形 +5 笑看人生1993 2024-09-22 11/550 2024-09-25 14:44 by 冰蓝夜游神
[论文投稿] SCI投稿状态 25+4 jorden8 2024-09-23 6/300 2024-09-25 10:24 by jorden8
[论文投稿] 爱斯维尔旗下Postharvest Biology and Technology,采用什么系统查查 20+4 sdsdsxh 2024-09-22 9/450 2024-09-25 07:33 by sdsdsxh
信息提示
请填处理意见