24小时热门版块排行榜    

查看: 3680  |  回复: 7

luo019922

铁虫 (初入文坛)

[求助] VASP高温高压AIMD后结构优化问题 已有2人参与

请教各位大神,我用VASP在高温高压(3000K,50GPa)下做AIMD模拟(NVT)熔体结构达到平衡之后,我从平衡的结构中选取一个snapshot去做结构优化之后得到的是我在高温高压下的稳定结构吗?我知道结构优化可以说是基于0K的,但我的理解是高温的AIMD其实每一个离子步中也是固定离子再对电子波函数迭代的,所以不太理解能否在高温AIMD之后再结构优化?希望得到指导,感谢!
回复此楼
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

zhangfrank

金虫 (正式写手)

【答案】应助回帖

★
感谢参与,应助指数 +1
luo019922: 金币+1, ★有帮助 2017-03-01 18:17:59
我也是新手 只说一点个人理解
我看到的类似工作,都是把高温下的一个结构,拿来做降温处理(有的降温到0,有的降温到300k),然后把降温后的做结构弛豫(DFT计算,在0k)

欢迎交流,互相帮助
计算顺利点吧
2楼2017-02-27 21:54:31
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

三叶虫的脚印

新虫 (初入文坛)

高温高压AIMD后,结构优化是可以的。只是可能会得到亚稳相的结构。
3楼2017-03-01 18:12:55
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

三叶虫的脚印

新虫 (初入文坛)

高温高压AIMD后,结构优化是可以的。只是可能会得到亚稳相的结构。
4楼2017-03-01 18:15:59
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

三叶虫的脚印

新虫 (初入文坛)

高温高压AIMD后,结构优化是可以的。只是可能会得到亚稳相的结构。
5楼2017-03-01 18:19:31
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

卡开发发

专家顾问 (著名写手)

Ab Initio Amateur

引用回帖:
2楼: Originally posted by zhangfrank at 2017-02-27 21:54:31
我也是新手 只说一点个人理解
我看到的类似工作,都是把高温下的一个结构,拿来做降温处理(有的降温到0,有的降温到300k),然后把降温后的做结构弛豫(DFT计算,在0k)

欢迎交流,互相帮助

那个目的大概是为了退火得到稳定结构,就是升温降温的循环少了点。
不一定挂在论坛,计算问题问题欢迎留言。
6楼2017-03-01 20:28:01
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

卡开发发

专家顾问 (著名写手)

Ab Initio Amateur

如果只是为了得到结构的话,你测量的情况显然是就应当是系综平均的结构,都高温高压下了,你很难保证随着温度增加系统平均的体积会如何变化,这些结构不太可能老老实实待在稳定状态上,倘若直接考虑优化显然不能描述体积膨胀的过程:

低温下或者只考虑简谐近似的话,平均结构显然是与直接优化没什么区别,因为简谐近似下势能面在平衡点附近是近似对称的。

随着温度增加,结构的分布可以越偏离平衡结构,但由于势能面本身就不是对称的,平均结构就会朝一个方向偏,所以导致晶格膨胀,这个问题如果只是进行普通优化就能描述的。这还不是稳态与亚稳态的问题。

其余性质其实也得建立在有限温度的过程,3000K和高压环境还是应当慎重处理。

这应该是对你之前的问题做了个详细的说明,还有楼上的竟然把别人的答案直接抄来,这样也可以?
不一定挂在论坛,计算问题问题欢迎留言。
7楼2017-03-01 20:44:04
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

dxcharlary

专家顾问 (职业作家)

【答案】应助回帖

★ ★ ★ ★ ★
感谢参与,应助指数 +1
luo019922: 金币+5, ★★★很有帮助 2017-03-22 23:40:24
你做relax不是相当于退火吗?
你是要单纯通过这种方式来找结构的话,relax是可行的。但是如果你想通过这种方式模拟实验的话不行。
8楼2017-03-02 09:30:21
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 luo019922 的主题更新
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[公派出国] 售SCI文章,我:8O5.5.1.O.54,科目齐全,可+急 +4 1uaZt2kAZj2S 2026-10-31 7/350 2026-11-01 12:40 by BCdvBAEpt3zm
[考研] 售SCI一区T0P文章,我:8.O.55.1.O54,科目全,可伽急 +4 1uaZt2kAZj2S 2026-10-31 8/400 2026-11-01 12:38 by BCdvBAEpt3zm
[找工作] 售SCI一区T0P文章,我:8.O.55.1.O54,科目全,可伽急 +3 1uaZt2kAZj2S 2026-10-30 9/450 2026-11-01 12:32 by BCdvBAEpt3zm
[考博] 售SCI文章,我:8O5.5.1.O.54,科目齐全,可+急 +3 1uaZt2kAZj2S 2026-10-30 10/500 2026-11-01 12:27 by BCdvBAEpt3zm
[论文投稿] 售SCI-T0P文章,我:8O.5.5.1.O.54,科目齐全,可+急 +3 1uaZt2kAZj2S 2026-10-30 12/600 2026-11-01 12:19 by BCdvBAEpt3zm
[硕博家园] 售一区SCI文章T0P,我:8O.551.O54,科目全,可十急 +3 1uaZt2kAZj2S 2026-10-30 13/650 2026-11-01 12:14 by BCdvBAEpt3zm
[硕博家园] 售SCI一区T0P文章,我:8O.55.1.O.54,科目全,可伽急 +4 X5bHxOULxFqb 2026-10-30 18/900 2026-11-01 12:11 by BCdvBAEpt3zm
[公派出国] 售SCI一区T0P文章,我:8.O.55.1.O54,科目全,可伽急 +3 hKe7U1kad5Sg 2026-10-31 5/250 2026-11-01 11:14 by BCdvBAEpt3zm
[考博] 售SCI一区文章,我:8.O.55.1.O.54,科目齐全,可伽急 +3 xvJS1BYDxvBS 2026-10-31 4/200 2026-11-01 11:01 by BCdvBAEpt3zm
[考研] 售SCI一区文章,我:8O5.5.1.O5.4,科目全,可伽急 +4 xvJS1BYDxvBS 2026-10-31 9/450 2026-11-01 10:56 by BCdvBAEpt3zm
[硕博家园] 售SCI一区T0P文章,我:8.O55.1.O.54,科目全,可十急 +4 xvJS1BYDxvBS 2026-10-31 9/450 2026-11-01 10:53 by BCdvBAEpt3zm
[考博] 售SCI一区文章,我:8.O.55.1.O.54,科目齐全,可伽急 +4 1uaZt2kAZj2S 2026-10-30 16/800 2026-11-01 10:27 by BCdvBAEpt3zm
[考博] 售SCI一区T0P文章,我:8.O.55.1.O.5.4,科目全,可+急 +4 X5bHxOULxFqb 2026-10-30 21/1050 2026-11-01 06:46 by YXM1FDCSE8Bn
[基金申请] 售SCI-T0P文章,我:8O.5.5.1.O.54,科目齐全,可+急 +5 X5bHxOULxFqb 2026-10-30 21/1050 2026-11-01 06:46 by YXM1FDCSE8Bn
[博后之家] 售SCI文章,我:8O5.5.1.O.54,科目齐全,可+急 +3 1uaZt2kAZj2S 2026-10-30 15/750 2026-11-01 05:22 by YXM1FDCSE8Bn
[论文投稿] 售一区SCI文章T0P,我:8O.551.O54,科目全,可十急 +3 1uaZt2kAZj2S 2026-10-31 8/400 2026-11-01 03:59 by YXM1FDCSE8Bn
[考研] 售SCI-T0P文章,我:8O.5.5.1.O.54,科目齐全,可+急 +3 1uaZt2kAZj2S 2026-10-30 12/600 2026-11-01 02:19 by YXM1FDCSE8Bn
[硕博家园] 越来越担心硕士毕业后没有竞争力了 +8 zuesli 2026-10-27 14/700 2026-10-30 22:44 by 好运连连公主
[文学芳草园] 闲聊 +3 myrtle 2026-10-25 6/300 2026-10-29 15:28 by myrtle
[考博] 申博求助 20+3 Katrinae 2026-10-26 4/200 2026-10-27 09:56 by 北京莱茵编辑
信息提示
请填处理意见