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

24小时热门版块排行榜    

小木虫论坛-学术科研互动平台 » 计算模拟区 » 第一性原理 » Vasp&MedeA » 关于Si光学性质的计算
291/1返回列表
查看: 8076  |  回复: 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返回列表
普通表情 龙 兔 虎 猫 高级回复(可上传附件)
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[基金申请] 省基金对国自然申请是有利还是有弊 +4 zycumt04 2024-06-26 5/250 2024-06-26 11:47 by msjy
[找工作] 老博士选高校 5+3 woodhead8386 2024-06-24 5/250 2024-06-26 11:13 by 516638798
[教师之家] 高校辞职，要求赔偿，这到底合不合理 +23 传动_海神 2024-06-23 35/1750 2024-06-26 10:30 by 阳光道路下
[有机交流] 做什么表征可以检测塑料中碳氢氧的含量 5+3 pzr的sci之路 2024-06-25 5/250 2024-06-26 10:13 by zqdsb
[硕博家园] 博士该不该读？ +6 L1009225316 2024-06-25 6/300 2024-06-26 08:41 by 陶渊明2024
[论文投稿] 我在写关于多维空间的论文，希望能与大家交流，准备在science上发表 +7 lgf519 2024-06-24 9/450 2024-06-26 08:29 by xxxxxxOdis
[考博] 一名额难求啊 +14 pinbo拼搏 2024-06-19 14/700 2024-06-25 16:55 by wjykycg
[基金申请] 基你太美 +13 lyfbangong 2024-06-24 14/700 2024-06-25 14:16 by 华氏911.
[硕博家园] 难道我真的要放弃吗？ +41 133456 2024-06-20 43/2150 2024-06-24 21:57 by sza999
[金属] EBSD的解析率只有10% +3 wallace6666 2024-06-20 7/350 2024-06-24 16:52 by wallace6666
[金属] 寻找钛合金热压缩代做 +4 liuyang358 2024-06-23 4/200 2024-06-24 13:50 by 搬砖狗不放弃
[基金申请] 国自然资助比率是不是要下降了？？ +8 今晚推荐22 2024-06-21 12/600 2024-06-24 11:15 by Pickfoot
[基金申请] 青年和面上，哪个上会难度更大 +12 今晚推荐22 2024-06-21 18/900 2024-06-24 11:08 by 半简体
[有机交流] 三颈瓶加热 5+4 wzjwx888 2024-06-20 4/200 2024-06-24 10:11 by w87437985
[基金申请] b口yq上会情况 +5 XUAAAAA 2024-06-19 6/300 2024-06-24 01:13 by sunrise3148
[基金申请] 江苏省333人才工程 出校后被刷的概率大不大？ +8 maxbirdzhang 2024-06-19 12/600 2024-06-23 20:39 by 查拉图斯特拉
[基金申请] 请问三个字的中文名字，在文中到底怎么表示？ 5+3 zhongyuan764 2024-06-22 6/300 2024-06-23 17:42 by 鱼翔浅底1
[论文投稿] ACS 编辑的意见 10+3 哈哈妞1993 2024-06-20 3/150 2024-06-21 17:06 by 投个论文
[有机交流] 怎么萃取出锡盐内包裹的化合物 +4 硕六过 2024-06-19 5/250 2024-06-21 09:50 by 光超嘟嘟
[论文投稿] 第一篇论文投稿进程记录 +4 慎独的小花卷 2024-06-20 9/450 2024-06-20 20:37 by 慎独的小花卷
信息提示
关闭
请填处理意见
关闭