搜档网
当前位置:搜档网 › VASP 中INCAR设置

VASP 中INCAR设置

VASP 中INCAR设置
VASP 中INCAR设置

INCAR(参数的首字母代表的参数性质:I初始化 L逻辑开关 E能量 N数目 T温度)一般要设置的关键词为:SYSTEM,ENCUT,ISTART,ICHARG,PREC,ISMEAR,SIGMA。

缺省值就是默认值。是指一个属性、参数在被修改前的初始值。

SYSTEM

任务的名字

ENCUT

平面波的切断能(energy cutoff in eV): default taken from POTCAR-file impor tant! 重要到几乎最好不要手工去设置除非文献告诉你要用多少,或者经过结果可靠性的验证当然,为了测试一下提交的任务,也不妨先设个较小的值附加说明:当且仅当POT CAR里头没有设置ENCUT时(其实貌似没有才是常态),才受PREC设置影响从POTCAR 里找出相应的ENMAX/ENMIN值来设置。PREC= Low Medium Accurate High ENC UT= ENMIN ENMAX ENMAX 130%ENMAX 对于多个元素的POTCAR不同的ENMA X/ENMIN,都取最大值

ISTART

决定是否读取WAVECAR文件。取值:0~2,缺省0/1 for无/有前次计算的WAVECAR(波函数)

0:begin 'from scratch',根据INIWAV初始化波函数

1:restart with constant energy cut-off,从WAVECAR读取波函数(重定义平面波集)2:restart with constant basis set,从WAVECAR读取波函数(平面波集不变)

在进行能带结构、电子态密度等性质计算时设置为1,其他情况一般都设置为0。ICHARG

决定如何建立初始电荷密度。取值:0~2,缺省值: if ISTART=0 2 else 0

0:由初始波函数计算电荷密度

1:从CHGCAR文件读取电荷密度

2:使用原子电荷密度的叠加

+10非自洽计算

PREC

precession: medium, high or low(VASP.4.5+ also: normal, accurate) Default: Medium VASP4.5+采用了优化的accurate来替代high,所以一般不推荐使用 high。不过high可以确保'绝对收敛',作为参考值有时也是必要的。同样受推荐的是normal,作为日常计算选项,可惜的是说明文档提供的信息不足。受PREC影响的参数有四类:ENCUT; NGX,NGY,NGZ; NGXF,NGYF,NGZF; ROPT; 如果设置了PREC,这些参数就都不需要出现了当然直接设置相应的参数也是同样效果的,这里不展开了,随后详释

ISMEAR

决定每个波函数的部分占位f nk如何设置。取值:-5 | -4 | -3 | -2 | 0 | N

-5-带有Blochl修正的四面体方法,

-4-不带Blochl修正的四面体方法,

-3-根据INCAR文件中提供的smearing参数执行循环,

-2-从WAVECAR文件中读取,

-1-Fermi-smearing,

0-Gaussian smearing,

>0-method of Methfessel-Paxton order N,缺省值:1采用部分占位波函数,用一个函数来平滑积分,尤其是对于金属体系可减少k点。采用所谓部分占有波函数,用一个函数来平滑积分,尤其是对于金属体系可减少k点

SIGMA

决定smearing的宽度(eV)。取值:实数,缺省值:0.2

其他参数:

NWRITE

输出内容详细程度。取值:0~4,缺省值:2如果是做长时间动力学计算的话,最好选0或1(首末步/每步核运动输出),短时运算用2,选3则会在出错的时候给出说明信息。

ISPIN

是否进行spin polarized calculation。取值:1,2(1-no,2-yes),缺省值:2 MAGMOM

在ICHARG=2或在CHGCAR中未包含磁化密度(ICHARG=1)时,指定每个原子的初始磁化时刻。取值:实数数组,缺省值:对ISPIN=2 NIONS*1.0,对非共线型磁化体系

3*NIONS*1.0

INIWAV

如何设置初始波函数,只在ISTART=0时使用。取值:0,1(0-最低动能的平面波,1-随机数),缺省值:1。

IDIPOL

控制计算单极、偶极和四极修正。取值:1~4。

1~3只计算第一/二/三晶矢方向,适于厚板(slab)的计算

4 所有方向都计算,适于计算孤立分子

LREAL

决定投射是在实空间还是倒易空间进行。取值:.TRUE.(实空间)/.FALSE.(倒易空间),缺省值:.FALSE.

用于求解赝势的非局域部分用到的一个积分,在倒格空间里采用平面波基组求解,在实空间里则采用积分球求解。

其他还有两个选项:O or On,A or Auto。

On和.TRUE.的差别在于是否使用King-Smith算法优化,设为Auto则进行自动选择,推荐使用。

ROPT

在LREAL=Auto or On时,优化控制每个核周围的积分球内的格点数。取值:实数数组For LREAL=On

PREC=Low, 700 points in the real space sphere (ROPT=0.67)

PREC=Med, 1000 points in the real space sphere (ROPT=1.0)

PREC=High, 1500 points in the real space sphere (ROPT=1.5)

For LREAL=Auto

PREC=Low, accuracy 10-2 (ROPT=0.01)

PREC=Med, accuracy 2*10-3 (ROPT=0.002)

PREC=High accuracy 2*10-4 (ROPT=2E-4)

EDIFF

指定电子自洽循环的全局中断条件,用于控制收敛精度。取值:实数,缺省值:10-4

注意,即使EDIFF=0,NELM步电子自洽循环也会执行。

是电子结构部分自洽迭代循环是,判断是否自洽了的条件,上次和当前两次迭代中总能和本征值的变化都小于EDIFF,则电子构件部分迭代循环停止。如果EDIFF=0,则进行NELM步迭代后停止迭代。

EDIFFG

指定离子弛豫循环的中断条件,用于控制核运动的收敛精度。取值,实数,缺省值:10*EDIFF EDIFFG>0在两个离子步的总自由能之差小于EDIFFG时停止

EDIFFG<0在所有的力都小于EDIFFG时停止。

EDIFFG=0在NSW步弛豫后停止

此参数不支持MD,仅用于弛豫。

ISIF = 3

calculate stress and what to relax

NSW

指定核运动步数(最大离子步数)缺省:0

NBLOCK,KBLOCK

取值:整数,缺省值:NBLOCK = 1,KBLOCK = NSW

在NBLOCK离子步后对成对相关函数和DOS进行计算,并且把离子配置写入XDATCAR 文件。

在KBLOCK*NBLOCK步主循环后平均的成对相关函数和DOS被写入PCDAT和DOSCAR 文件。

IBRION

决定离子怎样更新和运动。取值:-1~3,5~8(-1-无更新,0-MD,1-RMM-DIIS,2-共轭梯度算法,3-Damped MD,5,6:有限差分,7,8:密度函数扰动理论),

缺省值:if NSW=0/1,-1,else 0

这个参数是和ISIF, IALGO/ALGO一起决定怎么算的最重要的参数。

~3是三种弛豫的方法,根据ISIF决定是否固定离子位置、晶胞大小和形状,在INCAR中必须设置参数POTIM。

0是标准的ab-initio MD,不受ISIF影响,即不改变晶胞大小和形状。

5~8支持Hessian Matrix和phonon frequency的计算以及部分固定的MD。

NGX ;NGY ;NGZ

控制FFT网格在三个晶矢方向上的格点数量。

NGFX,NGFY,NGFZ

控制第二次更精确的FFT网格的格点数量。

也是两类重要的最好不要去动的参数。在未指定的情况下将根据PREC的设置从POTCAR 中自动读取。

PREC=High/Accurate,基组中向量的2倍值,用来避免wrap around errors,得到精确解。

PREC=Low/Medium/Normal,基组中向量的3/4倍值(已足够精确到1 meV/atom)。POTIM

BRION=0时,给出MD每步步长(fs),IBRION=1~3时,给出最小化的度量常量。取值:实数,缺省值:IBRION=0无缺省,必须指定,IBRION=1,2,3 0.5

RWIGS

Wigner-Seitz半径DOS计算用。用在计算分波态密度以及每根能带对应的波函数按spd和位置投影

NPAR

NPAR 并行计算band的节点数,每一个节点计算一个band当然可以提高并行效率,减少通讯量,不过貌似现在硬件的主要限制还是内存,而这个选项的使用可能会大幅增加内存的需求

ISYM = 0

决定是否计算应力张量以及弛豫中晶胞变化的自由度。取值:0~6,缺省值:if IBRION=0(MD) 0 else 2 ;0-nonsym 1-usesym 是否使用对称性

ISIF│calculate │ calculate │relax │ change │ change

│ force │stress tensor │ions │ cell shape │cell volume

──┼─────┼───────┼───┼──────┼──────

0 │ yes │ no │yes │ no │no

1 │ yes │ trace only │yes │ no │no

2 │ yes │ yes │yes │ no │no

3 │ yes │ yes │yes │ yes │yes

4 │ yes │ yes │yes │ yes │no

5 │ yes │ yes │no │ yes │no

6 │ yes │ yes │no │ yes │yes

7 │ yes │ yes │no │ no │yes

vasp参数详解

转载]VASP参数设置详解 已有1283 次阅读2011-9-15 14:26 |个人分类:vasp|系统分类:科研笔记 转自小木虫,略有增减 软件主要功能: 采用周期性边界条件(或超原胞模型)处理原子、分子、团簇、纳米线(或管)、薄膜、晶体、准晶和无定性材料,以及表面体系和固体 l 计算材料的结构参数(键长、键角、晶格常数、原子位置等)和构型 l 计算材料的状态方程和力学性质(体弹性模量和弹性常数) l 计算材料的电子结构(能级、电荷密度分布、能带、电子态密度和ELF) l 计算材料的光学性质 l 计算材料的磁学性质 l 计算材料的晶格动力学性质(声子谱等) l 表面体系的模拟(重构、表面态和STM模拟) l 从头分子动力学模拟 l 计算材料的激发态(GW准粒子修正) 计算主要的四个参数文件:INCAR ,POSCAR,POTCAR ,KPOINTS,下面简要介绍,详细权威的请参照手册 INCAR文件: 该文件控制VASP进行何种性质的计算,并设置了计算方法中一些重要的参数,这些参数主要包括以下几类: l 对所计算的体系进行注释:SYSTEM l 定义如何输入或构造初始的电荷密度和波函数:ISTART,ICHARG,INIWAV l 定义电子的优化 –平面波切断动能和缀加电荷时的切断值:ENCUT,ENAUG –电子部分优化的方法:ALGO,IALGO,LDIAG –电荷密度混合的方法:IMIX,AMIX,AMIN,BMIX,AMIX_MAG,BMIX_MAG,WC,INIMIX,MIXPRE,MAXMIX –自洽迭代步数和收敛标准:NELM,NELMIN,NELMDL,EDIFF l 定义离子或原子的优化 –原子位置优化的方法、移动的步长和步数:IBRION,NFREE,POTIM,NSW –分子动力学相关参数:SMASS,TEBEG,TEEND,POMASS,NBLOCK,KBLOCK,PSTRESS –离子弛豫收敛标准:EDIFFG l 定义态密度积分的方法和参数 – smearing方法和参数:ISMEAR,SIGMA –计算态密度时能量范围和点数:EMIN,EMAX,NEDOS –计算分波态密度的参数:RWIGS,LORBIT l 其它 –计算精度控制:PREC –磁性计算:ISPIN,MAGMOM,NUPDOWN –交换关联函数:GGA,VOSKOWN –计算ELF和总的局域势:LELF,LVTOT –结构优化参数:ISIF –等等。 主要参数说明如下:

VASP参数设置详解

VASP参数设置详解 计算材料2010-11-30 20:11:32 阅读197 评论0 字号:大中小订阅 转自小木虫,略有增减 软件主要功能: 采用周期性边界条件(或超原胞模型)处理原子、分子、团簇、纳米线(或管)、薄膜、晶体、准晶和无定性材料,以及表面体系和固体 l 计算材料的结构参数(键长、键角、晶格常数、原子位置等)和构型 l 计算材料的状态方程和力学性质(体弹性模量和弹性常数) l 计算材料的电子结构(能级、电荷密度分布、能带、电子态密度和ELF) l 计算材料的光学性质 l 计算材料的磁学性质 l 计算材料的晶格动力学性质(声子谱等) l 表面体系的模拟(重构、表面态和STM模拟) l 从头分子动力学模拟 l 计算材料的激发态(GW准粒子修正) 计算主要的四个参数文件:INCAR ,POSCAR,POTCAR ,KPOINTS,下面简要介绍,详细权威的请参照手册 INCAR文件: 该文件控制VASP进行何种性质的计算,并设置了计算方法中一些重要的参数,这些参数主要包括以下几类: 对所计算的体系进行注释:SYSTEM

●定义如何输入或构造初始的电荷密度和波函数:ISTART,ICHARG,INIWAV ●定义电子的优化 –平面波切断动能和缀加电荷时的切断值:ENCUT,ENAUG –电子部分优化的方法:ALGO,IALGO,LDIAG –电荷密度混合的方法:IMIX,AMIX,AMIN,BMIX,AMIX_MAG,BMIX_MAG,WC,INIMIX,MIXPRE,MAXMIX –自洽迭代步数和收敛标准:NELM,NELMIN,NELMDL,EDIFF ●定义离子或原子的优化 –原子位置优化的方法、移动的步长和步数:IBRION,NFREE,POTIM,NSW –分子动力学相关参数:SMASS,TEBEG,TEEND,POMASS,NBLOCK,KBLOCK,PSTRESS –离子弛豫收敛标准:EDIFFG ●定义态密度积分的方法和参数 –smearing方法和参数:ISMEAR,SIGMA –计算态密度时能量范围和点数:EMIN,EMAX,NEDOS –计算分波态密度的参数:RWIGS,LORBIT ●其它 –计算精度控制:PREC –磁性计算:ISPIN,MAGMOM,NUPDOWN –交换关联函数:GGA,VOSKOWN –计算ELF和总的局域势:LELF,LVTOT –结构优化参数:ISIF –等等。 主要参数说明如下: ?SYSTEM:该输入文件所要执行的任务的名字。取值:字符串,缺省值:SYSTEM ?NWRITE:输出内容详细程度。取值:0~4,缺省值:2

vasp并行总结 安装

一、安装vasp前的软件要求: ①C++编译器用intel的(l_cprof_p_11.1.07) ②Fortran编译器用intel的ifort11(l_cprof_p_11.1.072) ③l_mpi_p_3.2.011 ④MKL 有非商业版本可以免费下载,本来要用l_mkl_p_10.2.5.035的,但发现ifort11里/home/bjwang/intel/Compiler/11.1/072/mkl就有,这里免装,在.bashrc里把ifort11下mkl的路径包括进去。 附完整安装后的.bashrc source /home/bjwang/intel/Compiler/11.1/072/bin/intel64/ifortvars_intel64.sh---ifort source /home/bjwang/intel/Compiler/11.1/073/bin/intel64/iccvars_intel64.sh ---c++ export LD_LIBRARY_PATH=/home/bjwang/intel/Compiler/11.1/072/ifort/mkl/lib/em64t/:$LD_LIBRA RY_PATH ------ifort中包含的mkl source /home/bjwang/intel/impi/3.2.0.011/bin64/mpivars.sh --------l_mpi 请确认mpi、C++、Fortran编译器都已正确安装,并设定好相关的PATH路径和 LD_LIBRARY_PATH路径,具体参见新浪博客的其它相关内容。 示意如下:

VASP控制参数文件INCAR的简单介绍

限于能力,只对部分最基本的一些参数(>,没有这个标志的参数都是可以不出现的) 详细说明,在这里只是简单介绍这些参数的设置,详细的问题在后文具体示例中展开。 部分可能会干扰VASP运行的参数在这里被刻意隐去了,需要的同学还是请查看VASP自带的帮助文档原文。 参数列表如下: >SYSTEM name of System 任务的名字*** >NWRITE verbosity write-flag (how much is written) 输出内容详细程度0-3 缺省2 如果是做长时间动力学计算的话最好选0或1(首末步/每步核运动输出) 据说也可以结合shell的tail或grep命令手动输出 >ISTART startjob: restart选项0-3 缺省0/1 for 无/有前次计算的WAVECAR(波函数) 1 'restart with constant energy cut-off' 2 'restart with constant basis set' 3 'full restart including wave function and charge prediction' ICHARG charge: 1-file 2-atom 10-const Default:if ISTART=0 2 else 0 ISPIN spin polarized calculation (2-yes 1-no) default 2 MAGMOM initial mag moment / atom Default NIONS*1 INIWAV initial electr wf. : 0-lowe 1-rand Default 1 only used for start jobs (ISTART=0) IDIPOL calculate monopole/dipole and quadrupole corrections 1-3 只计算第一/二/三晶矢方向适于slab的计算 4 全部计算尤其适于就算孤立分子 >PREC precession: medium, high or low(VASP.4.5+ also: normal, accurate) Default: Medium VASP4.5+采用了优化的accurate来替代high,所以一般不推荐使用 high。不过high可以确保'绝对收敛',作为参考值有时也是必要的。 同样受推荐的是normal,作为日常计算选项,可惜的是说明文档提供的信息不足。 受PREC影响的参数有四类:ENCUT; NGX,NGY,NGZ; NGXF, NGYF, NGZF; ROPT 如果设置了PREC,这些参数就都不需要出现了 当然直接设置相应的参数也是同样效果的,这里不展开了,随后详释

VASP-INCAR参数设置

V A S P-I N C A R参数设置-CAL-FENGHAI-(2020YEAR-YICAI)_JINGBIAN

1. 结构优化 (Opt) SYSTEM = opt ISTART = 0 INIWAV = 1 ICHARG = 2 ISPIN = 2 LREAL = Auto ENCUT = 400 PREC = high NSW= 600 NELM = 60 IBRION = 2 ISIF = 2 POTIM = 0.1 ALGO= Fast LVDW = .TRUE. EDIFF = 1E-5 EDIFFG = 1E-4 or -0.05 # 体系需计算TS时,全部结构优化EDIFFG均设置为-0.05 ISMEAR = 0 SIGMA = 0.2 LCHARG = .FALSE. LWAVE = .FALSE.

2. 过渡态搜索 (TS): 计算时先进行低精度计算,再进行高精度计算 SYSTEM= TS ISTART = 0 INIWAV = 1 ICHARG = 2 ISPIN = 2 LREAL = Auto ENCUT = 400 PREC = high NSW = 600 NELMIN = 6 IBRION = 3 or 1 # 过渡态计算低精度为3,高精度为1 ISIF = 2 POTIM = 0.01 ALGO = Fast LVDW = .TRUE. EDIFF = 1E-5 EDIFFG = -1 or -0.05 # 过渡态计算低精度为-1,高精度为-0.05 ISMEAR = 0 SIGMA = 0.05 LCHARG= .FALSE. LWAVE= .FALSE. IMAGES=8 # TS专属设置 SPRING=-5 # TS专属设置 LCLIMB=.TRUE. # TS专属设置

vaspkit.014安装使用说明

vaspkit-0.14 修正版发布 安装下载指南 1、下载网址: https://www.sodocs.net/doc/e912665571.html,/downloads221/sourcecode/unix_linux/detail1038949.html (免费的) 2、利用xshell将软件包拷到系统中,建立文件夹vaspkit。 3、在gunzip vaspkit-0.14.tar.gz等到vaspkit-0.14.tar,再tar –xvf vaspkit-0.14.tar得到文件 夹vaspkit-0.14。 4、进入vaspkit-0.14文件夹,cat INSTALL,发现 Just do one thing: make Then you will get vaspkit program.然后退出,在vaspkit-0.14目录下输入make命令,得到vaspkit 可执行文件。 5、再将vaspkit文件分别拷到需要进行计算的文件下进行计算即可,输入./vasokit出现如下的运行页面。按照需要输入不同的数字即可。 相信很多人都为vasp计算结果的处理感到很头痛,网上虽然有很多处理vasp计算结果的软件或者脚本,不过要么只针对于某一个问题,要么使用不太友好。为此,我弄了一个软件包,把结构可视化,提取态密度,能带,电荷,自旋密度等数据等功能整合到一个软件包里,方便大家的使用。由于时间有限,我只花了两天时间做得,所以vaspkit工具包比较粗糙,而且有很多bug, 目前只实现最基本的功能。由于vaspkit 目前处于测试阶段,要不断的更新,所以不上传到论坛了。感兴趣的朋友可以到下载。希望大家能提出宝贵的意见以及需要增加的功能。同时,我们希望对vaspkit感兴趣的朋友一起加入并不断地增加功能,让vaspkit能成为一个很强大的数据处理工具。在此,特别感谢论坛版主csfn, foxhunter, 语过添情以及蓝等朋友的支持。 当前版本:0.14 . 1,修正若干bug; 2,增加输出xsf格式结构文件; 3,增加输出DOS积分; 4,增加提取指定原子指定轨道的能带数据. 主要功能有: 1,把POSCAR 或者CONTCAR转化为cif或xsf格式。 2,提取电荷,自旋密度,使用vesta软件可视化。 3,提取总,投影态密度数据; 4,提取总,投影能带数据。 期待增加功能: 1,增加计算差分电荷数据; 2,xsf格式电荷密度。 文件输出格式: 1,结构可视化文件输出文件为poscar.cif 或contcar.cif, 用MS, VESTA或其它可视化软件软件打开;2,总态密度输出文件为tdos.dat, 偏态密度文件为pdos.dat, 使用origin或gnuplot绘制。 3,能带数据输出文件名为band-x.dat . 同样使用origin或gnuplot绘制。 4,电荷和自旋密度输出文件名分别为chgden.dat 和spnden.dat, 用VESTA打开。注意同时把OUTCAR

vasp计算参数设置

软件主要功能: 采用周期性边界条件(或超原胞模型)处理原子、分子、团簇、纳米线(或管)、薄膜、晶体、准晶和无定性材料,以及表面体系和固体 l 计算材料的结构参数(键长、键角、晶格常数、原子位置等)和构型 l 计算材料的状态方程和力学性质(体弹性模量和弹性常数) l 计算材料的电子结构(能级、电荷密度分布、能带、电子态密度和ELF) l 计算材料的光学性质 l 计算材料的磁学性质 l 计算材料的晶格动力学性质(声子谱等) l 表面体系的模拟(重构、表面态和STM模拟) l 从头分子动力学模拟 l 计算材料的激发态(GW准粒子修正) 计算主要的四个参数文件:INCAR ,POSCAR,POTCAR ,KPOINTS,下面简要介绍,详细权威的请参照手册 INCAR文件: 该文件控制VASP进行何种性质的计算,并设置了计算方法中一些重要的参数,这些参数主要包括以下几类: l 对所计算的体系进行注释:SYSTEM l 定义如何输入或构造初始的电荷密度和波函数:ISTART,ICHARG,INIWA V l 定义电子的优化 –平面波切断动能和缀加电荷时的切断值:ENCUT,ENAUG –电子部分优化的方法:ALGO,IALGO,LDIAG –电荷密度混合的方法:IMIX,AMIX,AMIN,BMIX,AMIX_MAG,BMIX_MAG,WC,INIMIX,MIXPRE,MAXMIX –自洽迭代步数和收敛标准:NELM,NELMIN,NELMDL,EDIFF l 定义离子或原子的优化 –原子位置优化的方法、移动的步长和步数:IBRION,NFREE,POTIM,NSW –分子动力学相关参数:SMASS,TEBEG,TEEND,POMASS,NBLOCK,KBLOCK,PSTRESS –离子弛豫收敛标准:EDIFFG l 定义态密度积分的方法和参数 –smearing方法和参数:ISMEAR,SIGMA –计算态密度时能量范围和点数:EMIN,EMAX,NEDOS –计算分波态密度的参数:RWIGS,LORBIT l 其它 –计算精度控制:PREC –磁性计算:ISPIN,MAGMOM,NUPDOWN –交换关联函数:GGA,VOSKOWN –计算ELF和总的局域势:LELF,LVTOT –结构优化参数:ISIF –等等。 主要参数说明如下: ? SYSTEM:该输入文件所要执行的任务的名字。取值:字符串,缺省值:SYSTEM

vasp 安装心得

VASP5.2安装心得 2014-05-07 来源:小木虫作者: yysskk 花了五天时间终于学会怎么装VASP了,在此写下心得体会,供后人参考。个人觉得最难的一步就是makefile文件,网上流传着各种各样的版本,每个人都说自己编译成功了,却又各不相同,也说不清为什么,给新手极大的困扰。在此会详细介绍makefile的文件结构。其余大部分内容都是参考前人的,就不一一注明出处了。 一、系统、编译程序及准备工作 我用的是centos6.5+icc2011+ifort2011+openmpi1.6.5 1.1编译器安装 系统安装不说了,网上教程多得是。Icc和ifort可以申请免费非商业版本,icc和ifort都各自带了一个MKL,使用的时候别搞混了。装2011的时候会缺组件,用yum都可以免费下载。装编译器的时候会要求关闭selinux,按照给出的步骤关闭即可。之后会说系统不兼容,但是可以继续装,默认安装路径是在/opt下面。装完之后会有提示,把安装目录 /bin/ifortvars.sh 写到环境变量中,注意32/64位系统的参数不一样。C语言编译器建议用icc,毕竟是intel出品,针对自家cpu肯定有大量优化,效率上高于gcc是肯定的。这是装完之后的提示: For csh/tcsh: $ source install-dir/bin/compilervars.csh intel64 For bash: $ source install-dir/bin/compilervars.sh intel64 To invoke the installed compilers:

初学VASP中电子态密度计算设置参考

初学VASP中电子态密度计算基本设置参考主要分成三步:一、结构优化;二、静态自洽计算;三、非自洽计算以Al-FCC为例子 第一步结构优化 输入文件(INCAR, POTCAR, POSCAR, KPOINT) INCAR文件 System=Al ISTART=0 ISMEAR=1 SIGMA=0.2 ISPIN=2 GGA=91; VOSKOWN=1; EDIFF=0.1E-05; EDIFFG=-0.01 IBRION=2 NSW=50 ISIF=2 (OR 3) NPAR=10 POTCAR 文件直接在势库中拷贝 POSCAR文件 Al 4.05 1.0 0.0 0.0 0.0 1.0 0.0

0.0 0.0 1.0 4 Direct 0.0 0.0 0.0 0.5 0.5 0.0 0.5 0.0 0.5 0.0 0.5 0.5 KPOINT 文件 Automatic generation Mohkorst Pack 15 15 15 0.0 0.0 0.0 第二步静态自洽计算 INCAR:PREC = Medium,ISTART = 0,ICHARG = 2,ISMEAR = -5输入文件(INCAR, POTCAR, POSCAR, KPOINT) INCAR文件 System=Al ISTART=0 ISMEAR=1 SIGMA=0.2 ISPIN=2

GGA=91; VOSKOWN=1; EDIFF=0.1E-05; EDIFFG=-0.01 #IBRION=2 #NSW=50 #ISIF=2 (OR 3) NPAR=10 POTCAR 文件直接在势库中拷贝 POSCAR文件 Al 4.05 1.0 0.0 0.0 0.0 1.0 0.0 0.0 0.0 1.0 4 Selective Dynamic Direct 0.0 0.0 0.0 T T T 0.5 0.5 0.0 T T T 0.5 0.0 0.5 T T T 0.0 0.5 0.5 T T T KPOINT 文件 Automatic generation

VASP遇到小总结问题

VASP 计算的过程遇到的问题 01、第一原理计算的一些心得 (1)第一性原理其实是包括基于密度泛函的从头算和基于Hartree-Fock自洽计算的从头算,前者以电子密度作为基本变量(霍亨伯格-科洪定理),通过求解Kohn-Sham方程,迭代自洽得到体系的基态电子密度,然后求体系的基态性质;后者则通过自洽求解Hartree-Fock方程,获得体系的波函数,求基态性质; 评述:K-S方程的计算水平达到了H-F水平,同时还考虑了电子间的交换关联作用。 (2)关于DFT中密度泛函的Functional,其实是交换关联泛函 包括LDA,GGA,杂化泛函等等 一般LDA为局域密度近似,在空间某点用均匀电子气密度作为交换关联泛函的唯一变量,多数为参数化的CA-PZ方案; GGA为广义梯度近似,不仅将电子密度作为交换关联泛函的变量,也考虑了密度的梯度为变量,包括PBE,PW,RPBE等方案,BL YP泛函也属于GGA; 此外还有一些杂化泛函,B3L YP等。 (3)关于赝势 在处理计算体系中原子的电子态时,有两种方法,一种是考虑所有电子,叫做全电子法,比如WIEN2K中的FLAPW方法(线性缀加平面波);此外还有一种方法是只考虑价电子,而把芯电子和原子核构成离子实放在一起考虑,即赝势法,一般赝势法是选取一个截断半径,截断半径以内,波函数变化较平滑,和真实的不同,截断半径以外则和真实情况相同,而且赝势法得到的能量本征值和全电子法应该相同。 赝势包括模守恒和超软,模守恒较硬,一般需要较大的截断能,超软势则可以用较小的截断能即可。另外,模守恒势的散射特性和全电子相同,因此一般红外,拉曼等光谱的计算需要用模守恒势。 赝势的测试标准应是赝势与全电子法计算结果的匹配度,而不是赝势与实验结果的匹配度,因为和实验结果的匹配可能是偶然的。 (4)关于收敛测试 (a)Ecut,也就是截断能,一般情况下,总能相对于不同Ecut做计算,当Ecut增大时总能变化不明显了即可;然而,在需要考虑体系应力时,还需对应力进行收敛测试,而且应力相对于Ecut的收敛要比总能更为苛刻,也就是某个截断能下总能已经收敛了,但应力未必收敛。 (b)K-point,即K网格,一般金属需要较大的K网格,采用超晶胞时可以选用相对较小的K网格,但实际上还是要经过测试。 (5)关于磁性 一般何时考虑自旋呢?举例子,例如BaTiO3中,Ba、Ti和O分别为+2,+4和-2价,离子全部为各个轨道满壳层的结构,就不必考虑自旋了;对于BaMnO3中,由于Mn+3价时d 轨道还有电子,但未满,因此需考虑Mn的自旋,至于Ba和O则不必考虑。其实设定自旋就是给定一个原子磁矩的初始值,只在刚开始计算时作为初始值使用,具体的可参照磁性物理。 (6)关于几何优化 包括很多种了,比如晶格常数和原子位置同时优化,只优化原子位置,只优化晶格常数,还有晶格常数和原子位置分开优化等等。

VASP安装说明(简单易懂)

VASP安装说明(简单易懂) 1.安装linux 由于不熟悉linux的操作,所以很多问题无法解决,所以直接在windows下面安装了虚拟机然后安装linux操作系统。 (1) 虚拟机的安装 我所使用的是VMware,安装过程同普通的windows下软件的安装。下载地址:www.https://www.sodocs.net/doc/e912665571.html,或者www.https://www.sodocs.net/doc/e912665571.html, (2) VMware配置linux,过程如下: a)选择File菜单下的“New Virtual Machine”出现新虚拟机向导后单击“下一步” 选择“Typical”典型安装。 b)再单击“下一步”,在选择操作系统界面的“Guest Operation System”中选择 “Linux”,然后单击Version对应的下拉菜单选择具体的Linux版本,此处我选择的是“Red Hat Linux”。 c)单击“下一步”进入安装目录选择界面。该界面上面的文本框是系统的名字,保持默认值即 可,下面的文本框需要选择虚拟机操作系统的安装位置。 d)根据需要选择好后,单击“下一步”按钮,出现设置虚拟机内存大小的界面。Linux9.O对 内存的要求是:文本模式至少需要64MB;图形化模式至少需要

128MB,推荐使用192MB。我选 择的是192MB。 e)单击“下一步”按钮进入网络连接方式选择界面。VMware有四种网络设置方式,一般来说, Bridged方式使虚拟机就像网络内一台独立的计算机一样,最为方便好用(四种连网方式的区 别可参考VMware的有关资料)。此处我选择Brided方式。 f)单击“下一步”按钮进入虚拟磁盘的设置界面。这里有三种方式(Create a new virtual disk、Use an existing virtual disk、Use a physical disk)可供选 择、建议初学者选择“Create a new Virtual disk”,其含义是新建一个虚拟磁盘,该 虚拟磁盘只是主机下的一个独立文件。 g)在“下一步”中设置磁盘大小。在此、我们采用默认的4GB。 h)单击“下一步”进入文件存放路径选择界面。 在此界面可单击Browse按钮进行设置。此处我们使用默认值,单击“完成”按钮。 至此,完成一个虚拟机的配置。 (2) VMware下linux的安装 安装前需做好准备工作:购买一套Red Hat Linux 9.0的安装盘,共三张。或下载Red Hat Linux 9.0安装所需要的三个

VASP经典学习教程,有用

VASP 学习教程太原理工大学量子化学课题组 2012/5/25 太原

目录 第一章Linux命令 (1) 1.1 常用命令 (1) 1.1.1 浏览目录 (1) 1.1.2 浏览文件 (1) 1.1.3 目录操作 (1) 1.1.4 文件操作 (1) 1.1.5 系统信息 (1) 第二章SSH软件使用 (2) 2.1 软件界面 (2) 2.2 SSH transfer的应用 (3) 2.2.1 文件传输 (3) 2.2.2 简单应用 (3) 第三章V ASP的四个输入文件 (3) 3.1 INCAR (3) 3.2 KPOINTS (4) 3.3 POSCAR (4) 3.4 POTCAR (5) 第四章实例 (5) 4.1 模型的构建 (5) 4.2 V ASP计算 (8) 4.2.1 参数测试 (8) 4.2.2 晶胞优化(Cu) (13) 4.2.3 Cu(100)表面的能量 (2) 4.2.4 吸附分子CO、H、CHO的结构优化 (2) 4.2.5 CO吸附于Cu100表面H位 (4) 4.2.6 H吸附于Cu100表面H位 (5) 4.2.7 CHO吸附于Cu100表面B位 (6) 4.2.8 CO和H共吸附于Cu100表面 (7) 4.2.9 过渡态计算 (8)

第一章Linux命令 1.1 常用命令 1.1.1 浏览目录 cd: 进入某个目录。如:cd /home/songluzhi/vasp/CH4 cd .. 上一层目录;cd / 根目录; ls: 显示目录下的文件。 注:输入目录名时,可只输入前3个字母,按Tab键补全。1.1.2 浏览文件 cat:显示文件内容。如:cat INCAR 如果文件较大,可用:cat INCAR | more (可以按上下键查看) 合并文件:cat A B > C (A和B的内容合并,A在前,B在后) 1.1.3 目录操作 mkdir:建立目录;rmdir:删除目录。 如:mkdir T-CH3-Rh111 1.1.4 文件操作 rm:删除文件;vi:编辑文件;cp:拷贝文件 mv:移动文件;pwd:显示当前路径。 如:rm INCAR rm a* (删除以a开头的所有文件) rm -rf abc (强制删除文件abc) tar:解压缩文件。压缩文件??rar 1.1.5 系统信息 df:分区占用大小。如:df -h du:各级目录的大小。 top:运行的任务。 ps ax:查看详细任务。 kill:杀死任务。如:kill 12058 (杀死PID为12058的任务)注:PID为top命令的第一列数字。

VASP使用总结

VASP计算的理论及实践总结 一、赝势的选取 二、收敛测试 1、VASP测试截断能和K 点 2、MS测试 三、结构弛豫 四、VASP的使用流程(计算性质) 1、VASP的四个输入文件的设置 2、输出文件的查看及指令 3、计算单电能 (1) 测试截断能 (2) 测试K点 4、进行结构优化 5、计算弹性常数 6、一些常用指令

一、赝势的选取 VASP赝势库中分为:PP和PAW两种势,PP又分为SP(标准)和USPP(超软)。 交换关联函数分为:LDA(局域密度近似)和GGA(广义梯度近似)。GGA 又分为PW91和PBE。 在VASP中,其中pot ,pot-gga是属于超软势(使用较少)。Paw, paw-pbe ,和paw-gga是属于PAW。采用较多的是PAW-pbe 和PAW-gga。 此外vasp 中的赝势分为几种,包扩标准赝势(没有下标的)、还有硬(harder)赝势(_h)、软(softer)赝势(_s), 所谓的硬(难以赝化),就是指该元素原子的截断动能比较大,假想的势能与实际比较接近,计算得到的结果准确,但比较耗时,难以收敛。软(容易赝化),表示该元素原子的截断动能比较小,赝势模型比较粗糙,但相对简单,可以使计算很快收敛(比如VASP开发的超软赝势)。即硬的赝势精度高,但计算耗时。软的精度低,容易收敛,但节省计算时间。 另一种情况:如Gd_3,这是把f电子放入核内处理,对于Gd来说,f电子恰好半满。所以把f电子作为价电子处理的赝势还是蛮好的(类似还有Lu,全满)。(相对其他的4f元素来说,至于把f电子作为芯内处理,是以前对4f元素的通用做法。计算结果挺好) 常用的做法是:用两种赝势测试一下对自己所关心的问题的影响情况。在影响不大的情况下,选用不含4f电子的赝势(即后缀是3),一来减少计算量,二来避免DFT对4f电子的处理。 【1.赝势的选择: vasp的赝势文件放在目录~/vasp/potentials 下,可以看到该目录又包含五个子目录pot pot_GGA potpaw potpaw_GGA potpaw_PBE ,其中每一个子目录对应一种赝势形式。

初学VASP最重要的INCAR参数

初学VASP(六) 最重要的INCAR参数 初学VASP(六) 最重要的INCAR参数 INCAR是决定how to do 的文件 限于能力,只对部分最基本的一些参数(>,没有这个标志的参数都是可以不出现的) 详细说明,在这里只是简单介绍这些参数的设置,详细的问题在后文具体示例中展开。 部分可能会干扰VASP运行的参数在这里被刻意隐去了,需要的同学还是请查看VASP自带 的帮助文档原文。 参数列表如下: >SYSTEM name of System 任务的名字 *** >NWRITE verbosity write-flag (how much is written) 输出内容详细程度 0-3 缺省2 如果是做长时间动力学计算的话最好选0或1(首末步/每步核运动输出) 据说也可以结合shell的tail或grep命令手动输出 >ISTART startjob: restart选项 0-3 缺省0/1 for 无/有前次计算的WAVECAR(波函数) 1 'restart with constant energy cut-off' 2 'restart with constant basis set' 3 'full restart including wave function and charge prediction' ICHARG charge: 1-file 2-atom 10-const Default:if ISTART=0 2 else 0 ISPIN spin polarized calculation (2-yes 1-no) default 2 MAGMOM initial mag moment / atom Default NIONS*1 INIWAV initial electr wf. : 0-lowe 1-rand Default 1 only used for start jobs (ISTART=0) IDIPOL calculate monopole/dipole and quadrupole corrections 1-3 只计算第一/二/三晶矢方向适于slab的计算 4 全部计算尤其适于就算孤立分子 >PREC precession: medium, high or low(VASP.4.5+ also: normal, accurate)

CTI平台软件安装说明

CTI平台软件安装说明

目录 1 安装SFU (5) 1.1安装文件释放 (5) 1.2 安装过程 (6) 1.2.1 安装执行 (6) 1.2.2 填写用户信息 (6) 1.2.3 软件协议签订 (7) 1.2.4 安装方式选择 (7) 1.2.5 组件选择 (8) 1.2.6 安全设置 (8) 1.2.7 映射方式选择 (9) 1.2.8 映射文件选择 (9) 1.2.9 安装目录选择 (10) 1.2.10 安装过程 (10) 1.2.11 安装完成 (11) 1.3 配置过程 (11) 1.3.1 检查启动服务 (11) 1.3.2 Configuration配置 (11) 1.3.3 Maps配置 (12) 1.3.4 保存配置 (14) 1.4 安装完成 (15) 2配置SFU (15) 2.1设置SFU服务器共享目录 (15) 2.2设置windows访问权限 (16) 2.3创建座席录音子目录 (18) 2.3.1 检查子目录的windows权限 (19) 2.4设置rt-sc的nfs共享路径 (19) 2.5检查sfu服务器nfs共享是否生效 (20) 2.5.1 通过交换机来检查 (20) 2.5.2 通过vasp来检查 (21) 3 SFU服务器上安装JDK、TOMCAT (21) 3.1 安装JDK (21) 3.2 安装TOMCA T (21) 4安装ccmp (23) 4.1 SFU服务器上安装ccmp (23) 4.2 CCMP配置 (23) 4.3 配置数据库IP地址 (23) 4.4 配置录音路径 (23) 5 安装vasptools (24) 5.1 SFU服务器上安装vasptools (24) 5.2 V ASPTOOLS配置 (24) 5.3 系统管理配置 (25) 6安装PCMPlayer (25)

VASP-INCAR参数设置

1. 结构优化(Opt) SYSTEM = opt ISTART = 0 INIWAV = 1 ICHARG = 2 ISPIN = 2 LREAL = Auto ENCUT = 400 PREC = high NSW= 600 NELM = 60 IBRION = 2 ISIF = 2 POTIM = 0.1 ALGO= Fast LVDW = .TRUE. EDIFF = 1E-5 EDIFFG = 1E-4 or -0.05 # 体系需计算TS时,全部结构优化EDIFFG均设置为-0.05 ISMEAR = 0 SIGMA = 0.2 LCHARG = .FALSE. LWAVE = .FALSE.

2. 过渡态搜索(TS): 计算时先进行低精度计算,再进行高精度计算 SYSTEM= TS ISTART = 0 INIWAV = 1 ICHARG = 2 ISPIN = 2 LREAL = Auto ENCUT = 400 PREC = high NSW = 600 NELMIN = 6 IBRION = 3 or 1 # 过渡态计算低精度为3,高精度为1 ISIF = 2 POTIM = 0.01 ALGO = Fast LVDW = .TRUE. EDIFF = 1E-5 EDIFFG = -1 or -0.05 # 过渡态计算低精度为-1,高精度为-0.05 ISMEAR = 0 SIGMA = 0.05 LCHARG= .FALSE. LWAVE= .FALSE. IMAGES=8 # TS专属设置 SPRING=-5 # TS专属设置 LCLIMB=.TRUE. # TS专属设置

vasp编译

(1)准备 1. sudo apt-get install g++ sudo apt-get install build-essential sudo apt-get install gcc-multilib sudo apt-get install libstdc++5 sudo apt-get install openjdk-6-jre-headless sudo apt-get install ia32-libs sudo apt-get install lib32stdc++6 sudo apt-get install libc6-dev-i386 sudo apt-get install g++-multilib sudo apt-get install gfrotran 1 2 3 tar -2. xzvf xxx #xxx 指压缩包的名称 cd 解压后目录 ./install.sh ##########安装过程很简单,安装提示就行######## 安装完成后,修改环境变量 1 vi ~/.bashrc 在末尾添加 source /opt/intel/bin/compilervars.sh intel64 (64位版本) source /opt/intel/bin/compilervars.sh ia32 (32位版本) 1 source ~/.bashrc 4. 编译mpi 首先下载mpich2的安装文件 1 2 3 4 5 6 tar mpich2压缩包 cd 文件夹 ./configure --prefix=安装目录 f77=ifort f90=ifort fc=ifort make make install vi ~/.bashrc 笔者安装目录在/opt/mpich2,也可以放在其他文件夹中,注意执行make install 用户需要有写入权限 在文件最后写入如下代码 1 2 3 4 #mpi2 start PATH=/usr/local/mpich2-1.2.1p1/bin:$PATH export PATH #mpi2 end 再source 一下 1 source ~/.bashrc 5.编译fftw

vasp参数介绍

内容 描述体系:SYSTEM 设置如何输入或构造初始的电荷密度和波函数:ISTART, ICHARG, INIWA V 设置电子的优化: 平面波切断动能和缀加电荷时的切断值:ENCUT, ENAUG 电子部分优化的方法:ALGO, IALGO, LDIAG 电荷密度混合的方法:IMIX, AMIX, AMIN, BMIX, AMIX_MAG, BMIX_MAG, WC, INIMIX, MIXPRE, MAXMIX 自洽迭代步数和收敛标准:NELM, NELMIN, NELMDL, EDIFF 设置原子的驰豫: 原子如何移动以及步长和步数:IBRION, NFREE, POTIM, NSW 分子动力学相关参数:SMASS, TEBEG, TEEND, POMASS,NBLOCK, KBLOCK, PSTRESS 原子驰豫收敛标准:EDIFFG 定义态密度积分的方法和参数: smearing方法和参数:ISMEAR, SIGMA 计算态密度时能量范围和点数:EMIN, EMAX, NEDOS 计算分波态密度的参数:RWIGS, LORBIT 其他: 计算精度控制:PREC 磁性计算:ISPIN, MAGMOM, NUPDOWN 交换关联函数:GGA, VOSKOWN 计算ELF和总的局域势:LELF, LVTOT 结构优化参数:ISIF 初始化 SYSTEM: 注释所计算的体系,以示说明。 NWRITE: 默认值为2,可赋予值为0| 1 | 2 |3 |4。决定OUTCAR中输出内容的详细程度ISTART: 如果计算目录中有W A VECAR文件,则默认值为1,否则为0。可赋予值为0| 1 | 2 | 3 。决定是否读入W A VECAR: 0: 开始新的计算,按INIWA V初始化波函数 1: 接着计算,通常用在测试ENCUT的收敛性以及计算结合能曲线(也就是体积和总能的关系) 2: 接着计算,通常用在希望保持基矢不变的计算中 3: 接着计算,读入上一次计算得到的电荷密度和波函数,不推荐用 ICHARG: 如果ISTART=0,则默认值为2,否则为0。可赋予值为0 | 1 | 2。决定了如何构造初始的电荷密度 0: 从初始的波函数构造 1:读入CHGCAR读入,并同原子密度进行线性插值 2: 构造原子密度 11:读入自洽的CHGCAR,并进行能带计算或态密度的非自洽计算 12: 非自洽的原子密度计算 电子自洽收敛 INIWAV: 默认值为1,可赋予值为0 | 1。只在开始新的计算(也就是ISTART=0)中有效。决定了如何初始化波函数 0: 采用”jellium 波函数”

VASP5.3.5并行版详细安装教程

vasp5.3.5安装教程 制作者:甄江苏(湘潭大学材料科学与工程学院)1.安装Intel编译器: 。。。。。。。。 然后添加环境变量,在/etc/profile文件中最下面插入: source /opt/intel/composer_xe_2015.2.164/bin/ifortvars.sh intel64 source /opt/intel/composer_xe_2015.2.164/mkl/bin/mklvars.sh intel64 然后在终端执行: [root]# source /etc/profile 这里提供的安装包会把其他的Intel的组件包括并行编译器也安装上,但是后面不用它的mpirun这个命令,而是用mpich的mpirun命令,所以进入/opt/intel/composer_xe_2015.2.164/文件夹,把文件夹mpirt改名,让它的路径失效,比如改成mpirt11111。 此外,安装完后后,最好把licence文件放到安装目录里,Intel编译器指定的license 存放目录之一是/opt/intel/composer_xe_2015.2.164/license/这个文件夹,这个文件夹要自己新建,然后把license文件拷贝进去即可,不然的话,对于共享此软件的其他节点来说,会找不到license。 然后进入/opt/intel/composer_xe_2015.2.164/mkl/interfaces/fftw3xf/文件夹,执行:[root]# make libintel64 将会生成一个名为libfftw3xf_intel.a的静态库文件,用于后面vasp的链接。 2.安装openmpi:(openmpi总是出问题,无法多节点并行,原因不明,所以建 议用后面的mpich) 进入openmpi解压目录,执行: [root]# ./configure --prefix=/opt/openmpi CC=icc CXX=icpc F77=ifort FC=ifort [root]# make all install 然后添加环境变量,在/etc/profile中添加: export PATH=$PATH:/opt/openmpi/bin export LD_LIBRARY_PATH=$LD_LIBRARY_PATH:/opt/openmpi/lib export MANPATH=$MANPATH:/opt/openmpi/share/man

相关主题