# Electronic relaxation ENCUT=500 EDIFF=1E-8 ISMEAR=0 SIGMA=0.05 NELMIN=5 # Ionic relaxation NSW=100 EDIFFG=-1e-3 ISIF=2 IBRION=2 POTIM=0.20 # Other Tags #PSTRESS= #Write flags LWAVE=.FALSE. LCHARG=.FALSE. KPAR=2 ! cat > incar_stc<<! ISTART=0 ICHARG=2 PREC=Accurate GGA = ...
atoms.get_stress(),ASE的stress单位依赖于势函数,DP势时,单位为eV/Å^3. 正/负值对应拉/压应力 VASP 设置ISIF,计算stress,OUTCAR中搜关键字"in kB" in kB pxx pyy pzz pxy pxz pyz external pressure= (pxx+pyy+pzz)/3 kB VASP中stres的单位为kB,即kbar=1000 bar=0.1 GPa, 正/负值对应压/拉...
若不设置 vdw 校正,则输出 “--” 。 8、Pullay:OUTCAR 中最后一个 Pullay stress 对应的值。不加压计算,该值为 0。 9、sta_sym、dyn_sym、mag_sym:分别为 static configuration、dynamic configuration、magnetic configuration 对应的点群。 10、diff(%):三个格矢三个分量以及各个晶格常数 (POSCAR - CONTC...
sif=3对于2D系统没法使用，否则很容易遇到两种情形，要么真空层塌掉，要么就是stress根本收敛不了。除了...
PREC是控制计算精度最重要的一个参数,它决定了ENCUT、FFT网格大小、和ROPT的默认值。可能的取值为Low,Medium,High,Normal,Accurate(后两个只能在4.5以上版本中才起作用)。在一般的计算时推荐:4.5版本中用Normal,4.4版本中用Medium。当要提高力和Stresstensor的计算精度时,可以采用High或Accurate,并手动设置ENCUT的值。
其中,C_ij是第i个应力分量(i = 1,2,...,6)和第j个应变分量(j = 1,2,...,6)之间的弹性常数,stress_i是在第i个应变分量下计算得到的应力,stress_0_i是在未应变状态下计算得到的应力,strain_j是第j个应变分量。 使用以上计算方法,我们可以得到材料的弹性常数。然而,需要注意的是,计算弹性常数时,需要...
详情可以看这里 https://www.vasp.at/wiki/index.p ... ns_and_Pulay_stress
Trace only表示仅有总压力是正确的,总压力也是在OUTCAR文件中这一行“external pressure=.. kB”给出的。在对原胞的体积或形状进行优化时,ENCUT要略取的大一些(比如取为1.3*ENCUT的默认值或者设置PREC=High)以消除Pulay Stress导致的误差。 §5.7 ISMEAR和SIGMA ...
一般计算取NORMAL;当要提高Stress tensor计算精度时,HIGH 或ACCURATE,并手动设置ENCUT。 d. EDIFF & EDIFFG EDIFF 判断电子结构部分自恰迭代时自恰与否,一般取默认值=1E-4; EDIFFG 控制离子部分驰豫 e. ISTART & ICHARGE ISTART = 1, ICHARG = 11:能带结构、电子态密度计算时; ...
54、harge: 1-file 2-atom 10-co nstEDIFFG= -0.001NSW = 1000 # number of steps for IOMIBRION = 2 # ionic relax: 0-MD 1-quasi-New 2-CGISIF = 3 # stress and relaxation(这里变为3,是因为我们拟合的未必是最合适的,通过这个参数可以适当调节到更好)后面的可以不要了,如果算不动加上也无妨...