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

24小时热门版块排行榜    

小木虫论坛-学术科研互动平台 » 计算模拟区 » 第一性原理 » Vasp&MedeA » 关于Si光学性质的计算
291/1返回列表
查看: 8054  |  回复: 28
只看楼主@他人 存档 新回复提醒 (忽略) 收藏    在APP中查看

liliangfang

荣誉版主 (著名写手)

[交流] 关于Si光学性质的计算

总结一下计算光学性质，欢迎大家进来交流
1.优化后DFT计算
INCAR
SYSTEM = Si-Diamond
ENCUT = 400
ISTART = 0
ICHARG = 2
PREC = Accurate
NSW = 200
IBRION = 2
ISIF = 2
EDIFF = 0.00001
EDIFFG = -0.001
LRPA = FALSE
LCORR = TRUE
IALGO = 48
LWAVE = TRUE
ISPIN = 2

KPOINTS
Automatic mesh
0
M
9   9   9
0   0   0

POSCAR
Si-Diamond:                             
   5.46900000000000     
     0.0000000000000000    0.5000000000000000    0.5000000000000000
     0.5000000000000000    0.0000000000000000    0.5000000000000000
     0.5000000000000000    0.5000000000000000    0.0000000000000000
   Si
   2
Direct
  0.8750000000000000  0.8750000000000000  0.8750000000000000
  0.1250000000000000  0.1250000000000000  0.1250000000000000

  0.00000000E+00  0.00000000E+00  0.00000000E+00
  0.00000000E+00  0.00000000E+00  0.00000000E+00
计算
2.HSE06
INCAR
SYSTEM = Si-Diamond
ENCUT = 400
ISTART = 1
ICHARG = 2
PREC = Accurate
NSW = 0
IBRION = -1
ISIF = 2
EDIFF = 0.00001
EDIFFG = -0.01
LRPA = FALSE
NPAR = 4
IALGO = 48
LHFCALC = TRUE
HFSCREEN = 0.2
ALGO = All
TIME = 0.4
ENCUTFOCK = 0
PRECFOCK = Fast
ISPIN = 2
其他文件不变从DFT中全部拷贝
vasp计算
3.optical
INCAR
SYSTEM = Si-Diamond
ENCUT = 400
ISTART = 1
ICHARG = 2
ISMEAR = -5
PREC = Accurate
NSW = 0
IBRION = -1
ISIF = 2
EDIFF = 0.00001
EDIFFG = -0.01
LRPA = FALSE
IALGO = 48
NBANDS = 24
NPAR = 1
CSHIFT = 0.1
LNABLA = .FALSE.
LEPSILON = .FALSE.
LRPA = .TRUE.              
LREAL = .FALSE.
LWAVE = .FALSE.   
LCHARG = .FALSE.
LOPTICS = TRUE
NBANDS增加到2-3倍，k点增加
KPOINTS
Automatic mesh
0
M
11   11   11
0     0     0
其他不变从HSE中全部拷贝
计算
从COTCAR中对出介电函数作图如下





欢迎大家加入讨论，给意见
另外谁又计算过金属的光学性质，大家交流下

[ Last edited by liliangfang on 2011-9-19 at 16:37 ]
回复此楼

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

仿真建模与计算 first principle vasp光学 计算化学-光学性质
模拟，理论 vasp 第一性原理 DFT+U，HSE等

» 本帖已获得的红花（最新10朵）

mazuju028

» 猜你喜欢

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

» 抢金币啦！回帖就可以得到:

查看全部散金贴
1楼2011-09-19 16:33:20
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
回帖支持 ( 显示支持度最高的前 50 名 )

liliangfang

荣誉版主 (著名写手)
引用回帖:
19楼: Originally posted by mazuju028 at 2011-11-17 10:09:49:
很好的帖子啊。
楼主计算光学的过程基本上是：先结构优化，再用HSE06静态计算，最后用静态计算的波函数及电荷密度进行光学计算。
我有几处小疑问想问一下，
1，第二部分静态计算时，ALGO=All，和IALGO=48有没有 ...

不好意思，这是个问题，当初发帖的时候没检查，我已经在另外一个帖子中更正了，两种我都试过，还是HSE好多了，带隙有修正。另外我一直尝试GW，当时在计算能带的时候一直出问题，所以后面光学的也就没进行，不知道版主有没有用过GW

[ Last edited by liliangfang on 2011-11-17 at 11:12 ]
赞一下(2人) 回复此楼
20楼2011-11-17 11:10:41
 已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

valenhou001

至尊木虫 (职业作家)
★ ★ ★
liliangfang(金币+1):谢谢参与
zzy870720z(金币+1): 谢谢提示 2011-09-20 12:29:28
zzy870720z(金币+1): 谢谢提示 2011-09-20 12:29:34
liliangfang(金币+1): 2011-09-20 18:40:18
liliangfang(金币+1): 2011-09-27 09:14:16
liliangfang(金币+2): 2011-11-02 08:43:11
此论坛有人讨论过vasp对金属的光学性质的计算。vasp里面还没有考虑带内的跃迁。台大的郭光宇老师他的组在对vasp的光学性质进行了一些扩展，对CNT的光学性质做了广泛的计算。你可在PRB上找找他的文章，看一下里面的处理方法。
赞一下(4人) 回复此楼
2楼2011-09-20 07:29:44
 已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

guolianshun

金虫 (著名写手)
★
liliangfang(金币+1):谢谢参与
引用回帖:
1楼: Originally posted by liliangfang at 2011-09-19 16:33:20:
总结一下计算光学性质，欢迎大家进来交流
1.优化后DFT计算
INCAR
SYSTEM = Si-Diamond
ENCUT = 400
ISTART = 0
ICHARG = 2
PREC = Accurate
NSW = 200
IBRION = 2
ISIF = 2
EDIFF = 0.00001
EDIFFG  ...

这些东西怎么讨论呀，计算出来有什么用呀。
赞一下(2人) 回复此楼
5楼2011-09-26 16:48:42
 已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

gougou531

禁虫 (正式写手)
★
liliangfang(金币+1):谢谢参与
本帖内容被屏蔽
16楼2011-09-28 08:24:33
 已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

mazuju028

木虫 (正式写手)
★
小木虫(金币+0.5):给个红包，谢谢回帖
送鲜花一朵
很好的帖子啊。
楼主计算光学的过程基本上是：先结构优化，再用HSE06静态计算，最后用静态计算的波函数及电荷密度进行光学计算。
我有几处小疑问想问一下，
1，第二部分静态计算时，ALGO=All，和IALGO=48有没有冲突？在vasp说明中，ALGO=All和IALGO=58好像是一起使用的。
2，楼主不用HSE06只用PBE计算过吗？结果如何？
2，楼主最后光学计算中LRPA设置了两个不同的值；

[ Last edited by mazuju028 on 2011-11-17 at 10:13 ]
赞一下(2人) 回复此楼
19楼2011-11-17 10:09:49
 已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
普通回帖

liliangfang

荣誉版主 (著名写手)
引用回帖:
2楼: Originally posted by valenhou001 at 2011-09-20 07:29:44:
此论坛有人讨论过vasp对金属的光学性质的计算。vasp里面还没有考虑带内的跃迁。台大的郭光宇老师他的组在对vasp的光学性质进行了一些扩展，对CNT的光学性质做了广泛的计算。你可在PRB上找找他的文章，看一下里面的 ...

5.2也需要修正吗
赞一下 回复此楼

» 本帖已获得的红花（最新10朵）

liqian16
3楼2011-09-20 18:41:06
 已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

liqian16

银虫 (初入文坛)
★
liliangfang(金币+1):谢谢参与
送鲜花一朵
liliangfang(金币+1): 加了 2011-09-27 09:15:10
引用回帖:
3楼: Originally posted by liliangfang at 2011-09-20 18:41:06:
5.2也需要修正吗

楼主你好，我也想用vasp计算半导体的光学性质，但是没有经验，可以加我QQ597006602具体交流一下吗
赞一下(1人) 回复此楼
4楼2011-09-26 16:46:27
 已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

liliangfang

荣誉版主 (著名写手)
引用回帖:
5楼: Originally posted by guolianshun at 2011-09-26 16:48:42:
这些东西怎么讨论呀，计算出来有什么用呀。

这是算介电函数的，计算一些光学性质，一直存在问题，所以拿出来讨论下
赞一下 回复此楼
6楼2011-09-27 09:15:54
 已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

guolianshun

金虫 (著名写手)
★
youzhizhe(金币+1): 谢谢交流。 2011-09-27 10:29:11
引用回帖:
6楼: Originally posted by liliangfang at 2011-09-27 09:15:54:
这是算介电函数的，计算一些光学性质，一直存在问题，所以拿出来讨论下

光学性质的横坐标如果是波长wavelength，纵坐标是吸收absorption，这样的图能反应光吸收吗？
赞一下(1人) 回复此楼
7楼2011-09-27 09:27:40
 已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

beyondstar

木虫 (正式写手)
★
liliangfang(金币+1):谢谢参与
VASP如何得到吸收谱，反射率，似乎只有介电函数，
赞一下(1人) 回复此楼
8楼2011-09-27 13:00:27
 已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

liliangfang

荣誉版主 (著名写手)
引用回帖:
7楼: Originally posted by guolianshun at 2011-09-27 09:27:40:
光学性质的横坐标如果是波长wavelength，纵坐标是吸收absorption，这样的图能反应光吸收吗？

这你就的看固体物理了
赞一下 回复此楼
9楼2011-09-27 14:38:59
 已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

liliangfang

荣誉版主 (著名写手)
引用回帖:
8楼: Originally posted by beyondstar at 2011-09-27 13:00:27:
VASP如何得到吸收谱，反射率，似乎只有介电函数，

公式推导
回复此楼
10楼2011-09-27 14:39:47
 已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

guolianshun

金虫 (著名写手)
引用回帖:
9楼: Originally posted by liliangfang at 2011-09-27 14:38:59:
这你就的看固体物理了

咨询问题，本人又不是固体物理专业，能不能简单的说说，大家互相帮助呀。
赞一下 回复此楼
11楼2011-09-27 15:41:26
 已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

beyondstar

木虫 (正式写手)
引用回帖:
10楼: Originally posted by liliangfang at 2011-09-27 14:39:47:
公式推导

是自己编程吗
回复此楼
12楼2011-09-27 15:55:02
 已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

liliangfang

荣誉版主 (著名写手)
引用回帖:
11楼: Originally posted by guolianshun at 2011-09-27 15:41:26:
咨询问题，本人又不是固体物理专业，能不能简单的说说，大家互相帮助呀。

我也非物理，一些计算光学性质的文章都有推导，如果你不是做这方面的，就没有必要了
赞一下 回复此楼
13楼2011-09-27 20:55:04
 已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

liliangfang

荣誉版主 (著名写手)
引用回帖:
12楼: Originally posted by beyondstar at 2011-09-27 15:55:02:
是自己编程吗

不是，好多光学性质的文章都有的，你可以看看
赞一下 回复此楼
14楼2011-09-27 20:55:54
 已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

beyondstar

木虫 (正式写手)
引用回帖:
14楼: Originally posted by liliangfang at 2011-09-27 20:55:54:
不是，好多光学性质的文章都有的，你可以看看

我是说VASP能直接给出这些结果吗，还是你自己编程后处理的，公式我都知道
赞一下 回复此楼
15楼2011-09-27 23:03:27
 已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

liliangfang

荣誉版主 (著名写手)
★
贺仪(金币+1): 多谢指教！ 2011-09-28 13:48:20
引用回帖:
15楼: Originally posted by beyondstar at 2011-09-27 23:03:27:
我是说VASP能直接给出这些结果吗，还是你自己编程后处理的，公式我都知道

计算可以得到介电函数，其他的公式推导
赞一下(1人) 回复此楼
17楼2011-09-28 09:08:55
 已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

liliangfang

荣誉版主 (著名写手)
★
zzy870720z(金币+1): 鼓励交流 2011-09-28 15:48:41
引用回帖:
16楼: Originally posted by gougou531 at 2011-09-28 08:24:33:
我想问一下：计算光学性质的流程包括哪几部分，应该先优化，需要静态计算吗？在流程上和计算dos和band有什么区别吗？或者想要计算一个体系的光学性质优化的时候需要哪些文件，静态计算时需要哪些文件，计算光学的 ...

这个你得先看vasp说明了，我也是自己摸索的，你可以先看下别人怎么做的，然后自己计算总结
赞一下(1人) 回复此楼
18楼2011-09-28 09:10:39
 已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

mazuju028

木虫 (正式写手)
★
小木虫(金币+0.5):给个红包，谢谢回帖
引用回帖:
20楼: Originally posted by liliangfang at 2011-11-17 11:10:41:
不好意思，这是个问题，当初发帖的时候没检查，我已经在另外一个帖子中更正了，两种我都试过，还是HSE好多了，带隙有修正。另外我一直尝试GW，当时在计算能带的时候一直出问题，所以后面光学的也就没进行， ...

呵呵，我也刚开始接触用vasp计算光学性质，GW还没有用过。以后还要多多请教。
对了如果用GW来算光学性质的话，是否也是直接在光学计算中用GW得到的波函数和电荷密度文件就行了？
赞一下(1人) 回复此楼
21楼2011-11-17 15:10:54
 已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

liliangfang

荣誉版主 (著名写手)
引用回帖:
21楼: Originally posted by mazuju028 at 2011-11-17 15:10:54:
呵呵，我也刚开始接触用vasp计算光学性质，GW还没有用过。以后还要多多请教。
对了如果用GW来算光学性质的话，是否也是直接在光学计算中用GW得到的波函数和电荷密度文件就行了？

我在尝试呢，不过老出毛病，目前在摸索
赞一下 回复此楼
22楼2011-11-17 15:44:51
 已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

mengnan0501

铜虫 (小有名气)
★
小木虫(金币+0.5):给个红包，谢谢回帖
说实话，我看完帖子后思维有点混乱，我有几点疑问，
1，LZ究竟是用HSE06优化Si的结构还是普通的DFT优化晶格结构的？如果是HSE06，那你的第一步为何是DFT跑的。
2，如果是用HSE，为何光学性质（第3部分）那儿你不加上HSE那一部分设置，
3，你在结构优化时为何用的是ISIF=2，这种固定体积和形状的优化其实并不好，全优化或固定体积优化应用更多。
赞一下(2人) 回复此楼
23楼2011-11-17 17:49:18
 已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

tsglss

铁杆木虫 (著名写手)
★
小木虫: 金币+0.5, 给个红包，谢谢回帖
学习了，以后肯定能用上
赞一下 回复此楼
24楼2013-05-12 16:17:52
 已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

拨云见日no1

新虫 (初入文坛)
引用回帖:
23楼: Originally posted by mengnan0501 at 2011-11-17 17:49:18
说实话，我看完帖子后思维有点混乱，我有几点疑问，
1，LZ究竟是用HSE06优化Si的结构还是普通的DFT优化晶格结构的？如果是HSE06，那你的第一步为何是DFT跑的。
2，如果是用HSE，为何光学性质（第3部分）那儿你不加 ...

同问
回复此楼
25楼2013-09-16 16:02:36
 已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

liliangfang

荣誉版主 (著名写手)
引用回帖:
25楼: Originally posted by 拨云见日no1 at 2013-09-16 16:02:36
同问...

请仔细看帖
回复此楼
26楼2013-09-17 08:23:14
 已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

匿名

用户注销 (职业作家)
★
小木虫: 金币+0.5, 给个红包，谢谢回帖
本帖仅楼主可见
赞一下
27楼2016-11-23 11:55:30
 已阅   申请1ST强帖   回复此楼   编辑   查看我的主页

黑挖掘机

新虫 (初入文坛)
★
小木虫: 金币+0.5, 给个红包，谢谢回帖
引用回帖:
2楼: Originally posted by valenhou001 at 2011-09-20 07:29:44
此论坛有人讨论过vasp对金属的光学性质的计算。vasp里面还没有考虑带内的跃迁。台大的郭光宇老师他的组在对vasp的光学性质进行了一些扩展，对CNT的光学性质做了广泛的计算。你可在PRB上找找他的文章，看一下里面的处 ...

您好，我想请问，我再计算光学性质的那一步加了LOPTICS = .TRUE.；CSHIFT = 0.1。用于计算单层的小体系（单胞14个原子），能计算出来实部与虚部的矩阵，但是我计算比较大的体系（单胞64个原子），当然相应的NBANDS我也有改，而就是计算不出来介电函数实部与虚部的矩阵，这是为什么呢？
赞一下 回复此楼
28楼2019-07-22 13:00:06
 已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

黑挖掘机

新虫 (初入文坛)
★
小木虫: 金币+0.5, 给个红包，谢谢回帖
楼主您好，我用相同的办法计算较小的体系（14个原子的单胞）可计算出介电函数实部与虚部矩阵，但是计算较大的体系（64个原子的单胞），就计算不出来介电函数实部与虚部的矩阵，生成的日志文件只有direction 1，而小体系有direction1，direction2与direction3。请问这是什么问题呢（当然NBANDS我有根据自洽的结果做相应改动）
赞一下 回复此楼
29楼2019-07-22 13:04:24
 已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 liliangfang 的主题更新
291/1返回列表
普通表情 龙 兔 虎 猫 高级回复(可上传附件) 论文翻译-EnPapers美助理教授团队-SCI英文润色
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考博] 这个博士要读吗 +5 Sea Breeze 2024-06-16 7/350 2024-06-17 01:23 by Elio-
[教师之家] 关于2023的收入 +35 小龙虾2008 2024-06-10 36/1800 2024-06-16 19:30 by guoguo0119
[考博] 34岁读博士晚吗 +27 emitdne 2024-06-13 27/1350 2024-06-16 18:57 by yugege2009
[找工作] 初始合伙人来啦！（生物试剂耗材标准品） +8 欢快的小科研人 2024-06-15 14/700 2024-06-16 16:00 by 欢快的小科研人
[基金申请] 博后基金，博管会会提前知道消息吗？ +4 yuyiang 2024-06-13 4/200 2024-06-16 11:40 by yangyuzhong4
[论文投稿] 编辑是什么意思 15+3 s090604054 2024-06-15 3/150 2024-06-16 10:29 by bobvan
[文学芳草园] 累并快乐着 +13 MYHLD521 2024-06-14 13/650 2024-06-15 22:59 by zeolitess
[基金申请] 关于博后基金的bug问题 +6 lxr1991 2024-06-14 9/450 2024-06-15 21:17 by since—2010
[论文投稿] 投稿时忘记修改一作 +7 gll123456 2024-06-13 11/550 2024-06-15 11:49 by gll123456
[基金申请] 博后基金，以往的结果点不开，怎么回事呢？最后一次机会了，两次都没中前面。 +7 kyukitu 2024-06-14 13/650 2024-06-15 06:46 by 我是王小帅
[基金申请] 面上基金有一个（两个）C是不是就没戏了.... 5+3 zzzm116 2024-06-13 17/850 2024-06-14 22:26 by lzt8076
[基金申请] 工材E口JQ有消息了吗 +4 babyduck 2024-06-11 4/200 2024-06-14 17:23 by firepick
[基金申请] E12面上申请 +4 汉风之遗 2024-06-13 4/200 2024-06-14 15:28 by 天外飞去来
[硕博家园] 关于硕博连读的一些疑问？ +4 Lwenter 2024-06-14 4/200 2024-06-14 14:32 by ou0551
[论文投稿] 最近写了一篇控制优化领域的文章，可以投哪里啊？有没有水一些的期刊推荐 +7 香瓜木香 2024-06-12 13/650 2024-06-14 07:05 by 香瓜木香
[论文投稿] 投稿后发现其他作者的邮箱填错了该怎么办呀 10+4 在飞的猪 2024-06-13 6/300 2024-06-14 04:45 by 小虫子咔咔
[基金申请] 博士后面上项目状态还是专家评审吗 10+9 Thatcheremu 2024-06-13 55/2750 2024-06-13 21:23 by 乌合麒麟
[硕博家园] 科研求助 +5 杲www 2024-06-12 6/300 2024-06-13 16:16 by 姓李名明
[考博] 博导选择 +3 bing85977 2024-06-12 3/150 2024-06-13 15:34 by 我是邱尧
[论文投稿] 摩擦磨损论文投稿 +3 jmysan 2024-06-12 3/150 2024-06-13 08:36 by 莱茵润色
信息提示
关闭
请填处理意见
关闭