Bader分析
一、背景
理查德·贝德(Richard Bader)开发了一种将分子分解为原子的直观方法。他对原子的定义纯粹基于电子电荷密度。Bader使用所谓的零磁通表面来划分原子。零通量表面是2D表面,其上电荷密度垂直于表面。通常在分子系统中,电荷密度在原子之间达到最小值,这是将原子彼此分开的自然位置。除了作为分子中原子可视化的直观方案外,Bader的定义通常也可用于电荷分析。例如,Bader体积内的电荷与原子的总电子电荷很接近。电荷分布可用于确定相互作用的原子或分子的多极矩。Bader的分析也被用来定义原子的硬度,可以用来量化从原子中去除电荷的成本。
二、程序下载
1)下载网址:http://theoryNaN.utexas.edu/vtsttools/scripts.html (用到chgsum.pl脚本)
http://theoryNaN.utexas.edu/henkelman/code/bader/download/bader_lnx_64.tar.gz (下载bader处理脚本)
2)将前面下载的内容解压,然后复制到~/bin 文件夹中。
3)如果后面运行的时候遇到权限问题:执行:
chmod u+x ~/bin/chgsum.pl
chmod u+x ~/bin/bader
三、运行
1 优化自己的体系结构:
使用VASP计算,获得稳定的结构。
2 bader电荷分析的计算部分:
1) 将前面优化完的CONTCAR 重命名为POSCAR
2) 设置INCAR:
LAECHG=.TRUE.
LCHARG = .TRUE.
NSW = 0
IBRION = -1 (前面有了NSW = 0, 这个也可以不设置)
3) KPOINTS可以用之前优化计算的,也可设置的稍微大些。
chgsum.pl AECCAR0 AECCAR2
注意: AECCAR0和AECCAR2 是由LAECHG=.TRUE.这个参数控制输出的。结束后,会得到 CHGCAR_sum 这个文件。
四、分析
可以使用该命令运行该程序
bader CHGCAR -ref CHGCAR_sum
生成以下输出文件:ACF.dat,BCF.dat,AtomVolumes.dat。
ACF.dat包含每个原子的坐标,根据Bader分区与其相关的电荷,根据Bader分区的整体百分比以及到表面的最小距离。如果使用了伪电位,则应将该距离与核心区域的最大截止半径进行比较。
BCF.dat包含每个Bader最大值的坐标,该体积内的电荷,最近的原子以及到该原子的距离。
AtomVolumes.dat包含已分配给每个原子的每个卷的编号。
Bader分析的更多相关文章
- alias导致virtualenv异常的分析和解法
title: alias导致virtualenv异常的分析和解法 toc: true comments: true date: 2016-06-27 23:40:56 tags: [OS X, ZSH ...
- 火焰图分析openresty性能瓶颈
注:本文操作基于CentOS 系统 准备工作 用wget从https://sourceware.org/systemtap/ftp/releases/下载最新版的systemtap.tar.gz压缩包 ...
- 一起来玩echarts系列(一)------箱线图的分析与绘制
一.箱线图 Box-plot 箱线图一般被用作显示数据分散情况.具体是计算一组数据的中位数.25%分位数.75%分位数.上边界.下边界,来将数据从大到小排列,直观展示数据整体的分布情况. 大部分正常数 ...
- 应用工具 .NET Portability Analyzer 分析迁移dotnet core
大多数开发人员更喜欢一次性编写好业务逻辑代码,以后再重用这些代码.与构建不同的应用以面向多个平台相比,这种方法更加容易.如果您创建与 .NET Core 兼容的.NET 标准库,那么现在比以往任何时候 ...
- UWP中新加的数据绑定方式x:Bind分析总结
UWP中新加的数据绑定方式x:Bind分析总结 0x00 UWP中的x:Bind 由之前有过WPF开发经验,所以在学习UWP的时候直接省略了XAML.数据绑定等几个看着十分眼熟的主题.学习过程中倒是也 ...
- 查看w3wp进程占用的内存及.NET内存泄露,死锁分析
一 基础知识 在分析之前,先上一张图: 从上面可以看到,这个w3wp进程占用了376M内存,启动了54个线程. 在使用windbg查看之前,看到的进程含有 *32 字样,意思是在64位机器上已32位方 ...
- ZIP压缩算法详细分析及解压实例解释
最近自己实现了一个ZIP压缩数据的解压程序,觉得有必要把ZIP压缩格式进行一下详细总结,数据压缩是一门通信原理和计算机科学都会涉及到的学科,在通信原理中,一般称为信源编码,在计算机科学里,一般称为数据 ...
- ABP源码分析一:整体项目结构及目录
ABP是一套非常优秀的web应用程序架构,适合用来搭建集中式架构的web应用程序. 整个Abp的Infrastructure是以Abp这个package为核心模块(core)+15个模块(module ...
- HashMap与TreeMap源码分析
1. 引言 在红黑树--算法导论(15)中学习了红黑树的原理.本来打算自己来试着实现一下,然而在看了JDK(1.8.0)TreeMap的源码后恍然发现原来它就是利用红黑树实现的(很惭愧学了Ja ...
随机推荐
- css 超出部分显示省略号 汇总
单行: 加宽度 overflow: hidden; text-overflow:ellipsis; white-space: nowrap; 多行: display: -webkit-box; -we ...
- view视图总结
视图实质上存储的是一段sql. 创建方式: 1.create view 视图名 as 子查询语句 , 特点:与主表数据同步,对当前视图进行修改,会同时将根表一并修改. 2.create or ...
- 60. Permutation Sequence (JAVA)
The set [1,2,3,...,n] contains a total of n! unique permutations. By listing and labeling all of the ...
- Linux上安装JDK1.8,tomcat9,以及mysql8的步骤
(该篇是在centos7上安装JDK1.8.0_201 tomcat9.0.16 和 mysql8.0.15) 一.安装JDK 方式一 1.首先,下载JDK(链接http://www.oracle. ...
- java 如何读取 properties 配置文件
- Jmeter分布式测试dubbo接口1
最近工作中接到一个需求,需要对一个Dubbo接口进行压力测试,测试其性能,之前一直使用jmeter做压力测试,在踏了好多坑之后,决定把这些记录下来,顺便也希望能帮助到大家. 开始测试之前,我们需要先知 ...
- Vue基础第一章
Vue的简单示例 <!DOCTYPE html> <html> <head> <meta charset="utf-8"> < ...
- Python面向对象的三大特性之异常处理
一.错误与异常 程序中难免会出现错误,而错误分为两种 1.语法错误:(这种错误,根本过不了python解释器的语法检测,必须在程序执行前就改正) 2.逻辑错误:(逻辑错误),比如用户输入的不合适等一系 ...
- centos6.5下修改系统的roo用户/非root用户的密码
1.修改系统root用户的密码 [........~]# passwd然后输入新密码,若提示密码太简单,无需理会,直接敲回车: 然后再次输入新密码,即可修改成功. 2.修改系统非root用户的密码:e ...
- 初始Turtle
---恢复内容开始--- turtle库介绍 海龟作图最初源自20世纪60年代的Logo(创建于1967年,是一种教育编程语言)编程语言. 常用函数说明 画笔的属性 设置画笔的宽度 t.pensi ...