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

24小时热门版块排行榜    

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

405963512

铜虫 (正式写手)

[交流] 【求助】vasp优化，不动了

大家好，请问个优化的问题，我优化有机分子晶体，参数选择如下
SYSTEM=MAKXUS-CRYSTAL
ENCUT=400
ISTART=0;ICHARG=2
EDIFF=0.01;EDIFFG=-0.3
#GGA=PE;VOSKOWN=1!(EX FUNCTION)
PREC=Accurate
#ISPIN=2
ISMEAR=-5
#;SIGMA=0.2!(DOS=-5,BAND=0)
LREAL=AUTO;IALGO=48
NSW=100;IBRION=2;ISIF=2
#LORBIT=11;RWIGS=

#NBANDS=?(BAND STUCTURE)
   
#NSW=60;IBRION=5;NFREE=2;POTIM=0.005!(Freq)
优化的过程中，
vasp.4.6.28 25Jul05 complex
POSCAR found :  5 types and   92 ions
LDA part: xc-table for Pade appr. of Perdew
POSCAR, INCAR and KPOINTS ok, starting setup
FFT: planning ...          16
reading WAVECAR
entering main loop
       N       E                     dE             d eps       ncg     rms          rms(c)
RMM:   1     0.457607845383E+04    0.45761E+04   -0.14304E+05  6669   0.111E+03
RMM:   2     0.141410925763E+04   -0.31620E+04   -0.40686E+04  6669   0.313E+02
RMM:   3     0.106984143053E+03   -0.13071E+04   -0.14647E+04  6669   0.199E+02
RMM:   4    -0.350068053850E+03   -0.45705E+03   -0.45069E+03  6669   0.161E+02
RMM:   5    -0.529845862522E+03   -0.17978E+03   -0.13845E+03  6669   0.927E+01
RMM:   6    -0.594816965220E+03   -0.64971E+02   -0.47587E+02  6669   0.734E+01
RMM:   7    -0.625096069873E+03   -0.30279E+02   -0.23504E+02  6669   0.377E+01
RMM:   8    -0.638126707492E+03   -0.13031E+02   -0.92435E+01  6669   0.334E+01
RMM:   9    -0.651314270335E+03   -0.13188E+02   -0.11581E+02 18889   0.160E+01
RMM:  10    -0.652155507676E+03   -0.84124E+00   -0.80296E+00 17092   0.280E+00
RMM:  11    -0.652376124155E+03   -0.22062E+00   -0.15340E+00 12364   0.101E+00
RMM:  12    -0.652409163161E+03   -0.33039E-01   -0.30327E-01  9804   0.475E-01    0.833E+01
RMM:  13    -0.803751613690E+03   -0.15134E+03   -0.15536E+03 18253   0.832E+01    0.577E+01
RMM:  14    -0.636190294339E+03    0.16756E+03   -0.26571E+02 13413   0.301E+01    0.341E+01
RMM:  15    -0.597116211420E+03    0.39074E+02   -0.29988E+02 14574   0.235E+01    0.486E+01
RMM:  16    -0.584882735735E+03    0.12233E+02   -0.14395E+01 13729   0.789E+00    0.675E+00
RMM:  17    -0.583352343181E+03    0.15304E+01   -0.29749E+00 15161   0.262E+00    0.635E+00
RMM:  18    -0.581960799223E+03    0.13915E+01   -0.51541E+00 13685   0.509E+00    0.140E+00
在13步出现能量差突然变大，再之后，到18步就不动了，请大家指教是什么原因？
我的构型从数据库提取的。金币不多，大家帮忙！谢谢
回复此楼

» 猜你喜欢

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

fanchen021

铁杆木虫 (正式写手)
★
qasd(金币+1,VIP+0):xiexie~ 11-30 20:14
ISMEAR 改成0 试试呢。
SIGMA 前的注释号去掉
赞一下(10人) 回复此楼
2楼2009-11-29 19:52:47
 已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

gavinliu7390

木虫 (著名写手)

叶落鹰飞
vasp很难处理好分子的优化。
这个其实很正常。
如果你看一下CONTCAR，你会发现优化出来的结构很奇怪。
赞一下 回复此楼
真理是一点点接近的！
3楼2009-11-29 21:06:53
 已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

wuchenwf

荣誉版主 (职业作家)
★ ★
405963512(金币+1,VIP+0): 11-29 23:12
qasd(金币+1,VIP+0):xiexie~ 11-30 20:14
1 首先楼主要确认以前做其它体系的时候没有问题
2 如楼上所说优化分子体系很麻烦，你需要控制如真空层等因素（你是这个搞的吧），
3 你的体系有多大，我觉得收敛标准太低了，你可以提高点
赞一下(11人) 回复此楼
4楼2009-11-29 21:17:09
 已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

405963512

铜虫 (正式写手)
引用回帖:
Originally posted by wuchenwf at 2009-11-29 21:17:
1 首先楼主要确认以前做其它体系的时候没有问题
2 如楼上所说优化分子体系很麻烦，你需要控制如真空层等因素（你是这个搞的吧），
3 你的体系有多大，我觉得收敛标准太低了，你可以提高点

我是做能带计算，算别的没问题，可能这个原子数目有点多？90多个吧。
还有，精度提高了能有效？

[ Last edited by 405963512 on 2009-11-29 at 23:14 ]
赞一下 回复此楼
5楼2009-11-29 23:13:08
 已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

wuchenwf

荣誉版主 (职业作家)
像楼上说的分子优化的时候不是vasp强项，你可以用dmol3试试
赞一下 回复此楼
6楼2009-11-29 23:14:02
 已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

405963512

铜虫 (正式写手)
引用回帖:
Originally posted by wuchenwf at 2009-11-29 23:14:
像楼上说的分子优化的时候不是vasp强项，你可以用dmol3试试

你说的对，这个分子晶体我用dmol3优化过去了，还是挺快的，难道vasp不行，或者这两者能混用？
赞一下 回复此楼
7楼2009-11-29 23:15:06
 已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

wuchenwf

荣誉版主 (职业作家)
★
405963512(金币+1,VIP+0): 11-29 23:21
你用dmol得到的结果再用vasp搞一下如何
赞一下(1人) 回复此楼
8楼2009-11-29 23:18:19
 已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

405963512

铜虫 (正式写手)
引用回帖:
Originally posted by wuchenwf at 2009-11-29 23:18:
你用dmol得到的结果再用vasp搞一下如何

你是这个意思啊，我尝试一下，谢谢
赞一下 回复此楼
9楼2009-11-29 23:21:12
 已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

wangzg9292

木虫 (正式写手)
★
wuchenwf(金币+1,VIP+0):谢谢 11-30 20:50
引用回帖:
Originally posted by 405963512 at 2009-11-29 23:13:

我是做能带计算，算别的没问题，可能这个原子数目有点多？90多个吧。
还有，精度提高了能有效？

[ Last edited by 405963512 on 2009-11-29 at 23:14 ]

做能带计算最好用ISMEAR默认值，或者0来计算
赞一下(10人) 回复此楼
10楼2009-11-30 08:28:19
 已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 405963512 的主题更新
121/2返回列表12下一页
普通表情 龙 兔 虎 猫 高级回复(可上传附件)
信息提示
关闭
请填处理意见
关闭