应《网络安全法》要求，自2017年10月1日起，未进行实名认证将不得使用互联网跟帖服务。为保障您的帐号能够正常使用，请尽快对帐号进行手机号验证，感谢您的理解与支持！

24小时热门版块排行榜    

小木虫论坛-学术科研互动平台 » 计算模拟区 » 第一性原理 » Vasp&MedeA » 用收敛过的WAVECAR初始频率计算不收敛
51/1返回列表
查看: 1055  |  回复: 5
只看楼主 @他人 存档 新回复提醒 (忽略) 收藏    在APP中查看
当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖

zhangguangping

木虫 (著名写手)

[求助] 用收敛过的WAVECAR初始频率计算不收敛 已有1人参与

请教大家一个问题：我用10个核跑了一个电偶极矩矫正的slab体系的电子结构，收敛了，得到了WAVECAR，然后用它来初始了相同体系的做频率的计算（也做了偶极矫正），计算的时候用的是8个核，为什么不收敛了呢？
预期的情况是用收敛过的WAVECAR初始频率计算，会很快收敛开始的电子自洽，然后进入频率计算。现在是开始的电子自洽都不收敛了。不理解

附：
1.偶极矫正的INCAR
SYSTEM = M2-Au-E0
IALGO  = 48
NELM   = 500
EDIFF = 1E-4
EDIFFG = -0.02
NSW = 0
ISMEAR = 1
SIGMA  = 0.1
PREC   = Accurate
ENCUT  = 400
IBRION = 2
ISPIN  = 1
LCHARG = .FALSE.
LREAL  = Auto
LWAVE  = .TRUE.
LDIPOL = .TRUE.
IDIPOL = 3
EFIELD = 0


2.偶极矫正的频率计算INCAR
SYSTEM = M2-Au-freq
IALGO  = 48
NELM   = 500
EDIFF = 1E-4
EDIFFG = -0.02
NSW = 2
ISMEAR = 1
SIGMA  = 0.1
PREC   = Accurate
ENCUT  = 400
IBRION = 5
POTIM = 0.01
NFREE = 2
ISPIN  = 1
LCHARG = .FALSE.
LREAL  = Auto
LWAVE  = .FALSE.
LDIPOL = .TRUE.
IDIPOL = 3
NWRITE = 3
回复此楼

» 猜你喜欢
弘德明志博学笃行
1楼 2016-03-30 11:13:52
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

zhangguangping

木虫 (著名写手)
引用回帖:
2楼: Originally posted by fzuer at 2016-03-30 17:02:12
个人感觉是频率计算中会对结构做微扰，估计这种情况下有的构型就没有收敛

不知道适当的POTIM是否可以解决这个问题

现在发现问题的原因了，但是不明白为什么。

原因是采用偶极矫正的WAVECAR再次初始化相同的偶极矫正任务的时候，收敛的WAVECAR不是一个好的初始猜测。这个时候需要用没有偶极矫正的WAVECAR才行。

我用没有偶极矫正收敛的WAVECAR初始偶极矫正的计算或者频率计算（做了偶极矫正）
这个时候很快收敛：
entering main loop
       N       E                     dE             d eps       ncg     rms          rms(c)
RMM:   1    -0.617285334113E+03   -0.61729E+03   -0.28948E-05   958   0.329E-03    0.356E-01
RMM:   2    -0.633045700125E+03   -0.15760E+02   -0.22436E-01  2091   0.263E-01    0.131E+01
RMM:   3    -0.617613909416E+03    0.15432E+02   -0.14867E-01  2087   0.233E-01    0.325E+00
RMM:   4    -0.617282326516E+03    0.33158E+00   -0.10908E-02  2085   0.729E-02    0.156E-01
RMM:   5    -0.617282748456E+03   -0.42194E-03   -0.12005E-03  1997   0.282E-02    0.204E-01
RMM:   6    -0.617281688730E+03    0.10597E-02   -0.10445E-04  1208   0.850E-03    0.374E-02
RMM:   7    -0.617281678568E+03    0.10162E-04   -0.16863E-05  1042   0.361E-03
   1 F= -.61728168E+03 E0= -.61723961E+03  d E =-.126213E+00
Finite differences POTIM= 0.01000 DOF=  78
bond charge predicted
赞一下(1人) 回复此楼
弘德明志博学笃行
5楼2016-03-31 06:42:21
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
查看全部 6 个回答

fzuer

铜虫 (小有名气)

【答案】应助回帖

感谢参与，应助指数 +1
个人感觉是频率计算中会对结构做微扰，估计这种情况下有的构型就没有收敛

不知道适当的POTIM是否可以解决这个问题
赞一下 回复此楼
功能材料的DFT计算
2楼2016-03-30 17:02:12
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

zhangguangping

木虫 (著名写手)
引用回帖:
2楼: Originally posted by fzuer at 2016-03-30 17:02:12
个人感觉是频率计算中会对结构做微扰，估计这种情况下有的构型就没有收敛

不知道适当的POTIM是否可以解决这个问题

感谢您的回复。我觉得不是您说的这个问题。现在的自洽是在频率计算的时候一开始的电子结构计算。这个计算是对给出的体系得到电子结构。然后进行下面的“围绕”去计算。现在是对未“微扰”的体系的电子结构不收敛，我是用的已经收敛过的WAVECAR初始化的。非常奇怪的问题。
赞一下 回复此楼
弘德明志博学笃行
3楼2016-03-30 18:50:31
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

zhangguangping

木虫 (著名写手)
引用回帖:
2楼: Originally posted by fzuer at 2016-03-30 17:02:12
个人感觉是频率计算中会对结构做微扰，估计这种情况下有的构型就没有收敛

不知道适当的POTIM是否可以解决这个问题

更有趣的是，我用已经收敛的WAVECAR去初始化得到WAVECAR的电子结构计算，现在也不收敛了。不知道当时怎么收敛好的。
赞一下 回复此楼
弘德明志博学笃行
4楼2016-03-30 18:51:54
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
查看全部 6 个回答
普通表情 龙 兔 虎 猫
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[教师之家] 售SCI一区T0P文章，我:8.O55.1.O.54,科目全,可十急 +3 ero8OE6tv9cu 2026-09-02 8/400 2026-09-03 14:48 by T0rGB46095mJ
[硕博家园] 售一区SCI文章T0P，我:8O.551.O54,科目全,可十急 +3 YLHlRHNwYkce 2026-09-02 5/250 2026-09-03 13:55 by T0rGB46095mJ
[基金申请] 国社科又开始会评了，不知道这次命运如何 +10 雨打竹帘 2026-08-30 14/700 2026-09-03 13:17 by qsd10086
[找工作] 售SCI一区T0P文章，我:8O.55.1.O.54,科目全,可伽急 +3 YLHlRHNwYkce 2026-09-02 8/400 2026-09-03 12:05 by T0rGB46095mJ
[教师之家] 售SCI一区文章，我:8.O.55.1.O.54,科目齐全,可伽急 +3 YLHlRHNwYkce 2026-09-02 8/400 2026-09-03 12:03 by T0rGB46095mJ
[基金申请] 麻烦专家们看看评委们的意见（F口面上） +8 gdd2018 2026-08-28 13/650 2026-09-03 06:22 by 羽毛枫f
[公派出国] 售SCI-T0P文章，我:8O.5.5.1.O.54,科目齐全,可+急 +3 ero8OE6tv9cu 2026-09-02 4/200 2026-09-03 03:14 by rM1TE0WVDIIY
[博后之家] 售SCI文章，我:8O5.5.1.O.54,科目齐全,可+急 +3 ero8OE6tv9cu 2026-09-02 6/300 2026-09-02 23:59 by rM1TE0WVDIIY
[基金申请] 科研人应该花精力去思考如何解决问题，而不是去凝练问题 +7 瞬息宇宙 2026-09-01 14/700 2026-09-02 17:33 by ma0526
[基金申请] 要骂人了，新模版改版就是要淡化问题凝练这种虚的东西，结果有个评委还在说凝练得不够 +9 瞬息宇宙 2026-08-31 17/850 2026-09-02 14:04 by anjeeshine
[基金申请] 基金系统什么内容也没有 30+4 winsaint 2026-08-27 10/500 2026-09-02 11:35 by 大不刘6
[基金申请] 学科评审组评审是指会评吗？ +5 瞬息宇宙 2026-08-31 5/250 2026-09-02 10:13 by 雨冰共舞
[基金申请] 面上函评意见出来了，像什么等级？ 20+4 Tsingking1 2026-08-27 17/850 2026-09-01 19:51 by 超级无敌华子
[论文投稿] 小白求助 投论文要求的highlights应该如何写 5+3 l1963982152 2026-08-29 4/200 2026-09-01 09:04 by 北京莱茵编辑
[基金申请] 面上意见出来了 +12 黄鸟于飞Chao 2026-08-29 23/1150 2026-08-31 18:57 by 黄鸟于飞Chao
[基金申请] 能否申诉？ +7 echo8914667 2026-08-30 8/400 2026-08-31 17:00 by yihongxu
[基金申请] 29号明天会评吗 +4 笨笨唐 2026-08-28 4/200 2026-08-31 09:30 by huixian257
[基金申请] 为什么到现在没收到通知？ +5 tannykie 2026-08-29 5/250 2026-08-30 21:05 by purplejack
[考博] 找导师 +6 yuanjiabao 2026-08-29 7/350 2026-08-30 14:40 by 生科新手
[基金申请] 国自然评审意见 +13 wangmingqi 2026-08-28 19/950 2026-08-29 10:22 by Poppy1104
信息提示
关闭
请填处理意见
关闭