| 查看: 2637 | 回复: 11 | ||||
| 当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖 | ||||
[交流]
【求助】利用侯博的方法优化结构出现的问题（求验证）已有5人参与
|
||||
|
利用侯博的文章中优化的部分进行优化，照例对ENCUT，KPT，和SIGMA进行收敛测试，得到ENCUT=650eV，K=4*4*4，SIGMA=0.1。然后采用结构参数优化的复杂情况，因为体系为正交。 其中INCAR.relax： SYSTEM = NH3BH3-pmn21 ENCUT = 650 ISTART = 0; ICHARG = 2 ISMEAR = 1; SIGMA = 0.1 NSW = 60; IBRION = 2 ISIF = 5 POTIM = 0.2 EDIFF = 1E-5; EDIFFG = -1E-3 PREC = Accurate INCAR.static： SYSTEM = NH3BH3-pmn21 ENCUT = 650 ISTART = 0; ICHARG = 2 ISMEAR = -5 PREC = Accurate KPT： Auto generation 0 Gamma 8 8 8 0.0 0.0 0.0 POTCAR为： paw_GGA 但是计算得到的结果居然为： 5.041 98.6200 -62.534945 5.141 104.6000 -65.167658 5.241 110.8300 -67.385080 5.341 117.2900 -69.230747 5.441 124.0100 -70.744046 5.541 130.9700 -71.959857 5.641 138.1900 -72.909812 5.741 145.6700 -73.622417 5.841 153.4200 -74.123520 。。。 与理论期望不符。。。 望高手指点 |
回复此楼» 收录本帖的淘帖专辑推荐
 VASP |
» 猜你喜欢
博士读完未来一定会好吗
已经有24人回复
-
导师想让我从独立一作变成了共一第一
已经有9人回复
-
到新单位后，换了新的研究方向，没有团队，持续积累2区以上论文，能申请到面上吗
已经有11人回复
-
读博
已经有4人回复
-
JMPT 期刊投稿流程
已经有4人回复
-
心脉受损
已经有5人回复
-
Springer期刊投稿求助
已经有4人回复
-
小论文投稿
已经有3人回复
-
申请2026年博士
已经有6人回复
» 本主题相关价值贴推荐，对您同样有帮助:
重金请教采用准谐近似求热膨胀系数的方法
已经有16人回复- 【求助】计算机辅助药物设计
已经有9人回复
- 【求助成功】计算能带结构时出现错误，请教
已经有19人回复
ApolloYang
木虫 (正式写手)
- 应助: 15 (小学生)
- 金币: 2209.6
- 散金: 114
- 红花: 4
- 帖子: 390
- 在线: 236.9小时
- 虫号: 582657
- 注册: 2008-07-25
- 性别: GG
- 专业: 凝聚态物性 II ：电子结构
★ ★
小木虫(金币+0.5):给个红包，谢谢回帖交流
zzy870720z(金币+1):谢谢提示 2010-07-27 16:26:45
小木虫(金币+0.5):给个红包，谢谢回帖交流
zzy870720z(金币+1):谢谢提示 2010-07-27 16:26:45
|
引用侯博士的话 “对六角晶格, E是 a,c的函数，根据体积V与a,c的关系，可以得到E也是a, V的函数，或者E是V, c的函数或 E是V, c/a的函数。明白了这个几个关系之后就能大致明白那个例子里面所展示的内容。考虑的是计算不同V下面的总能，对给定一个V，初始化一个a, c/a，保持体积不变去优化优化a, c/a。计算了一组V下的总能，就有了E~V的数据，然后采用状态方程拟合，得到平衡体积V0，然后再初始化该体积下的a, c/a再进行固定体积的优化得到平衡态下的a0, c0/a0。” |
6楼2010-07-27 15:13:19
valenhou001
至尊木虫 (职业作家)
- 1ST强帖: 13
- 应助: 241 (大学生)
- 金币: 25701.7
- 散金: 602
- 红花: 166
- 帖子: 3782
- 在线: 873.8小时
- 虫号: 1007127
- 注册: 2010-04-27
- 专业: 凝聚态物性 II ：电子结构
2楼2010-07-26 18:07:30
nkleof
银虫 (正式写手)
- 应助: 0 (幼儿园)
- 金币: 2093.2
- 散金: 649
- 红花: 5
- 帖子: 813
- 在线: 165.9小时
- 虫号: 152581
- 注册: 2006-01-02
- 专业: 凝聚态物性 II ：电子结构
★
小木虫(金币+0.5):给个红包，谢谢回帖交流
小木虫(金币+0.5):给个红包，谢谢回帖交流
|
跟帖请问一下，在优化时遇到所得晶格常数随ENCUT值增大而增大的情况（从ENCUT=300到800eV，晶格常数从9.836涨到9.852，而实验值在9.814），问题可能出在哪里？采用的优化方法如下： INCAR --》 SYSTEM = CuCr2S4 cubic cell spinel PREC = High ENCUT = 900 LREAL = A IALGO = 48 NPAR = 8 NSIM = 8 ISTART = 0 ICHARG = 2 NELMIN = 4 LORBIT = 12 EDIFF = 1E-5 ISMEAR = -5 LWAVE = .False. NSW = 60 IBRION = 2 ISIF = 3 ISPIN = 2 MAGMOM = 2*0 4*4 8*0 KPOINTS --》 Auto 0 M 13 13 13 0 0 0 k点网格做过收敛测试，从此往后更密的网格造成的总能变化在0.00015eV/cell，每原胞14个原子。 体系是准金属性的，但是在费米面处的态密度值很小，且是磁性体系。是否有可能是ISMEAR=-5造成的问题？ |
3楼2010-07-26 19:20:38
valenhou001
至尊木虫 (职业作家)
- 1ST强帖: 13
- 应助: 241 (大学生)
- 金币: 25701.7
- 散金: 602
- 红花: 166
- 帖子: 3782
- 在线: 873.8小时
- 虫号: 1007127
- 注册: 2010-04-27
- 专业: 凝聚态物性 II ：电子结构
4楼2010-07-27 09:11:41