第一性计算原理供参习.doc

Vasp 我所用第一原理是基于密度泛函（DFT）的从头计算，是以电子密度作为基本变量（HK定理），通过求解kohn-sham方程，迭代自洽得到体系的基态电子密度，然后求体系的基态性质。还有一种是基于hartree-fock自洽计算，通过自洽求解HF方程，获得体系的波函数，求基态性质。KS方程的计算水平达到了HF水平，同时还考虑了电子间的交换关联作用。 关于DFT中密度泛函的Function其实是交换关联泛函，包括LDA,GGA,杂化泛函等等。一般LDA为局域密度近似，在空间某点用均匀电子气密度作为交换关联泛函的唯一变量，多数为参数化的CA-PZ方案；GGA为广义梯度近似，不仅将电子密度作为交换关联泛函的变量，也考虑了密度的梯度为变量，包括PBE,PE.RPBE等方案。 在处理计算体系中原子的电子态时有两种方法，一种是考虑所有电子叫做全电子法，比如WIEN2K中的FLAPW方法（线性缀加平面波）；另一种是只考虑价电子而把芯电子和原子核构成离子实放在一起考虑即赝势法，一般贋势法是选取一个截断半径，截断半径以内波函数变化较平滑，和真实的不同，截断半径以外则和真实情况相同，而且贋势法得到的本征值和全电子法应该相同。贋势的测试标准应是贋势与全电子法计算结果的匹配度，而不是贋势与实验结果的匹配度，因为和实验结果的匹配可能是偶然的。 关于Ecut的收敛测试。一般情况下，总能相对于不同Ecut做计算，当截断能增大时总能变化不明显即可。但是在需要考虑体系应力时，还需要对应力进行收敛测试，而且应力相对于截断能要比总能更为苛刻。也就是某个截断能下总能已经收敛了，但应力未必收敛。(力的计算是在能量的基础上进行的，能量对坐标的一阶导数得到力。计算量的增大和误差的传递导致力收敛慢。K点也是需要经过测试的。 何时需要考虑自旋？例如BaTiO3中，三个元素分别为=+2，+4，-2价，离子全部为各个轨道满壳层的结构，此时就不必考虑自旋了。对于BaMnO3中，由于Mn+4价时d轨道还有电子但未满，因此需要考虑Mn（4s23d5）的自旋，Ba和O就不必考虑。其实设定自旋就是给定一个原子磁矩的初始值，只在刚开始计算时作为初始值使用。 几何优化包括晶格常数和原子位置的优化，一般情况下也有不优化几何结构直接计算电子结构的，但是对于缺陷形成的计算则往往要优化。 软件大致分为基于平面波的软件，如CASTEP,PWSCF.ABINIT等，计算量大概和体系原子数目的三次方相关；还有基于原子轨道线性组合的软件，比如openmx等，计算量和体系原子数目相关，一般可模拟较多原子数目的体系。 VASP是使用贋势和平面波基组，进行从头量子力学分子动力学计算的软件包。VASP中的方法基于有限温度下的局域密度近似（用自由能作为变量）以及对每一MD步骤用有效矩阵对角方案和有效Pulay混合求解瞬时电子基态。这些技术可以避免元氏的Car-Parrinello方法存在的一切问题，而后者是基于电子、离子运动方程同时积分的方法。离子和电子的相互作用超缓Vanderbilt贋势（US-PP）或投影扩充波（PAW）方法描述。两种技术都可以相当程度地减少过度金属或第一行元素的每个原子所必须的平面波数量。VASP可以很容易地计算力与张力，用于把原子衰减到其瞬时基态中。！！！！！！！！！！！！！！！！！！！！！！！！！！！！！！！！！！！！！！！！！！！！！！！！！！！ VASP程序亮点： 使用PAW方法或超软贋势，因此基组尺寸非常小，描述材料一般需要原子不超过100个平面波，大多数情况下甚至每原子50个平面波就能得到可靠结果。 2. 在平面波程序中，某些部分代码的执行是三次标度。在VASP中，三次标度部分的前因子足可忽略，导致关于体系尺寸的高效标度。因此可以在实空间求解势的非局域贡献，并使正交化的次数最少。当体系具有大约2000个电子能带时，三次标度部分与其它部分可比，因此VASP可用于直到4000个价电子的体系。VASP使用传统的自洽场循环计算电子基态。这一方案与数值方法组合会实现有效、稳定、快速的Kohn-Sham方程自洽求解方案。程序使用的迭代矩阵对角化方案（RMM-DISS和分块Davidson）可能是目前最快的方案。VASP包含全功能的对称性代码，可以自动确定任意构型的对称性。5. 对称性代码还用于设定Monkhorst-Pack特殊点，可以有效计算体材料和对称的团簇。Brillouin区的积分使用模糊方法或四面体方法。四面体方法可以用Bl?chl校正去掉线性四面体方法的二次误差，实现更快的k点收敛速度。1、VASP能够进行哪些过程的计算？怎样设置？ 我们平时最常用的研究方法是做单点能计算，结构优化、从头计算的分子动力学和电子结构相关性质的计算。 一般我们的研究可以按照这样的过程来进行 如