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
相关内容
热门资讯
普联技术取得电源装置及电源系统...
国家知识产权局信息显示，普联技术有限公司取得一项名为“电源装置及电源系统”的专利，授权公告号CN22...
电荷转移与活化效率
一、为什么电荷转移决定“活化效率？ 在催化反应、光伏、光电化学、电池等领域，活化效率指的是反应物分...
主力榜丨同类流入金额第一！半导...
4月20日，半导体设备ETF国泰（159516）报收0.908元，收涨1.0%，成交金额10.97亿...
又一家功率半导体企业宣布涨价，...
4月17日，国内功率半导体企业捷捷微电发布《投资者关系活动记录表》，正式披露因原材料价格持续高位及产...
SK海力士：开始生产专为英伟达...
SK海力士：开始生产专为英伟达Vera Rubin芯片设计的SOCAMM2芯片。 来源：金融界AI...
华为Pura90 Pro搭载麒...
钛媒体App 4月20日消息，华为Pura 90系列今日发布。其中，华为Pura90 Pro 和Pr...
中信建投：霍尔木兹海峡开关状态...
每经AI快讯，中信建投发布研报称，4月17日，伊朗外长Araqchi宣布，霍尔木兹海峡已对商业船只恢...
集成电路ETF国泰（15954...
4月20日，集成电路ETF国泰（159546）涨超2%，AI硬件投资方兴未艾。 东吴证券指出，科技方...
富门子电气取得加固型电容柜专利...
国家知识产权局信息显示，襄阳市富门子电气有限公司取得一项名为“一种加固型电容柜”的专利，授权公告号C...