24小时热门版块排行榜    

查看: 4707  |  回复: 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 的主题更新
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[基金申请] 投票:  有多少人是今天查系统知道结果的? +14 爱看书的可乐 2026-08-26 16/800 2026-08-26 17:44 by skyjudy
[基金申请] 范进中举一文的中心思想 +9 炎黄贵胄 2026-08-22 10/500 2026-08-26 15:40 by semaglutide
[基金申请] 2026年的国家社科基金项目通讯评审的新规则与新动向、新挑战 +6 process2012 2026-08-23 9/450 2026-08-26 15:33 by process2012
[基金申请] 2026年8月25日国自然放榜前突然收到列入评审专家邮件,有关系吗? +25 木水思豆 2026-08-25 28/1400 2026-08-26 14:53 by draco1987
[基金申请] 能否退出参与的面上项目解除限项 +23 koalala 2026-08-24 26/1300 2026-08-26 14:29 by 宝贝虫子
[基金申请] 出来了 +9 trojank 2026-08-26 9/450 2026-08-26 14:25 by 宝贝虫子
[基金申请] 哪位高人中了,把查询到的截图贴出来让我看看,让我长长见识 +4 yuleib84 2026-08-26 5/250 2026-08-26 13:40 by yuleib84
[基金申请] 为什么国自然不能直接公布 +4 bjdxyxy 2026-08-26 4/200 2026-08-26 13:12 by qingmu1201
[基金申请] 科研孤儿太难了 +21 我4大白菜 2026-08-20 22/1100 2026-08-26 11:29 by nlgza
[基金申请] 国际合作可查了,中了面上 +18 Ldrop2023 2026-08-26 18/900 2026-08-26 11:15 by cmrandy
[基金申请] 系统进不去 +4 yanglien 2026-08-26 5/250 2026-08-26 11:10 by wenfengw83
[基金申请] 怎么看青基中了没有啊 +4 叶九微 2026-08-26 4/200 2026-08-26 10:52 by xiacongxu
[基金申请] 在坚冰还盖着北海的时候,我看到了怒放的梅花。 (金币+10) +6 ziyangfang 2026-08-25 9/450 2026-08-25 20:26 by huagongfeihu
[基金申请] 今日不放榜?网传国自然预计 8 月 27 日可查结果 +17 医学老男孩 2026-08-20 22/1100 2026-08-25 15:36 by 医学老男孩
[基金申请] 如果此刻你正在为国基感到焦虑,不妨来听听这首《基金之外》 +8 scalable 2026-08-24 8/400 2026-08-25 12:52 by jnhyjjm
[基金申请] 人气不行了 +11 fansofjerry 2026-08-21 11/550 2026-08-25 11:04 by 孤独的英雄6
[教师之家] 导师吐槽:我怎么摊上了这么个极品研究生! +3 苏东坡二世 2026-08-23 3/150 2026-08-25 10:35 by shisan1313
[基金申请] 建议基金发布提前给出明确的时间点 +13 kulium 2026-08-21 16/800 2026-08-24 16:27 by superceng
[基金申请] 让我中一个面上吧! +13 大萍1987 2026-08-20 16/800 2026-08-24 10:23 by 太傻了
[基金申请] 看来今天不会放榜了? +8 chengyan1220 2026-08-21 11/550 2026-08-21 17:52 by dcqxinyang
信息提示
请填处理意见