24小时热门版块排行榜    

CyRhmU.jpeg
查看: 4046  |  回复: 8
当前只显示满足指定条件的回帖,点击这里查看本话题的所有回帖

ljzhou86

金虫 (正式写手)

[交流] vasp 分子动力学NVT系综模拟,能量不断在增加,已有6人参与

最近模拟一个有机无机杂化的太阳能电池,MAPbI3阳离子MA掺杂FA的分子动力学模拟,建立了一个模型,由于涉及到小分子的振动无序与相互朝向及不同位点掺杂的情况,模型很多,所以想通过我想通过vasp的分子动力学NVT模拟之后,提取能量低的构型再做优化,最后算电子结构。现在遇到如下问题:
模拟过程中,能量不断在增加,很不正常。按理能量应该在一平衡位置波动才对。大家有知道怎么改进的吗?是因为初始模型的问题导致的,还是因为相关参数的特别设置?incar与相关的结果图 (由于机时限制,跑了两次2.5ps, 总共5ps,这是第二次的结果)如下: 
IBRION=0
NPAR =8

ISIF=0
#NBLOCK = 1
#KBLOCK =50
EDIFF = 1E-4
ISMEAR= 0;
SIGMA=0.05
PREC= Low

AMIX = 0.02 (对涉及到小分子,这个参数似乎很重要。试过用默认的值,能量向不断减小的发散;对于纯无机的材料,不设这个没有这个问题。)
#ALGO = F
IALGO=38
LREAL=AUTO
NELMIN=5
NELM = 70
NSW=15000
POTIM=1
ISYM=0

SMASS=1
TEBEG=400
TEEND=400

LWAVE= .FALSE.
LCHARG= .FALSE.
IVDW=11
GGA= PE

vasp 分子动力学NVT系综模拟,能量不断在增加,
2cell_energy.png
回复此楼

» 猜你喜欢

» 本主题相关商家推荐: (我也要在这里推广)

» 本主题相关价值贴推荐,对您同样有帮助:

已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

ljzhou86

金虫 (正式写手)

引用回帖:
2楼: Originally posted by luzihen at 2017-03-30 14:04:51
是不是因为PREC= Low,能量漂移太严重了.你仔细看下手册

动力学模拟,体系大,用PREC=low应该是可以的。我先前也试过用accurate,同样有偏离
3楼2017-03-31 02:49:09
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
查看全部 9 个回答

luzihen

金虫 (正式写手)

★ ★ ★
小木虫: 金币+0.5, 给个红包,谢谢回帖
漫天飘雪: 金币+2, 谢谢交流 2017-03-30 17:00:14
是不是因为PREC= Low,能量漂移太严重了.你仔细看下手册
2楼2017-03-30 14:04:51
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

羚轩惠子

新虫 (初入文坛)


小木虫: 金币+0.5, 给个红包,谢谢回帖
想问一下你的问题解决了么?  我最近做的时候也遇到了相同的问题,想向你学习一下。
4楼2017-04-12 21:42:00
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

largehoo

木虫 (小有名气)


小木虫: 金币+0.5, 给个红包,谢谢回帖
我的是往小漂,用了NVT系综,体系温度从300K减小到50K左右
5楼2017-05-15 18:50:29
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
普通表情 高级回复(可上传附件)
信息提示
请填处理意见