中科大MS--Wien2k(全势)简介.ppt

Wien2k简介 ftp 210 45 126 12 temp wien2k VASP赝势的代表 wien2k全势的代表 VASP是基于赝势的方法 即用一个平缓的势代替内部电子和原子的核的势 Wien2k是全电子势的方法 划分电子势的方法 即采用MT球的方法 球内采用球斜函数展开 间隙用平面波 Wien2k耗时 但是更精确 因此预测电子结构 化学键比较准确 当然做其他性质的计算也是很精确的 VASP相对来说 计算省时 适合大体系的计算 在有足够的资源支持下 而且体系不大 选Wien2k 反之 选Vasp 大体系也可采用Siesta AtomistixToolKit ATK Pwscf等可在小木虫获得相关的资料和帮助 球内部 原子间 Wien2k的安装 Wien2k只有Linux版本 安装比MaterialStudio要复杂 需要安装的数学工具 因为wien2k 08是用fortran90编写的 所以需要安装相应的编译器如intelfortran以及数学库intelmkl 配置好编译器和库的路径 etc ld so conf中添加 还需要看具体路径 opt intel fce 10 1 018 lib opt intel cmkl 9 1 lib em64t在 etc profile中添加 PATH opt intel fce 10 1 018 bin PATH PATH opt intel idbe 10 1 018 bin PATH exportPATHMANPATH opt intel idbe 10 1 018 man MANPATH MANPATH opt intel fce 10 1 018 man MANPATH exportMANPATHLD LIBRARY PATH opt intel cmkl 9 1 lib em64t LD LIBRARY PATH LD LIBRARY PATH opt intel fce 10 1 018 lib LD LIBRARY PATH exportLD LIBRARY PATHINCLUDE opt intel cmkl 9 1 include INCLUDE exportINCLUDEIFLPATH opt intel fce 10 1 018 libexportIFPATHMKLPATH opt intel cmkl 9 1 lib em64texportMKLPATHINTEL LICENSE FILE opt intel licenses exportINTEL LICENSE FILE 编译过程 gunzip gz 解压缩所有的文档chmod X expand lapw 改变属性使expand lapw成为可运行文档 expand lapw 解压缩所有的lapw文档 siteconfig lapw 根据自己的系统进行编译配置 Compileroptions FR w mp1 prec div pad ip O3 xPLDFLAGS设置 L opt intel cmkl 9 1 lib em64t Bstatic lguide Bdynamic lpthread并行过程的设置 sharedmemory一项选择y 接下来选择y和mpif90 如果为集群 可以选择n rsh y和mpif90 RP LIBS的设置 L opt intel cmkl 9 1 lib em64t libmkl scalapack aL opt intel cmkl 9 1 lib em64t libmkl lapack aL opt intel cmkl 9 1 lib em64t libmkl em64t aL opt intel cmkl 9 1 lib em64t libguide aL opt intel fce 10 1 018 lib libsvml aL opt intel cmkl 9 1 lib em64t libmkl blacs aL opt intel cmkl 9 1 lib em64t libmkl gfortran a 编译 运行接口界面 w2web 帐号密码 使用 请大家把工作目录设在 home msemsi wien2k work savefileandcleanup nn bondlengthfactor 2 Numberofk points 1000 文件管理 计算过程可查看TiC dayfile 性能分析 类似于materialstudio中的analysis 电子密度图 如果有XCRYSDEN 需要另外安装 选择2D plot设定分辨率为100点 firstline 选择3个原子定义一个面选择矩形区域 并放大计算电子密度 产生轮廓选择颜色晶体结构用球棒模型表示 没有XCRYSDEN 则需要编辑TiC in5 1 104 originofplot x y z denominator 1304 x endofplot3 104 y endofplot323 x y znumberofshells100100 x yplottingmesh chooseratio similartox ylengthRHOANGVALNODEBUGORTHO Calculateelectrondensity xlapw5 Plotoutput usingrhoplot afterthefirstpreviewselectarangezmin 0 5tozmax 2 分析态密度谱 Calculatepartialcharges xlapw2 qtl 分析能带结构 各轨道的共享 课后练习 输入结构并计算物性 NiAl3N4215 P 43m a 0 4436nmN