24小时热门版块排行榜

返回列表

当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖

huyaoqiao

铁虫 (小有名气)

应助: 23 (小学生)
金币: 442.9
散金: 50
红花: 12
帖子: 287
在线: 94.7小时
虫号: 3446117
注册: 2014-09-28
专业: 半导体材料

[求助] 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.

回复此楼

» 猜你喜欢

求个博导看看已经有17人回复
青基代表作，AAAI之类的A会的special track在国内认可度高吗？还是归为workshop之流？已经有3人回复
上海工程技术大学【激光智能制造】课题组招收硕士已经有6人回复
带资进组求博导收留已经有11人回复
自荐读博已经有5人回复
上海工程技术大学张培磊教授团队招收博士生已经有4人回复
求助院士们，这个如何合成呀已经有4人回复
临港实验室与上科大联培博士招生1名已经有9人回复
写了一篇“相变储能技术在冷库中应用”的论文，论文内容以实验为主，投什么期刊合适？已经有6人回复
最近几年招的学生写论文不引自己组发的文章已经有11人回复

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

1楼 2020-01-30 11:51:50

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

叶石纹

新虫 (初入文坛)

应助: 1 (幼儿园)
金币: 97.3
帖子: 45
在线: 79.5小时
虫号: 3945685
注册: 2015-06-29
性别: GG
专业: 化学工程及工业化学

引用回帖:

5楼: Originally posted by 叶石纹 at 2020-01-31 11:36:28
这个可以对模拟元胞进行缩放，计算不同体积对应的压力，拟合p-v曲线，得到你想要的压力对应的体积。我最近的文章计算方法里有这方面的内容，可以参考https://doi.org/10.1016/j.molliq.2019.112063
...

据我所知，楼主这个并不是npt，vasp现在可用的是郎之万npt，参数设置明显不是楼主这个。

发自小木虫Android客户端