24小时热门版块排行榜    

查看: 1972  |  回复: 6

高朋满座007

铁杆木虫 (正式写手)

[求助] VASP怎样在固定的温度例如1000K做结构优化? 已有2人参与

求助各位朋友,我想在1000K温度时做晶格结构的优化,能量和力为收敛标准,能否指导INCAR怎么写?感激不尽
回复此楼
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

gavinliu7390

木虫 (著名写手)

叶落鹰飞

【答案】应助回帖


感谢参与,应助指数 +1
ljw4010: 金币+1, 谢谢交流! 2014-11-08 10:00:33
SYSTEM = MD
ENCUT = 500
EDIFF = 1e-5
NSW = 10000
POTIM = 1.0
ISMEAR = 0
SIGMA=0.05
ISIF=2
LREAL=T
ISYM=0
LCHARG=F
LWAVE=F

#MD tag
IBRION=0
SMASS=5.0
TEBEG=10000

NBLOCK=5
KBLOCK=1
IALGO=48
真理是一点点接近的!
2楼2014-11-07 08:57:48
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

高朋满座007

铁杆木虫 (正式写手)

引用回帖:
2楼: Originally posted by gavinliu7390 at 2014-11-07 08:57:48
SYSTEM = MD
ENCUT = 500
EDIFF = 1e-5
NSW = 10000
POTIM = 1.0
ISMEAR = 0
SIGMA=0.05
ISIF=2
LREAL=T
ISYM=0
LCHARG=F
LWAVE=F

#MD tag
IBRION=0
SMASS=5.0
TEBEG=10000

NBLOCK=5
KBLOCK=1 ...

我试了一下,发现温度波动的范围特别大
3楼2014-11-10 20:25:41
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

1134684571

木虫 (小有名气)

【答案】应助回帖


liliangfang: 金币+1, 谢谢交流 2014-11-12 07:49:01
引用回帖:
3楼: Originally posted by 高朋满座007 at 2014-11-10 20:25:41
我试了一下,发现温度波动的范围特别大...

TEBEG=1000是第一性分子动力学的初始温度
Default:     
TEBEG = 0
TEEND = TEBEG
也就是说TEEND=1000怎么会温度变动大呢?
4楼2014-11-11 15:52:40
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

gavinliu7390

木虫 (著名写手)

叶落鹰飞

引用回帖:
3楼: Originally posted by 高朋满座007 at 2014-11-10 20:25:41
我试了一下,发现温度波动的范围特别大...

波动的大小,跟你原子数的多少有关系。原子越多,波动越小。
真理是一点点接近的!
5楼2014-11-12 08:54:07
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

高朋满座007

铁杆木虫 (正式写手)

引用回帖:
5楼: Originally posted by gavinliu7390 at 2014-11-12 08:54:07
波动的大小,跟你原子数的多少有关系。原子越多,波动越小。...

谢谢,我体系中只有5个原子,NSW=100可否?我看到步数越多温度偏离越大
6楼2014-11-13 00:12:30
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

elfyu

新虫 (初入文坛)

引用回帖:
5楼: Originally posted by gavinliu7390 at 2014-11-12 08:54:07
波动的大小,跟你原子数的多少有关系。原子越多,波动越小。...

您好 我也是非晶材料的计算方向 AIMD初学者一名,想请教您几个问题,希望您能在百忙之中回复,真心谢谢您!我看文献,里面说:故首先让体系在2000 K 下运行2000 个离子步, 然后依次降温, 在1700 K, 1500 K, 1300 K,1200 K, 1100 K, 1000 K, 900 K, 800 K, 700 K, 600K, 500 K, 400 K 和300 K 下分别运行2000 个离子步后获得非晶态结构,同时调整模拟盒子的大小保持压强为0 GPa. 每个温度下, 使用最后1000 个构型进行结构分析. 请问您,1. 使用NVT系综,设置AMASS=0,是否意味着T保持不变,TEBEG要等于TEEND?2.降温时候用如何设置?AMASS=-1?3. 里面的“调整模拟盒子的大小保持压强为0 GPa“如何设置?4. 此外,“每个温度下, 使用最后1000 个构型进行结构分析”这就话什么含义,我不是很明白;5. 最后,请问如何分析键对和维诺多面体?希望您百忙之中能指点指点,真心谢谢!
7楼2019-10-19 11:23:52
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 高朋满座007 的主题更新
信息提示
请填处理意见