计算化学公社

 找回密码 Forget password
 注册 Register
Views: 1537|回复 Reply: 4
打印 Print 上一主题 Last thread 下一主题 Next thread

[VASP] vasp+phonopy计算声子谱报错

[复制链接 Copy URL]

4

帖子

0

威望

167

eV
积分
171

Level 3 能力者

如题,计算声子谱就会出现下面的问题
我是在优化后(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

求助求助!!!



5万

帖子

99

威望

5万

eV
积分
112356

管理员

公社社长

2#
发表于 Post on 2023-9-7 19:52:36 | 只看该作者 Only view this author
发帖的时候看清楚板块再发,这在首页的公告栏里已经明确提醒了,求助帖别发到资源分享子版块,这次给你移动了,以后发错地方直接删帖+扣分处理。
北京科音自然科学研究中心http://www.keinsci.com)致力于计算化学的发展和传播,长期开办高质量的各种计算化学类培训:初级量子化学培训班中级量子化学培训班高级量子化学培训班量子化学波函数分析与Multiwfn程序培训班分子动力学与GROMACS培训班CP2K第一性原理计算培训班,内容介绍以及往届资料购买请点击相应链接查看。这些培训是计算化学从零快速入门以及进一步全面系统性提升研究水平的高速路!培训各种常见问题见《北京科音办的培训班FAQ》
欢迎加入“北京科音”微信公众号获取北京科音培训的最新消息、避免错过网上有价值的计算化学文章!
欢迎加入人气非常高、专业性特别强的综合性理论与计算化学交流QQ群“思想家公社QQ群”:1号:18616395,2号:466017436,3号:764390338,搜索群号能搜到哪个说明目前哪个能加,合计9000人。北京科音培训班的学员在群中可申请VIP头衔,提问将得到群主Sobereva的最优先解答。
思想家公社的门口Blog:http://sobereva.com(发布大量原创计算化学相关博文)
Multiwfn主页:http://sobereva.com/multiwfn(十分强大的量子化学波函数分析程序)
ResearchGate:https://www.researchgate.net/profile/Tian_Lu
Money and papers are rubbish, get a real life!

16

帖子

0

威望

1306

eV
积分
1322

Level 4 (黑子)

3#
发表于 Post on 2023-9-7 21:57:55 | 只看该作者 Only view this author
显然第一个离子步完全没收敛

1060

帖子

0

威望

3259

eV
积分
4319

Level 6 (一方通行)

4#
发表于 Post on 2023-9-8 10:10:32 | 只看该作者 Only view this author
你的结构已经优化好了,算声子谱为什么还要加 ISIF = 2, POTIM = 0.5 呢……

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

1

帖子

0

威望

7

eV
积分
8

Level 1 能力者

5#
发表于 Post on 2024-5-22 16:40:31 | 只看该作者 Only view this author
您好,最后这个问题是怎么解决的。LR internal error: the vector H(1)-e(1) S(1) |phi(0)> is not orthogona

本版积分规则 Credits rule

手机版 Mobile version|北京科音自然科学研究中心 Beijing Kein Research Center for Natural Sciences|京公网安备 11010502035419号|计算化学公社 — 北京科音旗下高水平计算化学交流论坛 ( 京ICP备14038949号-1 )|网站地图

GMT+8, 2024-11-24 16:02 , Processed in 0.189130 second(s), 27 queries , Gzip On.

快速回复 返回顶部 返回列表 Return to list