计算化学公社

标题: vasp+phonopy计算声子谱报错 [打印本页]

作者
Author: honghong 时间: 2023-9-7 17:29
标题: vasp+phonopy计算声子谱报错
如题，计算声子谱就会出现下面的问题
我是在优化后（1.0E-8）的晶胞基础上计算的。

错误提示如下：
DAV: 198 -0.630501814786E+03 0.75800E+01 -0.13223E+02  9780 0.270E+01 0.697E+00
DAV: 199 -0.638080200661E+03 -0.75784E+01 -0.14537E+02  9780 0.323E+01 0.886E+00
DAV: 200 -0.630511331345E+03 0.75689E+01 -0.13046E+02  9840 0.270E+01
1 F= -.63051133E+03 E0= -.63051133E+03  d E =-.218979E-15
Linear response reoptimize wavefunctions to high precision
DAV: 1 -0.637575821222E+03 -0.70645E+01 -0.14542E+02  9780 0.323E+01
DAV: 2 -0.637598035219E+03 -0.22214E-01 -0.22213E-01 12180 0.394E-01
DAV: 3 -0.637602591328E+03 -0.45561E-02 -0.45561E-02 15000 0.140E-01
Linear response G [H, r] |phi>, progress :
  Direction: 1
   N    E                   dE          d eps    ncg    rms
LRF_COMMUTATOR internal error: the vector H(1)-e(1) S(1) |phi(0)> is not orthog
onal to |phi(0)> (1.430865557188708E-004,-9.491780666027842E-012)
LRF_COMMUTATOR internal error: the vector H(1)-e(1) S(1) |phi(0)> is not orthog
onal to |phi(0)> (1.321338402906881E-004,-3.076494571249699E-012)
LRF_COMMUTATOR internal error: the vector H(1)-e(1) S(1) |phi(0)> is not orthog
onal to |phi(0)> (1.387371312654702E-004,1.902559975510483E-012)
LRF_COMMUTATOR internal error: the vector H(1)-e(1) S(1) |phi(0)> is not orthog
onal to |phi(0)> (1.085040326267353E-004,-5.290858357140674E-011)

我的INCAR：
ENCUT = 600
ISTART = 0
ICHARG = 2
ISMEAR = 0
SIGMA = 0.05
NELM = 200
NSW = 1
ALGO = N
IBRION = 8
POTIM = 0.5
ISIF = 2
ISYM = 2
EDIFF = 1.0E-8
EDIFFG = -0.005
LCHARG = F
LWAVE = F
LREAL = F
LPLANE= T
LEPSILON = T

求助求助！！！

作者
Author: sobereva 时间: 2023-9-7 19:52
发帖的时候看清楚板块再发，这在首页的公告栏里已经明确提醒了，求助帖别发到资源分享子版块，这次给你移动了，以后发错地方直接删帖+扣分处理。

作者
Author: Doria 时间: 2023-9-7 21:57
显然第一个离子步完全没收敛

作者
Author: 乐平 时间: 2023-9-8 10:10
你的结构已经优化好了，算声子谱为什么还要加 ISIF = 2, POTIM = 0.5 呢……

另外，如果你不需要计算介电性质，有效电荷，也不需要加 LEPSILON = T

作者
Author: mtian2022 时间: 2024-5-22 16:40
您好，最后这个问题是怎么解决的。LR internal error: the vector H(1)-e(1) S(1) |phi(0)> is not orthogona

欢迎光临计算化学公社 (http://bbs.keinsci.com/)