bader分析(简易版)

更新时间:2023-06-09 01:29:19 阅读: 评论:0

1.计算。在静态计算中,加入以下参数,其中NGX参数与POSCAR中的晶格参数成比例,太小bader计算不会收敛,所以,要收敛之后,结果才可靠。
#bader charge
LAECHG=.TRUE.
NGX=100
NGY=100
NGZ=200
NGXF=150
NGYF=150
NGZF=300
2.^
3.算完之后,生成文件AECCAR0 AECCAR1 AECCAR2
将chgsum脚本上传到计算文件夹
将其变成可执行脚本:chmod 777
然后使用命令:./ AECCAR0 AECCAR2
生成CHARGCAR_sum文件
4.将chgsum脚本上传到计算文件夹
然后将bader变成可执行文件: chmod 777 bader
输入命令:./bader CHGCAR -ref CHGCAR_sum
生成文件:   
就包含价电子的电荷信息:
5.
6.判断计算的可靠性
把文件拖到Origin中,CHARGE这一列的数字就是每个原子对应的电荷,而最下面那一行,NUMBER OF ELECTRONS:126,代表整个体系的总电量有126e。
我们可以通过对比总电量和所有原子加和电量之间的差值来判断是否收敛:
我们把CHGCAR这一列数据复制到Excel中,选中所有数据,按求和符号,或者在另一个格子中输入=sum(),Enter,得到总电荷是
126 -  =  (e)        / 126 =  %
我们发现:CHGCAR这一列数据这和+VACUUM CHARGE这个数值=NUMBER OF ELECTRONS
VACUUM CHARGE这个值是真空的电荷,如果这个值不为零,就意味这个有电子转移到真空中
标准的情况应该是:
#
所以,判断bader计算是否收敛,也可以通过判断VACUUM CHARGE这个值是否为0来判断。
7.应用:怎么用bader分析研究体系的电荷转移总量。
首先,我们可以从文件的最后一行NUMBER OF ELECTRONS知道这个体系的总电量。
我们计算一个体系在电荷转移前后的bader电荷,然后在ACF文件中读取其电量。如果两个文件中的原子数相等,就可以直接作差,得到电荷转移的量,如果原子数不同,就拿对应原子的电量作差。

本文发布于:2023-06-09 01:29:18,感谢您对本站的认可!

本文链接:https://www.wtabcd.cn/zhishi/a/168624535969828.html

版权声明:本站内容均来自互联网,仅供演示用,请勿用于商业和其他非法用途。如果侵犯了您的权益请与我们联系,我们将在24小时内删除。

本文word下载地址:bader分析(简易版).doc

本文 PDF 下载地址:bader分析(简易版).pdf

标签:计算   电量   文件   体系   收敛   电荷转移
留言与评论(共有 0 条评论)
   
验证码:
推荐文章
排行榜
Copyright ©2019-2022 Comsenz Inc.Powered by © 实用文体写作网旗下知识大全大全栏目是一个全百科类宝库! 优秀范文|法律文书|专利查询|