如果你说的是要运行Vasp的运行命令的话,直接键入VASP就可以运行了. 但是vasp必须先在linux下编译,才能执行.不知道你编译了没有.
如果你问是的计算弹性常数的步骤的话,量子化学网上,候拄峰写了一个关于用VASP计算六角结构弹性常数的步骤,你可以参考一下
如何建立所需文件夹来计算NaCl结构的弹性常数
你好!我是一名新手。但是我现在急需用vasp解决一个问题。我用castep计算得到的弹性常数和一篇引用文献中计算得到的结果正好相反(我计算为负值,则此结构在力学上是不稳定的;但文献中vasp计算得到的是正值,也就是此结构在力学上是稳定的。)。所以我想对NaCl结构进行几何优化,并且计算它的弹性常数, 来验证一下,并找出原因。我想请问一下在Linux系统下的计算步骤和命令。我们组刚买来vasp,可是我急需解决这个问题。所以恳请您的指导。
回复:tpp001
如果你说的是要运行Vasp的运行命令的话,直接键入VASP就可以运行了. 但是vasp必须先在linux下编译,才能执行.不知道你编译了没有.
如果你问是的计算弹性常数的步骤的话,量子化学网上,候拄峰写了一个关于用VASP计算六角结构弹性常数的步骤,你可以参考一下
楼上各位讨论同位素的兄弟好像弄错了一个问题: 同位素不是质子而是中子数不同吧? 如果是中子的话用赝势就不好处理了吧,中子应该与核外电子没啥作用,估计主要的影响就是原子核质量的变化
请问:算态密度时,出现warning:stress and forces are not correct??
计算TiC的态密度时,出现warning:stress and forces are not correct, 请问楼主,这是怎么回事?如何修复呢?非常感谢!
出现这种情况下,可能是由于你的K点的选择发生变化,或者是截断能的变化,引起体系内部受力发生改变,尽管官方网站上提供了一些可改性的参考,但是他们自己也认为这一点很难处理,不过值得庆幸的是,它对计算的结果影响不大.
用castep计算CeBr3和Ce3Al的结构优化,我用一般的手段不能收敛,用All Bands/EDFT也不收敛,加上自旋极化后还是无法收敛。提示:
*Warning* max. SCF cycles performed but system has not reached the groundstate. Warning: electronic minimisation did not converge during finite basis set correction.
不知道这些难以弛豫收敛的体系怎么对付?
如果是要使用一个大的晶胞,那么有没有必要去对它的原胞做一个结构优化来达到稳定,最
后用这个优化的晶格常数进行扩胞呢?还有你说DOS和设置的RWIGS有关,那么这个RWIGS应该怎样设才合理?
求助,bulk晶体体系的电导率如何计算?
有人推荐ATK,但是好像只能计算two-prob system
回复cuihuang,你的问题计算电导率,其实按照电导率的公式来看,它是与载流子的迁移率以及扩散系数是相关,但如何计算载流子的迁移率和扩散系数,目前还没发现很好的用第一性原理去计算这些东西的文章以及方法.
但是有篇文章可以参考一下:Computer Physics Communications 175 (2006) 67–71
比如对二维正方格子, 可以取象附件里那样的A和B两种基矢,这两种取法对计算能带有什么样的影响吗? 主要是两者的高对称点是等价的吗? 如果不等价在计算的时候在KPOINTS文件里面是不是要加上偏移量?
回复dwma:如果你的结构优化只是保持原子位置不变,针对体积,或者形变的做优化的话,可以通过选择ISIF的设置来实现.
或者这也是可行的,就是让你所说的两个原子静止,其他原子驰豫,那可以通过设置你的POSCAR跟相应的优化(驰豫)开关(ISIF,NSW,IBRION..)可以实现.
但是如果你的问题是保持一个体系中的两个原子的间距不变的情况下,仍然要求原子在晶体势场,或者力的的影响下运动的话,VASP的现在的模块是无法实现的,当然可能可以通过修改程序的方法达到你的目的
pdos的计算问题,新手求助
我计算的是Mg2Ni储氢材料
在DOSCAR中没有看到有s,p,d轨道的分别显示,但是在说明书中说好像能分别显示。。 PREC=Accurate ENCUT=250 EDIFF=1E-5 LORBIT=1 NSW=0 ISTART=1 ICHARG=11
RWIGS= 1.524 1.286 ISMEAR=-5
我哪里做错了或者设置的不对了将LORBIT=1改成LORBIT=11试试 谢谢
回复mumu03101::
你的第一个问题可以回答你,但是第二个问题我不太清楚你的目的
VASP计算过程中,可能由于断电或者其他原因导致VASP未计算完而中途推出.
要续算,我们首先得搞清楚它的原理,它在计算过程中是使用迭代的算法直至达到收敛精度,在每一个迭代的过程中,它都会产生一个新的电荷密度,和构造一个新的波函数,以备下一步的计算,因此当断电后,它也会产生相应的电荷密度(CHGCAR)和响应的波函数(WAVECAR),即便它没收敛,但它相比初试的电荷密度和波函数更接近收敛值,因此你要续算,只要设置好ICHARG和ISTART的开关,你应该可以达到你的目的,当然具体选择多少可以根据你的实际情况来
周期性扩展
怎样将得到的OUTCAR或者CONTCAR里面的结构进行周期行扩展,我因为计算的是表面,只计算了一个比较小的超胞,想验证这个超胞周期行扩展后是否还是正确的!
回复yhjiangshiyu:你所说的输出态密度的顺序是否跟POSCAR一致,这完全取决于你处理DOS的程序.
你的第二个问题:各个原子的DOS的和是不跟总的一致,因为各个原子的DOS是根据RWIGS的大小求出来,那么RWIGS的球并不能填充到整个空间,所以要偏小,但是这个差距应该不是很大.
当然引起不相等原因,还有很多,比如你的程序有问题.
怎样将得到的OUTCAR或者CONTCAR里面的结构进行周期行扩展,我因为计算的是表面,只计算了一个比较小的超胞,想验证这个超胞周期行扩展后是否还是正确的!
我用MS进行过手动平移,老板说手动平移不能验证周期性,请问有没有程序或者其他软件能让我得到的OUTCAR或者CONTCAR里面的原子坐标进行一个周期性的扩展!
你用ms进行手动平移是可行的 建议:
首先将contcar的格式转换成ms的格式,然后平移,这个是完全可行的。
具体ms的格式,你把一个ms自带的结构输出为一个.car格式的文件,读一读,参考下就会明白是什么格式,然后contcar转换成car格式即可。
我是无机化学专业,现在是研一,以前学的是土木工程,想来做量子化学计算,晶体的能带结构和态密度,具体是想看看一些黑色氧化物的吸光性质怎么样,我导师不是做计算的,她也不懂,所以就想知道这个事可以从什么地方开始,需要看什么书,可以用什么软件去计算
你用ms进行手动平移是可行的 建议:
首先将contcar的格式转换成ms的格式,然后平移,这个是完全可行的。
具体ms的格式,你把一个ms自带的结构输出为一个.car格式的文件,读一读,参考下就会明白是什么格式,然后contcar转换成car格式即可。
实例14
CASTEP,VASP,Quantum-ESPRESSO区别
CASTEP,VASP都是比较常用的平面波,赝势,周期性DFT计算的常用软件,我在晶体晶面计算的时候常看到,还经常被冠以First-principle calculation,一直很困惑这两个软件似乎都可以优化结构,都可以计算体系能量,但是区别在哪?价格上差那么多,难道就仅仅是因为一个图形界面好?今天又在JACS上看到了计算晶体,晶面结构还可以用另外一个软件,Quantum-ESPRESSO,查了下也是大名鼎鼎,里面的pwscf包也是平面波,赝势,周期性DFT,而且是免费的,实在弄不清楚这三个宣传都是平面波,赝势,周期性DFT的软件在优化结构,计算单点能方面有什么区别?为什么从CASTEP的10多万到VASP的4万再到Quantum-ESPRESSO开源,免费这么大的区别
CASTEP 贵,图形界面,容易上手
VASP 性价比最高,赝势最佳,算的很快
Quantum-ESPRESSO 免费,功能多,赝势上缺乏PAW。 个人推荐PW,算的很好,上手其实比vasp要快。
给Quantum-ESPRESSO 做个广告吧: Quantum-ESPRESSO 功能强大,含有 PW DFT计算 CPMD计算
PWCOND输运计算 PH 声子计算
后处理WANT wannier函数 SAX 光学性质
不过缺点就是赝势多为US PP
NC PP的赝势都很少,不用说PAW PP。
这点要做好心理准备,不然发现要算的体系你没有赝势就
我做的体系连JACS都发了 一般不会发生没有赝势的情况吧
感觉还是US PP更有效一点儿。ABINIT就是没有US PP,搞的我用NC PP算大晶胞那叫一个痛苦。
NC 赝势可以用ATOM之类的自己作,但是很费时间,而且还不一定好使,建议多交流现成赝势不一定,美国能源部实验室有个组,Pt赝势到现在还没有公开,只有他们组的结果算得和实验符合,一发就是牛文。 开源的还是小心赝势的问题。
PW注意测试赝势,连C的有几个PP,都是值得怀疑。
实例15
VASP算稀土永磁材料的磁学性能用哪种算法和赝势比较好?
用VASP计算稀土永磁材料(比如Sm-Co)的磁学性能用哪种算法和赝势比较好啊?
LDA GGA LSDA+U? PBE PW91?
不知哪位大侠能提供帮助/
PBE是比较好的交换关联能,但是对于磁性,最后加上U结果可能会好点。但是U的确定需要从文献和其他软件得到
我算的磁性没有f电子,这是为什么呢?是赝势的问题还是将f电子限制在芯内了?
电子的确是很深的,一般很难和其他原子的电子相互作用,这也是La系和Ac的元素的化学表现很相似的原因
那么怎样才能使磁性计算出现f电子呢?确实让人纠结啊!
请问使用PBEsol+U进行优化和性质计算,如何设置INCAR? 在vasp5.2手册上找不到PBEsol+U的说明,只有LSDA+U的
PBE是GGA类的交换关联能,LSDA的设置是可以同样用于GGA的
本人在用vasp计算遇到一个问题
用vasp进行结构优化过程中,在建立POTCAR赝势文件的时候,用到了K元素的paw gga ,但是软件里面出现了k pv和 k sv,本人不知道用哪一个,有位高手知道,请指点迷津,
X_sv表示把s电子作为半芯态处理,X_pv考虑把p电子作为半芯态来处理。这些元素一般很难赝化,特别是与电负性很强的元素(如F)结合时计算误差比较大。选择哪个,主要与计算精度有关。对于K,可选用pv这个。
谢谢你的回答,我的k点设置的是1.1.1.试问用pv这个可以吗?? 我再问你一个问题:我把优化好的Na结构中CONTCAR,直接用来做优化K元素的POSCAR可以吗?
简单的表示:cp CONTCAR (Na优化好的) POSCAR(作为k的) 只需要改变他的POTCAR.
他们两个最初的结构都一样,我想这样可以缩短时间,不知道可以吗? 谢谢你的指点, 本人急需
有哪位有钾(K)元素的paw gga的赝势文件,我的学校用的是盗版的vasp软件,但是赝势文件夹里面只出现了k pv和 k sv,本人需要正常的赝势, 有哪位k的赝势文件,麻烦你传一下。 非常的感谢 谢谢!!!!!!
你好,经查证,正版vasp里面也是没有K的potpaw_GGA的不加pv和sv的赝势 所以根本不需要其他人来传了