total drift太大的问题
在做一个结构的弛豫的时候遇到了这个问题
total drift: 19.934551 34.617080 2.916653
一般这个都是接近0的,之前也一直以为它是上面力的总和。
在我这个例子里,无论力和压力收敛的再好,这个total drift始终很大,非常困扰。
从vaspwiki上还有冲浪得到的结果来看,他们提出的方法包括
增加截断能
增大FFT网格(NGXF)
关掉实空间投影 LREAL
三个都试了一圈无果。后来在组里师弟的帮助下发现是由于ISMEAR的SIGMA取的太大的缘故。
我之前取得都是0.2,改成0.02之后就没有这个问题了。
VASP 加SOC后受力异常的问题
最近在计算WSe2中的SOC,本来想的是这么流行的体系应该很好算,结果VASP狠狠给我上了一课。弛豫好的结构放到SOC里受力都是个位数的 eV/Angstrom。怎么改参数都不行,搜也搜不着,而且别人都用VASP算这个发了文章了。没办法只能一个参数一个参数的测试,最后发现竟然是vasp_ncl
版本的并行的问题。去掉NCORE
参数后受力就正常了,目测是vasp_ncl
的并行的问题。可以考虑使用KPAR
并行,这个经测试没问题
顺便搜了一下vasp并行策略的问题,根据官方的PPT里说的,KPAR是k点并行是最有效率的,NPAR是能带并行次之。
计算SOC的参数
#SOC related input
ISYM = -1
# ISPIN = 2
LSORBIT = .TRUE.
# LNONCOLLINEAR = .TRUE. 打开LSORBIT后这个会自动打开
GGA_COMPAT = .FALSE.
LMAXMIX = 6 ##我这里是WSe2,涉及f轨道需要设为6
VASP优化结构中断问题
我这里是因为没有设置EDIFF,VASP在优化结构的时候还没达到力收敛标准的时候就能量收敛了(因为程序取的默认值比我要的偏大)。所以报这个错误。加上一个EDIFF参数就好了。
ZBRENT: fatal error in bracketing please rerun with smaller EDIFF, or copy CONTCAR to POSCAR and continue
结构“噪音”问题
计算TiSe2的时候算2x2x2结构的时候取的超胞,可能是程序检测到我算的结构可能不是最小元胞所以报这个错误,不会影响到后续的计算。
The determination of the symmetry of your systems shows a strong dependence on the tolerance parameter SYMPREC. This is often the result of "noise" in the structural parameters (lattice vectors and/or ionic positions). Set NWRITE=3 for more details.
VASP+phonopy生成力常数问题
今天遇到一个非常愚蠢的问题,虽然非常没脑子,但是还是记录在这里,以防再次遇到或遇到和我一样没脑的操作
在执行
phonopy -fc vasprun.xml
的时候遇到报错
(phonopy) [xbliu@login01 phonon]$ phonopy -fc vasprun.xml
_
_ __ | |__ ___ _ __ ___ _ __ _ _
| '_ \| '_ \ / _ \| '_ \ / _ \ | '_ \| | | |
| |_) | | | | (_) | | | | (_) || |_) | |_| |
| .__/|_| |_|\___/|_| |_|\___(_) .__/ \__, |
|_| |_| |___/
2.11.0
Python version 3.9.7
Spglib version 1.16.2
<_io.TextIOWrapper name='vasprun.xml' mode='r' encoding='UTF-8'>
Traceback (most recent call last):
File "/home/users/xbliu/anaconda3/envs/phonopy/bin/phonopy", line 43, in <module>
main(**argparse_control)
File "/home/users/xbliu/anaconda3/envs/phonopy/lib/python3.9/site-packages/phonopy/cui/phonopy_script.py", line 1552, in main
settings, confs, cell_filename = read_phonopy_settings(
File "/home/users/xbliu/anaconda3/envs/phonopy/lib/python3.9/site-packages/phonopy/cui/phonopy_script.py", line 1361, in read_phonopy_settings
phonopy_conf_parser = PhonopyConfParser(
File "/home/users/xbliu/anaconda3/envs/phonopy/lib/python3.9/site-packages/phonopy/cui/settings.py", line 1319, in __init__
self.read_file() # store .conf file setting in self._confs
File "/home/users/xbliu/anaconda3/envs/phonopy/lib/python3.9/site-packages/phonopy/cui/settings.py", line 309, in read_file
left, right = [x.strip() for x in line.split('=')]
ValueError: too many values to unpack (expected 2)
命令输错了,应该是
phonopy --fc vasprun.xml
两个横杠
VERY BAD NEWS! internal error in subroutine IBZKPT
应该是对称性的问题,把SYMPREC 设置为 1e-8也就是默认值后解决参考VERY BAD NEWS! internal error in subroutine IBZKPT - My Community (vasp.at)
计算超胞声子遇到问题
HNFORM: ERROR: k-point generating vectors and reciprocal lattice are incommensurate.
虽然提示是k点网格和晶格的问题,但是调整k点网格和晶格常数一点用也没有,网上查到可以通过ISYM=0
来解决,确实能解决,但是计算量增加也是非常多,最后发现可以改SYMPREC = 1.0E-2
解决