小木虫 --- 600万学术达人喜爱的学术科研平台
热门搜索：
&&【求助】VASP计算的KPOINTS 设置，求助
【求助】VASP计算的KPOINTS 设置，求助
1.请问关于Kpoints 设置，如果需要的只是 lattice relaxation, 找到稳定，低能的结构， 是不是可以选择auto-matic? auto-matic K点的设置是不是按照晶格常数a.b c的比来设置？
2. 如果需要算band structure 或者DOS 时候， 才需要手动选取高对称K点，right?
3.lattice relaxation步骤如何？
贴个INCAR过来，请行家帮忙看看。
SYSTEM=Mn 200GPa
PREC=Normal
EDIFF=1e-4
ISMEAR =1; SIGMA=0.1
POTIM=0.05
PSTRESS=2000
LWAVE = FALSE
LCHARG = FALSE
#SYMPREC=1e-9
这个INCAR，relaxation之后， 得到CONTCAR， 然后cp CONTCAR POSCAR,
继续relaxation, 直到一个循环就可以达到计算精度，就算relaxation完全，是这样吗？
3. “结构优化时每个原子受力”， 从哪里可以看到？
relaxation时Kpoints用auto就可以了，设置标准可以参考 http://emuch.net/bbs/viewthread.php?tid=2766048&fpage=1
band structure 的时候要自己设置高对称点， 在网上搜下Si能带的计算。dos计算是加密Kpoints就可以了。这个说得很详细http://emuch.net/bbs/viewthread.php?tid=2111169
计算晶格常数，不改变形状用ISIF=7就可以了，不过一般都是设置好多常数再做拟合比较准确。按你上面说的做relaxation后得到的CONTCAR就是最有结构了，计算其他性质前再做一次静态计算就可以了（改变ISIF=2）。
受力从OUTCAR里面找force就行。
有一句没看懂，“设置许多参数然后拟合”是指哪些参数？如何拟合？这样做的好处是什么？不这样做，直接像我的INCAR那样设置结果会有什么缺陷？
还有，我从OUTCAR里面发现了好几项force，都代表什么呢？整个outcar文件，除了能量等有限的几项外，好多项我不懂都是什么意思，有没有一个详细的说明关于outcar文件中各项的意义？
我从候老师的资料中知道，优化结构时，参数的选择由enthalpy 的精确度来确定。有人说，还要同时考虑“每个原子受力”，我不知道outcar中的哪个force对应的是原子受力。请指教。谢谢！
指教谈不上！我就懂一点点
-------------------------------- Iteration& & 2(& &1)&&---------------------
这个表示第二个离子步，第一次迭代（电子步），一个离子步中离子的坐标不改变，不断改变电荷密度，最终收敛，并计算每个离子的受力，若受力小于设定值，则停止得出最优几何构型。若大于则根据受力来改变离子（原子核）坐标，在开始下一次离子步，不断循环。
计算晶格常数时，多个参数，就是设置不同的晶格常数，分别计算能量，最后用体积与能量的方程拟合，我一般是用的ISIF=7做的优化，这样好像不好，你看看这个吧。
http://emuch.net/bbs/viewthread.php?tid=2267882
期待高手给出更详细的解答，
感谢感谢， 现在好像懂了，实际操练起来再回来问。
这位大哥能不能再给点指教啊&&我想做无定型结构中离子受力，不知道该从何做起，
学术必备与600万学术达人在线互动！
扫描下载送金币小木虫 --- 600万学术达人喜爱的学术科研平台
热门搜索：
&&【求助】VASP计算的KPOINTS 设置，求助
【求助】VASP计算的KPOINTS 设置，求助
1.请问关于Kpoints 设置，如果需要的只是 lattice relaxation, 找到稳定，低能的结构， 是不是可以选择auto-matic? auto-matic K点的设置是不是按照晶格常数a.b c的比来设置？
2. 如果需要算band structure 或者DOS 时候， 才需要手动选取高对称K点，right?
3.lattice relaxation步骤如何？
贴个INCAR过来，请行家帮忙看看。
SYSTEM=Mn 200GPa
PREC=Normal
EDIFF=1e-4
ISMEAR =1; SIGMA=0.1
POTIM=0.05
PSTRESS=2000
LWAVE = FALSE
LCHARG = FALSE
#SYMPREC=1e-9
这个INCAR，relaxation之后， 得到CONTCAR， 然后cp CONTCAR POSCAR,
继续relaxation, 直到一个循环就可以达到计算精度，就算relaxation完全，是这样吗？
3. “结构优化时每个原子受力”， 从哪里可以看到？
relaxation时Kpoints用auto就可以了，设置标准可以参考 http://emuch.net/bbs/viewthread.php?tid=2766048&fpage=1
band structure 的时候要自己设置高对称点， 在网上搜下Si能带的计算。dos计算是加密Kpoints就可以了。这个说得很详细http://emuch.net/bbs/viewthread.php?tid=2111169
计算晶格常数，不改变形状用ISIF=7就可以了，不过一般都是设置好多常数再做拟合比较准确。按你上面说的做relaxation后得到的CONTCAR就是最有结构了，计算其他性质前再做一次静态计算就可以了（改变ISIF=2）。
受力从OUTCAR里面找force就行。
有一句没看懂，“设置许多参数然后拟合”是指哪些参数？如何拟合？这样做的好处是什么？不这样做，直接像我的INCAR那样设置结果会有什么缺陷？
还有，我从OUTCAR里面发现了好几项force，都代表什么呢？整个outcar文件，除了能量等有限的几项外，好多项我不懂都是什么意思，有没有一个详细的说明关于outcar文件中各项的意义？
我从候老师的资料中知道，优化结构时，参数的选择由enthalpy 的精确度来确定。有人说，还要同时考虑“每个原子受力”，我不知道outcar中的哪个force对应的是原子受力。请指教。谢谢！
指教谈不上！我就懂一点点
-------------------------------- Iteration& & 2(& &1)&&---------------------
这个表示第二个离子步，第一次迭代（电子步），一个离子步中离子的坐标不改变，不断改变电荷密度，最终收敛，并计算每个离子的受力，若受力小于设定值，则停止得出最优几何构型。若大于则根据受力来改变离子（原子核）坐标，在开始下一次离子步，不断循环。
计算晶格常数时，多个参数，就是设置不同的晶格常数，分别计算能量，最后用体积与能量的方程拟合，我一般是用的ISIF=7做的优化，这样好像不好，你看看这个吧。
http://emuch.net/bbs/viewthread.php?tid=2267882
期待高手给出更详细的解答，
感谢感谢， 现在好像懂了，实际操练起来再回来问。
这位大哥能不能再给点指教啊&&我想做无定型结构中离子受力，不知道该从何做起，
学术必备与600万学术达人在线互动！
扫描下载送金币拒绝访问 |
| 百度云加速
请打开cookies.
此网站 () 的管理员禁止了您的访问。原因是您的访问包含了非浏览器特征(3b9b3-ua98).
重新安装浏览器，或使用别的浏览器文档分类：
下载后只包含 1 个 DOC 格式的文档，没有任何的图纸或源代码，
下载前请先预览，预览内容跟原文是一样的，在线预览图片经过高度压缩，下载原文更清晰。
您的浏览器不支持进度条
下载文档到电脑，查找使用更方便
还剩?页未读，继续阅读
播放器加载中，请稍候...
该用户其他文档
下载所得到的文件列表VASP计算AgGaS2能带及态密度及光学性质.doc
文档介绍：
VASP计算AgGaS能带及态密度及光学性质第一步:构型优化.准备四个输入文件POSCARINCARPOTCARKPOINTSPOSCAR:从ms中导入AgGaS结构,选择CASTEP,file,save,并保存成原包。这样,得到一隐藏文件.cell,将它用编辑器打开,从中的到vasp所需的POSCAR信息,修改得到POSCAR。AgGaSbulk.-....-....-.Direct...................POTCAR:用PBE-GGA的赝势,提取,AgGaO的赝势合并成一个赝势。(一般VASP有自带).下面是INCARSYSTEM=optimizationofAgGaSLPLANE=.TRUE.NPAR=ElecronicminimisationISTART=LREAL=.FALSE.PREC=LowEDIFF=e-EDIFFG=-.IALGO=NELMIN=ISYM=GGA=PBEISPIN=NBANDS=OUTPUTCONTROLLCHARG=.TRUE.LWAVE=.TRUE.LVTOT=.FALSE.IONICRELAXATIONNBLOCK=NSW=IBRION=-DOSRELATED(disabled)ISMEAR=(tetrahedron/gaussian/m-p)SIGMA=..再下面KPOINTSAM至此,四文件已准备好,进行计算,mpiexec-np~/bin/vasp..-mk-mp-pgi&/dev/null&vasp.out&(运行命令)(问题:计算前对某些参数的测试,比如截断能,晶格参数等,标准时什么?vasp运行后,不知道结果是否满足要求呢?)第二部:将构型优化后产生的CONTCAR文件拷贝为POSCAR文件mpiexec-np~/bin/vasp..-mk-mp-pgi&/dev/null&vasp.out&(运行命令)第三部:能带的计算.建立新的计算目录mkdirband创建目录cdband进入目录cp../INCAR.复制cp../PO*.cp../CHG*.对INCAR文件作如下修改:ICHARG=NBANDS=与上面构型优化时一致建立新的KPOINTS手动定义K点AgGaSRec(关键字要对齐第一行)..............K点坐标获得途径布里渊区各特殊点的坐标可利用CASTEP获得然后编辑viinplin特殊点lin特殊点之间的取点数lin第一个K点的坐标lin第二个K点的坐标....(同上..)....保存运行vasp_kpoints&inp就可以得到坐标手动复制到KPOINTS运行VASPmpiexec-np~/bin/vasp..-mk-mp-pgi&/dev/null&vasp.out&(运行命令)在运行VaspbandDrawthebandstructureofVASPTwofiles,EIGENVALandOUTCARareneeded----------------------------------------------------------------Spin-polarizedcalculation(T/F)?(是否自旋极化计算)fNumberoflinesalongtheBZ:(布里渊区特殊线数目)Kpointsequenceusedtoseparatelines:(各线段始终点序号)SettheFermileveltozero(T/F)?(是否费米设置为零点)t可以得到相应的band.dat用origin绘图然后标注各个特殊点Fermienergy=.vexp=.ev(可能是参数设计太粗糙…)不过放大还算清楚..(用MS作图K网格相同)第四部:DOS的计算及其绘图其POSCARPOTCARKPOINTS可以相同于第一步优化设置INCARSYSTEM=optimizationofAgGaSNPAR=LPLANE=.TRUE.ElecronicminimisationISTART=(从WAVECAR读入波函数)LREAL=.FALSE.PREC=MediumEDIFF=e-EDIFFG=-.IALGO=NELMIN=ISYM=ISIF=ISPIN=ISMEAR=SIGMA=.NBANDS=LORBIT=(投影到s.p.d)NEDOS=(DOS的点数)ICHARG=(从CHG读入电荷密度)EMIN=-EMAX=OUTPUTCONTROLLCHARG=.TRUE.LWAVE=.TRUE.LVTOT=.FALSE.IONICRELAXATIONNBLOCK=NSW=IBRION=-POTIM=.运行vasp命令同上然后编写DOSF(是否自旋极化计算)(轨道数)(要计算的DOS)-(为负数时,表示所考察的为原子的DOS)(原子编号)-(为正正数时,表示为原子轨道的DOS)(下面都同上)-保存运行vaspdos(如果没有设定路径的话。要自己指定路径)可以得到fort.文件运行cdos不过需要加一条DOS数就可以得到DOSDOSDOS......用origin分别绘图…..(可能是DOS编号出现错误了的原因导致Ga(p)和S(p)没有画出..)第五部:光学性质的计算(一)线性的计算.也是进行构型优化(同能带计算时结构优化)..只是计算中INCAR文件需包含关键词:ISIF=所有计算均需采用PAW型赝势.构型优化后将CONTCAR拷贝成POSCAR在INCAR里面NPAR=(计算光学性质时,必须按能带依次处理)LOPTICS=.TRUE.(计算光学性质)NSW=IBRION=-(无需进行构型优化)计算结束后(正常),会得到OPTIC(用于计算线性光学性质)和momentum_matrix(用于计算NLO即非线性光学性质).编辑OPTCRT文件常用关键词及其含义ISYMM=OMMIN=OMMAX=NEDOS=NBCON=LJDOS=.TRUE.LDOS=.TRUE.LKRAMERS=.TURE.GAMMA=.LSEARCH=.TRUE.EMINSEARCH=EMAXSEARCH=SCISSOR=..创建一个新目录,将OPTICOPTCTR,KPOINTS,POSCAR复制到其下命令cp../OPTIC.OPTCTR,KPOINTS,POSCAR.修改KPOINTS创建个临时目录复制vasp计算所需要的四个输入文件.运行单机版vasp。(由于机子没有装单机版,所以这一步COPY到其他的机子上做的)单机运行命令~/bin/vasp..-mk-sp-pgi&/dev/null&vasp.out&将其中的IBZKPT复制为KPOINTS运行vasp命令:mpiexec–np~/bin/vasp..-mk-mp-pgi&/dev/null&vasp.out&正常运行将产生个文件EPS----存放介电函数数据JDOS----存放JDOS数据在运行vasp_lo用origin软件根据EPS数据绘图在此基础上才可以得到其他线性光学性质:()折射率(refractiveindex)()吸收系数(adso1
内容来自淘豆网转载请标明出处.