| 查看: 3223 | 回复: 14 | ||
| 【悬赏金币】回答本帖问题，作者huyaoqiao将赠送您 30 个金币 | ||
| 当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖 | ||
[求助]
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. |
回复此楼» 猜你喜欢
来用亲身体会，一起说说年轻老师的辛苦
已经有25人回复
-
连续3年国社科青年本子都没中，今年最后一次青年了，写本子完全浪费时间
已经有18人回复
-
双非博士还是985科研助理，然后再读博
已经有5人回复
-
Nano-Micro Letters投稿可以用PDF上传吗
已经有4人回复
-
博三一直没文章怎么办
已经有30人回复
-
推荐转投（ transfer pending）是否有用？
已经有7人回复
-
评审规则突发奇想
已经有25人回复
-
“直接受聘正高专业技术职务”怎么理解
已经有9人回复
-
蹲一个男朋友
已经有29人回复
-
又见奇葩学生，可咋办
已经有26人回复
» 本主题相关商家推荐: (我也要在这里推广)
valenhou001
至尊木虫 (职业作家)
- 1ST强帖: 13
- 应助: 241 (大学生)
- 金币: 25632.7
- 散金: 602
- 红花: 166
- 帖子: 3768
- 在线: 873.7小时
- 虫号: 1007127
- 注册: 2010-04-27
- 专业: 凝聚态物性 II ：电子结构
3楼2020-01-30 12:00:18
10