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
相关内容
热门资讯
亚德诺半导体斥13.5亿美元现...
亚德诺 半导体表示，将以13.5亿美元现金收购Alif Semiconductor。该交易将帮助亚德...
半导体设备ETF国泰（1595...
9月9日，A股市场主要指数盘中震荡分化，上证指数飘红，涨幅0.34%，深证成指飘红，涨幅0.18%，...
埃科光电：公司持续关注人形机器...
证券之星消息，埃科光电(688610)09月09日在投资者关系平台上答复投资者关心的问题。 投资者提...
中瓷电子：约1.11亿股限售股...
每经AI快讯，中瓷电子（SZ 003031，收盘价：120.6元）9月9日晚间发布公告称，公司限售股...
龙胜电器取得超薄防尘开关面板专...
国家知识产权局信息显示，龙胜电器有限公司取得一项名为“一种超薄防尘开关面板”的专利，授权公告号CN2...
广州单片机培训哪家强，有真实项...
在广州，想要学习单片机开发的人越来越多，无论是电子、通信、自动化专业的在校生，还是希望从传统行业跨界...
A股PCB概念股集体上涨，华正...
A股市场PCB概念股集体上涨，其中，则成电子涨11%， 华工科技10CM涨停， 强达电路涨超8%， ...
星钥半导体完成超10亿元融资，...
企查查APP显示，9月4日，星钥半导体（武汉）有限公司完成超10亿元A轮融资，由宁德时代旗下溥泉资本...
险资1.57万亿持仓曝光！PC...
周末，金融圈被一则重磅消息彻底刷屏：财政部大手笔注资，8家机构集体“回血”！ 工商银行1000亿、农...