请选择 进入手机版 | 继续访问电脑版
第13届北京科音初级量子化学培训班将于10月5~8日于北京举办,请点击此链接查看详情。这是新人一次性正确、完整学习量子化学计算的最好、最快机会,能少走无数弯路,欢迎参加并相互转告!(已报满)

计算化学公社

 找回密码
 现在注册!
查看: 1046|回复: 15

[VASP] VASP计算发射光谱步骤求助

[复制链接]

26

帖子

0

威望

248

eV
积分
274

Level 3 能力者

发表于 2020-2-6 10:37:51 | 显示全部楼层 |阅读模式
大家好!最近再用VASP计算单层MoS2得荧光发射光谱,但是看了好多文献或者论坛上大家得讨论只有用VASP计算吸收光谱的,不知道用下面的步骤对不对。

VASP发射光谱计算流程

VASP发射光谱计算流程

谢谢大家~



964

帖子

0

威望

3377

eV
积分
4341

Level 6 (一方通行)

发表于 2020-2-6 11:10:13 | 显示全部楼层
垂直激发对应的是吸收光谱。
发射没戏。
欢迎关注“量子化学”公众号,也欢迎投稿。

122

帖子

0

威望

1683

eV
积分
1805

Level 5 (御坂)

发表于 2020-2-6 13:36:22 | 显示全部楼层
计算发射光谱需要优化出激发态构型

182

帖子

0

威望

2937

eV
积分
3119

Level 5 (御坂)

发表于 2020-2-6 17:39:42 | 显示全部楼层
您好,您介绍的主要是对于有机体系的计算,对于二维材料或者其他材料VASP能否计算有待考证,吸收光谱是通过计算导带价带的gap计算得到的吗?如果是的话,对于发射应该也是适用的,前提是直接带隙,且半高宽比较窄,类似于钙钛矿材料
If anything is impossible? It should be done!!

26

帖子

0

威望

248

eV
积分
274

Level 3 能力者

 楼主| 发表于 2020-2-6 18:09:39 | 显示全部楼层
本帖最后由 酱汁牛肉 于 2020-2-6 18:11 编辑
三石草祭 发表于 2020-2-6 17:39
您好,您介绍的主要是对于有机体系的计算,对于二维材料或者其他材料VASP能否计算有待考证,吸收光谱是通过 ...

您好!吸收光谱我是根据VASP官网介绍的TDDFT的BSE方法计算得来的
https://www.vasp.at/wiki/index.p ... ies_of_Si_using_BSE
他是通过计算optical response function从而得到的吸收光谱。

26

帖子

0

威望

248

eV
积分
274

Level 3 能力者

 楼主| 发表于 2020-2-6 18:41:37 | 显示全部楼层
pwzhou 发表于 2020-2-6 13:36
计算发射光谱需要优化出激发态构型

您好~我之前利用FERWE与FERDO固定电子态,尝试过进行激发态构型计算,我的体系是1X1的二维的MoS2体系,POSCAR如下:
Mo1 S2
   1.00000000000000
     3.1708696595022934    0.0000000045612923    0.0000000000000000
    -1.5854353218348112    2.7460533682335244    0.0000000000000000
     0.0000000000000000    0.0000000000000000   23.1245887822981189
   Mo   S
     1     2
Direct
  0.6666669999999968  0.3333330000000032  0.5000000000000000
  0.3333330000000032  0.6666669999999968  0.4321391015134495
  0.3333330000000032  0.6666669999999968  0.5678608984865505


然后进行计算的INCAR如下:
SYSTEM = MoS2

GGA = PE
IVDW =11

LREAL  = Auto
ENCUT  = 450
PREC   = Accurate
EDIFF  = 1e-5  
ALGO  = Normal
NELM=  200
     
ISMEAR =-2
ISPIN  = 2
FERWE = 11*1.0 1*0.0 1*1.0 15*0.0 #The syntax is BANDS*OCC
FERDO = 11*1.0 1*1.0 1*0.0 15*0.0

  EDIFFG =-0.02  
  NSW    = 500   
  IBRION = 2     
  NFREE=2


但是总是报错说:Error reading item 'FERWE' from file INCAR.
其中POTCAR使用的是Mo_pv(9个价电子)与S(6个价电子),总共有24个电子
不知道您知道是什么问题吗?

26

帖子

0

威望

248

eV
积分
274

Level 3 能力者

 楼主| 发表于 2020-2-6 18:42:33 | 显示全部楼层
zjxitcc 发表于 2020-2-6 11:10
垂直激发对应的是吸收光谱。
发射没戏。

您好~垂直发射能是不是可以计算发射光谱呀?

2016

帖子

3

威望

7616

eV
积分
9692

Level 6 (一方通行)

Ab Initio Amateur

发表于 2020-2-7 11:35:12 | 显示全部楼层
酱汁牛肉 发表于 2020-2-6 18:41
您好~我之前利用FERWE与FERDO固定电子态,尝试过进行激发态构型计算,我的体系是1X1的二维的MoS2体系,PO ...

1、MoS2,Mo_pv有9个价电子,S有6个价电子,该有9+6+6=21电子(最好确定一下POTCAR的顺序)。
2、周期性的电子态是由能带指标n和Bloch vector k两个来标记的,FERWE和FERDO使用的时候应该看一下到底是哪个能级激发到哪个能级,你可以参考https://cms.mpi.univie.ac.at/vas ... _SMEARINGS_tag.html,其中:
The band-index runs fastest.

可以先做SCF来看一下能级是怎么排的。
3、1x1的晶格来做是否合理我不太能够确定。

评分

参与人数 1eV +2 收起 理由
sobereva + 2

查看全部评分

满招损,谦受益。热衷于理论和方法研究水平不高但欢迎讨论。

26

帖子

0

威望

248

eV
积分
274

Level 3 能力者

 楼主| 发表于 2020-2-7 14:42:48 | 显示全部楼层
卡开发发 发表于 2020-2-7 11:35
1、MoS2,Mo_pv有9个价电子,S有6个价电子,该有9+6+6=21电子(最好确定一下POTCAR的顺序)。
2、周期性 ...

谢谢~~
不好意思我看错了Mo_pv是12个电子。
我是计算了SCF之后,从第13个能级开始时occupation为0,由此确定FERWE与FERDO的电子排布的。
谢谢老师的建议,我也觉得1X1有可能不太对,感觉有点小了,我再在4X4的表面上试一下~

91

帖子

0

威望

884

eV
积分
975

Level 4 (黑子)

发表于 2020-2-8 11:20:34 | 显示全部楼层
看到过用Castep计算固体荧光的,现在Castep学术免费了,可以尝试下!

26

帖子

0

威望

248

eV
积分
274

Level 3 能力者

 楼主| 发表于 2020-2-12 14:10:29 | 显示全部楼层
薛定谔的老猫 发表于 2020-2-8 11:20
看到过用Castep计算固体荧光的,现在Castep学术免费了,可以尝试下!

谢谢!!

27

帖子

0

威望

120

eV
积分
147

Level 2 能力者

发表于 2020-8-29 21:31:44 | 显示全部楼层
卡开发发 发表于 2020-2-7 11:35
1、MoS2,Mo_pv有9个价电子,S有6个价电子,该有9+6+6=21电子(最好确定一下POTCAR的顺序)。
2、周期性 ...

卡开发发前辈,关于用固定占据来近似激发态计算,我的理解是:
1. 计算基态能量,然后改变电子占据,不优化结构,直接算能量,那么这样的能量差可以认为垂直激发能;
2. 改变电子占据,然后优化结构,然后计算能量,然后减去基态能量,这样的能量差可以认为是发射的能量。
请问我的理解合适吗?
还有就是我看了你之前发的帖子,里面说要读取基态自洽的WAVECAR,那么不管是否需要做结构优化(改变电子占据后),都要读取基态自洽的WAVECAR吗?
期待回复,非常感谢。

2016

帖子

3

威望

7616

eV
积分
9692

Level 6 (一方通行)

Ab Initio Amateur

发表于 2020-8-29 22:13:40 | 显示全部楼层
Xian 发表于 2020-8-29 21:31
卡开发发前辈,关于用固定占据来近似激发态计算,我的理解是:
1. 计算基态能量,然后改变电子占据,不 ...

1、2按照deltaSCF的框架,确实如此,不过我尚不确定约束占据数的情况优化结构会不会有什么问题,最好做好详细的测试。
2、按道理说,读取应该不是必要的,只是可能会出现SCF比较难收敛的情况才这样做。
满招损,谦受益。热衷于理论和方法研究水平不高但欢迎讨论。

27

帖子

0

威望

120

eV
积分
147

Level 2 能力者

发表于 2020-8-31 19:07:43 | 显示全部楼层
卡开发发 发表于 2020-8-29 22:13
1、2按照deltaSCF的框架,确实如此,不过我尚不确定约束占据数的情况优化结构会不会有什么问题,最好做好 ...

非常感谢你的回复。我测试了一个小体系,发现电子步很难收敛,而且离子步跑了2步就有错误信息(Error EDDDAV: Call to ZHEGV failed. Returncode =   7 1   8),没法继续了

2016

帖子

3

威望

7616

eV
积分
9692

Level 6 (一方通行)

Ab Initio Amateur

发表于 2020-8-31 20:22:08 | 显示全部楼层
Xian 发表于 2020-8-31 19:07
非常感谢你的回复。我测试了一个小体系,发现电子步很难收敛,而且离子步跑了2步就有错误信息(Error EDD ...

试着减小初猜步长(POTIM)。
满招损,谦受益。热衷于理论和方法研究水平不高但欢迎讨论。
您需要登录后才可以回帖 登录 | 现在注册!

本版积分规则

手机版|北京科音自然科学研究中心|京公网安备 11010502035419号|计算化学公社 — 北京科音旗下高水平计算化学交流论坛 ( 京ICP备14038949-1号 )

GMT+8, 2020-9-19 00:09 , Processed in 0.155686 second(s), 28 queries .

快速回复 返回顶部 返回列表