24小时热门版块排行榜    

查看: 3667  |  回复: 13

Gina88

木虫 (正式写手)

[求助] VASP分子动力学优化的INCAR如下,怎么CONTCAR和POSCAR没有什么变化啊?

但是文献中结构退火后有很大的变化的,为什么我这里没有什么变化啊?
我的INCAR设置有什么错误吗?

ISTART=0;ICHARG=2

PREC=Accurate
ENCUT=300
GGA=91;VOSKOWN=1
ISMEAR=0;SIGMA=0.05

LREAL=A
ALGO=Veryfast;IALGO=48

ISIF=4
ISYM=0

NSW=1000
IBRION=0;POTIM=2
TEBEG=875;TEEND=875
SMASS=0

NBLOCK=50;KBLOCK=1
LCHARG=F
回复此楼

» 猜你喜欢

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

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

已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
回帖支持 ( 显示支持度最高的前 50 名 )

enola

捐助贵宾 (正式写手)

★ ★ ★
ben_ladeng(金币+3): 楼主反映很好,追加奖励 2011-08-19 09:44:25
引用回帖:
7楼: Originally posted by Gina88 at 2011-08-09 08:32:19:
麻烦帮看一下上贴呗。
多谢大侠了。
上贴中第一图是所建模型在退火之前的样子。
第二图是退火之后的样子。
好像好有动态的快照。
不知道这个动态的快照怎么得到的。

仅仅看那段文字,不是很清楚,我猜想建好结构后,作者先是结构优化(0K),然后为了验证结构稳定性,通过几千步的MD,把温度从0K升到875K(600度),来看结构的变化,模拟结果是只是结构热膨胀拉。。
timeflies..
8楼2011-08-09 11:39:12
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
普通回帖

enola

捐助贵宾 (正式写手)

★ ★
franch(金币+2): 谢谢回帖交流 2011-08-09 13:34:04
不知道楼主想算什么啊?退火嘛?但是你的输入文件不是退火啊:
开始温度和最后温度都一样的。
timeflies..
2楼2011-08-08 17:05:22
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

uuv2010

荣誉版主 (职业作家)

优秀版主

★ ★
franch(金币+2): 谢谢回帖交流 2011-08-09 13:34:13
呵呵,是啊
TEBEG=875;TEEND=875
为什么这样?
3楼2011-08-08 17:16:56
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

Gina88

木虫 (正式写手)

引用回帖:
3楼: Originally posted by uuv2010 at 2011-08-08 17:16:56:
呵呵,是啊
TEBEG=875;TEEND=875
为什么这样?

我想温度不变啊,就像经典分子动力学中的dynamics计算
4楼2011-08-08 17:57:44
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

uuv2010

荣誉版主 (职业作家)

优秀版主


zzy870720z(金币+1): 鼓励交流 2011-08-10 14:26:16
如果温度不变的话就不是模拟退火了,模拟退火应该温度降低,即TEBEG>TEEND
我没做过,希望做过模拟退火的虫子们给予指导。
5楼2011-08-08 18:20:43
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

Gina88

木虫 (正式写手)

引用回帖:
5楼: Originally posted by uuv2010 at 2011-08-08 18:20:43:
如果温度不变的话就不是模拟退火了,模拟退火应该温度降低,即TEBEG>TEEND
我没做过,希望做过模拟退火的虫子们给予指导。

但是文献中没有说明有两个温度啊,所以我不知道另一个温度是什么。

6楼2011-08-09 08:30:19
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

Gina88

木虫 (正式写手)

引用回帖:
2楼: Originally posted by enola at 2011-08-08 17:05:22:
不知道楼主想算什么啊?退火嘛?但是你的输入文件不是退火啊:
开始温度和最后温度都一样的。

麻烦帮看一下上贴呗。
多谢大侠了。
上贴中第一图是所建模型在退火之前的样子。
第二图是退火之后的样子。
好像好有动态的快照。
不知道这个动态的快照怎么得到的。
7楼2011-08-09 08:32:19
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

enola

捐助贵宾 (正式写手)

★ ★ ★
franch(金币+3): 鼓励交流 2011-08-09 13:34:59
引用回帖:
7楼: Originally posted by Gina88 at 2011-08-09 08:32:19:
麻烦帮看一下上贴呗。
多谢大侠了。
上贴中第一图是所建模型在退火之前的样子。
第二图是退火之后的样子。
好像好有动态的快照。
不知道这个动态的快照怎么得到的。

动态快照的话好弄,每一步会有一个坐标,把每一步坐标记录下来,些个Xcrysden的输入文件就好啦。。
timeflies..
9楼2011-08-09 11:40:30
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

Gina88

木虫 (正式写手)

ben_ladeng: 可以通知版主奖励应助人的。现在退金币给你哈。 2011-08-19 09:45:20
引用回帖:
8楼: Originally posted by enola at 2011-08-09 11:39:12:
仅仅看那段文字,不是很清楚,我猜想建好结构后,作者先是结构优化(0K),然后为了验证结构稳定性,通过几千步的MD,把温度从0K升到875K(600度),来看结构的变化,模拟结果是只是结构热膨胀拉。。

不是说温度和速度是成比例的吗,温度越高速度越大吗。
所以一般的优化是找到势能局域最小值。而分子动力学则可能找到势能全局最小值啊。
那我的这个设置就是需要将开始的温度改成0吗?
那么这样MD模拟后该怎么看动态图呢,我还是不清楚,能麻烦您详细告诉我一下吗?

还有,非常感谢您的回答,但是能不能麻烦您将帖子设成应助贴呢,呵呵,希望用金币表达一点谢意:)
10楼2011-08-10 14:06:55
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 Gina88 的主题更新
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[教师之家] 评正教授需要两个国家级项目,有人用子课题糊弄,结果在评审前资格公示时被举报拿下了 +22 瞬息宇宙 2024-11-12 31/1550 2024-11-15 22:02 by 潇湘之迷
[考博] 中南大学 粉末冶金国家重点实验室 陈超教授课题组拟招收 1~2位博士研究生通知 +7 中南大学-金属材 2024-11-14 18/900 2024-11-15 20:22 by 中南大学-金属材
[论文投稿] 职称论文 +10 xiaoqi_66 2024-11-09 13/650 2024-11-15 16:45 by sugarmei
[教师之家] 大学老师 +11 考研一路顺风 2024-11-13 14/700 2024-11-15 15:24 by LNP@mRNA
[基金申请] 去年七月底入站的还能申请下一批吗? +4 brightwo 2024-11-14 4/200 2024-11-15 14:18 by gazi1111
[硕博家园] 实验进度 +9 麻辣烫。 2024-11-11 14/700 2024-11-15 13:14 by hanweiguo
[有机交流] 二甲胺的使用 20+5 太阳谷 2024-11-14 13/650 2024-11-15 13:11 by xiaomei1031
[考博] 联系的博导说有时间聊一聊 我回复了自己可以的时间后 博导就没再回复 +4 evan0489 2024-11-15 4/200 2024-11-15 13:10 by hanweiguo
[硕博家园] 课题和专业一点都不相关怎么办 +7 imt4n9 2024-11-11 10/500 2024-11-15 01:55 by 成长虫
[论文投稿] OE返修遇到expired +3 隔壁老王来了 2024-11-14 7/350 2024-11-14 20:46 by 隔壁老王来了
[论文投稿] 投稿系统中的通讯作者和文章中的通讯作者不一样,文章目前被录用了? +4 babybabygo 2024-11-12 5/250 2024-11-14 19:13 by 走了002
[教师之家] 处在人生职业的分水岭 +4 otani 2024-11-13 4/200 2024-11-14 14:17 by mddzwo
[有机交流] 硫谜用双氧水氧化出现的困惑(老板说一定要用双氧水氧化做好)求指教 10+4 luowenwu 2024-11-10 15/750 2024-11-14 10:52 by luolinfeng
[基金申请] 博后面上出结果了 +5 wl87925139 2024-11-14 5/250 2024-11-14 10:26 by kingmax18996
[论文投稿] 要不要撤稿另投 10+4 wangzhesd 2024-11-09 8/400 2024-11-14 09:09 by LIU_V
[论文投稿] 求JACS和ACIE的最新文章模板 10+4 Jane060390 2024-11-10 4/200 2024-11-13 09:13 by 北京莱茵润色
[硕博家园] 如何看清华北大硕博的毕业论文 +6 2025lucky 2024-11-10 6/300 2024-11-12 22:46 by 楚扬长流
[论文投稿] 论文接受后,编辑部发邮件提醒缺少附件 15+4 上善若水明泪 2024-11-10 4/200 2024-11-12 14:24 by 北京莱茵润色
[论文投稿] 爱思唯尔投稿系统里的通讯作者可以和文章里的通讯作者标注不同吗 +7 Omnissiah 2024-11-10 7/350 2024-11-12 14:07 by holypower
[基金申请] 要求论文发表日期在项目执行期内,论文发表日期是在线日期还是见刊日期 +6 Jiangnanyu1 2024-11-09 7/350 2024-11-10 17:49 by jurkat.1640
信息提示
请填处理意见