搜索
第一原理计算方法简介及MaterialsStudio中Cestep使用第一原理计算方法简介第一性原理方法（First-principlesmethod），有时候也称为从头计算（abinitio），其基本思想是将多原子构成的体系当作电子和原子核（或原子实）组成的多粒子系统，从量子力学第一性原理（多电子体系的Schrödinger方程）出发，对材料进行“非经验性”的模拟。原则上，第一性原理方法无可调经验参数，只用到几个基本物理常数，如光速c、Planck常数h、电子电量e、电子质量me以及原子的各种同位素的质量，因此处理不同体系时候具有较好的可移植性（transferability）。但是，在具体实行时，仍依赖于具体近似方法的选取，从而带来系统误差。多粒子体系（电子＋核）的薛定谔方程三个近似a.非相对论近似(忽略了电子运动的相对论效应)vec,ve~108cm/s3×1010cm/s),me=m0求解非相对论的薛定谔方程，而不是相对论的狄拉克方程b.Born-Oppenheimer近似，核固定近似中子/质子的质量是电子质量的约1835倍，即电子的运动速率比核的运动速率要高3个数量级，因此可以实现电子运动方程和核运动方程的近似脱耦。这样，电子可以看作是在一组准静态原子核的平均势场下运动。c.单电子近似把体系中的电子运动看成是每个电子在其余电子的平均势场作用中运动，从而把多电子的薛定谔方程简化单电子方程。HartreeFock方程薛定谔方程简化为：将总Hamilton分解成单电子贡献H0和电子－电子相互作用U。应用变分法计算多电子波函数方程，可得Hartree-Fock方程。量子化学分子轨道方法分子轨道方法：在Hartree-Fock框架下，将单电子波函数用原子轨道（Slater型－STO，Gaussian型－GTO）的线性叠加表示来求解。组态相互作用方法（采用多个Slater行列式考虑电子关联）Mφller-Plesset（MP）修正（将关联作用作为微扰修正）Hartree-Fock方法（忽略交换作用，严格计算电子积分）半经验方法，如CNDO，MNDO，MINDO，AM1，PM3等（同样忽略交换作用，近似计算电子积分）精度，计算量密度泛函理论Hartree-Fock方法的主要缺限：(1)完全忽略电子关联效应；(2)计算量偏大，随系统尺度4次方关系增长。20世纪60年代，Hohenberg，Kohn和Sham(沈吕九)提出了密度泛函理论(DFT)。DFT理论奠定了将多电子问题转化为单电子方程的理论基础，给出了单电子有效势计算的可行方法，DFT在计算物理、计算化学、计算材料学等领域取得巨大成功。1998年，W.Kohn与分子轨道方法的奠基人Pople分享了诺贝尔化学奖。密度泛函理论的主要目标就是用电子密度取代波函数做为研究的基本量。用电子密度更方便处理。密度泛函理论Hohenberg-Kohn第一定理指出体系的基态能量仅仅是电子密度的泛函。Hohenberg-Kohn第二定理证明了以基态密度为变量，将体系能量最小化之后就得到了基态能量。根据以上两定理，将薛定谔方程转变为Kohn-Sham方程电子与原子核间的库仑势电子间的库仑势交换关联势（未知）密度函数密度泛函理论LDA和GGA近似Kohn-Sham方程原则是精确的，但遗憾的是交换关联势是未知的。要进行具体计算，就必须使用近似方法求出交换关联势。常用的近似方法有局域密度近似(LocalDensityApproximation)和广义梯度近似(GeneralizedGradientApproximation),在某些情况下，广义梯度近似改善了局域密度近似的计算结果，但它并不总是优于局域密度近似。密度泛函理论基组(basisset)求解Kohn-Sham方程，选取适当的基组，将波函数对其展开，将方程求解转化为线性代数问题。一般选用如下基组展开:(Linearized)augmentedplanewaves-(L)APW’s(Linearized)muffin-tinorbitals-(L)MTO’sProjectoraugmentedwaves-PAW’s密度泛函理论赝势(pseudopotential)赝势就是把离子实的内部势能用假想的势能取代真实的势能，但在求解波动方程时，不改变能量本征值和离子实之间区域的波函数。模守恒赝势NCP(NormConservingPseudopotential)和超软赝势USPP(UltrasoftPseudoptential)第一原理常用计算软件根据对势函数及内层电子的处理方法不同主要分为两大类，一种是波函数中包含了高能态和内层电子，而势函数只是原子核的贡献，这称为全电子（allelectroncalculation）法，另一种处理方法是势函数为原子核和内层电子联合产生的势，称为离子赝势，波函数只是高能态电子的函数，这称为赝势（pseudo-potential）法。第一原理计算软件CodeNameBasisSetPotentialsPlaneWavePseudopotentialCodes操作系统WebSiteABINITPlanewavePseudo,PAWLinux概述MaterialStudio的特点：•采用服务器/客户机模式的软件环境，Microsoft标准用户界面，不需要登录服务器。HttpGatewayFtpXP,2000,2003,Vista,2008MaterialStudio简介ModuleparallelWindowsLinuxIA32LinuxIA64ModuleparallelWindowsLinuxIA32LinuxIA64MaterialsVisualizer√ONETEP√√√√AdsorptionLocator√√Polymorph√√AmorphousCell√√√QMERA√√Blends√√QSARandQSARPlus√√CASTEPandNMRCASTEP√√√√Reflex-PatternProcessingandPowderDiffraction√COMPASS√√√√Reflex-PowderIndexing√√√CCDC√√Reflex-PowderRefinement√√Conformers√√ReflexPlus√√√Discover√√√√ReflexQPA√√DMol3√√√√Sorption√√DPD√√√Synthia√Equilibria√√√VAMP√√√Forcite√√X-Cell√√√Gaussian√√√√Mesotek√√√GULP√√√√Morphology√√MesoDyn√√√√2.能够容易地创建并研究分子模型或材料结构，使用极好的制图能力来显示结果。3.与其它标准PC软件整合的工具使得容易共享这些数据Origin,Matlab…。4.采用材料模拟中领先的十分有效并广泛应用的模拟方法（LDA,GGA）。5.可模拟的内容：催化剂、聚合物、固体化学、结晶学、晶粉衍射以及材料特性等。主要模块：Visualizer建模模块AmorphousCellBlendsCASTEPConformersDMol3DPDDiscoverEquilibriaForciteGULPMesoDynMorphologyOnetepPolymorphQMERAReflexSynthiaVAMPGaussian计算和分析模块分子力学与分子动力学MS.DISCOVERMS.COMPASSMS.AmorphousCellMS.ForciteMS.ForcitePlusMS.GULPMS.EquilibriaMS.Sorption晶体、结晶与X射线衍射MS.PolymorphPredictorMS.MorphologyMS.X-CellMS.ReflexMS.ReflexPlusMS.ReflexQPA量子力学MS.Dmol3MS.CASTEPMS.NMRCASTEPMS.VAMP高分子与介观模拟MS.SynthiaMS.BlendsMS.DPDMS.MesoDynMS.MesoPro定量结构-性质关系MS.QSARMS.QSARPlusMS.Dmol3DescriptorVisualizer:图形化建模模块可构建计算的模型：晶胞，分子，晶体表面，纳米结构，聚合物等锐钛矿TiO2TiO2(111)Pt(110)-CO(2x1)碳纳米管TiO2纳米棒MaterialStudio晶体结构模型建立建立方铅矿PbS晶体结构模型（实验15）步骤（1）查晶体结构数据（如ICSD、PDF数据库）（2）通过软件建模（如MaterialStudio中模块Visualizer、Diamond）Castep使用CASTEP是特别为固体材料学而设计的一个现代的量子力学基本程序，其使用了密度泛函(DFT）平面波赝势方法，进行第一原理量子力学计算，以探索如半导体，陶瓷，金属，矿物和沸石等材料的晶体和表面性质。典型的应用包括表面化学，键结构，态密度和光学性质等研究，CASTEP也可用于研究体系的电荷密度和波函数的3D形式。此外，CASTEP可用于有效研究点缺陷（空位，间隙和置换杂质）和扩展缺陷（如晶界和位错）的性质。MaterialStudio使用组件对话框中的CASTEP选项来准备，启动，分析和监测CASTEP计算工作。计算：允许选择计算选项（如基集，交换关联势和收敛判据），作业控制和文档控制。分析：允许处理和演示CASTEP计算结果。这一工具提供加速整体直观化以及键结构图，态密度图形和光学性质图形。CASTEP模块CambridgeSerialTotalEnergyPackage)CASTEP的任务CASTEP计算包括单点的能量计算，几何优化或分子动力学。可提供这些计算中的每一个以便产生特定的物理性能。在CASTAP计算中有很多运行步骤，可分为如下几组：结构定义：必须规定包含所感兴趣结构的周期性的3D模型文件，有大量方法规定一种结构：可使用构建晶体（BuildCrystal)或构建真空板(BuildVacuumStab)来构建，也可从已经存在的结构文档中引入，还可修正已存在的结构。注意：CASTEP仅能在3D周期模型文件基础上进行计算，必须构建超单胞，以便研究分子体系。提示：CASTAP计算所需时间随原子数平方的增加而增加。因此，建议是用最小的初晶胞来描述体系，可使用Build\Symmetry\PrimitiveCell菜单选项来转换成初晶胞。CASTEP的任务计算设置：合适的3D模型文件一旦确定，必须选择计算类型和相关参数，例如，对于动力学计算必须确定系综和参数，包括温度，时间步长和步数。选择运行计算的磁盘并开始CASTEP作业。结果分析：计算完成后，相关的CASTEP作业的文档返回用户，在项目面板适当位置显示。这些文档进一步处理能获得所需的观察量如光学性质。CASTAP中选择一项任务1从模块面板（ModuleExplorer)选择CASTAP\Calculation2选择设置表3从任务列表中选择所要求的任务CASTEP能量任务CASTEP能量任务允许计算特定体系的总能量以及物理性质。除了总能量之外，在计算之后还可报告作用于原子上的力；也能创建电荷密度文件；利用材料可视化（MaterialVisualizer)允许目测电荷密度的立体分布；还能报告计算中使用的Monkhorst-Park的k点的电子能量，因此在CASTEP分析中可生成态密度图。对于能够得到可靠结构信息的体系