| 查看: 1872 | 回复: 4 | |||
泽润东方金虫 (正式写手)
|
[求助]
关于缩放系数的优化 已有1人参与
|
|
关于VASP计算我些问题想请教各位。 结构弛豫：空间群是Pnma，三个晶格常数互不相等。 实验方面的结构数据如下： **************************** a = 0.3846 nm, b = 1.1268 nm, c = 1.1321 nm, α = 90 °, β = 90 °, γ = 90 ° 。 **************************** 以实验结果为基础开始优化，采用GGA-PBE泛函。步骤如下： 1.从0.98-1.8之间选取了二十多个数值作为缩放系数进行了一系列的结构弛豫。 2.选择步骤1中能量最低的结果，作为初始结构，首先固定a方向晶格常数，保持b/a和c/a不变（ISIF=4），进行晶格常数优化，但是保持a，b，c之间的夹角不变；其次再调整a方向的常数，重复上面的步骤。反复进行了十多个测试。 3.最后选择步骤2中得到的结构作为最稳定的结构。如下： **************************** direct lattice vectors 3.915449903 0.000000000 0.000000000 0.000000000 11.390399806 0.000000000 0.000000000 0.000000000 11.492100194 **************************** 4.在步骤3的基础上重新做了结构弛豫，测试了K点、Ecut。最后确定下来结构和Ecut。 之后的静待自洽、能带以及态密度的计算都是步骤的结构为基础进行的。 目前我比较了一下我的结构跟之前发表的理论值与实验值。误差稍微有些大： present EV-GGA（已发表的理论值） Exp1 Exp2 Exp3 a(angstr) 3.87 3.84 3.85 3.83 3.84 b(angstr) 11.21 11.22 11.27 11.23 11.22 c(angstr) 11.96 11.30 11.32 11.29 11.31 Vol(angstr3) 519.8 - 490.62 - - 我的问题如下： a. 在此基础上所做的能带和态密度可信度高么 ？ b. 如果不可信，应如何从新开始结构参数的优化 ？ c. 出现计算误差的原因有哪些 ？ |
回复此楼» 猜你喜欢
278求调剂
已经有3人回复
-
301求调剂
已经有3人回复
-
275求调剂
已经有3人回复
-
一志愿，福州大学材料专硕339分求调剂
已经有3人回复
-
本人考085602 化学工程 专硕
已经有15人回复
-
283求调剂
已经有3人回复
-
274求调剂
已经有5人回复
-
302求调剂
已经有6人回复
-
11408 一志愿西电，277分求调剂
已经有3人回复
-
化学工程321分求调剂
已经有13人回复
» 本主题相关价值贴推荐，对您同样有帮助:
- 关于matlab的参数估计
已经有15人回复
- VASP 计算E-V曲线求助。。
已经有10人回复
- 求一份详细的产品设计开发项目策划书
已经有7人回复
- 【求助】怎么用晶胞和原子坐标等信息写POSCAR
已经有38人回复
- 【求助】结构优化中，如何选择POTIM？
已经有7人回复
- 【原创】晶格优化中ISIF取值的探究（六方MgO为例）
已经有9人回复
★ ★ ★ ★ ★ ★ ★
感谢参与，应助指数 +1
liliangfang: 金币+2, 谢谢交流 2014-10-26 19:13:30
泽润东方: 金币+5, ★有帮助 2014-11-18 13:38:20
感谢参与，应助指数 +1
liliangfang: 金币+2, 谢谢交流 2014-10-26 19:13:30
泽润东方: 金币+5, ★有帮助 2014-11-18 13:38:20
 本帖仅楼主可见 |
2楼2014-10-23 21:12:30
泽润东方
金虫 (正式写手)
- 应助: 2 (幼儿园)
- 金币: 1702.8
- 散金: 406
- 帖子: 470
- 在线: 222.7小时
- 虫号: 684187
- 注册: 2008-12-29
- 性别: GG
- 专业: 基础物理学
3楼2014-11-18 13:41:31
4楼2014-11-19 10:40:28
对马拉二胡
至尊木虫 (职业作家)
- 应助: 27 (小学生)
- 金币: 25915.6
- 散金: 244
- 红花: 13
- 帖子: 4065
- 在线: 938小时
- 虫号: 2471734
- 注册: 2013-05-19
- 性别: GG
- 专业: 无机非金属材料
5楼2014-11-21 19:11:05
