24小时热门版块排行榜    

查看: 3259  |  回复: 14
【悬赏金币】回答本帖问题,作者huyaoqiao将赠送您 30 个金币
当前只显示满足指定条件的回帖,点击这里查看本话题的所有回帖

huyaoqiao

铁虫 (小有名气)

[求助] VASP MD 得到 amorphous结构已有1人参与

大家好。我想利用VASP跑分子动力学模拟得到amorphous结构。思路大概是从crystalline phase起先melt再quench。有一个问题是,amorphous phase的密度一般比crystalline phase密度要低,所以其晶格常数要大一些。如果始终用NVT ensemble的话晶格常数不会改变,这样得到的结构不太合理,因为得到的amorphous 结构密度和crystalline phase一样,会存在很大 residual stress。所以我试着用NpT ensemble (ISIF=3),保持压强为一个大气压,期望晶格常数会有所增加,得到合理的结构。现在问题是,NpT  MD时候晶格常数一直不变,请问怎么回事?

INCAR file:
# General parameters for this run
System= SnO2
NCORE = 12
# Electronic specification

PREC  = Normal    ! standard precision, we recomend the use of PREC = Normal
ENCUT = 150       ! cutoff should be set manually, smaller if ultra-large supercell
ISMEAR= 0
SIGMA = 0.01
LREAL = Auto      ! Real space to spped up for large molecular dynamics runs
ALGO  = Very Fast ! For large molecular dynamics runs
NELMIN= 4         ! Set minimum steps of electronic optimization
EDIFF = 1E-3      ! Do not use default, too much drift

# MD related
# NpT

ISYM  = 0         ! strongly recommened for MD, no symmertry
IBRION= 0         ! molecular dynamics
NSW   = 7500      ! 7500 steps
POTIM = 2         ! timestep 2 fs

# SMASS = -1        ! Nose Hoover thermostat
ISIF  = 3
MDALGO= 3         !
TEBEG = 3000
TEEND = 300       ! temperature
PSTRESS= 0.001

# Output control

NBLOCK= 50        ! XDATCAR updated every 50 steps
NWRITE = 0        ! For MD, use NWRITE = 0 or NWRITE = 1
LCHARG = .FALSE.
LWAVE  = .FALSE.
回复此楼
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

叶石纹

新虫 (初入文坛)

这个可以对模拟元胞进行缩放,计算不同体积对应的压力,拟合p-v曲线,得到你想要的压力对应的体积。我最近的文章计算方法里有这方面的内容,可以参考https://doi.org/10.1016/j.molliq.2019.112063

发自小木虫Android客户端
5楼2020-01-31 11:36:28
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
查看全部 15 个回答

Nano2016

新虫 (著名写手)

你们是做什么材料的SnO2?

发自小木虫IOS客户端
2楼2020-01-30 12:00:03
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

valenhou001

至尊木虫 (职业作家)

对不同密度的超原胞分别进行升温熔化,退火降温的模拟,取样统计分析。

发自小木虫Android客户端
3楼2020-01-30 12:00:18
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

时光!

新虫 (正式写手)

4楼2020-01-30 12:36:02
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
不应助 确定回帖应助 (注意:应助才可能被奖励,但不允许灌水,必须填写15个字符以上)
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[教师之家] 博士高校求职 安建大vs西科大 +5 chengmy19 2024-06-01 12/600 2024-06-02 19:35 by kanyechris
[论文投稿] 希望嫡长子顺利出场! +5 C刊霸王 2024-05-31 6/300 2024-06-02 19:26 by 子非愚
[基金申请] +5 河马の史诗 2024-06-02 5/250 2024-06-02 19:05 by 328838485
[论文投稿] 没收到邮件 10+4 荣小撇 2024-05-31 9/450 2024-06-02 13:25 by bobvan
[硕博家园] 论大家对6070后普通教授导师的看法 +7 SNaiL1995 2024-05-28 11/550 2024-06-02 13:16 by 小虫子咔咔
[硕博家园] 每到中夜,情难自抑 +34 sioc-sunj 2024-05-28 58/2900 2024-06-02 02:53 by csyky2007
[考博] 申请2024或2025年博士研究生 +5 嘟噜嘟1 2024-05-29 12/600 2024-06-01 22:36 by 嘟噜嘟1
[考博] 24or25材料专业申博 +4 农夫三拳有点痛 2024-05-30 11/550 2024-06-01 14:45 by Napoleonsky
[论文投稿] 求Sci期刊推荐 10+4 甄小鱼 2024-05-30 7/350 2024-06-01 10:41 by bobvan
[基金申请] 数理的人才答辩通知发了吗 +7 yzy3327 2024-05-30 7/350 2024-06-01 09:15 by 撒哈拉西部
[硕博家园] 哈工大硕博招生!博士每月入学! +4 nailooo 2024-05-30 5/250 2024-06-01 06:47 by anevay
[论文投稿] 纠结选哪一个期刊,电化学领域 50+8 Freya163 2024-05-28 10/500 2024-05-31 15:09 by wzykobe
[文学芳草园] 对对子啊 +5 天若孤独 2024-05-29 7/350 2024-05-31 09:00 by wjykycg
[材料综合] 真空封石英管 北京 +4 dessha 2024-05-29 5/250 2024-05-30 16:40 by mpdfwxgui
[有机交流] 液相纯度高,但产品析不太出来 10+4 cui19236 2024-05-27 9/450 2024-05-30 07:45 by yuanjijoy
[博后之家] 2024公派博后申请 +4 326lhpqk 2024-05-27 5/250 2024-05-29 20:03 by @古月胡
[论文投稿] 有没有老师需要发表论文 +4 金老师论文助理- 2024-05-29 4/200 2024-05-29 16:51 by liuyupu132
[论文投稿] 真急着毕业,CPB主编终审17天了,邮件催稿了两次,就是一点动静没有 5+3 kkkk夏 2024-05-28 6/300 2024-05-29 11:18 by hitsdu
[论文投稿] 核心初审被拒,理由是“选题的意义不明确,文章写得不像是科技论文”,怎么改 5+3 工藤雷花樱 2024-05-27 8/400 2024-05-29 10:09 by topedit
[有机交流] 奇怪的物质 100+4 桃桃PXS 2024-05-27 7/350 2024-05-28 10:22 by 091602
信息提示
请填处理意见