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

24小时热门版块排行榜    

CyRhmU.jpeg
小木虫论坛-学术科研互动平台 » 计算模拟区 » 第一性原理 » Vasp&MedeA » 【求助】有机含金属晶体优化不收敛
41/1返回列表
查看: 405  |  回复: 3
只看楼主@他人 存档 新回复提醒 (忽略) 收藏    在APP中查看
当前主题已经存档。
【有奖交流】积极回复本帖子，参与交流，就有机会分得作者 405963512 的 8 个金币

405963512

铜虫 (正式写手)

[交流] 【求助】有机含金属晶体优化不收敛已有3人参与

大家好，请教个vasp优化的问题。
我的体系是有机含金属的晶体。分子中含有57原子，晶胞中两个分子。
在优化时输入文件
INCAR:
SYSTEM=F16CuPc-ribbon
ENCUT=500
ISTART=0;ICHARG=2
EDIFF=0.001;EDIFFG=-0.01
PREC=Accurate
ISMEAR=-5
LREAL=AUTO;IALGO=48
NSW=300;IBRION=2;ISIF=2
POTIM=0.3
但是优化过程不收敛，显式能量在震荡：

RMM:  16    -0.852798608629E+03    0.22539E+02   -0.10543E+01  3316   0.586E+00    0.681E+01
RMM:  17    -0.839190701776E+03    0.13608E+02   -0.43329E+01  2886   0.106E+01    0.133E+01
RMM:  18    -0.834180303109E+03    0.50104E+01   -0.76078E+00  3087   0.424E+00    0.867E+00
RMM:  19    -0.834193859965E+03   -0.13557E-01   -0.68726E-01  3104   0.133E+00    0.900E+00
RMM:  20    -0.833648909951E+03    0.54495E+00   -0.10395E+00  2950   0.179E+00    0.641E+00
RMM:  21    -0.833873805380E+03   -0.22490E+00   -0.55067E-01  2870   0.132E+00    0.758E+00
RMM:  22    -0.833127702597E+03    0.74610E+00   -0.48444E-01  2892   0.106E+00    0.257E+00
RMM:  23    -0.833229782385E+03   -0.10208E+00   -0.66226E-02  2783   0.373E-01    0.721E+00
RMM:  24    -0.833322878288E+03   -0.93096E-01   -0.34482E-02  2614   0.333E-01    0.718E+00
RMM:  25    -0.833131040765E+03    0.19184E+00   -0.36504E-02  2500   0.383E-01    0.677E+00
RMM:  26    -0.833445974919E+03   -0.31493E+00   -0.34157E-02  2419   0.362E-01    0.752E+00
RMM:  27    -0.833183354916E+03    0.26262E+00   -0.28995E-02  2262   0.354E-01    0.728E+00
RMM:  28    -0.833081141141E+03    0.10221E+00   -0.56150E-03  1902   0.135E-01    0.721E+00
RMM:  29    -0.833085097342E+03   -0.39562E-02   -0.21071E-03  1698   0.776E-02    0.688E+00
RMM:  30    -0.833352184105E+03   -0.26709E+00   -0.28333E-02  2298   0.393E-01    0.737E+00
RMM:  31    -0.833083572117E+03    0.26861E+00   -0.24805E-02  2120   0.326E-01    0.686E+00
RMM:  32    -0.833364680698E+03   -0.28111E+00   -0.25323E-02  2148   0.335E-01    0.737E+00
RMM:  33    -0.833072378214E+03    0.29230E+00   -0.25370E-02  2132   0.331E-01    0.685E+00
RMM:  34    -0.833367304470E+03   -0.29493E+00   -0.25167E-02  2127   0.329E-01    0.736E+00
RMM:  35    -0.833069849200E+03    0.29746E+00   -0.25507E-02  2121   0.329E-01    0.168E+00
RMM:  36    -0.833152832427E+03   -0.82983E-01   -0.10457E-02  1805   0.189E-01    0.720E+00
RMM:  37    -0.833089416410E+03    0.63416E-01   -0.10508E-02  1858   0.211E-01    0.691E+00
RMM:  38    -0.833344024313E+03   -0.25461E+00   -0.26024E-02  2135   0.349E-01    0.732E+00
RMM:  39    -0.833083460867E+03    0.26056E+00   -0.25620E-02  2101   0.329E-01    0.689E+00
RMM:  40    -0.833351117922E+03   -0.26766E+00   -0.25394E-02  2120   0.328E-01    0.732E+00
RMM:  41    -0.833078464599E+03    0.27265E+00   -0.25813E-02  2105   0.329E-01    0.688E+00
RMM:  42    -0.833352065423E+03   -0.27360E+00   -0.25390E-02  2115   0.327E-01    0.731E+00
RMM:  43    -0.833076433871E+03    0.27563E+00   -0.25864E-02  2107   0.329E-01    0.688E+00
RMM:  44    -0.833351437194E+03   -0.27500E+00   -0.25411E-02  2115   0.327E-01    0.730E+00
RMM:  45    -0.833075646604E+03    0.27579E+00   -0.25908E-02  2102   0.329E-01    0.688E+00
RMM:  46    -0.833350223385E+03   -0.27458E+00   -0.25431E-02  2112   0.326E-01    0.730E+00
RMM:  47    -0.833075546165E+03    0.27468E+00   -0.25946E-02  2103   0.329E-01    0.688E+00
RMM:  48    -0.833348789729E+03   -0.27324E+00   -0.25452E-02  2113   0.326E-01    0.729E+00
RMM:  49    -0.833075807235E+03    0.27298E+00   -0.25981E-02  2103   0.328E-01    0.688E+00
RMM:  50    -0.833347339404E+03   -0.27153E+00   -0.25476E-02  2112   0.326E-01    0.729E+00
RMM:  51    -0.833076255387E+03    0.27108E+00   -0.26011E-02  2103   0.329E-01    0.688E+00
RMM:  52    -0.833345943099E+03   -0.26969E+00   -0.25501E-02  2113   0.326E-01    0.729E+00
RMM:  53    -0.833076782018E+03    0.26916E+00   -0.26036E-02  2102   0.329E-01    0.688E+00
RMM:  54    -0.833344631231E+03   -0.26785E+00   -0.25521E-02  2111   0.326E-01    0.728E+00
RMM:  55    -0.833077342018E+03    0.26729E+00   -0.26040E-02  2102   0.330E-01    0.689E+00
以前也遇到过，但是不知道该怎么调整。
希望大家给个建议，谢谢大家。十金币鼓励。解决问题多加感谢
回复此楼

» 猜你喜欢

» 本主题相关商家推荐: (我也要在这里推广)
1楼2010-04-08 21:40:11
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

hn19870519

木虫 (正式写手)
★
qasd(金币+1):xiexie 2010-04-10 20:04
要不你改一下算法。。。IALGO=38试试？
赞一下(1人) 回复此楼
2楼2010-04-08 22:26:50
 已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

wvsw

木虫 (小有名气)

老虎大王
★
qasd(金币+1):xiexie 2010-04-10 20:05
ISMEAR改成0试试
赞一下(1人) 回复此楼
生于忧患，死于安乐
3楼2010-04-09 16:40:25
 已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

alick1102

铜虫 (小有名气)
★ ★
405963512(金币+2):谢谢，我先试试啊！ 2010-04-09 19:34
qasd(金币+2):xiexie 2010-04-10 20:05
引用回帖:
Originally posted by 405963512 at 2010-04-08 21:40:11:
大家好，请教个vasp优化的问题。
我的体系是有机含金属的晶体。分子中含有57原子，晶胞中两个分子。
在优化时输入文件
INCAR:
SYSTEM=F16CuPc-ribbon
ENCUT=500
ISTART=0;ICHARG=2
EDIFF=0.001;EDIFFG=-0. ...

这个可能是NBANDS取少了的原因。
NBANDS=(电子数)/2 + 2*(原子数)  ,如果这样取SCF还不收敛，建议再加50.

另外RMM-DIIS的迭代有时候会不收敛，你可以采取DAVISION和RMM联合的方法。
就是用 ALGO = very_FAST 取代你的IALGO=48

[ Last edited by alick1102 on 2010-4-9 at 17:04 ]
赞一下(1人) 回复此楼
到了博士的时候才知道，科学也不过是为了吃饭的。
4楼2010-04-09 17:01:35
 已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 405963512 的主题更新
41/1返回列表
普通表情 龙 兔 虎 猫 高级回复(可上传附件)
信息提示
关闭
请填处理意见
关闭