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
相关内容
热门资讯
原创 禁...
文/王新喜 中国对日本的所有两用物项的出口禁令已经即时生效，禁令执行的第四天，有日本客户联系了中国的...
华正新材取得复合层和光伏组件专...
国家知识产权局信息显示，浙江华正新材料股份有限公司取得一项名为“复合层和光伏组件”的专利，授权公告号...
国星光电：LED封装及组件产品...
证券之星消息，国星光电(002449)01月13日在投资者关系平台上答复投资者关心的问题。 投资者提...
我国自研原位传感器完成验证
近日，由中国科学院西安光学精密机械研究所吴国俊团队牵头研制的多型国产海洋生物地球化学原位传感器完成多...
2025 EMI滤波器生产厂家...
根据《2025年电磁兼容产业蓝皮书》数据，全球EMI滤波器市场规模预计以8.5%的年复合增长率增至2...
婚礼电子请柬的访客记录，让谁尴...
在婚礼前制作一份电子请柬，私发给亲朋好友，或分享到微信群、朋友圈，接收各方祝福，如今已成为婚礼流程中...
重庆A股23家上涨 电科芯片、...
1月13日，78家重庆A股上市公司中有23家上涨。 同花顺iFinD数据显示，电科芯片以上涨10.0...
锂威电子取得隔膜涂布装置专利，...
国家知识产权局信息显示，惠州锂威电子科技有限公司取得一项名为“一种隔膜涂布装置”的专利，授权公告号C...
徐工矿业申请新能源矿用挖掘机与...
国家知识产权局信息显示，徐州徐工矿业机械有限公司申请一项名为“一种新能源矿用挖掘机与电源车的安全通断...
比芯片荒的影响差远了，内存涨价...
2026一开年，“内存涨价比金条还快”多次冲上热搜。同样是给成本带来的压力，汽车的体感远没有消费电子...