VASP计算实例
VASP表面
Next: 采用VASP如何计算晶体的弹性常数 Up: VASP计算实例 Previous: 用VASP 计算Pd金属的晶格常数 Contents用VASP 计算表面能做表面计算时，第一步我们需要测试K 点的收敛性。
通常，在垂直表面方向用1个K 点就可以了，在平行表面方向，可以用和体材料类似的K 点密度。
其次，我们要测试真空厚度(vacuum thickness)的收敛性。
我们构造完一个slab 后，将真空厚度逐渐从增加到，体系的总能量改变不超过10meV 的时候，可以初步认为真空厚度达到标准。
以下是一个3层的(fcc) Pd slab 的能量随着真空厚度的变化。
其INCAR 文件如下：SYSTEM = undeformed fcc Pd (111) surface calculation Startparameter for this run: PREC = AccurateISTART = 0 job : 0-new 1-cont 2-samecut ICHARG = 2 charge: 1-file 2-atom 10-const ISPIN = 1 spin polarized calculation?Electronic Relaxation 1NELM = 90; NELMIN= 8; # of ELM steps EDIFF = 0.1E-03 stopping-criterion for ELM LREAL = .FALSE. real-space projection NBANDS = 40Ionic relaxationEDIFFG = 0.1E-2 stopping-criterion for IOM NSW = 0 number of steps for IOMIBRION = 2 ionic relax: 0-MD 1-quasi-New 2-CG ISIF = 2 stress and relaxationPOTIM = 0.10 time-step for ionic-motion TEIN = 0.0 initial temperatureTEBEG = 0.0; TEEND = 0.0 temperature during runDOS related values:ISMEAR = 1 ; SIGMA = 0.20 broadening in eV -4-tet -1-fermi 0-gausElectronic relaxation 2 (details)Write flagsLWAVE = F write WAVECARLCHARG = F write CHGCARLVTOT = .TRUE.其中因为Pd是金属，ISMEAR设置为method of Methfessel-Paxton。
vasp在计算磁性的实例
兄弟，问3个问题1，vasp在计算磁性的时候，oszicar中得到的磁矩和outcar中得到各原子磁矩之和不一致，在投稿的是否曾碰到有审稿人质疑，对于这个不一致你们一般是怎么解释的了？2，另外，磁性计算应该比较负责。
你应该还使用别的程序计算过磁性，与vasp结果比较是否一致，对磁性计算采用的程序有什么推荐。
ps：由于曾使用vasp和dmol算过非周期体系磁性，结构对磁性影响非常大，因此使用这两个程序计算的磁性要一致很麻烦。
还不敢确定到底是哪个程序可能不可靠。
3，如果采用vasp计算磁性，对采用的方法和设置有什么推荐。
1，OSZICAR中得到的磁矩是OUTCAR中最后一步得到的总磁矩是相等的。
总磁矩和各原子的磁矩(RMT球内的磁矩)之和之差就是间隙区的磁矩。
因为有间隙区存在,不一致是正常的。
2，如果算磁性，全电子的结果更精确，我的一些计算结果显示磁性原子对在最近邻的位置时，PAW与FPLAW给出的能量差不一致，在长程时符合的很好。
虽然并没有改变定性结论。
感觉PAW似乎不能很好地描述较强耦合。
我试图在找出原因，主要使用exciting和vasp做比较。
计算磁性推荐使用FP-LAPW, FP-LMTO, FPLO很吸引人(不过是商业的)，后者是O(N)算法。
3，使用vasp计算磁性，注意不同的初始磁矩是否收敛为同一个磁矩。
倒没有特别要注意的地方，个人认为。
归根结底，需要一个优秀的交换关联形式出现VASP计算是否也是像计算DOS和能带一样要进行三步（结构优化，静态自洽计算，非自洽计算），然后看最后一步的出的磁矩呢？一直想计算固体中某个原子的磁矩，根据OUTCAR的结果似乎不能分析，因为它里面总磁矩跟OSZICAR的值有一定的差别，据说是OUTCAR中只考虑WS半径内磁矩造成的。
最近看到一个帖子说是可以用bader电荷分析方法分析原子磁矩。
如法炮制之后发现给出的总磁矩与OSZICAR的结果符合的甚好，可是觉得没有根据，有谁知道这样做的依据吗，欢迎讨论！设置ISPIN=2计算得到的态密度成为自旋态密度。
vasp弹性常数计算过程[精彩]
![vasp弹性常数计算过程[精彩]](https://img.taocdn.com/s3/m/3c7ba60ca9114431b90d6c85ec3a87c240288adc.png)
按照HfN的NaCl结构，弛豫结果：K点最优是21*21*21，ENCUT取550eV，sigma取0.2，EV优化的结果也很理想，最优体积为23.28854892，换算之后，晶格常数是4.533457483，与参考文献一致。
对原来的晶格基矢矩阵做形变[e=(r,r,r,0,0,0)],把形变后的矩阵(矩阵换算规则：R’=R*(δ+e))代替原矩阵放到POSCAR中：000#!/bin/sh#PBS -N vasp#PBS -j oe#PBS -l nodes=1:ppn=12cd ${PBS_O_WORKDIR}source /public/software/profile.d/intel-env.shsource /public/software/profile.d/openmpi-intel-env.shechoecho "Starting V ASP run at" `date`echomaster=`hostname`echo "The job submission node is $master"echo "The working directory is " ${PBS_O_WORKDIR}echo "V ASP input file is" ${PBS_O_WORKDIR}/${inputfile}echoecho "V ASP execution start at" `date`echohostname > grep 'Linux' /etc/issue >> grep 'model name' /proc/cpuinfo |cut -d: -f2 |uniq -c >> grep 'cpu M' /proc/cpuinfo >> grep 'MemTotal' /proc/meminfo >> free -g >> ulimit -a >> cat $PBS_NODEFILE >> NP=`cat $PBS_NODEFILE | wc -l`#################################################################rm W A VECAR 2>/dev/nullechofor i in -0.02 -0.018 -0.016 -0.014 -0.012 -0.01 -0.008 -0.006 -0.004 -0.002 0 0.002 0.004 0.006 0.008 0.01 0.012 0.014 0.016 0.018 0.02; dok=`echo $i | awk '{printf"%.7f",$i*0.5+0.5}'`cat > POSCAR <<!Cubic HfN4.533357680.000000000000000 $k $k 原始矩阵：0 0.5 0.5$k 0.000000000000000 $k 0.5 0 0.5$k $k 0.000000000000 0.5 0.5 01 1Direct0.00000000 0.00000000 0.000000000.50000000 0.50000000 0.50000000!cd /home/yjhao/chenlongqing/r/r1mpirun -np $NP-machinefile $PBS_NODEFILE --mca btl self,sm,tcp -bind-to-core numactl --localalloc /public/software/vasp5.2/vasp >outE=`grep F= OSZICAR` ; echo "r=$i $E" >> TOTALcp CONTCAR POSCAR.$icp OUTCAR OUTCAr.$idonecp INCAR.static INCARcp POSCAR.-0.02 POSCARmpirun -np $NP-machinefile $PBS_NODEFILE --mca btl self,sm,tcp -bind-to-core numactl --localalloc /public/software/vasp5.2/vasp >outE=`grep F= OSZICAR` ; echo "r=-0.02 $E" >> TOTANcp POSCAR.-0.018 POSCARmpirun -np $NP-machinefile $PBS_NODEFILE --mca btl self,sm,tcp -bind-to-core numactl --localalloc /public/software/vasp5.2/vasp >outE=`grep F= OSZICAR` ; echo "r=-0.018 $E" >> TOTANcp POSCAR.-0.016 POSCARmpirun -np $NP-machinefile $PBS_NODEFILE --mca btl self,sm,tcp -bind-to-core numactl --localalloc /public/software/vasp5.2/vasp >outE=`grep F= OSZICAR` ; echo "r=-0.016 $E" >> TOTANcp POSCAR.-0.014 POSCARmpirun -np $NP-machinefile $PBS_NODEFILE --mca btl self,sm,tcp -bind-to-core numactl --localalloc /public/software/vasp5.2/vasp >outE=`grep F= OSZICAR` ; echo "r=-0.014 $E" >> TOTANcp POSCAR.-0.012 POSCARmpirun -np $NP-machinefile $PBS_NODEFILE --mca btl self,sm,tcp -bind-to-core numactl --localalloc /public/software/vasp5.2/vasp >outE=`grep F= OSZICAR` ; echo "r=-0.012 $E" >> TOTANcp POSCAR.-0.01 POSCARmpirun -np $NP-machinefile $PBS_NODEFILE --mca btl self,sm,tcp -bind-to-core numactl --localalloc /public/software/vasp5.2/vasp >outcp POSCAR.-0.008 POSCARmpirun -np $NP-machinefile $PBS_NODEFILE --mca btl self,sm,tcp -bind-to-core numactl --localalloc /public/software/vasp5.2/vasp >outE=`grep F= OSZICAR` ; echo "r=-0.008 $E" >> TOTANcp POSCAR.-0.006 POSCARmpirun -np $NP-machinefile $PBS_NODEFILE --mca btl self,sm,tcp -bind-to-core numactl --localalloc /public/software/vasp5.2/vasp >outE=`grep F= OSZICAR` ; echo "r=-0.006 $E" >> TOTANcp POSCAR.-0.004 POSCARmpirun -np $NP-machinefile $PBS_NODEFILE --mca btl self,sm,tcp -bind-to-core numactl --localalloc /public/software/vasp5.2/vasp >outE=`grep F= OSZICAR` ; echo "r=-0.004 $E" >> TOTANcp POSCAR.-0.002 POSCARmpirun -np $NP-machinefile $PBS_NODEFILE --mca btl self,sm,tcp -bind-to-core numactl --localalloc /public/software/vasp5.2/vasp >outE=`grep F= OSZICAR` ; echo "r=-0.002 $E" >> TOTANcp POSCAR.0 POSCARmpirun -np $NP-machinefile $PBS_NODEFILE --mca btl self,sm,tcp -bind-to-core numactl --localalloc /public/software/vasp5.2/vasp >outE=`grep F= OSZICAR` ; echo "r=0 $E" >> TOTANcp POSCAR.0.002 POSCARmpirun -np $NP-machinefile $PBS_NODEFILE --mca btl self,sm,tcp -bind-to-core numactl --localalloc /public/software/vasp5.2/vasp >outE=`grep F= OSZICAR` ; echo "r=0.002 $E" >> TOTANcp POSCAR.0.004 POSCARmpirun -np $NP-machinefile $PBS_NODEFILE --mca btl self,sm,tcp -bind-to-core numactl --localalloc /public/software/vasp5.2/vasp >outE=`grep F= OSZICAR` ; echo "r=0.004 $E" >> TOTANcp POSCAR.0.006 POSCARmpirun -np $NP-machinefile $PBS_NODEFILE --mca btl self,sm,tcp -bind-to-core numactl --localalloc /public/software/vasp5.2/vasp >outE=`grep F= OSZICAR` ; echo "r=0.006 $E" >> TOTANcp POSCAR.0.008 POSCARmpirun -np $NP-machinefile $PBS_NODEFILE --mca btl self,sm,tcp -bind-to-core numactl --localalloc /public/software/vasp5.2/vasp >outE=`grep F= OSZICAR` ; echo "r=0.008 $E" >> TOTANcp POSCAR.0.01 POSCARmpirun -np $NP-machinefile $PBS_NODEFILE --mca btl self,sm,tcp -bind-to-core numactl --localalloc /public/software/vasp5.2/vasp >outE=`grep F= OSZICAR` ; echo "r=0.01 $E" >> TOTANcp POSCAR.0.012 POSCARmpirun -np $NP-machinefile $PBS_NODEFILE --mca btl self,sm,tcp -bind-to-core numactl --localalloc /public/software/vasp5.2/vasp >outcp POSCAR.0.014 POSCARmpirun -np $NP-machinefile $PBS_NODEFILE --mca btl self,sm,tcp -bind-to-core numactl --localalloc /public/software/vasp5.2/vasp >outE=`grep F= OSZICAR` ; echo "r=0.014 $E" >> TOTANcp POSCAR.0.016 POSCARmpirun -np $NP-machinefile $PBS_NODEFILE --mca btl self,sm,tcp -bind-to-core numactl --localalloc /public/software/vasp5.2/vasp >outE=`grep F= OSZICAR` ; echo "r=0.016 $E" >> TOTANcp POSCAR.0.018 POSCARmpirun -np $NP-machinefile $PBS_NODEFILE --mca btl self,sm,tcp -bind-to-core numactl --localalloc /public/software/vasp5.2/vasp >outE=`grep F= OSZICAR` ; echo "r=0.018 $E" >> TOTANcp POSCAR.0.02 POSCARmpirun -np $NP-machinefile $PBS_NODEFILE --mca btl self,sm,tcp -bind-to-core numactl --localalloc /public/software/vasp5.2/vasp >outE=`grep F= OSZICAR` ; echo "r=0.02 $E" >> TOTAN############################################################################### ##########3exit 0相对应的，TOTAN文件的内容如下：r=-0.02 1 F= -.21835920E+02 E0= -.21835920E+02 d E =0.000000E+00r=-0.018 1 F= -.21837975E+02 E0= -.21837975E+02 d E =0.000000E+00r=-0.016 1 F= -.21839808E+02 E0= -.21839808E+02 d E =0.000000E+00r=-0.014 1 F= -.21841417E+02 E0= -.21841417E+02 d E =0.000000E+00r= -0.012 1 F= -.21842791E+02 E0= -.21842791E+02 d E =0.000000E+00r= -0.01 1 F= -.21843936E+02 E0= -.21843936E+02 d E =0.000000E+00r=-0.008 1 F= -.21844865E+02 E0= -.21844865E+02 d E =0.000000E+00r=-0.006 1 F= -.21845578E+02 E0= -.21845578E+02 d E =0.000000E+00r=-0.004 1 F= -.21846085E+02 E0= -.21846085E+02 d E =0.000000E+00r=-0.002 1 F= -.21846387E+02 E0= -.21846387E+02 d E =0.000000E+00r=0 1 F= -.21846497E+02 E0= -.21846497E+02 d E =0.000000E+00r=0.002 1 F= -.21846390E+02 E0= -.21846390E+02 d E =0.000000E+00r=0.004 1 F= -.21846095E+02 E0= -.21846095E+02 d E =0.000000E+00r=0.006 1 F= -.21845602E+02 E0= -.21845602E+02 d E =0.000000E+00r=0.008 1 F= -.21844905E+02 E0= -.21844905E+02 d E =0.000000E+00r=0.01 1 F= -.21843996E+02 E0= -.21843996E+02 d E =0.000000E+00r=0.012 1 F= -.21842878E+02 E0= -.21842878E+02 d E =0.000000E+00r=0.014 1 F= -.21841544E+02 E0= -.21841544E+02 d E =0.000000E+00r=0.016 1 F= -.21840003E+02 E0= -.21840003E+02 d E =0.000000E+00r=0.018 1 F= -.21838261E+02 E0= -.21838261E+02 d E =0.000000E+00r=0.02 1 F= -.21836314E+02 E0= -.21836314E+02 d E =0.000000E+00将TOTAN中的r和E0提取出来，再用Origin来做拟合（阶数选4，提高精度）可以看到，B2结果是25.228367816，带入下式：V0＝23.29163453C44=2B2*160.2/3V计算结果是C44＝115，与文献数据相符同理，施加形变e=(0,0,0,r,r,r)和e=(r,r,0,0,0)，大体类似，改变的只是变换后的矩阵，分别是：0.00 $k $k$k 0.0 $k$k $k 0.00和0.00 $k 0.50$k 0.0 0.50$k $k 0.00分别将TOTAN文件里面的能量提出，并且转换成，代入Origin里面进行拟合，得到的B2分别是186.4783989和106.79144这样计算之后：C11+2C12＝855（这个数据除以3应等于B0，根据EOS拟合，B0＝285，完全符合）, C11+C12＝734最后得出：C11＝613，C12＝121，C44＝115与参考文献一致：。
vasp在计算磁性的实例和讨论
兄弟，问3个问题1，vasp在计算磁性的时候，oszicar中得到的磁矩和outcar中得到各原子磁矩之和不一致，在投稿的是否曾碰到有审稿人质疑，对于这个不一致你们一般是怎么解释的了？2，另外，磁性计算应该比较负责。
你应该还使用别的程序计算过磁性，与vasp结果比较是否一致，对磁性计算采用的程序有什么推荐。
ps：由于曾使用vasp和dmol算过非周期体系磁性，结构对磁性影响非常大，因此使用这两个程序计算的磁性要一致很麻烦。
还不敢确定到底是哪个程序可能不可靠。
3，如果采用vasp计算磁性，对采用的方法和设置有什么推荐。
1，OSZICAR中得到的磁矩是OUTCAR中最后一步得到的总磁矩是相等的。
总磁矩和各原子的磁矩(RMT球内的磁矩)之和之差就是间隙区的磁矩。
因为有间隙区存在,不一致是正常的。
2，如果算磁性，全电子的结果更精确，我的一些计算结果显示磁性原子对在最近邻的位置时，PAW与FPLAW给出的能量差不一致，在长程时符合的很好。
虽然并没有改变定性结论。
感觉PAW似乎不能很好地描述较强耦合。
我试图在找出原因，主要使用exciting和vasp做比较。
计算磁性推荐使用FP-LAPW, FP-LMTO, FPLO很吸引人(不过是商业的)，后者是O(N)算法。
3，使用vasp计算磁性，注意不同的初始磁矩是否收敛为同一个磁矩。
倒没有特别要注意的地方，个人认为。
归根结底，需要一个优秀的交换关联形式出现VASP计算是否也是像计算DOS和能带一样要进行三步（结构优化，静态自洽计算，非自洽计算），然后看最后一步的出的磁矩呢？一直想计算固体中某个原子的磁矩，根据OUTCAR的结果似乎不能分析，因为它里面总磁矩跟OSZICAR的值有一定的差别，据说是OUTCAR中只考虑WS半径内磁矩造成的。
最近看到一个帖子说是可以用bader电荷分析方法分析原子磁矩。
如法炮制之后发现给出的总磁矩与OSZICAR的结果符合的甚好，可是觉得没有根据，有谁知道这样做的依据吗，欢迎讨论！设置ISPIN=2计算得到的态密度成为自旋态密度。
Vasp入门+实例
0.6 0.5
DOS
0.4 0.3 0.2 0.1 -6 -4 -2 0 2 4 6 8 10
Energy
(4). 做非自洽计算, 求电子结构
• 修改INCAR文件: 将参数ICHARG设为 11 • 修改KPOINTS输入文件
• 运行VASP程序,从输出文件EIGENVAL中提出电子结构
� 1� a1 � a ( i � 2 � 1� a2 � a ( i � 2 � � a3 � ck
第一原理电子结构计算程序：VASP
• 程序原理
• 输入文件
• 输出文件 • 应用
输入文件
POTCAR KPOINTS POSCAR INCAR
Choosing POTCAR file LDA GGA PAW_LDA PAW_GGA PAW_PBE(VASP4.5)
pseudopotentail file Brillouin zone sampling structural data steering parameters
POSCAR输入文件: 原胞中的原子位置
Diamond Si 3.9 0.0 0.5 0.5 0.5 0.0 0.5 0.5 0.5 0.0 1 Direct 0.0 0.0 0.0
基矢的公因子
基矢a1 基矢a2
基矢a3 原胞中的原子个数 坐标系选为基矢构成的坐标系
基矢坐标系下原子的位置
� 1 � � a1 � a( j � k ) 2 � 1 � � a2 � a(i � k ) 2 � 1 � � a3 � a(i � j ) 2
1
0
-1
0 0.07 0.14 0.21 0.28 0.34 0.41 0.48 0.55
vasp计算实例
USPP-GGA/ USPP-LDA/[yzhang@sugon test]$ cp ~/POTCAR-backup/USPP-GGA/AlAl/ Al_h/[yzhang@sugon test]$ cp ~/POTCAR-backup/USPP-GGA/AlAl/ Al_h/[yzhang@sugon test]$ cp ~/POTCAR-backup/USPP-GGA/AAg/ Al/ Al_h/ Ar/ As/ Au/[yzhang@sugon test]$ cp ~/POTCAR-backup/USPP-GGA/AAg/ Al/ Al_h/ Ar/ As/ Au/[yzhang@sugon test]$ cp ~/POTCAR-backup/USPP-GGA/AlAl/ Al_h/[yzhang@sugon test]$ cp ~/POTCAR-backup/USPP-GGA/AlAl/ Al_h/[yzhang@sugon test]$ cp ~/POTCAR-backup/USPP-USPP-GGA/ USPP-LDA/[yzhang@sugon test]$ cp ~/POTCAR-backup/USPP-USPP-GGA/ USPP-LDA/[yzhang@sugon test]$ lsINCAR KPOINTS POSCAR[yzhang@sugon test]$ cd[yzhang@sugon ~]$ lsbin POTCAR-backup yzhang[yzhang@sugon ~]$ cd POTCAR-backup/[yzhang@sugon POTCAR-backup]$ lsPA W-91 PA W-LDA PA W-PBE USPP-GGA USPP-LDA [yzhang@sugon POTCAR-backup]$ cd PA W-PBE/Display all 168 possibilities? (y or n)[yzhang@sugon POTCAR-backup]$ cd PA W-PBE/AlAl/ Al_h/[yzhang@sugon POTCAR-backup]$ cd PA W-PBE/AlAl/ Al_h/[yzhang@sugon POTCAR-backup]$ cd PA W-PBE/Al/DDE.Z POTCAR PSCTR.Z V_RHFIN.Z V_TABIN.Z [yzhang@sugon POTCAR-backup]$ cd PA W-PBE/Al/DDE.Z POTCAR PSCTR.Z V_RHFIN.Z V_TABIN.Z [yzhang@sugon POTCAR-backup]$ cd PA W-PBE/Al/DDE.Z POTCAR PSCTR.Z V_RHFIN.Z V_TABIN.Z [yzhang@sugon POTCAR-backup]$[yzhang@sugon POTCAR-backup]$ ;ls-bash: syntax error near unexpected token `;'[yzhang@sugon POTCAR-backup]$ lsPA W-91 PA W-LDA PA W-PBE USPP-GGA USPP-LDA [yzhang@sugon POTCAR-backup]$ cd[yzhang@sugon ~]$ lsbin POTCAR-backup yzhang[yzhang@sugon ~]$ cd yzhang/[yzhang@sugon yzhang]$ lstest[yzhang@sugon yzhang]$ cd test/[yzhang@sugon test]$ lsINCAR KPOINTS POSCAR[yzhang@sugon test]$ cp ~/POTCAR-backup/PA W-PA W-91/ PA W-LDA/ PA W-PBE/[yzhang@sugon test]$ cp ~/POTCAR-backup/PA W-PBE/AlAl/ Al_h/[yzhang@sugon test]$ cp ~/POTCAR-backup/PA W-PBE/AlAl/ Al_h/[yzhang@sugon test]$ cp ~/POTCAR-backup/PA W-PBE/Al/POTCAR .[yzhang@sugon test]$ lsINCAR KPOINTS POSCAR POTCAR[yzhang@sugon test]$ lsINCAR KPOINTS POSCAR POTCAR[yzhang@sugon test]$ lsINCAR KPOINTS POSCAR POTCAR[yzhang@sugon test]$ lsINCAR KPOINTS POSCAR POTCAR[yzhang@sugon test]$ vaspvasp.4.6.34 5Dec07 complexPOSCAR found : 1 types and 1 ionsLDA part: xc-table for Pade appr. of PerdewPOSCAR, INCAR and KPOINTS ok, starting setupW ARNING: wrap around errors must be expectedFFT: planning (1)reading W A VECARentering main loopN E dE d eps ncg r ms rms(c)DA V: 1 -0.308753898526E+01 -0.30875E+01 -0.59037E+02 430 0.159E+02 DA V: 2 -0.372313582553E+01 -0.63560E+00 -0.58860E+00 485 0.166E+01 DA V: 3 -0.372697132889E+01 -0.38355E-02 -0.38336E-02 440 0.126E+00 DA V: 4 -0.372697751054E+01 -0.61817E-05 -0.61813E-05 470 0.546E-02 DA V: 5 -0.372697751835E+01 -0.78117E-08 -0.79994E-08 415 0.158E-03 0.933E-01DA V: 6 -0.371445979461E+01 0.12518E-01 -0.14507E-03 415 0.247E-01 0.580E-01DA V: 7 -0.370504809115E+01 0.94117E-02 -0.47985E-03 410 0.428E-01 0.252E-02DA V: 8 -0.370505518255E+01 -0.70914E-05 -0.17041E-05 285 0.435E-021 F= -.37050552E+01 E0= -.37047594E+01 d E =-.887456E-03writing wavefunctions[yzhang@sugon test]$ lsCHG CONTCAR EIGENV AL INCAR OSZICAR PCDA T POTCAR WA VECARCHGCAR DOSCAR IBZKPT KPOINTS OUTCAR POSCAR vasprun.xml XDA TCAR[yzhang@sugon test]$ lsCHG CONTCAR EIGENV AL INCAR OSZICAR PCDA T POTCAR WA VECARCHGCAR DOSCAR IBZKPT KPOINTS OUTCAR POSCAR vasprun.xml XDA TCAR[yzhang@sugon test]$ vi OUTCAR[yzhang@sugon test]$ vi KPOINTS[yzhang@sugon test]$ lsCHG CONTCAR EIGENV AL INCAR OSZICAR PCDA T POTCAR WA VECARCHGCAR DOSCAR IBZKPT KPOINTS OUTCAR POSCAR vasprun.xml XDA TCAR[yzhang@sugon test]$ vaspvasp.4.6.34 5Dec07 complexPOSCAR found : 1 types and 1 ionsLDA part: xc-table for Pade appr. of Perdewfound W A VECAR, reading the headernumber of k-points has changed, file: 35 present: 220trying to continue reading WA VECAR, but it might failPOSCAR, INCAR and KPOINTS ok, starting setupW ARNING: wrap around errors must be expectedFFT: planning (1)reading W A VECARthe W A VECAR file was read sucessfullyinitial charge from wavefunctionentering main loopN E dE d eps ncg rms rms(c)DA V: 1 -0.338903546871E+01 -0.33890E+01 -0.38722E+02 3085 0.566E+01BRMIX: very serious problemsthe old and the new charge density differold charge density: 3.95554 new 3.000000.289E+00DA V: 2 -0.372640269291E+01 -0.33737E+00 -0.66811E-01 2640 0.444E+00 0.995E-01DA V: 3 -0.373013183172E+01 -0.37291E-02 -0.10958E-02 2705 0.823E-01 0.151E-01DA V: 4 -0.373004294057E+01 0.88891E-04 -0.84537E-05 2610 0.807E-021 F= -.37300429E+01 E0= -.37296573E+01 d E =-.115689E-02writing wavefunctions[yzhang@sugon test]$ lsCHG CONTCAR EIGENV AL INCAR OSZICAR PCDA T POTCAR WA VECARCHGCAR DOSCAR IBZKPT KPOINTS OUTCAR POSCAR vasprun.xml XDA TCAR[yzhang@sugon test]$ lsCHG CONTCAR EIGENV AL INCAR OSZICAR PCDA T POTCAR WA VECARCHGCAR DOSCAR IBZKPT KPOINTS OUTCAR POSCAR vasprun.xml XDA TCAR[yzhang@sugon test]$ lsCHG CONTCAR EIGENV AL INCAR OSZICAR PCDA T POTCAR WA VECARCHGCAR DOSCAR IBZKPT KPOINTS OUTCAR POSCAR vasprun.xml XDA TCAR[yzhang@sugon test]$ vi In[yzhang@sugon test]$ ;ls-bash: syntax error near unexpected token `;'[yzhang@sugon test]$ lsCHG CONTCAR EIGENV AL INCAR OSZICAR PCDA T POTCAR WA VECARCHGCAR DOSCAR IBZKPT KPOINTS OUTCAR POSCAR vasprun.xml XDA TCAR[yzhang@sugon test]$ lsCHG CONTCAR EIGENV AL INCAR OSZICAR PCDA T POTCAR WA VECARCHGCAR DOSCAR IBZKPT KPOINTS OUTCAR POSCAR vasprun.xml XDA TCAR[yzhang@sugon test]$ lsCHG CONTCAR EIGENV AL INCAR OSZICAR PCDA T POTCAR WA VECARCHGCAR DOSCAR IBZKPT KPOINTS OUTCAR POSCAR vasprun.xml XDA TCAR[yzhang@sugon test]$ vi INCAR[yzhang@sugon test]$ lsCHG CONTCAR EIGENV AL INCAR OSZICAR PCDA T POTCAR WA VECARCHGCAR DOSCAR IBZKPT KPOINTS OUTCAR POSCAR vasprun.xml XDA TCAR[yzhang@sugon test]$ lsCHG CONTCAR EIGENV AL INCAR OSZICAR PCDA T POTCARWA VECARCHGCAR DOSCAR IBZKPT KPOINTS OUTCAR POSCAR vasprun.xml XDA TCAR[yzhang@sugon test]$ lsCHG CONTCAR EIGENV AL INCAR OSZICAR PCDA T POTCAR WA VECARCHGCAR DOSCAR IBZKPT KPOINTS OUTCAR POSCAR vasprun.xml XDA TCAR[yzhang@sugon test]$ lsCHG CONTCAR EIGENV AL INCAR OSZICAR PCDA T POTCAR WA VECARCHGCAR DOSCAR IBZKPT KPOINTS OUTCAR POSCAR vasprun.xml XDA TCAR[yzhang@sugon test]$ lsCHG CONTCAR EIGENV AL INCAR OSZICAR PCDA T POTCAR WA VECARCHGCAR DOSCAR IBZKPT KPOINTS OUTCAR POSCAR vasprun.xml XDA TCAR[yzhang@sugon test]$ changed, file: 35 present: 220-bash: changed,: command not found[yzhang@sugon test]$ trying to continue reading W A VECAR, but it might fail-bash: trying: command not found[yzhang@sugon test]$ POSCAR, INCAR and KPOINTS ok, starting setup-bash: POSCAR,: command not found[yzhang@sugon test]$ W ARNING: wrap around errors must be expected-bash: W ARNING:: command not found[yzhang@sugon test]$ FFT: planning (1)-bash: FFT:: command not found[yzhang@sugon test]$ reading W A VECAR-bash: reading: command not found[yzhang@sugon test]$ changed, file: 35 present: 220-bash: changed,: command not found[yzhang@sugon test]$ trying to continue reading W A VECAR, but it might fail-bash: trying: command not found[yzhang@sugon test]$ POSCAR, INCAR and KPOINTS ok, starting setup-bash: POSCAR,: command not found[yzhang@sugon test]$ W ARNING: wrap around errors must be expected-bash: W ARNING:: command not found[yzhang@sugon test]$ FFT: planning (1)-bash: FFT:: command not found[yzhang@sugon test]$ reading W A VECAR-bash: reading: command not found[yzhang@sugon test]$ POINTS OUTCAR POSCAR vasprun.xml XDA TCAR-bash: POINTS: command not found[yzhang@sugon test]$ [yzhang@sugon test]$ ls-bash: [yzhang@sugon: command not found[yzhang@sugon test]$ CHG CONTCAR EIGENV AL INCAR OSZICAR PCDA T POTCAR W A VECAR-bash: CHG: command not found[yzhang@sugon test]$ CHGCAR DOSCAR IBZKPT KPOINTPOINTS OUTCAR POSCAR vasprun.xml XDA TCAR-bash: CHGCAR: command not found[yzhang@sugon test]$ [yzhang@sugon test]$ ls-bash: [yzhang@sugon: command not found[yzhang@sugon test]$ CHG CONTCAR EIGENV AL INCAR OSZICAR PCDA T POTCAR W A VECAR-bash: CHG: command not found[yzhang@sugon test]$ CHGCAR DOSCAR IBZKPT KPOINTg: command not found-bash: CHGCAR: command not found[yzhang@sugon test]$ [yzhang@sugon test]$ POSCAR, INCAR and KPOINTS ok, starting setup-bash: W ARNING:: command not found-bash: [yzhang@sugon: command not found[yzhang@sugon test]$ -bash: POSCAR,: command not found[yzhang@sugon test]$ FFT: planning (1)-bash: -bash:: command not found[yzhang@sugon test]$ [yzhang@sugon test]$ W ARNING: wrap around errors must be expected -bash: [yzhang@sugon: command not found[yzhang@sugon test]$ -bash: W ARNING:: command not found-bash: -bash:: command not found[yzhang@sugon test]$ [yzhang@sugon test]$ FFT: planning (1)-bash: [yzhang@sugon: command not found[yzhang@sugon test]$ -bash: FFT:: command not found-bash: changed,: command not found-bash: -bash:: command not found[yzhang@sugon test]$ [yzhang@sugon test]$ reading WA VECAR-bash: [yzhang@sugon: command not found[yzhang@sugon test]$ -bash: reading: command not found-bash: -bash:: command not found[yzhang@sugon test]$ [yzhang@sugon test]$ changed, file: 35 present: 220-bash: [yzhang@sugon: command not found[yzhang@sugon test]$ -bash: changed,: command not found-bash: -bash:: command not found[yzhang@sugon test]$ [yzhang@sugon test]$ trying to continue reading WA VECAR, but it might fail-bash: [yzhang@sugon: command not found[yzhang@sugon test]$ -bash: trying: command not found-bash: -bash:: command not found[yzhang@sugon test]$ [yzhang@sugon test]$ POSCAR, INCAR and KPOINTS ok, starting setup-bash: [yzhang@sugon: command not found[yzhang@sugon test]$ -bash: POSCAR,: command not found-bash: W ARNING:: command not found-bash: -bash:: command not found[yzhang@sugon test]$ FFT: planni[yzhang@sugon test]$ [yzhang@sugon test]$ W ARNING: wrap around errors must be expected-bash: [yzhang@sugon: command not found[yzhang@sugon test]$ -bash: W ARNING:: command not found-bash: -bash:: command not found[yzhang@sugon test]$ [yzhang@sugon test]$ FFT: planni[yzhang@sugon test]$。
VASP几个计算实例
用VASP计算Ｈ原子的能量氢原子的能量为。
在这一节中，我们用VASP计算H原子的能量。
对于原子计算，我们可以采用如下的INCAR文件PREC=ACCURATENELMDL = 5 make five delays till charge mixingISMEAR = 0; SIGMA=0.05 use smearing method采用如下的KPOINTS文件。
由于增加K点的数目只能改进描述原子间的相互作用，而在单原子计算中并不需要。
所以我们只需要一个K点。
Monkhorst Pack 0 Monkhorst Pack1 1 10 0 0采用如下的POSCAR文件atom 115.00000 .00000 .00000.00000 15.00000 .00000.00000 .00000 15.000001cart0 0 0采用标准的H的POTCAR得到结果如下：k-point 1 : 0.0000 0.0000 0.0000band No. band energies occupation1 -6.3145 1.000002 -0.0527 0.000003 0.4829 0.000004 0.4829 0.00000我们可以看到，电子的能级不为。
Free energy of the ion-electron system (eV)---------------------------------------------------alpha Z PSCENC = 0.00060791Ewald energy TEWEN = -1.36188267-1/2 Hartree DENC = -6.27429270-V(xc)+E(xc) XCENC = 1.90099128PAW double counting = 0.00000000 0.00000000entropy T*S EENTRO = -0.02820948eigenvalues EBANDS = -6.31447362atomic energy EATOM = 12.04670449---------------------------------------------------free energy TOTEN = -0.03055478 eVenergy without entropy = -0.00234530 energy(sigma->0) = -0.01645004我们可以看到也不等于。
VASP计算稀土[教学]
![VASP计算稀土[教学]](https://img.taocdn.com/s3/m/23653bd15ebfc77da26925c52cc58bd6318693f3.png)
实例1VASP算稀土永磁材料的磁学性能用哪种算法和赝势比较好？用VASP计算稀土永磁材料（比如Sm-Co）的磁学性能用哪种算法和赝势比较好啊？LDA GGALSDA+U?PBEPW91?PBE是比较好的交换关联能，但是对于磁性，最后加上U结果可能会好点。
但是U的确定需要从文献和其他软件得到我算的磁性没有f电子，这是为什么呢？是赝势的问题还是将f电子限制在芯内了？f电子的确是很深的，一般很难和其他原子的电子相互作用，这也是La系和Ac的元素的化学表现很相似的原因那么怎样才能使磁性计算出现f电子呢？确实让人纠结啊！请问使用PBEsol+U进行优化和性质计算,如何设置INCAR？在vasp5.2手册上找不到PBEsol+U的说明，只有LSDA+U的PBE是GGA类的交换关联能，LSDA的设置是可以同样用于GGA的实例2vasp计算中sigmma值稀土金属怎么取vasp计算中sigmma值稀土金属怎么取啊？计算出来老感觉不对。
取不同的sigma测试，然后，计算结果中取能量的哪一行，sigma-->0跟不趋近0的时候的比较，差别满足你的精度需求就是实例3关于VASP计算用不同赝势产生的能量差异！为啥同一个结构，用不同的赝势文件POTCAR，如PAW_PBE赝势和用US赝势来计算，为啥能量不一样？连初始第一步的能量的差别就很大啊？本人理论知识很浅，各位大侠能说说其中原理吗？这不只是精度的问题，因为能量就不在一个层次上！两套赝势的能量没有可比性.Generally the PAW potentials are more accurate than the ultra-soft pseudopotentials. There are two reasons forthis: first, the radial cutoffs (core radii) are smaller than the radii used for the USpseudopotentials, and second thePAW potentials reconstruct the exact valence wave function with a ll nodes in the core region.能量绝对值没有任何意义的，不同赝势能量参考态不一样，只有能量之差才有意义。
- 1、下载文档前请自行甄别文档内容的完整性，平台不提供额外的编辑、内容补充、找答案等附加服务。
- 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
- 3、如文档侵犯您的权益，请联系客服反馈,我们会尽快为您处理(人工客服工作时间：9:00-18:30)。
VASP计算实例目录一、氢气分子H2键长的计算 (3)1.基本文件 (3)2.赝势类型的选择 (3)3.截断能ENCUT参数的选择 (4)4.KPOINTS参数选择 (5)5.对晶格常数进行优化 (6)二、Si晶体晶格常数计算 (8)1.赝势类型选择 (8)2.截断能（ENCUT）参数的选定 (9)3.KPOINTS参数的选定 (11)4.SIGMA参数的选定 (12)5.晶格常数计算结果 (13)三、Si元素单原子能量计算 (14)1.由内聚能倒推单原子能量 (14)2.基本文件 (15)3.单原子能量计算 (15)四、Si的VASP力学常数计算 (16)1.计算所需文件 (16)2.计算与数据处理 (17)3.计算所用到的公式： (18)五、SI晶体的电子结构 (19)1.采用VASP计算能带的步骤 (19)2.电荷分布计算结果 (20)能带计算和结果 (21)3.态密度计算和结果 (21)六、Si晶体介电函数和光学性质的计算 (22)1.计算步骤 (22)2.用到的文件 (23)3.计算结果 (26)七、VASP的声子谱计算 (29)1.计算步骤 (29)2.基本文件 (30)3.声子谱、声子态密度计算和结果 (33)4.热学性质计算和结果 (34)八、化合物co2键长计算 (35)1.计算步骤 (35)2.基本文件 (35)一、氢气分子H2键长的计算1.基本文件准备基本文件INCAR、POTCAR、POSCAR、KPOINT以及脚本文件encut、k、optimize2.赝势类型的选择输入文件如下其中参数要靠经验初选INCAR:System = F2ISTART = 0ICHARG = 2NELMDL = 5ISMEAR = 0SIGMA = 0.1PREC = AccurateKPOINTS：Automatic meshM1 1 10 0 0POSCAR：O115.0 0.00 0.000.00 14.0 0.000.00 0.00 13.01D0.00 0.00 0.00分别选用五个贋势文件进行计算。
计算方法：把INCAR、POSCAR、POTCAR、KPOINTS四个文件分别放入同一个文件夹中，在终端中输入“VASP”得到五个能量。
GGA-91 E0= -.36620018E+01GGA-CA E0= -.35466314E+01PAW E0= -.34415634E+01PAW-GGA E0= -.34948621E+01PBE E0= -.35753416E+01选择单原子能量最低的，晶体最稳定，故选择GGA-91赝势文件POTCAR。
3.截断能ENCUT参数的选择准备文件KPOINTS、POSCAR、POTCAR以及脚本文件encut并放入同一个文件夹中encut：#!/bin/shrm WAVECARfor i in 550 600 650 700 750 800 850 900docat > INCAR <<!System = H2ISTART=0ISMEAR=0;SIGMA=0.1ICHARG = 2ENCUT=$iPREC=Accurateecho "ENCUT = $i eV" ; vaspE=`grep "TOTEN" OUTCAR | tail -1 | awk '{printf "%12.6f\n",$5}'` KP=`grep "irreducible" OUTCAR | tail -1 | awk '{printf "%5i \n", $2 }'`echo $i $KP $E >> commentdoneKPOINTS：Automatic meshM1 1 10 0 0POSCAR：F21.2094.000.00 0.000.00 4.00 0.000.000.00 4.002Direct0.000.00 0.000.000.00 0.25选用上次计算是最低能量的POTCAR运行脚本文件encut，得到comment文件内容如下：550 1 -3.654087600 1 -3.654701650 1 -3.651275700 1 -3.652794750 1 -3.652555800 1 -3.652600850 1 -3.653120900 1 -3.652554950 1 -3.652669总能变化在0.001eV左右就足够了，因此我们选择encut=700eV4.KPOINTS参数选择准备文件POSCA、INCAR、POTCAR以及脚本文件k，POSCAR、POTCAR与上面相同脚本文件k：#!/bin/shrm WAVECARfor i in 3 5 7 9 11 13 15docat > KPOINTS <<!Automatic meshM$i $i $i0 0 0!echo "k mesh = $i x $i x $i" ; vaspE=`grep "TOTEN" OUTCAR | tail -1 | awk '{printf "%12.6f \n", $5 }'`KP=`grep "irreducible" OUTCAR | tail -1 | awk '{printf "%5i \n", $2 }'`echo $i $KP $E >>commentdoneINCAR：System = O2ENCUT = 700ISTART = 0ICHARG = 2ISMEAR = 0SIGMA = 0.1PREC = Accurate运行脚本文件k，得到comment文件：3 6 -3.6523825 18 -3.6524117 40 -3.6523889 75 -3.652386总能变化在0.001eV左右就足够了，因此我们选择5作为k点值。
5.对晶格常数进行优化准备文件INCAR、KPOINTS、POTCAR以及脚本文件optimize，其余文件与上一样Optimize：#! /bin/shBIN=/usr/bin/vasprm -rf WAVECARfor i in 1.16 1.21 1.26 1.31 1.36 1.41 1.46 1.51 1.56 1.61 1.66;do cat >POSCAR <<!O2$i4.00 0.00 0.000.00 4.00 0.000.00 0.00 4.002Direct0.00 0.00 0.000.00 0.00 0.25!echo "a= $i" ;$BINE=`tail -1 OSZICAR`;echo $i $E >>RESULT.fccdonecat RESULT.fcccat RESULT.fcc运行脚本文件optimize，得到RESULT文件如下：1.16 1 F= -.15006150E+01 E0= -.15006150E+01 d E =-.179658E-311.21 1 F= -.24084278E+01 E0= -.24084278E+01 d E =-.145599E-161.26 1 F= -.30094870E+01 E0= -.30094870E+01 d E =-.247133E-211.31 1 F= -.33821346E+01 E0= -.33821346E+01 d E =-.403048E-161.36 1 F= -.35847564E+01 E0= -.35847564E+01 d E =-.800387E-691.41 1 F= -.36607958E+01 E0= -.36607958E+01 d E =-.270236E-211.46 1 F= -.36424044E+01 E0= -.36424044E+01 d E =-.866906E-241.51 1 F= -.35553122E+01 E0= -.35553122E+01 d E =-.289397E-271.56 1 F= -.34180694E+01 E0= -.34180694E+01 d E =-.998984E-161.61 1 F= -.32473112E+01 E0= -.32473112E+01 d E =-.696853E-231.66 1 F= -.30547759E+01 E0= -.30547759E+01 d E =-.540502E-15利用计算数据拟合二次曲线最低点就是计算出的晶格常数。
拟合后的曲线如下图：计算出晶格常数是1.41实验室是同样是1.41故参数选择合适计算准确。
二、Si晶体晶格常数计算1.赝势类型选择分别带入五个赝势库文件带入计算去能量最小值的POTCAR。
输入文件：INCAR:System = SiISTART = 0ENCUT = 240NELM = 100EDIFF = 1E-04EDIFFG = 1E-03NSW = 5IBRION = 2KPOINTS:Automatic meshMonkhorst Pack9 9 90 0 0POSCAR:Si5.430.00 0.50 0.500.50 0.00 0.500.50 0.50 0.002Direct0.00 0.00 0.000.25 0.25 0.25计算结果为：GGA-CA赝势文件POTCAR计算能量最低所以。
E0=-.11947589E+02eV即优化后Si晶体的总能量。
2.截断能（ENCUT）参数的选定利用encut脚本文件计算ENCUT参数，准备POSCAR,POTCAR,KPOINT文件。
双击encut进行计算。
结果文件在comment文件中。
选取变化在0.001的ENCUT值。
所用文件：POSCAR：Si5.430.00 0.50 0.500.50 0.00 0.500.50 0.50 0.002Direct0.00 0.00 0.000.25 0.25 0.25KPOINTS：Automatic meshMonkhorst Pack9 9 90 0 0POTCAR从赝势库GGA-CA里选择所对应的文件脚本文件encut：#!/bin/shrm WAVECARfor i in 150 200 250 300 350docat> INCAR <<!System = SiENCUT = $iISTART = 0ICHARG = 2ISMEAR = 0SIGMA = 0.1PREC = Accurate!echo "ENCUT = $i eV" ; vaspE=`grep "TOTEN" OUTCAR | tail -1 | awk '{printf "%12.6f\n",$5}'` KP=`grep "irreducible" OUTCAR | tail -1 | awk '{printf "%5i \n", $2 }'`echo $i $KP $E >> commentdone运行脚本文件encut，计算结果位于comment文件，内容如下：150 35 -11.900756200 35 -11.937845250 35 -11.945420300 35 -11.945545350 35 -11.945502选择总能变化在0.001eV左右就足够了。