VASP 计算Bader电荷应用案例？
创始人
2025-11-07 15:35:24
0次
Bader电荷分析基于Richard Bader提出的一种原子电荷分析方法，通过电子密度的拓扑结构来划分原子区域，从而计算每个原子的电荷分布。深圳华算科技有限公司将详细介绍晶体氧化镁Bader电荷计算流程和数据处理方法。
结构优化
首先通过晶体结构数据库下载MgO的cif文件，并生成POSCAR文件，如下图所示。
然后编写结构优化的INCAR和KPOINTS文件。INCAR文件：
ISTART=0 #随机产生初始波函数
ICHARG=2 #从原子电荷密度产生体系初始电荷密度
PREC=M #计算精度为Medium
ISPIN=1 #关闭自旋极化
ALGO=N #DAV算法优化电子波函数
NELM=60 #电子波函数最多计算60步
EDIFF=1E-5 #电子波函数能量收敛标准1E-5 eV
ENCUT=400 #平面波截断能400 eV
IVDW=11 #考虑范德华力修正
IBRION=2 #共轭梯度法优化晶体结构和原子坐标
NSW=100 #晶体结构和原子坐标优化步数最大100步
ISIF=3 #优化原子坐标和晶体结构
EDIFFG=-0.1 #原子残余力小于0.1 eV/A
ISMEAR=0 #费米能级附近电子占据数为高斯分布
SIGMA=0.1 #高斯分布展宽0.1 eV
KPOINTS文件：
Automatic generation #注释行
0 #自动产生K点网格
G #布里渊区K点网格以
Gamma点为中心
5 5 5 #K点网格密度
0 0 0 #K点网格中心平移矢量
自洽计算
完成结构优化计算后，保持KPOINTS，POTCAR文件不变，将CONTCAR文件复制成POSCAR文件，并对结构优化的INCAR文件作如下修改：IBRION=-1 #固定结构和原子位置保持不变
NSW=0 #仅作一步自洽计算
LAECHG=.T. #计算Bader电荷
数据处理
完成自洽计算后，运行以下两个脚本，从而可以得到每个原子Bader电荷，它位于ACF.dat文件中CHARGE一列，其中1-4为Mg原子，5-8为O原子，排序和POSCAR一致。chgsum.pl AECCAR0 AECCAR2bader CHGCAR -ref CHGCAR_sum
相关内容
热门资讯
半导体板块11月7日跌1.13...
证券之星消息，11月7日半导体板块较上一交易日下跌1.13%，大为股份领跌。当日上证指数报收于399...
主力榜丨同类流入金额第一！科创...
11月7日，科创芯片ETF富国（588810）报收1.677元，收跌1.47%，成交金额5660.2...
原创 1...
操心“太空霸权”主导位数十年的老美，终于栽跟头了。 前段时间，伊朗领导人乘坐的直升机因天气等原因坠毁...
2025年北京退休人员养老金计...
新京报讯 11月7日，北京市人力资源和社会保障局、北京市财政局、国家税务总局北京市税务局发布关于20...
福建舰官宣交付，美媒向华盛顿发...
据新华社7日消息，11月5日中国海军迎来了一个历史性时刻——我国自主设计建造的第一艘电磁弹射型航空母...
中航光电：与富士达无同业竞争实...
证券之星消息，中航光电(002179)11月07日在投资者关系平台上答复投资者关心的问题。 投资者提...
激光二极管在干涉仪上面的作用
我们都知道激光二极管在日常生活中应用的十分广泛，那么你知道市面上常见的激光干涉仪的工作原理吗？其里面...
著名表演艺术家刘大刚因病去世！...
搜狐娱乐讯 7日消息，《西游记》“沙僧”扮演者刘大刚因病去世！戏曲演员、豫剧大师马金凤弟子刘冰在社交...
11月7日国企改革（39997...
证券之星消息，11月7日，国企改革（399974）指数报收于1916.43点，跌0.1%，成交135...
namur电容式接近开关
这是(xiexinke79)整理的信息，希望能帮助到大家 在工业自动化控制领域，检测物体是否存在或位...