VASPKIT论坛

 找回密码
 立即注册

QQ登录

只需一步,快速开始

搜索
热搜: VASPKIT
查看: 27420|回复: 48

VASP+VASPKIT基于玻尔兹曼输运理论计算输运性质

[复制链接]

35

主题

629

帖子

3794

积分

管理员

Rank: 9Rank: 9Rank: 9

积分
3794
QQ
发表于 2021-10-12 11:33:39 | 显示全部楼层 |阅读模式
VASPKIT 1.3.1及以后版本支持基于半经验玻尔兹曼输运理论计算输运性质,当前版本可输出electronic conductivity, Seebeck coefficient, power factor 和electronic thermal conductivity。实现算法见文献【Computer Physics Communications 175 (2006) 67–71】。详细教程见http://vaspkit.cn/index.php/67.html
采用VASP+VASPKIT计算输运性质步骤如下:######################################################################################################
1) 准备优化好的POSCAR结构,POTCAR文件,静态计算INCAR;
2) 调用vaspkit -task 681产生KPOINTS,注意输运性质计算需要非常密集的K-mesh;
3) 执行VASP计算;
4) 准备VPKIT.in文件(VPKIT.in文件在未来版本中将被重名为INPUT.in),包含以下参数

1           # 1 for gradient method with constant relaxation time approximation         
3D          # 2D for slab, 3D for bulk                                                                     
-2.0  2.0   # Minimum and maximum values of chemical potential (float type, in units of eV), referenced to the Fermi energy
1000        # Number of intervals between the minimum and maximum chemical potential values
300.0       # Temperature (float type, in units of K)
1.0          # Relaxation time (float type, in units of s), Set 1.0 if you don't know the exact relaxation time
0.0          # Desired gap value in scissor correction (float type, in units of V), Zero means make no correction, available only non-spin-polarization calculation
1             # Finite difference method (integer type)


5) 调用vaspkit -task 682计算输运性质。
######################################################################################################

作为验证我们选择Si和GaCuS2两个体系,计算结果与Openmx+BoltzTrap得到的结果符合很好(见下图)。欢迎大家跟帖反馈bug及提出改进建议,该功能今后将不断完善和扩充。




boltzman.png
若使用VASPKIT软件请引用文献 Computer Physics Communications 267, 108033, (2021) ,您的支持是对我们付出的最大肯定。
VASPKIT英文网址https://vaspkit.com, 中文网址http://vaspkit.cn, 中文论坛http://bbs.vaspkit.cn, QQ学术交流群:331895604和498683602, 微信公众号:VASPKIT
回复

使用道具 举报

1

主题

7

帖子

31

积分

新手上路

Rank: 1

积分
31
发表于 2021-10-12 16:37:09 | 显示全部楼层
本帖最后由 binshao1991 于 2021-10-12 16:41 编辑

对比了下数量级,vaspkit处理出来的数据是自动乘上了弛豫时间了?  在哪里指定的?还是说自动一起计算出来了,所以目前只能处理体材料的?
回复

使用道具 举报

35

主题

629

帖子

3794

积分

管理员

Rank: 9Rank: 9Rank: 9

积分
3794
QQ
 楼主| 发表于 2021-10-12 16:44:39 | 显示全部楼层
binshao1991 发表于 2021-10-12 16:37
对比了下数量级,vaspkit处理出来的数据是自动乘上了弛豫时间了?  在哪里指定的?还是说自动一起计算出来 ...

VPKIT.in文件中第6行可以指定弛豫时间,默认为1,目前还不能计算弛豫时间,今后版本会考虑增加。
若使用VASPKIT软件请引用文献 Computer Physics Communications 267, 108033, (2021) ,您的支持是对我们付出的最大肯定。
VASPKIT英文网址https://vaspkit.com, 中文网址http://vaspkit.cn, 中文论坛http://bbs.vaspkit.cn, QQ学术交流群:331895604和498683602, 微信公众号:VASPKIT
回复

使用道具 举报

0

主题

1

帖子

6

积分

新手上路

Rank: 1

积分
6
发表于 2021-10-13 10:29:04 | 显示全部楼层
vaspkit 发表于 2021-10-12 16:44
VPKIT.in文件中第5行可以指定弛豫时间,默认为1,目前还不能计算弛豫时间,今后版本会考虑增加。 ...

请问这个计算是基于哪种模型?对于各向异性较强的结构准确性如何
回复

使用道具 举报

35

主题

629

帖子

3794

积分

管理员

Rank: 9Rank: 9Rank: 9

积分
3794
QQ
 楼主| 发表于 2021-10-13 10:34:03 | 显示全部楼层
1126191527 发表于 2021-10-13 10:29
请问这个计算是基于哪种模型?对于各向异性较强的结构准确性如何

实现算法见文献Computer Physics Communications 175 (2006) 67–71。该功能将在后期版本中不断改进和增强。
若使用VASPKIT软件请引用文献 Computer Physics Communications 267, 108033, (2021) ,您的支持是对我们付出的最大肯定。
VASPKIT英文网址https://vaspkit.com, 中文网址http://vaspkit.cn, 中文论坛http://bbs.vaspkit.cn, QQ学术交流群:331895604和498683602, 微信公众号:VASPKIT
回复

使用道具 举报

0

主题

1

帖子

6

积分

新手上路

Rank: 1

积分
6
发表于 2021-10-14 09:52:40 | 显示全部楼层
Carrier Mobility这个怎么处理?
回复

使用道具 举报

35

主题

629

帖子

3794

积分

管理员

Rank: 9Rank: 9Rank: 9

积分
3794
QQ
 楼主| 发表于 2021-10-14 11:54:43 | 显示全部楼层
1769915666 发表于 2021-10-14 09:52
Carrier Mobility这个怎么处理?

玻尔兹曼输运理论不涉及到Carrier Mobility计算。
若使用VASPKIT软件请引用文献 Computer Physics Communications 267, 108033, (2021) ,您的支持是对我们付出的最大肯定。
VASPKIT英文网址https://vaspkit.com, 中文网址http://vaspkit.cn, 中文论坛http://bbs.vaspkit.cn, QQ学术交流群:331895604和498683602, 微信公众号:VASPKIT
回复

使用道具 举报

0

主题

1

帖子

18

积分

新手上路

Rank: 1

积分
18
发表于 2021-10-14 13:01:51 | 显示全部楼层
我想请问,VASP计算得到的费米能不是一定的吗,为什么横坐标的化学势要取一个范围?
回复

使用道具 举报

35

主题

629

帖子

3794

积分

管理员

Rank: 9Rank: 9Rank: 9

积分
3794
QQ
 楼主| 发表于 2021-10-14 15:21:15 | 显示全部楼层
qwsy007 发表于 2021-10-14 13:01
我想请问,VASP计算得到的费米能不是一定的吗,为什么横坐标的化学势要取一个范围? ...

对半导体材料进行p型掺杂费米能可低于价带顶,进行n型费米能级可高于导带底,因此化学势的变化范围可以反映掺杂浓度。
若使用VASPKIT软件请引用文献 Computer Physics Communications 267, 108033, (2021) ,您的支持是对我们付出的最大肯定。
VASPKIT英文网址https://vaspkit.com, 中文网址http://vaspkit.cn, 中文论坛http://bbs.vaspkit.cn, QQ学术交流群:331895604和498683602, 微信公众号:VASPKIT
回复

使用道具 举报

0

主题

7

帖子

30

积分

新手上路

Rank: 1

积分
30
发表于 2021-10-20 19:40:39 | 显示全部楼层
非常感谢这么好用的功能,但是发现个小问题,还请楼主修正下,期待能够改正这个bug。
Carrier concentration 貌似有问题,concentratoin为负数,并且是一个恒定的值。
回复

使用道具 举报

您需要登录后才可以回帖 登录 | 立即注册

本版积分规则

QQ|Archiver|手机版|小黑屋|VASPKIT论坛 ( 陕ICP备2020017291号 )

GMT+8, 2024-11-13 01:03 , Processed in 0.088994 second(s), 25 queries .

Powered by Discuz! X3.4

Copyright © 2001-2021, Tencent Cloud.

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