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

24小时热门版块排行榜    

小木虫论坛-学术科研互动平台 » 计算模拟区 » 第一性原理 » Vasp&MedeA » vasp5.2光学计算中怎么加入scissor operation?
141/2返回列表12下一页
查看: 4110  |  回复: 13
只看楼主@他人 存档 新回复提醒 (忽略) 收藏    在APP中查看
本帖产生 1 个 1ST强帖 ，点击这里进行查看

mazuju028

木虫 (正式写手)

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

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

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

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

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

» 猜你喜欢

» 本主题相关价值贴推荐，对您同样有帮助:
锐气藏于胸，和气浮于脸，才气见于事，义气施与人
1楼2011-11-17 21:51:33
已阅   回复此楼   关注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把所有的能带重算了一遍。
赞一下(3人) 回复此楼
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。注意剪刀算符的大小是用能带结构确定的，吸收带的带边一般会比禁带宽度大一些（光学带宽大于禁带宽）。
赞一下(1人) 回复此楼
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的方法。可能还得回到用剪刀差的方法上来。我还有一个疑问，如果我用剪刀差的话，对于介电函数实部需不需要 ...

你所有的数据只要把对应的频率平移就可以了，不用再处理介电常数的实部和虚部
赞一下(2人) 回复此楼
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算的光吸收谱。
赞一下(1人) 回复此楼
10楼2011-12-15 09:30:38
 已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 mazuju028 的主题更新
141/2返回列表12下一页
普通表情 龙 兔 虎 猫
信息提示
关闭
请填处理意见
关闭