产品
ABINIT的主程序使用赝势和平面波，用密度泛函理论计算总能量，电荷密度，分子和周期性固体的电子结构，进行几何优化和分子动力学模拟，用TDDFT（对分子）或GW近似（多体微扰理论）计算激发态。此外还提供了大量的工具程序。程序的基组库包括了元素周期表1-109号所有元素。ABINIT适于固体物理，材料科学，化学和材料工程的研究，包括固体，分子，材料的表面，以及界面，如导体、半导体、绝缘体和金属。
1 材料设计软件包的整体功能描述
材料科学是现代科技发展的支柱之一，随着高科技的发展，对材料性能的要求也越来越高。传统的材料设计常常是通过‘炒菜’方法，建立在大量反复的实验基础之上。这种方法由于缺少必要的理论指导，不可避免地具有一定的盲目性，从而造成人力物力的巨大浪费，且材料设计的周期很长。即使用炒菜方法获得了一些经验规律，往往也是知其然，而不知其所以然。不能满足快速发展的高科技对新材料的需求。材料科学家们希望借助于计算机，获取材料性能方面的信息，预测材料的性能，从而指导材料设计，这也是材料设计工作发展的必然方向。
本课题所建立的软件平台包括三个部分：第一原理平面波赝势总能计算软件包ABINIT；基于第一原理计算的材料设计辅助软件；以及平台界面。下面对这三部分的功能简单说明。
1.1 第一原理平面波赝势总能计算软件包ABINIT
第一原理计算辅助材料设计是目前国际上的研究热点之一。一方面，第一原理方法可以精确地计算出材料中与宏观性能密切相关的量，如弹性模量，缺陷形成能等，通过这些物理量的计算，间接预测材料的力学性能；另一方面，第一原理计算可以为更高层次的材料模拟如有限元方法，Kinetic Monte Carlo方法等提供基本参数，从而直接得到可与实验比较的结果。更为重要的是，第一原理方法不需要任何经验或拟合的参数，是一种真正意义上的预测。
第一原理方法是在用户给定体系的几何结构和化学成分后，对体系做一些适当的简化近似（如Bonopenheimer绝热近似、局域密度近似、平面波近似等），自洽求解体系的Schrodinger方程，得到其电子分布密度，并给出系统能量等与电子结构密切相关的物理量。
该软件包（包括源程序和可执行程序）在深腾6800上的安装路径为/home4/ qmhu/ abinit。程序的编译过程及使用指南，请见软件包内的说明（/home4/ qmhu/ abinit/Info/）。
1.2 材料设计附件
第一原理计算涉及到固体物理，量子力学等复杂艰深的物理理论，软件编写者大多为物理等基础学科领域的科学家，这使得现有第一原理软件包一般注重在物理，化学等基础科学中的应用，给出的是一些基本的物理结果如系统总能，电荷分布，原子力等等。从这些基本物理量很难联系材料的宏观性能特别是力学性能。为了让第一原理计算能更有效地为实际材料设计服务，根据我们长期材料计算设计的经验，自行编程并收集整理了一些材料设计辅助软件，这些软件能以第一原理方法的计算结果为输入，得到与材料宏观性能密切相关的量。
各材料设计附件的功能如下：
1.2.1 PDConcen
PDConcen利用第一原理计算得到的“粗”点缺陷，自洽计算在不同的化学计量比或温度下金属间化合物中的点缺陷浓度；根据点缺陷浓度计算化合物的长程序参数；由长程序参数随温度的变化确定化合物的有序无序相变温度。
该程序（包括源程序和可执行程序）在深腾6800上的安装路径为/home4 /qmhu/accessory/pdconcen。
1.2.2 ANN
ANN可以用来寻找复杂数据间的关系，并利用这些关系进行预测与优化。
该软件源程序在深腾6800上的保存路径为/home4/qmhu/accessory/ann。因需针对6800进行优化，目前尚未编译安装。
1.2.3 ISOTROPY
ISOTROPY是美国Birgham Young大学Stokes 和 Hatch所编写的利用群论方法分析晶体中的固态相变的软件。该程序共包括8个功能模块，其中ISOTROPY为主模块，提供详细的群论知识；ISODISPLACE确定与空间群不可约表征相关的晶体结构的位移型畸变；FINDSYM根据晶胞中的原子位置确定晶体的空间群对称性；COPL根据母相和子相的空间群对称性，确定完整的相变序参量列表；INVARIANTS产生序参量各组分的不变多相式；SMODES寻找晶体中将动态矩阵（最小）块对角化的位移模式；RELATION确定两个晶体是否具有群－子群关系；FROZSL用冻声子法计算声子频率与位移模式；COMSUBS寻找重构相变中两个结构的共有子群。
该程序在深腾6800上的安装路径为/home4/qmhu/accessory/isobyu。作者不提供该软件的源程序。使用指南请见该目录下的isoman.pdf。
1.2.4 PHON
PHON是Alfe开发的一个计算声子谱及热力学的程序。它以用户所给出的晶体结构模型为输入，产生用第一原理计算该模型的晶格声子行为所应该施加的位移；通过分析这些位移所引起的作用在原子上的力(可以用第一原理计算得到)，该程序可以给出晶体的振动频率，声子谱等。更进一步，它还可以计算出晶体的一些热力学性能，如Helmholtz自由能，振动熵，比热，及内能等。
该程序在深腾上的安装路径为/home4/qmhu/accessory/phon。程序的编译及使用指南，请参见该软件内的说明（manual.ps）。
1.2.5 ATAT
ATAT是Walle及其合作者开发的一个利用第一原理结合MonteCarlo模拟计算相图及热力学性能的软件包。该软件包首先采用团簇展开(cluster expansion)，用第一原理方法计算一系列构型的形成能，然后用这些能量作为Monte Carlo模拟的哈密顿量，计算材料的热力学性能。
该程序在深腾上的安装路径为/home4/qmhu/accessory/atat。程序的编译及使用指南，请参见该软件包内的详细说明(/home4/qmhu/accessory/atat/doc)。
1.3 用户界面
因为第一原理计算涉及到固体物理，量子力学等复杂艰深的物理理论，使用时参数的输入及结果处理往往要求用户具有相当丰富的背景知识，门槛很高，这也限制了材料设计者的使用。基于此，结合远程计算的要求，我们开发了一个用户界面，并提供了丰富且简单易用的在线帮助。利用该界面，可以通过网页为超级机上的第一原理计算程序提供输入并计算，观察输出结果，还可以提供适时方便的在线帮助。
2 内核程序ABINIT与材料设计附件间的协作关系
2.1 ABINIT与PDConcen
PDConcen需要ABINIT提供的“粗”点缺陷形成能含缺陷化合物与不含缺陷化合物见的总能量差），及参考能量。
2.2 ABINIT与PHON
PHON根据晶体对称性，为ABINIT力场计算确定进行声子谱计算所需的原子位移；ABINIT对PHON确定的一系列结构结构计算，提供PHON声子谱计算所需的力场。
2.3 ABINIT与ATAT
ATAT首先选取一种团簇展开的形式，用ABINIT计算得到相应构型的总能量，ATAT根据这些总能量，进行拟合，得到有效相互作用能，作为Kinetic Monte Carlo模拟，计算相图及热力学参数。
3 输入输出方式说明
我们为ABINIT软件编写了用户界面，用户可以通过调用相关网页，远程输入计算参数，提交作业，并观察输出结果。
材料设计附件的输入输出尚未集成到该用户界面中，采用命令行及文本方式输入输出。
各输入输出参数的含义，在相关软件包中有详细说明。
4 用户指南及算例
用户指南请参见各软件包内的详细说明。
Abinit算例见/userdata3/qmhu/test/supercell。我们对含一个空位的Al超晶胞进行了并行效率测试，使用的CPU数为1到16个。
ABINIT的主程序使用赝势和平面波，用密度泛函理论计算总能量，电荷密度，分子和周期性固体的电子结构，进行几何优化和分子动力学模拟，用TDDFT（对分子）或GW近似（多体微扰理论）计算激发态。此外还提供了大量的工具程序。程序的基组库包括了元素周期表1-109号所有元素。ABINIT适于固体物理，材料科学，化学和材料工程的研究，包括固体，分子，材料的表面，以及界面，如导体、半导体、绝缘体和金属。
可以计算很多物理属性：
A. 计算倒格子中核与电子的总能量。
A.1. 计算使用平面波和赝势。
A.2. 总能量的计算使用密度泛函理论（DFT）。可以使用大多数重要的局域密度近似 (LDA)，包括Perdew-Zunger近似。可以使用两种不同的局域自旋密度(LSD)，包括Perdew Wang 92和M. Teter的LSD。还可以使用Perdew-Burke-Ernzerhof，revPBE，RPBE和HCTH等GGA (自旋极化和非极化)。
A.3. 自恰场计算生成DFT基态，以及相关的能量和密度。此后的非自恰计算可以对能带结构的大量k-点产生本征能量。态密度的计算即可以用四面体方法，也可以用模糊技术。
A.4. 程序可以使用多种不同的赝势。对整个周期表适用的有两种：Troullier-Martins型和Goedecker型（这种类型包括自旋-轨道耦合）。如果需要的话，有四个代码可以产生新的赝势。
A.5. 程序本身可以处理金属和绝缘体系。
A.6. 晶胞可以是正交或者非正交。计算可以输入任何对称性及相应的k-点集。
A.7. 电子体系可以用自旋极化和自旋非极化计算。一个特殊的选项可以有效地处理反铁磁性。可以对总能量计算非共线的磁性（不能用于力，张量，相应函数…）。可以禁止晶胞的总磁矩。
A.8. 总能量，力，张量和电子结构的计算可以考虑自旋-轨道耦合。
A.9. 能量可分解为不同的成分（局域势，XC，Hartree…）。
A.10. 计算内部电子本征值。
A.11. 230个空间群和1191个Shubnikov磁群的对称性分析。
B. 计算总能量和本征能量
B.1. 用解析公式计算Hellman-Feynman力。
B.2. 计算应力。
B.3. 极化的计算。
B.4. 响应的计算。
B.5. 计算近似的和准确的磁化系数矩阵和介电矩阵。
B.6. 解析计算电子本征能量的导数。
B.7. 计算光学传导性。
B.8. Born有效电荷的能带分解，以及局域化张量的计算。
C. 激发态
C.1. 用GW近似计算电离能和亲和能。
C.2. 用TDDFT计算原子和分子的（单重、三重）激发态和振荡强度。
D. 移动原子，改变晶胞参数
D.1. 用不同的方法寻找平衡构型。可以同时优化晶胞参数。优化过程中如果需要的话，可以固定指定的晶胞参数，角度，或原子位置。
D.2. 有两种算法进行分子动力学计算。
D.3. 自动分析键长键角。原子坐标的格式支持用可视化软件XMOL显示。
E. 分析和图形工具
E.1. 后期处理程序cut3d用于分析密度和势文件。它还可以改变文件格式，提取2D明面或者1D线。此外还可以分析波函文件。
E.2. 另一个后期处理程序aim，用于进行Bader的“原子中的分子”（AIM）密度分析。
E.3. 对可视化程序产生格式化数据：键结构（用XMGR显示），不同参数的总能量（用XMGR显示），电荷密度（3D轮廓线，先用cut3d，再用商业程序matlab；cut3d也可以产生2D密度图）。
E.4. 后期处理程序band2eps自动画出eps格式的声子散射曲线。
免责声明：本文系网络转载或改编，未找到原创作者，版权归原作者所有。如涉及版权，请联系删
155-2731-8020
